如何用fortran输出wrf土壤温度变量

在linux下我通常是用ncdump和ncview先看查看netcdf文件，知含搏确定要读取的变量及其维数后，再用fortran读，只要三条语句就可以搞定了（下面程序中的红色部分）老腔。

下面是一个简单的Fortran程序读取wrfout文件，搭祥程序名为read_netcdf.f90

! This program is for reading wrfout

implicit none

integer, parameter :: LLm=144, MMm=96, Ndas=25 !定义维数

integer ierr,ncid,varid,len_file

real h(LLm,MMm,Ndas) !变量名：地形高度

character*299 file_wrfout!文件名路径

include 'netcdf.inc'

file_wrfout='/data/zlzang/fly_wrfout/wrfout_d03_2010-06-03_00:00:00'

len_file=len_trim(file_wrfout)

1、Vtable.GFS是fnl数据的头文件，用于WPS读取fnl数据用的。

2、wrfout和fnl文件的变量芦缺不一样，wrfout中有很多变量是预报过程中新生成的。

3、WRF通过fnl获得初始场和边界场，然后经过转WRF程序运算得出预报，（WRF程序，说白了宴哗氏就是一大堆大气动力数学物理公式的化身）。

注：我从事WRF数晌散值模拟研究两年了，也发表了一些相关的论文，希望可以帮到你~

WRF模式垂直方向的层数是可以设置的，垂直方向采用地形跟随质量坐标（即值为从0-1的eta分层），输出数据为不同预报时效各层次的气象要素，包括气压P、位势高度H、温度T、纬向风祥返U、经向风V等等。

不知道你皮昌获取的WRF模式输出数据是模式直接谨握饥输出的数据，还是经过后处理软件ARWpost处理后的数据？ARWpost软件可以将模式输出的垂直eta分层数据通过interp插值到气压层或者高度层上（具体的气压层和高度层可以根据需要设置，如气压层设置为1000hPa、925hPa、850hPa、700hPa、500hPa、300hPa等标准气压层）。

因此，如果你获取的是WRF模式直接输出的数据，需要采用ARWpost软件通过interp进行垂直插值处理，得到你所需的垂直高度层或气压层的数据；如果你获取的是经过ARWpost软件处理后的数据，那么你需要问提供数据的人，插值时设置的是气压层还是高度层，具体的层数是如何设置的。

---