求助！不同点不同时次的剖面作法

jlmjlm2007

想做不同点不同时次的风场，结果老是显示Cannot contour grid - all undefined values  

以下是我写的gs
'reinit'
'open D:\fnl\fnl.ctl'
'set gxout fwrite'
'set fwrite D:\date\uv.dat'
'set grads off'
'set grid off'

*******************24日08时 21.3 121.9
'set t 1'
j=1
while(j<22)
'set z 'j
'set lat 21.3'
'set lon 121.9'
'd UGRDprs'
j=j+1
endwhile

j=1
while(j<22)
'set z 'j
'set lat 21.3'
'set lon 121.9'
'd VGRDprs'
j=j+1
endwhile
*******************24日14时 22.1 120.6
'set t 2'
j=1
while(j<22)
'set z 'j
'set lat 22.1'
'set lon 120.6'
'd UGRDprs'
j=j+1
endwhile

j=1
while(j<22)
'set z 'j
'set lat 22.1'
'set lon 120.6'
'd VGRDprs'
j=j+1
endwhile
*******************24日20时 23.1 119.2
'set t 3'
'set lat 23.1'
'set lon 119.2'
j=1
while(j<22)
'set z 'j
'set lat 23.1'
'set lon 119.2'
'd UGRDprs'
j=j+1
endwhile

j=1
while(j<22)
'set z 'j
'set lat 23.1'
'set lon 119.2'
'd VGRDprs'
j=j+1
endwhile

*******************25日02时 23.0 118.1
'set t 4'
j=1
while(j<22)
'set z 'j
'set lat 23.0'
'set lon 118.1'
'd UGRDprs'
j=j+1
endwhile

j=1
while(j<22)
'set z 'j
'set lat 23.0'
'set lon 118.1'
'd VGRDprs'
j=j+1
endwhile
*******************25日08时 23.7 117.4
'set t 5'
j=1
while(j<22)
'set z 'j
'set lat 23.7'
'set lon 117.4'
'd UGRDprs'
j=j+1
endwhile

j=1
while(j<22)
'set z 'j
'set lat 23.7'
'set lon 117.4'
'd VGRDprs'
j=j+1
endwhile
*******************25日14时 24.5 115.9
'set t 6'
j=1
while(j<22)
'set z 'j
'set lat 24.5'
'set lon 115.9'
'd UGRDprs'
j=j+1
endwhile

j=1
while(j<22)
'set z 'j
'set lat 24.5'
'set lon 115.9'
'd VGRDprs'
j=j+1
endwhile
'disable fwrite'
'reinit'
;



运行后,写的ctl，个人怀疑ctl没写对
DSET D:\date\uv.dat
*options template
index D:\date\uv.dat.idx
UNDEF -9999999.0
TITLE section
DTYPE grib2
ydef 1 linear -90.000000 1
xdef 1 linear 0.000000 1.000000
tdef 6 linear 00Z24aug2019 6hr
*  z has 21 levels -prs
zdef 21 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000
options pascals
vars 2
UGRDprs   21,100  0,2,2 ** (1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100) none U-Component of Wind [m/s]
VGRDprs   21,100  0,2,3 ** (1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100) none V-Component of Wind [m/s]
ENDVARS



然后画图时
set t 1 6
set z 1 21
d UGRDprs;VGRDprs
就显示Cannot contour grid - all undefined values  

求助各位大佬

nunu18

你的ctl描述的是纬向，纵向都只有一个格点，问题出在这儿估计

jlmjlm2007

nunu18 发表于 2020-4-6 16:11
你的ctl描述的是纬向，纵向都只有一个格点，问题出在这儿估计

我只要一个点的。要是就是不同点不同时次的气象要素。
看来得换写数据的思路。