WRF3.6版本中生成的ctl文件中计算大气顶和地面辐射通量的方法？

bingxue

请问各位大神，WRF3.6版本中生成的ctl文件中计算大气顶和地面辐射通量的方法？是用向下的短波辐射通量和长波辐射通量之和，减去向上的两者之和吗？那大气顶是指lev=100?,地面是指lev=1000?本人刚接触模式，还不太了解，希望各位指导一下！
dset e:\wrf2\test.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.6.1 MODEL
pdef  149 149 lcc  40.200  115.600   75.000   75.000  40.20000  40.20000  115.60000   3000.000   3000.000
xdef  418 linear  112.77216   0.01351351
ydef  317 linear   38.04006   0.01351351
zdef   19 levels  
1000.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef    2 linear 00Z04NOV2014     180MN      
VARS  151
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
LAP_HGT        1  0  Laplacian of orography (m)
U             19  0  x-wind component (m s-1)
V             19  0  y-wind component (m s-1)
W             19  0  z-wind component (m s-1)
PH            19  0  perturbation geopotential (m2 s-2)
PHB           19  0  base-state geopotential (m2 s-2)
T             19  0  perturbation potential temperature (theta-t0) (K)
MU             1  0  perturbation dry air mass in column (Pa)
MUB            1  0  base state dry air mass in column (Pa)
NEST_POS       1  0  - (-)
P             19  0  perturbation pressure (Pa)
PB            19  0  BASE STATE PRESSURE (Pa)
P_HYD         19  0  hydrostatic pressure (Pa)
Q2             1  0  QV at 2 M (kg kg-1)
T2             1  0  TEMP at 2 M (K)
TH2            1  0  POT TEMP at 2 M (K)
PSFC           1  0  SFC PRESSURE (Pa)
U10            1  0  U at 10 M (m s-1)
V10            1  0  V at 10 M (m s-1)
QVAPOR        19  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        19  0  Cloud water mixing ratio (kg kg-1)
QRAIN         19  0  Rain water mixing ratio (kg kg-1)
QICE          19  0  Ice mixing ratio (kg kg-1)
QSNOW         19  0  Snow mixing ratio (kg kg-1)
QGRAUP        19  0  Graupel mixing ratio (kg kg-1)
QNICE         19  0  Ice Number concentration (  kg-1)
QNRAIN        19  0  Rain Number concentration (  kg(-1))
SHDMAX         1  0  ANNUAL MAX VEG FRACTION (-)
SHDMIN         1  0  ANNUAL MIN VEG FRACTION (-)
SNOALB         1  0  ANNUAL MAX SNOW ALBEDO IN FRACTION (-)
TSLB           4  0  SOIL TEMPERATURE (K)
SMOIS          4  0  SOIL MOISTURE (m3 m-3)
SH2O           4  0  SOIL LIQUID WATER (m3 m-3)
SMCREL         4  0  RELATIVE SOIL MOISTURE (-)
SEAICE         1  0  SEA ICE FLAG (-)
XICEM          1  0  SEA ICE FLAG (PREVIOUS STEP) (-)
SFROFF         1  0  SURFACE RUNOFF (mm)
UDROFF         1  0  UNDERGROUND RUNOFF (mm)
IVGTYP         1  0  DOMINANT VEGETATION CATEGORY (-)
ISLTYP         1  0  DOMINANT SOIL CATEGORY (-)
VEGFRA         1  0  VEGETATION FRACTION (-)
GRDFLX         1  0  GROUND HEAT FLUX (W m-2)
ACGRDFLX       1  0  ACCUMULATED GROUND HEAT FLUX (J m-2)
ACSNOM         1  0  ACCUMULATED MELTED SNOW (kg m-2)
SNOW           1  0  SNOW WATER EQUIVALENT (kg m-2)
SNOWH          1  0  PHYSICAL SNOW DEPTH (m)
CANWAT         1  0  CANOPY WATER (kg m-2)
SSTSK          1  0  SKIN SEA SURFACE TEMPERATURE (K)
COSZEN         1  0  COS of SOLAR ZENITH ANGLE (dimensionless)
LAI            1  0  LEAF AREA INDEX (m-2/m-2)
VAR            1  0  OROGRAPHIC VARIANCE (-)
TKE_PBL       19  0  TKE from PBL (m2 s-2)
EL_PBL        19  0  Length scale from PBL (m)
MAPFAC_M       1  0  Map scale factor on mass grid (-)
MAPFAC_MX      1  0  Map scale factor on mass grid, x direction (-)
MAPFAC_MY      1  0  Map scale factor on mass grid, y direction (-)
MF_VX_INV      1  0  Inverse map scale factor on v-grid, x direction (-)
F              1  0  Coriolis sine latitude term (s-1)
E              1  0  Coriolis cosine latitude term (s-1)
SINALPHA       1  0  Local sine of map rotation (-)
COSALPHA       1  0  Local cosine of map rotation (-)
HGT            1  0  Terrain Height (m)
TSK            1  0  SURFACE SKIN TEMPERATURE (K)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINSH         1  0  ACCUMULATED SHALLOW CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
SNOWNC         1  0  ACCUMULATED TOTAL GRID SCALE SNOW AND ICE (mm)
GRAUPELNC      1  0  ACCUMULATED TOTAL GRID SCALE GRAUPEL (mm)
HAILNC         1  0  ACCUMULATED TOTAL GRID SCALE HAIL (mm)
CLDFRA        19  0  CLOUD FRACTION (-)
SWDOWN         1  0  DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE (W m-2)
GLW            1  0  DOWNWARD LONG WAVE FLUX AT GROUND SURFACE (W m-2)
SWNORM         1  0  NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT) (W m-2)
ACSWUPT        1  0  ACCUMULATED UPWELLING SHORTWAVE FLUX AT TOP (J m-2)
ACSWUPTC       1  0  ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (J m-2)
ACSWDNT        1  0  ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT TOP (J m-2)
ACSWDNTC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (J m-2)
ACSWUPB        1  0  ACCUMULATED UPWELLING SHORTWAVE FLUX AT BOTTOM (J m-2)
ACSWUPBC       1  0  ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (J m-2)
ACSWDNB        1  0  ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT BOTTOM (J m-2)
ACSWDNBC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (J m-2)
ACLWUPT        1  0  ACCUMULATED UPWELLING LONGWAVE FLUX AT TOP (J m-2)
ACLWUPTC       1  0  ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP (J m-2)
ACLWDNT        1  0  ACCUMULATED DOWNWELLING LONGWAVE FLUX AT TOP (J m-2)
ACLWDNTC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP (J m-2)
ACLWUPB        1  0  ACCUMULATED UPWELLING LONGWAVE FLUX AT BOTTOM (J m-2)
ACLWUPBC       1  0  ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (J m-2)
ACLWDNB        1  0  ACCUMULATED DOWNWELLING LONGWAVE FLUX AT BOTTOM (J m-2)
ACLWDNBC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (J m-2)
I_ACSWUPT      1  0  BUCKET FOR UPWELLING SHORTWAVE FLUX AT TOP (J m-2)
I_ACSWUPTC     1  0  BUCKET FOR UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (J m-2)
I_ACSWDNT      1  0  BUCKET FOR DOWNWELLING SHORTWAVE FLUX AT TOP (J m-2)
I_ACSWDNTC     1  0  BUCKET FOR DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (J m-2)
I_ACSWUPB      1  0  BUCKET FOR UPWELLING SHORTWAVE FLUX AT BOTTOM (J m-2)
I_ACSWUPBC     1  0  BUCKET FOR UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (J m-2)
I_ACSWDNB      1  0  BUCKET FOR DOWNWELLING SHORTWAVE FLUX AT BOTTOM (J m-2)
I_ACSWDNBC     1  0  BUCKET FOR DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (J m-2)
I_ACLWUPT      1  0  BUCKET FOR UPWELLING LONGWAVE FLUX AT TOP (J m-2)
I_ACLWUPTC     1  0  BUCKET FOR UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP (J m-2)
I_ACLWDNT      1  0  BUCKET FOR DOWNWELLING LONGWAVE FLUX AT TOP (J m-2)
I_ACLWDNTC     1  0  BUCKET FOR DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP (J m-2)
I_ACLWUPB      1  0  BUCKET FOR UPWELLING LONGWAVE FLUX AT BOTTOM (J m-2)
I_ACLWUPBC     1  0  BUCKET FOR UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (J m-2)
I_ACLWDNB      1  0  BUCKET FOR DOWNWELLING LONGWAVE FLUX AT BOTTOM (J m-2)
I_ACLWDNBC     1  0  BUCKET FOR DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (J m-2)
SWUPT          1  0  INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT TOP (W m-2)
SWUPTC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (W m-2)
SWDNT          1  0  INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT TOP (W m-2)
SWDNTC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (W m-2)
SWUPB          1  0  INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT BOTTOM (W m-2)
SWUPBC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (W m-2)
SWDNB          1  0  INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT BOTTOM (W m-2)
SWDNBC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (W m-2)
LWUPT          1  0  INSTANTANEOUS UPWELLING LONGWAVE FLUX AT TOP (W m-2)
LWUPTC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP (W m-2)
LWDNT          1  0  INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT TOP (W m-2)
LWDNTC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP (W m-2)
LWUPB          1  0  INSTANTANEOUS UPWELLING LONGWAVE FLUX AT BOTTOM (W m-2)
LWUPBC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (W m-2)
LWDNB          1  0  INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT BOTTOM (W m-2)
LWDNBC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (W m-2)
OLR            1  0  TOA OUTGOING LONG WAVE (W m-2)
ALBEDO         1  0  ALBEDO (-)
CLAT           1  0  COMPUTATIONAL GRID LATITUDE, SOUTH IS NEGATIVE (degree_north)
ALBBCK         1  0  BACKGROUND ALBEDO (-)
EMISS          1  0  SURFACE EMISSIVITY (-)
NOAHRES        1  0  RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET (W m{-2})
TMN            1  0  SOIL TEMPERATURE AT LOWER BOUNDARY (K)
XLAND          1  0  LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
UST            1  0  U* IN SIMILARITY THEORY (m s-1)
PBLH           1  0  PBL HEIGHT (m)
HFX            1  0  UPWARD HEAT FLUX AT THE SURFACE (W m-2)
QFX            1  0  UPWARD MOISTURE FLUX AT THE SURFACE (kg m-2 s-1)
LH             1  0  LATENT HEAT FLUX AT THE SURFACE (W m-2)
ACHFX          1  0  ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE (J m-2)
ACLHF          1  0  ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE (J m-2)
SNOWC          1  0  FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER) (-)
SR             1  0  fraction of frozen precipitation (-)
AFWA_RAIN      1  0  AFWA Diagnostic: Rain fall (mm)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
LAKEMASK       1  0  LAKE MASK (1 FOR LAKE, 0 FOR NON-LAKE) (-)
SST            1  0  SEA SURFACE TEMPERATURE (K)
SST_INPUT      1  0  SEA SURFACE TEMPERATURE FROM WRFLOWINPUT FILE (K)
pressure      19  0  Model pressure (hPa)
height        19  0  Model height (km)
tk            19  0  Temperature (K)
tc            19  0  Temperature (C)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.6.1 MODEL
@ global String comment START_DATE = 2014-11-04_00:00:00
@ global String comment SIMULATION_START_DATE = 2014-11-04_00:00:00
@ global String comment WEST-EAST_GRID_DIMENSION =   150
@ global String comment SOUTH-NORTH_GRID_DIMENSION =   150
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    35
@ global String comment DX =      3000.00
@ global String comment DY =      3000.00
@ global String comment STOCH_FORCE_OPT =     0
@ global String comment GRIDTYPE = C
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment DAMPCOEF =         0.01
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =     8
@ global String comment RA_LW_PHYSICS =     4
@ global String comment RA_SW_PHYSICS =     4
@ 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 =     1
@ global String comment SF_LAKE_PHYSICS =     0
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment GFDDA_INTERVAL_M =     0
@ global String comment GFDDA_END_H =     0
@ global String comment GRID_SFDDA =     0
@ global String comment SGFDDA_INTERVAL_M =     0
@ global String comment SGFDDA_END_H =     0
@ global String comment HYPSOMETRIC_OPT =     2
@ global String comment SF_URBAN_PHYSICS =     1
@ global String comment SHCU_PHYSICS =     0
@ global String comment MFSHCONV =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment SWRAD_SCAT =         1.00
@ global String comment W_DAMPING =     0
@ global String comment DT =        18.00
@ global String comment RADT =         3.00
@ global String comment BLDT =         0.00
@ global String comment CUDT =         5.00
@ global String comment AER_OPT =     0
@ global String comment SWINT_OPT =     0
@ global String comment AER_TYPE =     1
@ global String comment AER_AOD550_OPT =     1
@ global String comment AER_ANGEXP_OPT =     1
@ global String comment AER_SSA_OPT =     1
@ global String comment AER_ASY_OPT =     1
@ global String comment AER_AOD550_VAL =         0.12
@ global String comment AER_ANGEXP_VAL =         1.30
@ global String comment AER_SSA_VAL =         0.00
@ global String comment AER_ASY_VAL =         0.00
@ global String comment MOIST_ADV_OPT =     1
@ global String comment SCALAR_ADV_OPT =     1
@ global String comment TKE_ADV_OPT =     1
@ global String comment DIFF_6TH_OPT =     0
@ global String comment DIFF_6TH_FACTOR =         0.12
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment BUCKET_MM =        -1.00
@ global String comment BUCKET_J =        -1.00
@ global String comment PREC_ACC_DT =         0.00
@ global String comment SF_OCEAN_PHYSICS =     0
@ global String comment ISFTCFLX =     0
@ global String comment ISHALLOW =     0
@ global String comment ISFFLX =     1
@ global String comment ICLOUD =     1
@ global String comment ICLOUD_CU =     0
@ global String comment TRACER_PBLMIX =     1
@ global String comment SCALAR_PBLMIX =     0
@ global String comment GRAV_SETTLING =     0
@ global String comment DFI_OPT =     0
@ global String comment WEST-EAST_PATCH_START_UNSTAG =     1
@ global String comment WEST-EAST_PATCH_END_UNSTAG =   149
@ global String comment WEST-EAST_PATCH_START_STAG =     1
@ global String comment WEST-EAST_PATCH_END_STAG =   150
@ global String comment SOUTH-NORTH_PATCH_START_UNSTAG =     1
@ global String comment SOUTH-NORTH_PATCH_END_UNSTAG =   149
@ global String comment SOUTH-NORTH_PATCH_START_STAG =     1
@ global String comment SOUTH-NORTH_PATCH_END_STAG =   150
@ global String comment BOTTOM-TOP_PATCH_START_UNSTAG =     1
@ global String comment BOTTOM-TOP_PATCH_END_UNSTAG =    34
@ global String comment BOTTOM-TOP_PATCH_START_STAG =     1
@ global String comment BOTTOM-TOP_PATCH_END_STAG =    35
@ global String comment GRID_ID =     1
@ global String comment PARENT_ID =     0
@ global String comment I_PARENT_START =     1
@ global String comment J_PARENT_START =     1
@ global String comment PARENT_GRID_RATIO =     1
@ global String comment CEN_LAT =        40.20
@ global String comment CEN_LON =       115.60
@ global String comment TRUELAT1 =        40.20
@ global String comment TRUELAT2 =        40.20
@ global String comment MOAD_CEN_LAT =        40.20
@ global String comment STAND_LON =       115.60
@ global String comment POLE_LAT =        90.00
@ global String comment POLE_LON =         0.00
@ global String comment GMT =         0.00
@ global String comment JULYR =  2014
@ global String comment JULDAY =   308
@ global String comment MAP_PROJ =     1
@ global String comment MAP_PROJ_CHAR = Lambert Conformal
@ global String comment MMINLU = USGS
@ global String comment NUM_LAND_CAT =    24
@ global String comment ISWATER =    16
@ global String comment ISLAKE =    -1
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14

WMQ

请问楼主问题解决了吗，想要请教一下