fortran如何读取wrf输出文件

wzf171

    在使用fortran读取wrf输出文件(.dat)时，想读取某时次的数据，总共33个时次，假如想读取第20个时次的rainc数据，fortran程序该怎么写？求大神指导。。。ctl文件如下：
dset f:/2004/wrfout-d3_040630.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.2.1 MODEL
pdef   66  48 lcc  21.640  111.922   33.500   24.500  60.00000  30.00000  111.92208   3000.000   3000.000
xdef  154 linear  110.87812   0.01351351
ydef  110 linear   20.89588   0.01351351
zdef    2 levels  
850.00000
500.00000
tdef   33 linear 12Z30JUN2004     720MN      
VARS   10
PH             2  0  perturbation geopotential (m2 s-2)
PHB            2  0  base-state geopotential (m2 s-2)
T              2  0  perturbation potential temperature (theta-t0) (K)
T2             1  0  TEMP at 2 M (K)
HGT            1  0  Terrain Height (m)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
height         2  0  Model height (KM)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.2.1 MODEL
@ global String comment SIMULATION_START_DATE = 2004-06-30_12:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION =    67
@ global String comment SOUTH-NORTH_GRID_DIMENSION =    49
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    28
@ global String comment MAP_PROJ =     1
@ global String comment DX =      3000.00
@ global String comment DY =      3000.00
@ global String comment CEN_LAT =        21.64
@ global String comment CEN_LON =       111.92
@ global String comment TRUELAT1 =        30.00
@ global String comment TRUELAT2 =        60.00
@ global String comment MOAD_CEN_LAT =        21.64
@ global String comment STAND_LON =       111.92
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =     1
@ global String comment RA_LW_PHYSICS =     1
@ global String comment RA_SW_PHYSICS =     1
@ global String comment SF_SFCLAY_PHYSICS =     2
@ global String comment SF_SURFACE_PHYSICS =     2
@ global String comment BL_PBL_PHYSICS =     2
@ global String comment CU_PHYSICS =     0
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment W_DAMPING =     0
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment GRID_ID =     3
@ global String comment PARENT_ID =     2
@ global String comment I_PARENT_START =    23
@ global String comment J_PARENT_START =    17
@ global String comment PARENT_GRID_RATIO =     3
@ global String comment DT =        15.00
@ global String comment ISWATER =    16
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14

vera

我也遇到这个问题。解决了吗？请赐教。

longlong

请赐教！谢谢

同路人..

工作量不大的话，可以用fwrite写出来，再用fortran读就好读了。

colorlingling

在此飘过，楼主加油

whw334455

很简单   先用grads的fwrite写出来    再用Fortran程序读写即可

zhimin0122

厉害啊。。。。。。。。。

vivian9624

期待好的答案，飘过