FNL资料的变量说明

查了一顿也没找到fnl数据到底包含多少气象数据，这是用wgrib从fnl数据里面读出来的，自己整理的。希望哪位大侠给指导下，fnl数据。 1 PRES=Pressure 气压 2 HGT=Geopotential height [gpm] 位势高度 3 TMP=Temp.[K] 温度 4 VVEL=Pressure vertical velocit

-

 

 

查了一顿也没找到fnl数据到底包含多少气象数据，这是用wgrib从fnl数据里面读出来的，自己整理的。希望哪位大侠给指导下，fnl数据。

 

 1 PRES=Pressure 气压

  2 HGT=Geopotential height [gpm]  位势高度

  3 TMP=Temp. [K]  温度

  4 VVEL=Pressure vertical velocity [Pa/s]  气压垂直速度

  5 RH=Relative humidity [%]  相对湿度

  6 ABSV=Absolute vorticity [/s]  绝对涡度

  7 O3MR=Ozone mixing ratio [kg/kg] 臭氧混合比

  8 CLWMR=Cloud water [kg/kg]  云水

  9 5WAVH=5-wavegeopotential height [gpm]

 10 UGRD=u wind [m/s]   纬向风

 11 VGRD=v wind [m/s]   经向风

 12 SPFH=Specific humidity [kg/kg] 比湿

 13 PWAT=Precipitable water [kg/m^2] 降雨量

 14 VWSH=Vertical speed shear [1/s]  垂向风切变

 15 LFTX=Surface lifted index [K]  地面举升指数

 16 CAPE=Convective Avail. Pot. Energy [J/kg]  对流可用位能

 17 CIN=Convective inhibition [J/kg]  对流抑制能

 18 4LFTX=Best (4-layer) lifted index [K]

 19 PRMSL=Pressure reduced to MSL [Pa]

 20 POT=Potential temp. [K   位温

 21 TOZNE=Total ozone [Dobson]  总臭氧

 22 CWAT=Cloud water [kg/m^2]   云水

 23 SOILW=Volumetric soil moisture [fraction]

 24 WEASD=Accum. snow [kg/m^2] 雪深

 25 LAND=Land cover (land=1;sea=0) [fraction]土地标示

 26 ICEC=Ice concentration (ice=1;no ice=0) [fraction] 冰密集度

 27 HPBL=Planetary boundary layer height

 

­ 行星边界层高度

 28TCDC=Total cloud cover [%]   总云水

 29 GPA=Geopotential height anomaly [gpm] 位势高度异常

 30 5WAVA=5-wave geopot. height anomaly [gpm]

 

 

ｇｒａｄｓ处理多个ｃｔｌ

(2010-07-29 15:02:41)

标签：

杂谈 

 

先生成一个ctl

grib1:

C:\ProgramFiles\PCGrADS\win32>grib2ctl grib_file > grib_file.ctl

or

gs>!grib2ctl grib_file >grib_file.ctl

grib2:

 

For analyses:

 

$ g2ctl -O grib2_file>grib2_file.ctl

$ gribmap -O -i grib2_file.ctl

$ grads

Landscape mode? (no for portrait):

ga-> open grib2_file.ctl

 

For forecasts:

 

 

$ g2ctl grib2_file>grib2_file.ctl

$ gribmap -i grib2_file.ctl

$ grads

Landscape mode? (no forportrait):

ga->open grib2_file.ctl

 

 

用grads处理多个相同格式的数据时若单个单个处理非常麻烦，当文件非常多的时候是单个处理是不实际的。下面介绍一种方法;

 

   第一步，在这种情况下可以重新写一个ctl描述文件，其文件变量都和已知的ctl相同，若原来的n文件只是时间不同，那么新描述文件的时间维数是所有原文件的时间的和。同样，若其他维数不同时也用同样的方法处理。

   第二步，在第一行之后添加一行： options template 表示多个时间序列原始数据文件想用一个描述文件统一地描述。这些原数据的原文件名由dset定义的形势命名文件名。

   第三步，修改dset 的文件名。原路径不变，把文件名用％表示。其中：

   ％y2  代表两位数年

   ％y4  代表四位数年

   ％m1  代表一位或者两位数的月

   ％m2  代表两位数月（用0补齐1位数）

   ％mc  3个字符月份的缩写

   ％d1  1或2位天

   ％d2  两位天

   ％h1  1或者2位时

   ％h2  2位时

例如：

   原文件其中之一的文件名为gdas2006050812f00，且所有文件只有天和时的变化

   那么新描述文件的文件名为：gdas200605%d2%h2f00

 

   另外如果源文件里有index项的话，需要修改其idx的文件名，假设改成fnl.idx。并用在dos下用gribmap函数生成一个新的idx文件。gribmap-e-i fnl.ctl（加绝对路径）

openfnl.ctl就可以打开所有文件。