Slected variables from wrfout at "virtual tower" grid columns
The timeseries output does not include all of the variables available in the standard wrfout files. To supplement the timeseires output, selected variables are extracted for the virtual tower grid columns from one-minute frequency wrfout (d02) files
domains
- only the inner domain (d02)
- Frequency: 1 minutes
-
file location
/global/cfs/cdirs/m3314/ksa/simulations/{date}/{casename}/wrfout_tower
2D (Vertically integrated, surface, and soil variables)
file info
- Each file contains one-hour data, i.e., 60 samples for all the virtual towers
- File name format:
wrfout_{domain}_{casename}_vtowers_2D.all.{date-time period}.v[process version].nc
e.g.,
wrfout_d02_0518_03_vtowers_2D.all.2016-05-18_1400-1459.v1_0_0.nc
- Contains 7 x 7 grid columns centered at the virtual tower to help address sampling uncertainty and the model effective resolution.
dimensions:
soil_layers_stag = 4 ;
south_north = 7 ;
west_east = 7 ;
site = 89 ;
time = UNLIMITED ; // (60 currently)
nchars = 5 ;
variables:
char st_pfx(site, nchars) ;
st_pfx:_Encoding = "ascii" ;
st_pfx:longname = "name of virtual towers" ;
double time(time) ;
time:units = "seconds since 2016-05-18 00:00:00" ;
float ZS(soil_layers_stag) ;
ZS:units = "m" ;
ZS:longname = "DEPTHS OF CENTERS OF SOIL LAYERS" ;
float DZS(soil_layers_stag) ;
DZS:units = "m" ;
DZS:longname = "THICKNESSES OF SOIL LAYERS" ;
float XLAT(site, south_north, west_east) ;
XLAT:units = "degree_north" ;
XLAT:longname = "LATITUDE, SOUTH IS NEGATIVE" ;
float XLONG(site, south_north, west_east) ;
XLONG:units = "degree_east" ;
XLONG:longname = "LONGITUDE, WEST IS NEGATIVE" ;
float HGT(site, south_north, west_east) ;
HGT:units = "m" ;
HGT:longname = "Terrain Height" ;
int IVGTYP(site, south_north, west_east) ;
IVGTYP:units = "" ;
IVGTYP:longname = "DOMINANT VEGETATION CATEGORY" ;
int ISLTYP(site, south_north, west_east) ;
ISLTYP:units = "" ;
ISLTYP:longname = "DOMINANT SOIL CATEGORY" ;
float VEGFRA(site, south_north, west_east) ;
VEGFRA:units = "" ;
VEGFRA:longname = "VEGETATION FRACTION" ;
float Q2(time, south_north, west_east, site) ;
Q2:units = "kg kg-1" ;
Q2:longname = "QV at 2 M" ;
float T2(time, south_north, west_east, site) ;
T2:units = "K" ;
T2:longname = "TEMP at 2 M" ;
float TH2(time, south_north, west_east, site) ;
TH2:units = "K" ;
TH2:longname = "POT TEMP at 2 M" ;
float PSFC(time, south_north, west_east, site) ;
PSFC:units = "Pa" ;
PSFC:longname = "SFC PRESSURE" ;
float U10(time, south_north, west_east, site) ;
U10:units = "m s-1" ;
U10:longname = "U at 10 M" ;
float V10(time, south_north, west_east, site) ;
V10:units = "m s-1" ;
V10:longname = "V at 10 M" ;
float GRDFLX(time, south_north, west_east, site) ;
GRDFLX:units = "W m-2" ;
GRDFLX:longname = "GROUND HEAT FLUX" ;
float TSK(time, south_north, west_east, site) ;
TSK:units = "K" ;
TSK:longname = "SURFACE SKIN TEMPERATURE" ;
float RAINC(time, south_north, west_east, site) ;
RAINC:units = "mm" ;
RAINC:longname = "ACCUMULATED TOTAL CUMULUS PRECIPITATION" ;
float RAINSH(time, south_north, west_east, site) ;
RAINSH:units = "mm" ;
RAINSH:longname = "ACCUMULATED SHALLOW CUMULUS PRECIPITATION" ;
float RAINNC(time, south_north, west_east, site) ;
RAINNC:units = "mm" ;
RAINNC:longname = "ACCUMULATED TOTAL GRID SCALE PRECIPITATION" ;
float SNOWNC(time, south_north, west_east, site) ;
SNOWNC:units = "mm" ;
SNOWNC:longname = "ACCUMULATED TOTAL GRID SCALE SNOW AND ICE" ;
float GRAUPELNC(time, south_north, west_east, site) ;
GRAUPELNC:units = "mm" ;
GRAUPELNC:longname = "ACCUMULATED TOTAL GRID SCALE GRAUPEL" ;
float HAILNC(time, south_north, west_east, site) ;
HAILNC:units = "mm" ;
HAILNC:longname = "ACCUMULATED TOTAL GRID SCALE HAIL" ;
float SWDOWN(time, south_north, west_east, site) ;
SWDOWN:units = "W m-2" ;
SWDOWN:longname = "DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE" ;
float GLW(time, south_north, west_east, site) ;
GLW:units = "W m-2" ;
GLW:longname = "DOWNWARD LONG WAVE FLUX AT GROUND SURFACE" ;
float SWUPT(time, south_north, west_east, site) ;
SWUPT:units = "W m-2" ;
SWUPT:longname = "INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT TOP" ;
float SWUPTC(time, south_north, west_east, site) ;
SWUPTC:units = "W m-2" ;
SWUPTC:longname = "INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP" ;
float SWDNT(time, south_north, west_east, site) ;
SWDNT:units = "W m-2" ;
SWDNT:longname = "INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT TOP" ;
float SWDNTC(time, south_north, west_east, site) ;
SWDNTC:units = "W m-2" ;
SWDNTC:longname = "INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP" ;
float SWUPB(time, south_north, west_east, site) ;
SWUPB:units = "W m-2" ;
SWUPB:longname = "INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT BOTTOM" ;
float SWUPBC(time, south_north, west_east, site) ;
SWUPBC:units = "W m-2" ;
SWUPBC:longname = "INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM" ;
float SWDNB(time, south_north, west_east, site) ;
SWDNB:units = "W m-2" ;
SWDNB:longname = "INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT BOTTOM" ;
float SWDNBC(time, south_north, west_east, site) ;
SWDNBC:units = "W m-2" ;
SWDNBC:longname = "INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM" ;
float LWUPT(time, south_north, west_east, site) ;
LWUPT:units = "W m-2" ;
LWUPT:longname = "INSTANTANEOUS UPWELLING LONGWAVE FLUX AT TOP" ;
float LWUPTC(time, south_north, west_east, site) ;
LWUPTC:units = "W m-2" ;
LWUPTC:longname = "INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP" ;
float LWDNT(time, south_north, west_east, site) ;
LWDNT:units = "W m-2" ;
LWDNT:longname = "INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT TOP" ;
float LWDNTC(time, south_north, west_east, site) ;
LWDNTC:units = "W m-2" ;
LWDNTC:longname = "INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP" ;
float LWUPB(time, south_north, west_east, site) ;
LWUPB:units = "W m-2" ;
LWUPB:longname = "INSTANTANEOUS UPWELLING LONGWAVE FLUX AT BOTTOM" ;
float LWUPBC(time, south_north, west_east, site) ;
LWUPBC:units = "W m-2" ;
LWUPBC:longname = "INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM" ;
float LWDNB(time, south_north, west_east, site) ;
LWDNB:units = "W m-2" ;
LWDNB:longname = "INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT BOTTOM" ;
float LWDNBC(time, south_north, west_east, site) ;
LWDNBC:units = "W m-2" ;
LWDNBC:longname = "INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM" ;
float OLR(time, south_north, west_east, site) ;
OLR:units = "W m-2" ;
OLR:longname = "TOA OUTGOING LONG WAVE" ;
float ALBEDO(time, south_north, west_east, site) ;
ALBEDO:units = "-" ;
ALBEDO:longname = "ALBEDO" ;
float UST(time, south_north, west_east, site) ;
UST:units = "m s-1" ;
UST:longname = "U* IN SIMILARITY THEORY" ;
float HFX(time, south_north, west_east, site) ;
HFX:units = "W m-2" ;
HFX:longname = "UPWARD HEAT FLUX AT THE SURFACE" ;
float QFX(time, south_north, west_east, site) ;
QFX:units = "kg m-2 s-1" ;
QFX:longname = "UPWARD MOISTURE FLUX AT THE SURFACE" ;
float LH(time, south_north, west_east, site) ;
LH:units = "W m-2" ;
LH:longname = "LATENT HEAT FLUX AT THE SURFACE" ;
float CLDFRAC2D(time, south_north, west_east, site) ;
CLDFRAC2D:units = "%" ;
CLDFRAC2D:longname = "2-D MAX CLOUD FRACTION" ;
float WVP(time, south_north, west_east, site) ;
WVP:units = "kg m-2" ;
WVP:longname = "WATER VAPOR PATH" ;
float LWP(time, south_north, west_east, site) ;
LWP:units = "kg m-2" ;
LWP:longname = "LIQUID CLOUD WATER PATH" ;
float IWP(time, south_north, west_east, site) ;
IWP:units = "kg m-2" ;
IWP:longname = "ICE CLOUD WATER PATH" ;
float SWP(time, south_north, west_east, site) ;
SWP:units = "kg m-2" ;
SWP:longname = "SNOW CLOUD WATER PATH" ;
float WP_SUM(time, south_north, west_east, site) ;
WP_SUM:units = "kg m-2" ;
WP_SUM:longname = "SUM OF LWP+IWP+SWP" ;
float RE_QC(time, south_north, west_east, site) ;
RE_QC:units = "m" ;
RE_QC:longname = "MASS-WEIGHTED LIQUID CLOUD EFFECTIVE RADIUS" ;
float RE_QI(time, south_north, west_east, site) ;
RE_QI:units = "m" ;
RE_QI:longname = "MASS-WEIGHTED ICE EFFECTIVE RADIUS" ;
float RE_QS(time, south_north, west_east, site) ;
RE_QS:units = "m" ;
RE_QS:longname = "MASS-WEIGHTED SNOW EFFECTIVE RADIUS" ;
float TAU_QC(time, south_north, west_east, site) ;
TAU_QC:units = "" ;
TAU_QC:longname = "MASS-WEIGHTED LIQUID CLOUD OPTICAL THICKNESS" ;
float TAU_QI(time, south_north, west_east, site) ;
TAU_QI:units = "" ;
TAU_QI:longname = "MASS-WEIGHTED ICE OPTICAL THICKNESS" ;
float TAU_QS(time, south_north, west_east, site) ;
TAU_QS:units = "" ;
TAU_QS:longname = "MASS-WEIGHTED SNOW OPTICAL THICKNESS" ;
float CBASEHT(time, south_north, west_east, site) ;
CBASEHT:units = "m agl" ;
CBASEHT:longname = "CLOUD BASE HEIGHT" ;
float CTOPHT(time, south_north, west_east, site) ;
CTOPHT:units = "m agl" ;
CTOPHT:longname = "CLOUD TOP HEIGHT" ;
float TSLB(time, soil_layers_stag, south_north, west_east, site) ;
TSLB:units = "K" ;
TSLB:longname = "SOIL TEMPERATURE" ;
float SMOIS(time, soil_layers_stag, south_north, west_east, site) ;
SMOIS:units = "m3 m-3" ;
SMOIS:longname = "SOIL MOISTURE" ;
3D (height-dependent atmospheric variables)
still being processed (May 18 case done)
file info
- Each file contains one-hour data, i.e., 60 samples for one virtual tower
- File name format:
wrfout_{domain}_{casename}_vtowers_3D.[site name].{date-time period}.v1_0_0.nc
e.g.,
wrfout_d02_0518_03_vtowers_3D.VT032.2016-05-18_1600-1659.v1_0_0.nc
- Staggered variables are destaggered to the mass points
- Contains 7 x 7 grid columns centered at the virtual tower to help address sampling uncertainty and the model effective resolution.
dimensions:
bottom_top = 304 ;
bottom_top_stag = 305 ;
south_north = 7 ;
west_east = 7 ;
time = UNLIMITED ; // (60 currently)
variables:
double time(time) ;
time:units = "seconds since 2016-05-18 00:00:00" ;
float ZNU(bottom_top) ;
ZNU:units = "" ;
ZNU:longname = "eta values on half (mass) levels" ;
float ZNW(bottom_top_stag) ;
ZNW:units = "" ;
ZNW:longname = "eta values on full (w) levels" ;
float XLAT(south_north, west_east) ;
XLAT:units = "degree_north" ;
XLAT:longname = "LATITUDE, SOUTH IS NEGATIVE" ;
float XLONG(south_north, west_east) ;
XLONG:units = "degree_east" ;
XLONG:longname = "LONGITUDE, WEST IS NEGATIVE" ;
float HGT(south_north, west_east) ;
HGT:units = "m" ;
HGT:longname = "Terrain Height" ;
int IVGTYP(south_north, west_east) ;
IVGTYP:units = "" ;
IVGTYP:longname = "DOMINANT VEGETATION CATEGORY" ;
int ISLTYP(south_north, west_east) ;
ISLTYP:units = "" ;
ISLTYP:longname = "DOMINANT SOIL CATEGORY" ;
float VEGFRA(south_north, west_east) ;
VEGFRA:units = "" ;
VEGFRA:longname = "VEGETATION FRACTION" ;
float U(time, bottom_top, south_north, west_east) ;
U:units = "m s-1" ;
U:longname = "x-wind component" ;
float V(time, bottom_top, south_north, west_east) ;
V:units = "m s-1" ;
V:longname = "y-wind component" ;
float W(time, bottom_top, south_north, west_east) ;
W:units = "m s-1" ;
W:longname = "z-wind component" ;
float PH(time, bottom_top, south_north, west_east) ;
PH:units = "m2 s-2" ;
PH:longname = "perturbation geopotential" ;
float T(time, bottom_top, south_north, west_east) ;
T:units = "K" ;
T:longname = "perturbation potential temperature theta-t0" ;
float THM(time, bottom_top, south_north, west_east) ;
THM:units = "K" ;
THM:longname = "either 1) pert moist pot temp=(1+Rv/Rd Qv)*(theta)-T0, or 2) pert dry pot temp=t" ;
float P(time, bottom_top, south_north, west_east) ;
P:units = "Pa" ;
P:longname = "perturbation pressure" ;
float P_HYD(time, bottom_top, south_north, west_east) ;
P_HYD:units = "Pa" ;
P_HYD:longname = "hydrostatic pressure" ;
float CLDFRA(time, bottom_top, south_north, west_east) ;
CLDFRA:units = "" ;
CLDFRA:longname = "CLOUD FRACTION" ;
float QVAPOR(time, bottom_top, south_north, west_east) ;
QVAPOR:units = "kg kg-1" ;
QVAPOR:longname = "Water vapor mixing ratio" ;
float QCLOUD(time, bottom_top, south_north, west_east) ;
QCLOUD:units = "kg kg-1" ;
QCLOUD:longname = "Cloud water mixing ratio" ;
float QRAIN(time, bottom_top, south_north, west_east) ;
QRAIN:units = "kg kg-1" ;
QRAIN:longname = "Rain water mixing ratio" ;
float QICE(time, bottom_top, south_north, west_east) ;
QICE:units = "kg kg-1" ;
QICE:longname = "Ice mixing ratio" ;
float QSNOW(time, bottom_top, south_north, west_east) ;
QSNOW:units = "kg kg-1" ;
QSNOW:longname = "Snow mixing ratio" ;
float QGRAUP(time, bottom_top, south_north, west_east) ;
QGRAUP:units = "kg kg-1" ;
QGRAUP:longname = "Graupel mixing ratio" ;
float QNICE(time, bottom_top, south_north, west_east) ;
QNICE:units = " kg-1" ;
QNICE:longname = "Ice Number concentration" ;
float QNRAIN(time, bottom_top, south_north, west_east) ;
QNRAIN:units = " kg(-1)" ;
QNRAIN:longname = "Rain Number concentration" ;
float QNDROP(time, bottom_top, south_north, west_east) ;
QNDROP:units = " kg-1" ;
QNDROP:longname = "Droplet number mixing ratio" ;
float H_DIABATIC(time, bottom_top, south_north, west_east) ;
H_DIABATIC:units = "K s-1" ;
H_DIABATIC:longname = "MICROPHYSICS LATENT HEATING" ;
float QV_DIABATIC(time, bottom_top, south_north, west_east) ;
QV_DIABATIC:units = "g g-1 s-1" ;
QV_DIABATIC:longname = "MICROPHYSICS QV TENDENCY" ;
float QC_DIABATIC(time, bottom_top, south_north, west_east) ;
QC_DIABATIC:units = "g g-1 s-1" ;
QC_DIABATIC:longname = "MICROPHYSICS QC TENDENCY" ;
float m11(time, bottom_top, south_north, west_east) ;
m11:units = "m2 s-2" ;
m11:longname = "11 component of NBA subgrid stress tensor" ;
float m22(time, bottom_top, south_north, west_east) ;
m22:units = "m2 s-2" ;
m22:longname = "22 component of NBA subgrid stress tensor" ;
float m33(time, bottom_top, south_north, west_east) ;
m33:units = "m2 s-2" ;
m33:longname = "33 component of NBA subgrid stress tensor" ;
float m12(time, bottom_top, south_north, west_east) ;
m12:units = "m2 s-2" ;
m12:longname = "12 component of NBA subgrid stress tensor" ;
float m13(time, bottom_top, south_north, west_east) ;
m13:units = "m2 s-2" ;
m13:longname = "13 component of NBA subgrid stress tensor" ;
float m23(time, bottom_top, south_north, west_east) ;
m23:units = "m2 s-2" ;
m23:longname = "23 component of NBA subgrid stress tensor" ;