t639垂直速度场的CTL怎么写

acmilan

639中的垂直速度厂，怎么写ctl文件啊

做个霸气的木头

ctl怎么写是基础的说，你看看手册，对应你的数据情况，自己写，遇到问题了再来问~

edwardli

@echo dset D:\WPSP\T639\%e%.grib1>>T639.ctl
@echo index D:\WPSP\T639\%e%.grib1.idx>>T639.ctl
@echo undef 9.999E+20>>T639.ctl
@echo title %e%.grib1>>T639.ctl
@echo *  produced by grib2ctl v0.9.12.5p39c>>T639.ctl
@echo dtype grib 255>>T639.ctl
@echo options yrev>>T639.ctl
@echo ydef 91 linear 0.000000 1.000000>>T639.ctl
@echo xdef 181 linear 0.000000 1.000000>>T639.ctl
@echo tdef 41 linear 12Z%e:~6,2%JUL%e:~0,4% 6hr>>T639.ctl
@echo *  z has 9 levels, for prs>>T639.ctl
@echo zdef 9 levels>>T639.ctl
@echo 1000 925 850 700 600 500 400 300 200>>T639.ctl
@echo vars 23>>T639.ctl
@echo APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]>>T639.ctl
@echo DEPRprs 3 18,100,0 ** (profile) Dew point depression [K]>>T639.ctl
@echo EPOTprs 3 14,100,0 ** (profile) Pseudo-adiabatic pot. temp. [K]>>T639.ctl
@echo HGTprs 9 7,100,0 ** (profile) Geopotential height [gpm]>>T639.ctl
@echo KXmsl  0 133,102,0  ** unknown level K index [K]>>T639.ctl
@echo LRGHRprs 3 241,100,0 ** (profile) Large scale condensation heating [K/s]>>T639.ctl
@echo MCONVprs 3 135,100,0 ** (profile) Horizontal moisture divergence [kg/kg/s]>>T639.ctl
@echo PRESsfc  0 1,1,0  ** surface Pressure [Pa]>>T639.ctl
@echo PRMSLmsl  0 2,102,0  ** unknown level Pressure reduced to MSL [Pa]>>T639.ctl
@echo RELDprs 4 44,100,0 ** (profile) Relative divergence >>T639.ctl
@echo RELVprs 4 43,100,0 ** (profile) Relative vorticity >>T639.ctl
@echo RHprs 9 52,100,0 ** (profile) Relative humidity [%]>>T639.ctl
@echo RH2m  0 52,105,2 ** 2 m above ground Relative humidity [%]>>T639.ctl
@echo SPFHprs 9 51,100,0 ** (profile) Specific humidity [kg/kg]>>T639.ctl
@echo SXprs 3 134,100,0 ** (profile) Sweat index [K]>>T639.ctl
@echo TMPprs 9 11,100,0 ** (profile) Temp. [K]>>T639.ctl
@echo TMP2m  0 11,105,2 ** 2 m above ground Temp. [K]>>T639.ctl
@echo TTRADprs 3 216,100,0 ** (profile) Temp. tendency by all radiation [K/s]>>T639.ctl
@echo UGRDprs 9 33,100,0 ** (profile) u wind [m/s]>>T639.ctl
@echo UGRD10m  0 33,105,10 ** 10 m above ground u wind [m/s]>>T639.ctl
@echo VGRDprs 9 34,100,0 ** (profile) v wind [m/s]>>T639.ctl
@echo VGRD10m  0 34,105,10 ** 10 m above ground v wind [m/s]>>T639.ctl
@echo VVELprs 9 39,100,0 ** (profile) Pressure vertical velocity [Pa/s]>>T639.ctl
@echo ENDVARS>>T639.ctl



把前面的@echo 和后面的》T639.ctl去掉就行了  VVELprs就是垂直速度；
当然你也可以拷贝一个T639一个时次的改个名字，然后用grib2ctl -verf 2012060600.grib 来产生ctl文件

mofangbao

楼主下次发帖时，麻烦按照我给改的标题格式来写标题，这样才能让大家明白。

qxtjls

非常好，成功