图层大小怎么不一样，还有这个是涡度，那平流怎么算

caojunwei

'reinit'
'sdfopen E:\ncep-data\data\2016\uwnd.2016.nc'
'sdfopen E:\ncep-data\data\2016\vwnd.2016.nc'
'sdfopen E:\ncep-data\data\2016\hgt.2016.nc'
'E:/ncep-data/red_blue3.gs'
'enable print E:\pic\vort\500hpa.gmf'
year='2016'
mon='Apr'
lev='500'
t=1
while(t<3)
hr=0
while(hr<=23)
'c'
'set time 'hr'z't''mon''year''
'set lon 80 150'
'set lat 30 70'
'set z 6'
'set grid off'
'set gxout shaded'
'define dv=cdiff(vwnd.2,x)'
'define dx=cdiff(lon,x)*3.1416/180'
'define du=cdiff(uwnd.1*cos(lat*3.1416/180),y)'
'define dy=cdiff(lat,y)*3.1416/180'
'd (dv/dx-du/dy)/(6.37E6*cos(lat*3.1416/180))'
cbarn.gs
'set gxout contour'
'set cint 4'
'set ccolor 1'
'set cstyle 1'
'set cthick 5'
'set clskip 2'
'd hgt.3/10'
'draw title vort, 'lev'hPa,  'hr'z't''mon''year''
*****************************
hr=hr+12
'print'
endwhile
t=t+1
endwhile
'disable print'
*-----------------------------

a阿克

图层大小不同是指填色和地图范围不同吗？因为算涡度用的是cdiff这个中央差分，四侧的格点肯定不加入运算。你可以在计算涡度的时候把范围设置大一点，要填色的时候再将范围缩小。
比如还是要画这个涡度场。
'set lon 70 160'
'set lat 20 80'
'set lev 500'
'define vor=hcurl(uwnd.1,vwnd.2)'
'set lon 80 150'
'set lat 30 70'
'set gxout shaded'
'd vor'

caojunwei

a阿克 发表于 2017-3-13 16:26
图层大小不同是指填色和地图范围不同吗？因为算涡度用的是cdiff这个中央差分，四侧的格点肯定不加入运算。 ...

感谢感谢

karenlk

出图的范围要小于计算的范围

caojunwei

karenlk 发表于 2017-3-14 11:00
出图的范围要小于计算的范围

嗯嗯 弄明白了