用fwrite转多变量的nc资料时出现的问题

朕小通

我想把3个nc资料里的变量转化放到一个grd里面，如下是我的gs程序：
'reinit'
'set gxout fwrite'
'set fwrite d:\yjh\data_required\grd\1.grd'
'sdfopen d:\yjh\data_required\6\cl.chemsteply.1.nc'
'sdfopen d:\yjh\data_required\6\O3.chemsteply.1.nc'
'sdfopen d:\yjh\data_required\6\T.chemsteply.1.nc'
'set lon 75 105'
'set lat 28 40'
ti=1
while(ti<=47)
'set t ' ti
zi=1
while(zi<=23)
'set z ' zi
'd cl'
zi=zi+1
endwhile
zi=1
while(zi<=23)
'set z ' zi
'd o3.2'
zi=zi+1
endwhile
zi=1
while(zi<=23)
'set z ' zi
'd t.3'
zi=zi+1
endwhile
ti=ti+1
endwhile
'disable fwrite'
'reinit'
;

ctl文件如下：
dset  D:\yjh\data_required\grd\1.grd
title WACCM3
undef 9.96921e+36
xdef 7 linear 75 5
ydef 4 linear 28 4
zdef 23 levels 1000 925 850 775 700 600 500 400 300 250 200 150 100 70 50 30 20 10 7 5 3 2 1
tdef 47 linear 00:30Z01jun1995 30mn
vars 3
cl 23 -999 cl
o3 23 -999 o3
t 23 -999 t
endvars

用生成的grd作图时发现，物理量之间完全错位了，请教各位。



用grd画出来的图和原图不同，尤其是画t变量的时候，单位明显不对。（ps）

mofangbao

楼主能否描述一下3个nc之间的关系，维数是如何对应的

朕小通

mofangbao 发表于 2012-5-1 09:50

                               
登录/注册后可看大图

楼主能否描述一下3个nc之间的关系，维数是如何对应的

你好！3个nc是同一个模式转出来的资料，水平分辨率：5*4，格点数：73*46，共17层，47个时次，包括缺测值，都是完全相同的

river

朕小通 发表于 2012-5-1 10:13

                               
登录/注册后可看大图

你好！3个nc是同一个模式转出来的资料，水平分辨率：5*4，格点数：73*46，共17层，47个时次，包括缺测值， ...

一共十七层，那楼主的循环怎么是23呢，而且三个资料完全相同的话，根本不需要写那么多循环，时间外循环，层次内循环，然后直接d那三个变量就行了

nx花儿-LJQ

river 发表于 2013-4-18 09:23

                               
登录/注册后可看大图

一共十七层，那楼主的循环怎么是23呢，而且三个资料完全相同的话，根本不需要写那么多循环，时间外循环， ...

不能那样直接d三个变量 貌似出不来 我也在弄这个 还是搞不懂

安清静静

你的问题解决了么  我现在也遇上了这样的问题 想问问你

Pygmalion

我用FNL资料的时候也出现过类似问题，转成grd文件后画图会出现变量不对应，且二维面上变量值也会发生错位

沐谌Linnea

Pygmalion 发表于 2014-4-20 11:56
我用FNL资料的时候也出现过类似问题，转成grd文件后画图会出现变量不对应，且二维面上变量值也会发生错位

我想问一下，我用的也是fnl资料，想要提取两个变量（u 和v）,这样写fwrite 是不是不对 总是不能成功转成grd文件
'reinit'
'open h:\tc\200007\fnl_200007.ctl'
'open h:\tc\200008\fnl_200008.ctl'
'open h:\tc\200009\fnl_200009.ctl'
'open h:\tc\200010\fnl_200010.ctl'
'set fwrite h:\tc\filter\uvlpf.grd'
'set gxout fwrite'
'set lon 0 360'
'set lat -90 90'
'set t 1'
'set lev 1000'
'd ugrdprs.1;vgrdprs.1'
'd ugrdprs.2;vgrdprs.2'
'd ugrdprs.3;vgrdprs.3'
'd ugrdprs.4;vgrdprs.4'
'disable fwrite'

Pygmalion

沐谌Linnea 发表于 2014-4-21 14:46
我想问一下，我用的也是fnl资料，想要提取两个变量（u 和v）,这样写fwrite 是不是不对 总是不能成功转成g ...

我刚开始也是用的这个办法，画出图后就出现我在7楼描述的情况。。。后来我用wgrib2命令：
wgrib2 inputfilename | grep  "var" | grep "lev" | wgrib2 -i inputfilename -no_header -bin outputfilename.grd

把var替换成你需要的变量名，lev替换成你需要的层次（需要所有高度的话就写“mb”），出来以后再用cat命令把输出的grd文件整合。。。这是我目前找到的办法。。。