gs文件中define的一个问题，求大神指教啊

小胖子

mofangbao 发表于 2013-3-28 08:57

                               
登录/注册后可看大图

说句题外话,为啥你这里
‘set z 1’
第一个引号是中文的?

额，这个是我复制进来的时候加上去的，因为数据只有一层，而且试了的，不管有没有这一句，结果都一样，所以我原来的gs文件里就没加这句

lizhiboys

mofangbao 发表于 2013-3-28 08:57

                               
登录/注册后可看大图

说句题外话,为啥你这里
‘set z 1’
第一个引号是中文的?

我觉得是打印错误  为什么是中文的呢？

river

小胖子 发表于 2013-3-27 23:11

                               
登录/注册后可看大图

又有另外一个问题，我按照版主所说的方法，计算大气净辐射通量，发现在青藏高原北部附近有明显的 ...

楼主看看是不是计算的时候加减号和正负号弄错了

小胖子

river 发表于 2013-3-28 11:13

                               
登录/注册后可看大图

楼主看看是不是计算的时候加减号和正负号弄错了

应该没有啊，我都检查过至少5遍了，还发现一个问题：我先用grads把nc文件转成grd文件，再编ctl文件画图，和直接编gs文件画出来的图不一样啊，中高纬一致，但到了低纬，用grd数据画出的图就发生严重错误，难道是我gs转格式编的不对？
nc转grd的gs文件如下：
'reinit'
'sdfopen e:\TP-Q1Q2\1981-2010mon.mean-fluxdata\1981.1-2010.12monshtfl-sfc.nc'
'set fwrite e:\TP-Q1Q2\1981-2010mon.mean-fluxdata\result\1981.1-2010.12monshtfl-sfc.grd'
'set gxout fwrite'
'set x 1 192'
'set y 1 94'
i=1
while(i<=360)
'set t 'i''
'd shtfl'
i=i+1
endwhile
'disable fwrite'
'reinit'
;

小胖子

问题困扰我很久了

小胖子

blui 发表于 2013-3-28 13:57

                               
登录/注册后可看大图

有没有计算大气热源的程序可以分享一下吗？

算的过程中遇到了问题

小胖子

river 发表于 2013-3-28 12:16

                               
登录/注册后可看大图

gs看着没什么问题，把set z加上试试，再看看是不是ctl里缺测值设置的有问题

试过set z了，还是一样；我用q ctlinfo把所有nc数据都查了一遍，缺测值都是-9.99e+33，我的ctl里面也是这么编的，实在是想不通原因啊，还是缺测值的处理另有玄妙？

小胖子

Aires 发表于 2013-3-28 12:22

                               
登录/注册后可看大图

对不起，我不是学气象的

，还是很谢谢Aires版主啊

river

小胖子 发表于 2013-3-28 11:29

                               
登录/注册后可看大图

应该没有啊，我都检查过至少5遍了，还发现一个问题：我先用grads把nc文件转成grd文件，再编ctl文件画图， ...

gs看着没什么问题，把set z加上试试，再看看是不是ctl里缺测值设置的有问题

Aires

小胖子 发表于 2013-3-27 23:11

                               
登录/注册后可看大图

又有另外一个问题，我按照版主所说的方法，计算大气净辐射通量，发现在青藏高原北部附近有明显的 ...

对不起，我不是学气象的