求助：关于grads 中fwrite和ctl问题

bearwoo

我用grads将模式结果中的PM52.5作了时间平均以后提取出来，用Meteoinfo可以将ctl文件打开并画图，可是grads中打不开ctl,报错如下：

1.  原始ctl：
dset ^aero_norad2016121300_11d
options sequential big_endian
title model output from grapes
undef -9.99E+33
xdef  300 linear  100.0000    0.1500
ydef  200 linear   20.0000    0.1500
zdef   33 linear 1 1
tdef   792 linear 00z13DEC2016    60mn
vars   8
rh    3 99 relative humidity
aero_norad_ext   30 99 ext
pm25  3 99 pm2.5
vis_vert   0 99 vertical visibility
vis  0 99 old visibility
aod  0 0 aod
pblh 0 0 pblh
fktm 15 99 diffusion coefficient
endvars


2.  重新写入变量的gs文件：
'reinit'
'set gxout fwrite'
'set fwrite E:/py/work1/aqi/case1/model/myj_mean_pm25.dat'
*打开数据
'open H:/model/output/myj/case1/aero_hour/aero.ctl'
'set z 1'
'define pmmean=ave(pm25,t=1,t=792)'
'd pmmean'
'disable fwrite'


3. 新的ctl文件：
dset ^myj_mean_pm25.dat
undef -9.99E+33
xdef  300 linear  100.0000    0.1500
ydef  200 linear   20.0000    0.1500
zdef   1 level 1
tdef   1 level 1
vars   1
pm25  0 99 pm2.5
endvars


对grads不是很精通，请教一下这种问题应该怎么解决？

轻轻巧巧

新的ctl文件里zdef 和tdef有错，应该写为类似zdef 1 linear 1 1，tdef 1 linear  00z13DEC2016    60mn

ppioo

很好，支持

bearwoo

轻轻巧巧 发表于 2017-11-7 15:41
新的ctl文件里zdef 和tdef有错，应该写为类似zdef 1 linear 1 1，tdef 1 linear  00z13DEC2016    60mn

问题解决啦，感谢感谢！！

bearwoo

轻轻巧巧 发表于 2017-11-7 15:41
新的ctl文件里zdef 和tdef有错，应该写为类似zdef 1 linear 1 1，tdef 1 linear  00z13DEC2016    60mn

你好，我在处理气象要素时又遇见同样的问题，我按照原ctl写的。仍然打不开，能帮我看看是怎么回事吗？
原ctl:

1. dset ^posvar_norad2016121300_11d
   
2. options sequential big_endian
   
3. title post output from grapes
   
4. undef 9.999E+20
   
5. xdef   300  linear   100.0000    0.1500
   
6. ydef   200  linear    20.0000    0.1500
   
7. zdef   31 levels
   
8. 1000.000000
   
9. 975.0000000
   
10. 950.0000000
    
11. 925.0000000
    
12. 900.0000000
    
13. 850.0000000
    
14. 800.0000000
    
15. 750.0000000
    
16. 700.0000000
    
17. 650.0000000
    
18. 600.0000000
    
19. 550.0000000
    
20. 500.0000000
    
21. 450.0000000
    
22. 400.0000000
    
23. 350.0000000
    
24. 300.0000000
    
25. 250.0000000
    
26. 200.0000000
    
27. 150.0000000
    
28. 100.0000000
    
29. 70.00000000
    
30. 50.00000000
    
31. 30.00000000
    
32. 20.00000000
    
33. 10.00000000
    
34. 7.000000000
    
35. 5.000000000
    
36. 3.000000000
    
37. 2.000000000
    
38. 1.000000000
    
39. tdef   792 linear 00z13DEC2016    60mn
    
40. vars 13
    
41. ps 0 0 surface pressure
    
42. psl 0 0 sea level pressure
    
43. ts 0 0 surface temperature
    
44. t2m 0 0 t at 2m
    
45. pblh 0 0 v at 10m
    
46. ust 0 0 v at 10m
    
47. mol 0 0 mol at 10m
    
48. xland 0 0 xland 1 or 2
    
49. br 0 0 rib at 10m
    
50. znt 0 0 roughness lenth
    
51. mavail 0 0 surface moisture availability(0~1)
    
52. hfx 0 0 surface heat flux
    
53. qfx 0 0 surface vapour flux
    
54. endvars

复制代码

gs:

1. 'reinit'
   
2. 'set gxout fwrite'
   
3. 'set fwrite E:/py/work1/case1/met/model/mrf/mrf_mean_met.dat'
   
4. 'open H:/model/output/mrf/case1/met_hour/post.ctl'
   
5. 'set z 1'
   
6. 'define v1=ave(t2m,t=1,t=792)'
   
7. 'define v2=ave(pblh,t=1,t=792)'
   
8. 'define v3=ave(ust,t=1,t=792)'
   
9. 'define v4=ave(mol,t=1,t=792)'
   
10. 'd v1'
    
11. 'd v2'
    
12. 'd v3'
    
13. 'd v4'
    
14. 'disable fwrite'

复制代码

重新写的描述性ctl:

1. dset ^mrf_mean_met.dat
   
2. undef 9.999E+20
   
3. xdef  300 linear  100.0000    0.1500
   
4. ydef  200 linear   20.0000    0.1500
   
5. zdef   1 levels 1000.000000
   
6. tdef   1 linear 00z13DEC2016    60mn
   
7. vars   4
   
8. t2m 0 0 t at 2m
   
9. pblh 0 0 v at 10m
   
10. ust 0 0 v at 10m
    
11. mol 0 0 mol at 10m
    
12. endvars

复制代码

轻轻巧巧

错误提示是什么？

bearwoo

轻轻巧巧 发表于 2017-11-9 19:51
错误提示是什么？

原来是路径问题，之前竟然没注意到。。。不好意思哈!