FNL变量的提取

白兔糖0217

先附一个具体的FNL的控制文件：
dset d:/taifu/2013/fnl_20130627_00_00_c                                      [给定二进制原始数据集的文件名（包括路径）]
index d:/taifu/2013/fnl_20130627_00_00_c.idx                                  (对应的索引文件)
undef 9.999E+20                                          （定义缺测值或者缺值值，Grads在运算操作和图形操作时将忽略这些值点）
title d:/taifu/2013/fnl_20130627_00_00_c                                      （标题名称）
*  produced by grib2ctl v0.9.12.5p39c                                        （生成来源）
dtype grib 3                                                                (设置背景网格参数为3)
options yrev                                                               （y维与YDEF定义相反方式存放）
ydef 181 linear -90.000000 1        （设置网格点值与纬度的对应关系，181是Y方向格点数，取LEVELS表格距均匀，-90为开始，间隔为1）
xdef 360 linear 0.000000 1.000000   （设置网格点值与经度的对应关系，360是X方向格点数，取LEVELS表格距均匀，0为开始，间隔为1）
tdef 1 linear 00Z27jun2013 1mo   （设置网格值与时间的映射关系，1是数据集中的时次数，00Z27jun2013表开始时间，1mo为时间的增量值）
*  z has 26 levels, for prs                                                    （表示有26层等压面）
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10      （26层高度对应的气压值）
vars 137                                                                     （vars表示变量描述开始， 137代表数据集中的变量数)
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]                                       【 近地表四层等压面的抬升指数】no5WAVA500mb  0 230,100,500 ** 500 mb 5-wave geopot. height anomaly [gpm]            【500 hPa 等压面位势高度距平】no5WAVH500mb  0 222,100,500 ** 500 mb 5-wave geopotential height [gpm]                   【500 hPa 等压面位势高度 】
ABSVprs 26 41,100,0 ** (profile) Absolute vorticity                                                      【绝对涡度】
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]                                          【地表对流降水量】
ALBDOsfc  0 84,1,0  ** surface Albedo [%]                                                              【地表反照率】
APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]                                                 【地表总降水量】
CAPEsfc  0 157,1,0  ** surface Convective Avail. Pot. Energy [J/kg]                                         【对流有效位能】
CAPE180_0mb  0 157,116,46080 ** 180-0 mb above gnd Convective Avail. Pot. Energy [J/kg]
CFRZRsfc  0 141,1,0  ** surface Categorical freezing rain [yes=1;no=0]
CICEPsfc  0 142,1,0  ** surface Categorical ice pellets [yes=1;no=0]
CINsfc  0 156,1,0  ** surface Convective inhibition [J/kg]
CIN180_0mb  0 156,116,46080 ** 180-0 mb above gnd Convective inhibition [J/kg]
CLWMRprs 21 153,100,0 ** (profile) Cloud water [kg/kg]              【云水含量】
CPRATsfc  0 214,1,0  ** surface Convective precip. rate [kg/m^2/s]
CRAINsfc  0 140,1,0  ** surface Categorical rain [yes=1;no=0]
CSNOWsfc  0 143,1,0  ** surface Categorical snow [yes=1;no=0]
CWATclm  0 76,200,0 ** atmos column Cloud water [kg/m^2]         【气柱云水含量】
CWORKclm  0 146,200,0 ** atmos column Cloud work function [J/kg]
DLWRFsfc  0 205,1,0  ** surface Downward long wave flux [W/m^2]
DSWRFsfc  0 204,1,0  ** surface Downward short wave flux [W/m^2]
GFLUXsfc  0 155,1,0  ** surface Ground heat flux [W/m^2]
GPAprs 2 27,100,0 ** (profile) Geopotential height anomaly [gpm]
HGTsfc  0 7,1,0  ** surface Geopotential height [gpm]
HGTprs 26 7,100,0 ** (profile) Geopotential height [gpm]
HGTpv2  0 7,117,2000 ** pot vorticity = 2000 units level Geopotential height [gpm]
HGTpvneg2  0 7,117,34768 ** pot vorticity = -2000 units level Geopotential height [gpm]
HGThtfl  0 7,204,0 ** highest trop freezing level Geopotential height [gpm]
HGT0deg  0 7,4,0 ** 0C isotherm level Geopotential height [gpm]
HGTmwl  0 7,6,0 ** max wind level Geopotential height [gpm]
HGTtrp  0 7,7,0 ** tropopause Geopotential height [gpm]
HPBLsfc  0 221,1,0  ** surface Planetary boundary layer height [m]
ICECsfc  0 91,1,0  ** surface Ice concentration (ice=1;no ice=0) [fraction]
LANDsfc  0 81,1,0  ** surface Land cover (land=1;sea=0) [fraction]
LFTXsfc  0 131,1,0  ** surface Surface lifted index [K]
LHTFLsfc  0 121,1,0  ** surface Latent heat flux [W/m^2]
O3MRprs 6 154,100,0 ** (profile) Ozone mixing ratio [kg/kg]
PEVPRsfc  0 145,1,0  ** surface Potential evaporation rate [W/m^2]
POTsig995   0 13,107,9950 ** sigma=.995  Potential temp. [K]
PRATEsfc  0 59,1,0  ** surface Precipitation rate [kg/m^2/s]
PRESsfc  0 1,1,0  ** surface Pressure [Pa]
PRESmsl  0 1,102,0  ** unknown level Pressure [Pa]
PRESpv2  0 1,117,2000 ** pot vorticity = 2000 units level Pressure [Pa]
PRESpvneg2  0 1,117,34768 ** pot vorticity = -2000 units level Pressure [Pa]
PRESlcb  0 1,212,0 ** low cloud base Pressure [Pa]
PRESlct  0 1,213,0 ** low cloud top Pressure [Pa]
PRESmcb  0 1,222,0 ** mid-cloud base Pressure [Pa]
PRESmct  0 1,223,0 ** mid-cloud top Pressure [Pa]
PREShcb  0 1,232,0 ** high cloud base Pressure [Pa]
PREShct  0 1,233,0 ** high cloud top Pressure [Pa]
PRESmwl  0 1,6,0 ** max wind level Pressure [Pa]
PREStrp  0 1,7,0 ** tropopause Pressure [Pa]
PRMSLmsl  0 2,102,0  ** unknown level Pressure reduced to MSL [Pa]
PWATclm  0 54,200,0 ** atmos column Precipitable water [kg/m^2]
RHprs 21 52,100,0 ** (profile) Relative humidity [%]                             【相对湿度】
RH2m  0 52,105,2 ** 2 m above ground Relative humidity [%]
RHsig995   0 52,107,9950 ** sigma=.995  Relative humidity [%]
RHsg33_100   0 52,108,8548 ** sigma=0.33-1 layer Relative humidity [%]
RHsg44_72   0 52,108,11336 ** sigma=0.44-0.72 layer Relative humidity [%]
RHsg44_100   0 52,108,11364 ** sigma=0.44-1 layer Relative humidity [%]
RHsg72_94   0 52,108,18526 ** sigma=0.72-0.94 layer Relative humidity [%]
RH30_0mb  0 52,116,7680 ** 30-0 mb above gnd Relative humidity [%]
RHclm  0 52,200,0 ** atmos column Relative humidity [%]
RHhtfl  0 52,204,0 ** highest trop freezing level Relative humidity [%]
RH0deg  0 52,4,0 ** 0C isotherm level Relative humidity [%]
SHTFLsfc  0 122,1,0  ** surface Sensible heat flux [W/m^2]
SOILW0_10cm  0 144,112,10 ** 0-10 cm underground Volumetric soil moisture [fraction]
SOILW10_40cm  0 144,112,2600 ** 10-40 cm underground Volumetric soil moisture [fraction]
SOILW40_100cm  0 144,112,10340 ** 40-100 cm underground Volumetric soil moisture [fraction]
SOILW100_200cm  0 144,112,25800 ** 100-200 cm underground Volumetric soil moisture [fraction]
SPFH2m  0 51,105,2 ** 2 m above ground Specific humidity [kg/kg]
SPFH30_0mb  0 51,116,7680 ** 30-0 mb above gnd Specific humidity [kg/kg]
TCDCclm  0 71,200,0 ** atmos column Total cloud cover [%]
TCDCbcl  0 71,211,0 ** boundary cld layer Total cloud cover [%]
TCDClcl  0 71,214,0 ** low cloud level Total cloud cover [%]
TCDCmcl  0 71,224,0 ** mid-cloud level Total cloud cover [%]
TCDChcl  0 71,234,0 ** high cloud level Total cloud cover [%]
TCDCcvl  0 71,244,0 ** convective cld layer Total cloud cover [%]
TMAX2m  0 15,105,2 ** 2 m above ground Max. temp. [K]
TMIN2m  0 16,105,2 ** 2 m above ground Min. temp. [K]
TMPsfc  0 11,1,0  ** surface Temp. [K]
TMPprs 26 11,100,0 ** (profile) Temp. [K]
TMP1524m  0 11,103,1524 ** 1524 m above msl Temp. [K]
TMP914m  0 11,103,914 ** 914 m above msl Temp. [K]
TMP457m  0 11,103,457 ** 457 m above msl Temp. [K]
TMP2m  0 11,105,2 ** 2 m above ground Temp. [K]
TMPsig995   0 11,107,9950 ** sigma=.995  Temp. [K]
TMP0_10cm  0 11,112,10 ** 0-10 cm underground Temp. [K]
TMP10_40cm  0 11,112,2600 ** 10-40 cm underground Temp. [K]
TMP40_100cm  0 11,112,10340 ** 40-100 cm underground Temp. [K]
TMP100_200cm  0 11,112,25800 ** 100-200 cm underground Temp. [K]
TMP30_0mb  0 11,116,7680 ** 30-0 mb above gnd Temp. [K]
TMPpv2  0 11,117,2000 ** pot vorticity = 2000 units level Temp. [K]
TMPpvneg2  0 11,117,34768 ** pot vorticity = -2000 units level Temp. [K]
TMPlct  0 11,213,0 ** low cloud top Temp. [K]
TMPmct  0 11,223,0 ** mid-cloud top Temp. [K]
TMPhct  0 11,233,0 ** high cloud top Temp. [K]
TMPmwl  0 11,6,0 ** max wind level Temp. [K]
TMPtrp  0 11,7,0 ** tropopause Temp. [K]
TOZNEclm  0 10,200,0 ** atmos column Total ozone [Dobson]
UGWDsfc  0 147,1,0  ** surface Zonal gravity wave stress [N/m^2]
UFLXsfc  0 124,1,0  ** surface Zonal momentum flux [N/m^2]
UGRDprs 26 33,100,0 ** (profile) u wind [m/s]
UGRD1524m  0 33,103,1524 ** 1524 m above msl u wind [m/s]
UGRD914m  0 33,103,914 ** 914 m above msl u wind [m/s]
UGRD457m  0 33,103,457 ** 457 m above msl u wind [m/s]
UGRD10m  0 33,105,10 ** 10 m above ground u wind [m/s]
UGRDsig995   0 33,107,9950 ** sigma=.995  u wind [m/s]
UGRD30_0mb  0 33,116,7680 ** 30-0 mb above gnd u wind [m/s]
UGRDpv2  0 33,117,2000 ** pot vorticity = 2000 units level u wind [m/s]
UGRDpvneg2  0 33,117,34768 ** pot vorticity = -2000 units level u wind [m/s]
UGRDmwl  0 33,6,0 ** max wind level u wind [m/s]
UGRDtrp  0 33,7,0 ** tropopause u wind [m/s]
ULWRFsfc  0 212,1,0  ** surface Upward long wave flux [W/m^2]
ULWRFtoa  0 212,8,0 ** top of atmos Upward long wave flux [W/m^2]
USWRFsfc  0 211,1,0  ** surface Upward short wave flux [W/m^2]
USWRFtoa  0 211,8,0 ** top of atmos Upward short wave flux [W/m^2]
VGWDsfc  0 148,1,0  ** surface Meridional gravity wave stress [N/m^2]
VFLXsfc  0 125,1,0  ** surface Meridional momentum flux [N/m^2]
VGRDprs 26 34,100,0 ** (profile) v wind [m/s]                                                 【v  分量】
VGRD1524m  0 34,103,1524 ** 1524 m above msl v wind [m/s]
VGRD914m  0 34,103,914 ** 914 m above msl v wind [m/s]
VGRD457m  0 34,103,457 ** 457 m above msl v wind [m/s]
VGRD10m  0 34,105,10 ** 10 m above ground v wind [m/s]
VGRDsig995   0 34,107,9950 ** sigma=.995  v wind [m/s]
VGRD30_0mb  0 34,116,7680 ** 30-0 mb above gnd v wind [m/s]
VGRDpv2  0 34,117,2000 ** pot vorticity = 2000 units level v wind [m/s]
VGRDpvneg2  0 34,117,34768 ** pot vorticity = -2000 units level v wind [m/s]
VGRDmwl  0 34,6,0 ** max wind level v wind [m/s]
VGRDtrp  0 34,7,0 ** tropopause v wind [m/s]
VVELprs 21 39,100,0 ** (profile) Pressure vertical velocity [Pa/s]                                                          【垂直速度】
VVELsig995   0 39,107,9950 ** sigma=.995  Pressure vertical velocity [Pa/s]
VWSHpv2  0 136,117,2000 ** pot vorticity = 2000 units level Vertical speed shear [1/s]
VWSHpvneg2  0 136,117,34768 ** pot vorticity = -2000 units level Vertical speed shear [1/s]
VWSHtrp  0 136,7,0 ** tropopause Vertical speed shear [1/s]                                                【垂直风切变】
WATRsfc  0 90,1,0  ** surface Water runoff [kg/m^2]
WEASDsfc  0 65,1,0  ** surface Accum. snow [kg/m^2]
ENDVARS





我想问问，如果现在想提取【云水含量】和【气柱云水含量】这两项，该从哪里着手？
修改了原来的GS文件和读取的FORTRAN程序，结果却都是0？
恳请各位指导指导

river

去看grads实用手册

传说中的谁

请问什么是“原来的GS文件和读取的FORTRAN程序”

白兔糖0217

传说中的谁 发表于 2013-12-4 21:05

                               
登录/注册后可看大图

请问什么是“原来的GS文件和读取的FORTRAN程序”

论坛里面有关于FNL数据提取的帖子，我就是根据那个改的gs.FORTRAN程序是用来转换二进制文件的，是别人给的。与拿来可以用的，现在我自己加了变量后，输出不了。

半个瞎子

请问楼主，你那有你说的“GS文件和读取的FORTRAN程序”么？我在论坛里没找到，你是否可以邮件给我呢？万分感谢。307494208@qq.com

buzaiyouyu

太感谢了！

buzaiyouyu

很明确。。

buzaiyouyu

好好研究下

hikingbaby

为什么我下载的时候就可以直接选择变量？那样下载下来就是var只有一个了吗？FNL的资料以前没用过。很烦恼啊

沐谌Linnea

楼主，我想问一下你说的哪个“原来的GS文件和读取的FORTRAN程序” 在哪里找到啊？能不能发给我看看  xiamochen@sina.cn