MeteoInfo脚本示例：GRIB to ARL

良辰 · 发表于 2015-6-27 10:15:07

良辰发表于 2015-6-27 09:50
王老师，刚才我查看了一下原始的nc数据，里面数据没有缺失，但是垂直高度上，U\V一致，W比UV多。如果我 ...

File Name: D:\Temp\wrfout.nc
File type: NetCDF-3/CDM (netCDF)
Dimensions: 9
Time = 85;
DateStrLen = 19;
west_east = 201;
south_north = 171;
bottom_top = 38;
bottom_top_stag = 39;
soil_layers_stag = 4;
west_east_stag = 202;
south_north_stag = 172;
X Dimension: Xmin = -902263.9756635734; Xmax = 897736.0243364266; Xsize = 201; Xdelta = 9000.0
Y Dimension: Ymin = -763839.0042780734; Ymax = 766160.9957219266; Ysize = 171; Ydelta = 9000.0
Global Attributes: [TITLE = " OUTPUT FROM WRF V3.3.1 MODEL", START_DATE = "2015-06-17_12:00:00", SIMULATION_START_DATE = "2015-06-17_12:00:00", WEST-EAST_GRID_DIMENSION = 202, SOUTH-NORTH_GRID_DIMENSION = 172, BOTTOM-TOP_GRID_DIMENSION = 39, DX = 9000.0f, DY = 9000.0f, GRIDTYPE = "C", DIFF_OPT = 1, KM_OPT = 4, DAMP_OPT = 1, DAMPCOEF = 0.09f, KHDIF = 0.0f, KVDIF = 0.0f, MP_PHYSICS = 4, RA_LW_PHYSICS = 4, RA_SW_PHYSICS = 4, SF_SFCLAY_PHYSICS = 5, SF_SURFACE_PHYSICS = 2, BL_PBL_PHYSICS = 5, CU_PHYSICS = 5, SURFACE_INPUT_SOURCE = 1, SST_UPDATE = 0, GRID_FDDA = 1, GFDDA_INTERVAL_M = 360, GFDDA_END_H = 60, GRID_SFDDA = 0, SGFDDA_INTERVAL_M = 0, SGFDDA_END_H = 0, SF_URBAN_PHYSICS = 0, FEEDBACK = 0, SMOOTH_OPTION = 2, SWRAD_SCAT = 1.0f, W_DAMPING = 1, MOIST_ADV_OPT = 1, SCALAR_ADV_OPT = 1, TKE_ADV_OPT = 1, DIFF_6TH_OPT = 2, DIFF_6TH_FACTOR = 0.06f, FGDT = 0.0f, GUV = 6.0E-4f, GT = 6.0E-4f, GQ = 6.0E-4f, IF_RAMPING = 0, DTRAMP_MIN = 60.0f, OBS_NUDGE_OPT = 0, BUCKET_MM = -1.0f, BUCKET_J = -1.0f, PREC_ACC_DT = 0.0f, OMLCALL = 0, ISFTCFLX = 0, ISHALLOW = 0, DFI_OPT = 0, SHCU_PHYSICS = 0, WEST-EAST_PATCH_START_UNSTAG = 1, WEST-EAST_PATCH_END_UNSTAG = 201, WEST-EAST_PATCH_START_STAG = 1, WEST-EAST_PATCH_END_STAG = 202, SOUTH-NORTH_PATCH_START_UNSTAG = 1, SOUTH-NORTH_PATCH_END_UNSTAG = 171, SOUTH-NORTH_PATCH_START_STAG = 1, SOUTH-NORTH_PATCH_END_STAG = 172, BOTTOM-TOP_PATCH_START_UNSTAG = 1, BOTTOM-TOP_PATCH_END_UNSTAG = 38, BOTTOM-TOP_PATCH_START_STAG = 1, BOTTOM-TOP_PATCH_END_STAG = 39, GRID_ID = 1, PARENT_ID = 1, I_PARENT_START = 1, J_PARENT_START = 1, PARENT_GRID_RATIO = 1, DT = 27.0f, CEN_LAT = 43.500004f, CEN_LON = 126.5f, TRUELAT1 = 30.0f, TRUELAT2 = 60.0f, MOAD_CEN_LAT = 43.500004f, STAND_LON = 126.5f, POLE_LAT = 90.0f, POLE_LON = 0.0f, GMT = 12.0f, JULYR = 2015, JULDAY = 168, MAP_PROJ = 1, MMINLU = "MODIFIED_IGBP_MODIS_NOAH", NUM_LAND_CAT = 20, ISWATER = 17, ISLAKE = -1, ISICE = 15, ISURBAN = 13, ISOILWATER = 14]
: TITLE = " OUTPUT FROM WRF V3.3.1 MODEL"
: START_DATE = "2015-06-17_12:00:00"
: SIMULATION_START_DATE = "2015-06-17_12:00:00"
: WEST-EAST_GRID_DIMENSION = 202
: SOUTH-NORTH_GRID_DIMENSION = 172
: BOTTOM-TOP_GRID_DIMENSION = 39
: DX = 9000.0f
: DY = 9000.0f
: GRIDTYPE = "C"
: DIFF_OPT = 1
: KM_OPT = 4
: DAMP_OPT = 1
: DAMPCOEF = 0.09f
: KHDIF = 0.0f
: KVDIF = 0.0f
: MP_PHYSICS = 4
: RA_LW_PHYSICS = 4
: RA_SW_PHYSICS = 4
: SF_SFCLAY_PHYSICS = 5
: SF_SURFACE_PHYSICS = 2
: BL_PBL_PHYSICS = 5
: CU_PHYSICS = 5
: SURFACE_INPUT_SOURCE = 1
: SST_UPDATE = 0
: GRID_FDDA = 1
: GFDDA_INTERVAL_M = 360
: GFDDA_END_H = 60
: GRID_SFDDA = 0
: SGFDDA_INTERVAL_M = 0
: SGFDDA_END_H = 0
: SF_URBAN_PHYSICS = 0
: FEEDBACK = 0
: SMOOTH_OPTION = 2
: SWRAD_SCAT = 1.0f
: W_DAMPING = 1
: MOIST_ADV_OPT = 1
: SCALAR_ADV_OPT = 1
: TKE_ADV_OPT = 1
: DIFF_6TH_OPT = 2
: DIFF_6TH_FACTOR = 0.06f
: FGDT = 0.0f
: GUV = 6.0E-4f
: GT = 6.0E-4f
: GQ = 6.0E-4f
: IF_RAMPING = 0
: DTRAMP_MIN = 60.0f
: OBS_NUDGE_OPT = 0
: BUCKET_MM = -1.0f
: BUCKET_J = -1.0f
: PREC_ACC_DT = 0.0f
: OMLCALL = 0
: ISFTCFLX = 0
: ISHALLOW = 0
: DFI_OPT = 0
: SHCU_PHYSICS = 0
: WEST-EAST_PATCH_START_UNSTAG = 1
: WEST-EAST_PATCH_END_UNSTAG = 201
: WEST-EAST_PATCH_START_STAG = 1
: WEST-EAST_PATCH_END_STAG = 202
: SOUTH-NORTH_PATCH_START_UNSTAG = 1
: SOUTH-NORTH_PATCH_END_UNSTAG = 171
: SOUTH-NORTH_PATCH_START_STAG = 1
: SOUTH-NORTH_PATCH_END_STAG = 172
: BOTTOM-TOP_PATCH_START_UNSTAG = 1
: BOTTOM-TOP_PATCH_END_UNSTAG = 38
: BOTTOM-TOP_PATCH_START_STAG = 1
: BOTTOM-TOP_PATCH_END_STAG = 39
: GRID_ID = 1
: PARENT_ID = 1
: I_PARENT_START = 1
: J_PARENT_START = 1
: PARENT_GRID_RATIO = 1
: DT = 27.0f
: CEN_LAT = 43.500004f
: CEN_LON = 126.5f
: TRUELAT1 = 30.0f
: TRUELAT2 = 60.0f
: MOAD_CEN_LAT = 43.500004f
: STAND_LON = 126.5f
: POLE_LAT = 90.0f
: POLE_LON = 0.0f
: GMT = 12.0f
: JULYR = 2015
: JULDAY = 168
: MAP_PROJ = 1
: MMINLU = "MODIFIED_IGBP_MODIS_NOAH"
: NUM_LAND_CAT = 20
: ISWATER = 17
: ISLAKE = -1
: ISICE = 15
: ISURBAN = 13
: ISOILWATER = 14
Variations: 205
char Times(Time,DateStrLen);
float LU_INDEX(Time,south_north,west_east);
LU_INDEX: FieldType = 104
LU_INDEX: MemoryOrder = "XY "
LU_INDEX: description = "LAND USE CATEGORY"
LU_INDEX: units = ""
LU_INDEX: stagger = ""
LU_INDEX: coordinates = "XLONG XLAT"
float ZNU(Time,bottom_top);
ZNU: FieldType = 104
ZNU: MemoryOrder = "Z "
ZNU: description = "eta values on half (mass) levels"
ZNU: units = ""
ZNU: stagger = ""
float ZNW(Time,bottom_top_stag);
ZNW: FieldType = 104
ZNW: MemoryOrder = "Z "
ZNW: description = "eta values on full (w) levels"
ZNW: units = ""
ZNW: stagger = "Z"
float ZS(Time,soil_layers_stag);
ZS: FieldType = 104
ZS: MemoryOrder = "Z "
ZS: description = "DEPTHS OF CENTERS OF SOIL LAYERS"
ZS: units = "m"
ZS: stagger = "Z"
float DZS(Time,soil_layers_stag);
DZS: FieldType = 104
DZS: MemoryOrder = "Z "
DZS: description = "THICKNESSES OF SOIL LAYERS"
DZS: units = "m"
DZS: stagger = "Z"
float U(Time,bottom_top,south_north,west_east_stag);
U: FieldType = 104
U: MemoryOrder = "XYZ"
U: description = "x-wind component"
U: units = "m s-1"
U: stagger = "X"
U: coordinates = "XLONG_U XLAT_U"
float V(Time,bottom_top,south_north_stag,west_east);
V: FieldType = 104
V: MemoryOrder = "XYZ"
V: description = "y-wind component"
V: units = "m s-1"
V: stagger = "Y"
V: coordinates = "XLONG_V XLAT_V"
float W(Time,bottom_top_stag,south_north,west_east);
W: FieldType = 104
W: MemoryOrder = "XYZ"
W: description = "z-wind component"
W: units = "m s-1"
W: stagger = "Z"
W: coordinates = "XLONG XLAT"
float PH(Time,bottom_top_stag,south_north,west_east);
PH: FieldType = 104
PH: MemoryOrder = "XYZ"
PH: description = "perturbation geopotential"
PH: units = "m2 s-2"
PH: stagger = "Z"
PH: coordinates = "XLONG XLAT"
float PHB(Time,bottom_top_stag,south_north,west_east);
PHB: FieldType = 104
PHB: MemoryOrder = "XYZ"
PHB: description = "base-state geopotential"
PHB: units = "m2 s-2"
PHB: stagger = "Z"
PHB: coordinates = "XLONG XLAT"
float T(Time,bottom_top,south_north,west_east);
T: FieldType = 104
T: MemoryOrder = "XYZ"
T: description = "perturbation potential temperature (theta-t0)"
T: units = "K"
T: stagger = ""
T: coordinates = "XLONG XLAT"
float HFX_FORCE(Time);
HFX_FORCE: FieldType = 104
HFX_FORCE: MemoryOrder = "0 "
HFX_FORCE: description = "SCM ideal surface sensible heat flux"
HFX_FORCE: units = "W m-2"
HFX_FORCE: stagger = ""
float LH_FORCE(Time);
LH_FORCE: FieldType = 104
LH_FORCE: MemoryOrder = "0 "
LH_FORCE: description = "SCM ideal surface latent heat flux"
LH_FORCE: units = "W m-2"
LH_FORCE: stagger = ""
float TSK_FORCE(Time);
TSK_FORCE: FieldType = 104
TSK_FORCE: MemoryOrder = "0 "
TSK_FORCE: description = "SCM ideal surface skin temperature"
TSK_FORCE: units = "W m-2"
TSK_FORCE: stagger = ""
float HFX_FORCE_TEND(Time);
HFX_FORCE_TEND: FieldType = 104
HFX_FORCE_TEND: MemoryOrder = "0 "
HFX_FORCE_TEND: description = "SCM ideal surface sensible heat flux tendency"
HFX_FORCE_TEND: units = "W m-2 s-1"
HFX_FORCE_TEND: stagger = ""
float LH_FORCE_TEND(Time);
LH_FORCE_TEND: FieldType = 104
LH_FORCE_TEND: MemoryOrder = "0 "
LH_FORCE_TEND: description = "SCM ideal surface latent heat flux tendency"
LH_FORCE_TEND: units = "W m-2 s-1"
LH_FORCE_TEND: stagger = ""
float TSK_FORCE_TEND(Time);
TSK_FORCE_TEND: FieldType = 104
TSK_FORCE_TEND: MemoryOrder = "0 "
TSK_FORCE_TEND: description = "SCM ideal surface skin temperature tendency"
TSK_FORCE_TEND: units = "W m-2 s-1"
TSK_FORCE_TEND: stagger = ""
float MU(Time,south_north,west_east);
MU: FieldType = 104
MU: MemoryOrder = "XY "
MU: description = "perturbation dry air mass in column"
MU: units = "Pa"
MU: stagger = ""
MU: coordinates = "XLONG XLAT"
float MUB(Time,south_north,west_east);
MUB: FieldType = 104
MUB: MemoryOrder = "XY "
MUB: description = "base state dry air mass in column"
MUB: units = "Pa"
MUB: stagger = ""
MUB: coordinates = "XLONG XLAT"
float NEST_POS(Time,south_north,west_east);
NEST_POS: FieldType = 104
NEST_POS: MemoryOrder = "XY "
NEST_POS: description = "-"
NEST_POS: units = "-"
NEST_POS: stagger = ""
NEST_POS: coordinates = "XLONG XLAT"
float P(Time,bottom_top,south_north,west_east);
P: FieldType = 104
P: MemoryOrder = "XYZ"
P: description = "perturbation pressure"
P: units = "Pa"
P: stagger = ""
P: coordinates = "XLONG XLAT"
float PB(Time,bottom_top,south_north,west_east);
PB: FieldType = 104
PB: MemoryOrder = "XYZ"
PB: description = "BASE STATE PRESSURE"
PB: units = "Pa"
PB: stagger = ""
PB: coordinates = "XLONG XLAT"
float FNM(Time,bottom_top);
FNM: FieldType = 104
FNM: MemoryOrder = "Z "
FNM: description = "upper weight for vertical stretching"
FNM: units = ""
FNM: stagger = ""
float FNP(Time,bottom_top);
FNP: FieldType = 104
FNP: MemoryOrder = "Z "
FNP: description = "lower weight for vertical stretching"
FNP: units = ""
FNP: stagger = ""
float RDNW(Time,bottom_top);
RDNW: FieldType = 104
RDNW: MemoryOrder = "Z "
RDNW: description = "inverse d(eta) values between full (w) levels"
RDNW: units = ""
RDNW: stagger = ""
float RDN(Time,bottom_top);
RDN: FieldType = 104
RDN: MemoryOrder = "Z "
RDN: description = "inverse d(eta) values between half (mass) levels"
RDN: units = ""
RDN: stagger = ""
float DNW(Time,bottom_top);
DNW: FieldType = 104
DNW: MemoryOrder = "Z "
DNW: description = "d(eta) values between full (w) levels"
DNW: units = ""
DNW: stagger = ""
float DN(Time,bottom_top);
DN: FieldType = 104
DN: MemoryOrder = "Z "
DN: description = "d(eta) values between half (mass) levels"
DN: units = ""
DN: stagger = ""
float CFN(Time);
CFN: FieldType = 104
CFN: MemoryOrder = "0 "
CFN: description = "extrapolation constant"
CFN: units = ""
CFN: stagger = ""
float CFN1(Time);
CFN1: FieldType = 104
CFN1: MemoryOrder = "0 "
CFN1: description = "extrapolation constant"
CFN1: units = ""
CFN1: stagger = ""
float P_HYD(Time,bottom_top,south_north,west_east);
P_HYD: FieldType = 104
P_HYD: MemoryOrder = "XYZ"
P_HYD: description = "hydrostatic pressure"
P_HYD: units = "Pa"
P_HYD: stagger = ""
P_HYD: coordinates = "XLONG XLAT"
float Q2(Time,south_north,west_east);
Q2: FieldType = 104
Q2: MemoryOrder = "XY "
Q2: description = "QV at 2 M"
Q2: units = "kg kg-1"
Q2: stagger = ""
Q2: coordinates = "XLONG XLAT"
float T2(Time,south_north,west_east);
T2: FieldType = 104
T2: MemoryOrder = "XY "
T2: description = "TEMP at 2 M"
T2: units = "K"
T2: stagger = ""
T2: coordinates = "XLONG XLAT"
float TH2(Time,south_north,west_east);
TH2: FieldType = 104
TH2: MemoryOrder = "XY "
TH2: description = "POT TEMP at 2 M"
TH2: units = "K"
TH2: stagger = ""
TH2: coordinates = "XLONG XLAT"
float PSFC(Time,south_north,west_east);
PSFC: FieldType = 104
PSFC: MemoryOrder = "XY "
PSFC: description = "SFC PRESSURE"
PSFC: units = "Pa"
PSFC: stagger = ""
PSFC: coordinates = "XLONG XLAT"
float U10(Time,south_north,west_east);
U10: FieldType = 104
U10: MemoryOrder = "XY "
U10: description = "U at 10 M"
U10: units = "m s-1"
U10: stagger = ""
U10: coordinates = "XLONG XLAT"
float V10(Time,south_north,west_east);
V10: FieldType = 104
V10: MemoryOrder = "XY "
V10: description = "V at 10 M"
V10: units = "m s-1"
V10: stagger = ""
V10: coordinates = "XLONG XLAT"
float RDX(Time);
RDX: FieldType = 104
RDX: MemoryOrder = "0 "
RDX: description = "INVERSE X GRID LENGTH"
RDX: units = ""
RDX: stagger = ""
float RDY(Time);
RDY: FieldType = 104
RDY: MemoryOrder = "0 "
RDY: description = "INVERSE Y GRID LENGTH"
RDY: units = ""
RDY: stagger = ""
float RESM(Time);
RESM: FieldType = 104
RESM: MemoryOrder = "0 "
RESM: description = "TIME WEIGHT CONSTANT FOR SMALL STEPS"
RESM: units = ""
RESM: stagger = ""
float ZETATOP(Time);
ZETATOP: FieldType = 104
ZETATOP: MemoryOrder = "0 "
ZETATOP: description = "ZETA AT MODEL TOP"
ZETATOP: units = ""
ZETATOP: stagger = ""
float CF1(Time);
CF1: FieldType = 104
CF1: MemoryOrder = "0 "
CF1: description = "2nd order extrapolation constant"
CF1: units = ""
CF1: stagger = ""
float CF2(Time);
CF2: FieldType = 104
CF2: MemoryOrder = "0 "
CF2: description = "2nd order extrapolation constant"
CF2: units = ""
CF2: stagger = ""
float CF3(Time);
CF3: FieldType = 104
CF3: MemoryOrder = "0 "
CF3: description = "2nd order extrapolation constant"
CF3: units = ""
CF3: stagger = ""
int ITIMESTEP(Time);
ITIMESTEP: FieldType = 106
ITIMESTEP: MemoryOrder = "0 "
ITIMESTEP: description = ""
ITIMESTEP: units = ""
ITIMESTEP: stagger = ""
float XTIME(Time);
XTIME: FieldType = 104
XTIME: MemoryOrder = "0 "
XTIME: description = "minutes since simulation start"
XTIME: units = ""
XTIME: stagger = ""
float QVAPOR(Time,bottom_top,south_north,west_east);
QVAPOR: FieldType = 104
QVAPOR: MemoryOrder = "XYZ"
QVAPOR: description = "Water vapor mixing ratio"
QVAPOR: units = "kg kg-1"
QVAPOR: stagger = ""
QVAPOR: coordinates = "XLONG XLAT"
float QCLOUD(Time,bottom_top,south_north,west_east);
QCLOUD: FieldType = 104
QCLOUD: MemoryOrder = "XYZ"
QCLOUD: description = "Cloud water mixing ratio"
QCLOUD: units = "kg kg-1"
QCLOUD: stagger = ""
QCLOUD: coordinates = "XLONG XLAT"
float QRAIN(Time,bottom_top,south_north,west_east);
QRAIN: FieldType = 104
QRAIN: MemoryOrder = "XYZ"
QRAIN: description = "Rain water mixing ratio"
QRAIN: units = "kg kg-1"
QRAIN: stagger = ""
QRAIN: coordinates = "XLONG XLAT"
float QICE(Time,bottom_top,south_north,west_east);
QICE: FieldType = 104
QICE: MemoryOrder = "XYZ"
QICE: description = "Ice mixing ratio"
QICE: units = "kg kg-1"
QICE: stagger = ""
QICE: coordinates = "XLONG XLAT"
float QSNOW(Time,bottom_top,south_north,west_east);
QSNOW: FieldType = 104
QSNOW: MemoryOrder = "XYZ"
QSNOW: description = "Snow mixing ratio"
QSNOW: units = "kg kg-1"
QSNOW: stagger = ""
QSNOW: coordinates = "XLONG XLAT"
float qke(Time,bottom_top,south_north,west_east);
qke: FieldType = 104
qke: MemoryOrder = "XYZ"
qke: description = "twice TKE FROM MYNN"
qke: units = "m2 s-2"
qke: stagger = ""
qke: coordinates = "XLONG XLAT"
float LANDMASK(Time,south_north,west_east);
LANDMASK: FieldType = 104
LANDMASK: MemoryOrder = "XY "
LANDMASK: description = "LAND MASK (1 FOR LAND, 0 FOR WATER)"
LANDMASK: units = ""
LANDMASK: stagger = ""
LANDMASK: coordinates = "XLONG XLAT"
float TSLB(Time,soil_layers_stag,south_north,west_east);
TSLB: FieldType = 104
TSLB: MemoryOrder = "XYZ"
TSLB: description = "SOIL TEMPERATURE"
TSLB: units = "K"
TSLB: stagger = "Z"
TSLB: coordinates = "XLONG XLAT"
float SMOIS(Time,soil_layers_stag,south_north,west_east);
SMOIS: FieldType = 104
SMOIS: MemoryOrder = "XYZ"
SMOIS: description = "SOIL MOISTURE"
SMOIS: units = "m3 m-3"
SMOIS: stagger = "Z"
SMOIS: coordinates = "XLONG XLAT"
float SH2O(Time,soil_layers_stag,south_north,west_east);
SH2O: FieldType = 104
SH2O: MemoryOrder = "XYZ"
SH2O: description = "SOIL LIQUID WATER"
SH2O: units = "m3 m-3"
SH2O: stagger = "Z"
SH2O: coordinates = "XLONG XLAT"
float SMCREL(Time,soil_layers_stag,south_north,west_east);
SMCREL: FieldType = 104
SMCREL: MemoryOrder = "XYZ"
SMCREL: description = "RELATIVE SOIL MOISTURE"
SMCREL: units = ""
SMCREL: stagger = "Z"
SMCREL: coordinates = "XLONG XLAT"
float SEAICE(Time,south_north,west_east);
SEAICE: FieldType = 104
SEAICE: MemoryOrder = "XY "
SEAICE: description = "SEA ICE FLAG"
SEAICE: units = ""
SEAICE: stagger = ""
SEAICE: coordinates = "XLONG XLAT"
float XICEM(Time,south_north,west_east);
XICEM: FieldType = 104
XICEM: MemoryOrder = "XY "
XICEM: description = "SEA ICE FLAG (PREVIOUS STEP)"
XICEM: units = ""
XICEM: stagger = ""
XICEM: coordinates = "XLONG XLAT"
float SFROFF(Time,south_north,west_east);
SFROFF: FieldType = 104
SFROFF: MemoryOrder = "XY "
SFROFF: description = "SURFACE RUNOFF"
SFROFF: units = "mm"
SFROFF: stagger = ""
SFROFF: coordinates = "XLONG XLAT"
float UDROFF(Time,south_north,west_east);
UDROFF: FieldType = 104
UDROFF: MemoryOrder = "XY "
UDROFF: description = "UNDERGROUND RUNOFF"
UDROFF: units = "mm"
UDROFF: stagger = ""
UDROFF: coordinates = "XLONG XLAT"
int IVGTYP(Time,south_north,west_east);
IVGTYP: FieldType = 106
IVGTYP: MemoryOrder = "XY "
IVGTYP: description = "DOMINANT VEGETATION CATEGORY"
IVGTYP: units = ""
IVGTYP: stagger = ""
IVGTYP: coordinates = "XLONG XLAT"
int ISLTYP(Time,south_north,west_east);
ISLTYP: FieldType = 106
ISLTYP: MemoryOrder = "XY "
ISLTYP: description = "DOMINANT SOIL CATEGORY"
ISLTYP: units = ""
ISLTYP: stagger = ""
ISLTYP: coordinates = "XLONG XLAT"
float VEGFRA(Time,south_north,west_east);
VEGFRA: FieldType = 104
VEGFRA: MemoryOrder = "XY "
VEGFRA: description = "VEGETATION FRACTION"
VEGFRA: units = ""
VEGFRA: stagger = ""
VEGFRA: coordinates = "XLONG XLAT"
float GRDFLX(Time,south_north,west_east);
GRDFLX: FieldType = 104
GRDFLX: MemoryOrder = "XY "
GRDFLX: description = "GROUND HEAT FLUX"
GRDFLX: units = "W m-2"
GRDFLX: stagger = ""
GRDFLX: coordinates = "XLONG XLAT"
float ACGRDFLX(Time,south_north,west_east);
ACGRDFLX: FieldType = 104
ACGRDFLX: MemoryOrder = "XY "
ACGRDFLX: description = "ACCUMULATED GROUND HEAT FLUX"
ACGRDFLX: units = "J m-2"
ACGRDFLX: stagger = ""
ACGRDFLX: coordinates = "XLONG XLAT"
float SNOW(Time,south_north,west_east);
SNOW: FieldType = 104
SNOW: MemoryOrder = "XY "
SNOW: description = "SNOW WATER EQUIVALENT"
SNOW: units = "kg m-2"
SNOW: stagger = ""
SNOW: coordinates = "XLONG XLAT"
float SNOWH(Time,south_north,west_east);
SNOWH: FieldType = 104
SNOWH: MemoryOrder = "XY "
SNOWH: description = "PHYSICAL SNOW DEPTH"
SNOWH: units = "m"
SNOWH: stagger = ""
SNOWH: coordinates = "XLONG XLAT"
float CANWAT(Time,south_north,west_east);
CANWAT: FieldType = 104
CANWAT: MemoryOrder = "XY "
CANWAT: description = "CANOPY WATER"
CANWAT: units = "kg m-2"
CANWAT: stagger = ""
CANWAT: coordinates = "XLONG XLAT"
float SST(Time,south_north,west_east);
SST: FieldType = 104
SST: MemoryOrder = "XY "
SST: description = "SEA SURFACE TEMPERATURE"
SST: units = "K"
SST: stagger = ""
SST: coordinates = "XLONG XLAT"
float SSTSK(Time,south_north,west_east);
SSTSK: FieldType = 104
SSTSK: MemoryOrder = "XY "
SSTSK: description = "SKIN SEA SURFACE TEMPERATURE"
SSTSK: units = "K"
SSTSK: stagger = ""
SSTSK: coordinates = "XLONG XLAT"
float LAI(Time,south_north,west_east);
LAI: FieldType = 104
LAI: MemoryOrder = "XY "
LAI: description = "Leaf area index"
LAI: units = "area/area"
LAI: stagger = ""
LAI: coordinates = "XLONG XLAT"
float DTAUX3D(Time,bottom_top,south_north,west_east);
DTAUX3D: FieldType = 104
DTAUX3D: MemoryOrder = "XYZ"
DTAUX3D: description = "LOCAL U GWDO STRESS"
DTAUX3D: units = "m s-1"
DTAUX3D: stagger = ""
DTAUX3D: coordinates = "XLONG XLAT"
float DTAUY3D(Time,bottom_top,south_north,west_east);
DTAUY3D: FieldType = 104
DTAUY3D: MemoryOrder = "XYZ"
DTAUY3D: description = "LOCAL V GWDO STRESS"
DTAUY3D: units = "m s-1"
DTAUY3D: stagger = ""
DTAUY3D: coordinates = "XLONG XLAT"
float DUSFCG(Time,south_north,west_east);
DUSFCG: FieldType = 104
DUSFCG: MemoryOrder = "XY "
DUSFCG: description = "COLUMN-INTEGRATED U GWDO STRESS"
DUSFCG: units = "Pa m s-1"
DUSFCG: stagger = ""
DUSFCG: coordinates = "XLONG XLAT"
float DVSFCG(Time,south_north,west_east);
DVSFCG: FieldType = 104
DVSFCG: MemoryOrder = "XY "
DVSFCG: description = "COLUMN-INTEGRATED V GWDO STRESS"
DVSFCG: units = "Pa m s-1"
DVSFCG: stagger = ""
DVSFCG: coordinates = "XLONG XLAT"
float VAR(Time,south_north,west_east);
VAR: FieldType = 104
VAR: MemoryOrder = "XY "
VAR: description = "OROGRAPHIC VARIANCE"
VAR: units = ""
VAR: stagger = ""
VAR: coordinates = "XLONG XLAT"
float CON(Time,south_north,west_east);
CON: FieldType = 104
CON: MemoryOrder = "XY "
CON: description = "OROGRAPHIC CONVEXITY"
CON: units = ""
CON: stagger = ""
CON: coordinates = "XLONG XLAT"
float OA1(Time,south_north,west_east);
OA1: FieldType = 104
OA1: MemoryOrder = "XY "
OA1: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OA1: units = ""
OA1: stagger = ""
OA1: coordinates = "XLONG XLAT"
float OA2(Time,south_north,west_east);
OA2: FieldType = 104
OA2: MemoryOrder = "XY "
OA2: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OA2: units = ""
OA2: stagger = ""
OA2: coordinates = "XLONG XLAT"
float OA3(Time,south_north,west_east);
OA3: FieldType = 104
OA3: MemoryOrder = "XY "
OA3: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OA3: units = ""
OA3: stagger = ""
OA3: coordinates = "XLONG XLAT"
float OA4(Time,south_north,west_east);
OA4: FieldType = 104
OA4: MemoryOrder = "XY "
OA4: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OA4: units = ""
OA4: stagger = ""
OA4: coordinates = "XLONG XLAT"
float OL1(Time,south_north,west_east);
OL1: FieldType = 104
OL1: MemoryOrder = "XY "
OL1: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OL1: units = ""
OL1: stagger = ""
OL1: coordinates = "XLONG XLAT"
float OL2(Time,south_north,west_east);
OL2: FieldType = 104
OL2: MemoryOrder = "XY "
OL2: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OL2: units = ""
OL2: stagger = ""
OL2: coordinates = "XLONG XLAT"
float OL3(Time,south_north,west_east);
OL3: FieldType = 104
OL3: MemoryOrder = "XY "
OL3: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OL3: units = ""
OL3: stagger = ""
OL3: coordinates = "XLONG XLAT"
float OL4(Time,south_north,west_east);
OL4: FieldType = 104
OL4: MemoryOrder = "XY "
OL4: description = "OROGRAPHIC DIRECTION ASYMMETRY FUNCTION"
OL4: units = ""
OL4: stagger = ""
OL4: coordinates = "XLONG XLAT"
float H_DIABATIC(Time,bottom_top,south_north,west_east);
H_DIABATIC: FieldType = 104
H_DIABATIC: MemoryOrder = "XYZ"
H_DIABATIC: description = "MICROPHYSICS LATENT HEATING"
H_DIABATIC: units = "K s-1"
H_DIABATIC: stagger = ""
H_DIABATIC: coordinates = "XLONG XLAT"
float MAPFAC_M(Time,south_north,west_east);
MAPFAC_M: FieldType = 104
MAPFAC_M: MemoryOrder = "XY "
MAPFAC_M: description = "Map scale factor on mass grid"
MAPFAC_M: units = ""
MAPFAC_M: stagger = ""
MAPFAC_M: coordinates = "XLONG XLAT"
float MAPFAC_U(Time,south_north,west_east_stag);
MAPFAC_U: FieldType = 104
MAPFAC_U: MemoryOrder = "XY "
MAPFAC_U: description = "Map scale factor on u-grid"
MAPFAC_U: units = ""
MAPFAC_U: stagger = "X"
MAPFAC_U: coordinates = "XLONG_U XLAT_U"
float MAPFAC_V(Time,south_north_stag,west_east);
MAPFAC_V: FieldType = 104
MAPFAC_V: MemoryOrder = "XY "
MAPFAC_V: description = "Map scale factor on v-grid"
MAPFAC_V: units = ""
MAPFAC_V: stagger = "Y"
MAPFAC_V: coordinates = "XLONG_V XLAT_V"
float MAPFAC_MX(Time,south_north,west_east);
MAPFAC_MX: FieldType = 104
MAPFAC_MX: MemoryOrder = "XY "
MAPFAC_MX: description = "Map scale factor on mass grid, x direction"
MAPFAC_MX: units = ""
MAPFAC_MX: stagger = ""
MAPFAC_MX: coordinates = "XLONG XLAT"
float MAPFAC_MY(Time,south_north,west_east);
MAPFAC_MY: FieldType = 104
MAPFAC_MY: MemoryOrder = "XY "
MAPFAC_MY: description = "Map scale factor on mass grid, y direction"
MAPFAC_MY: units = ""
MAPFAC_MY: stagger = ""
MAPFAC_MY: coordinates = "XLONG XLAT"
float MAPFAC_UX(Time,south_north,west_east_stag);
MAPFAC_UX: FieldType = 104
MAPFAC_UX: MemoryOrder = "XY "
MAPFAC_UX: description = "Map scale factor on u-grid, x direction"
MAPFAC_UX: units = ""
MAPFAC_UX: stagger = "X"
MAPFAC_UX: coordinates = "XLONG_U XLAT_U"
float MAPFAC_UY(Time,south_north,west_east_stag);
MAPFAC_UY: FieldType = 104
MAPFAC_UY: MemoryOrder = "XY "
MAPFAC_UY: description = "Map scale factor on u-grid, y direction"
MAPFAC_UY: units = ""
MAPFAC_UY: stagger = "X"
MAPFAC_UY: coordinates = "XLONG_U XLAT_U"
float MAPFAC_VX(Time,south_north_stag,west_east);
MAPFAC_VX: FieldType = 104
MAPFAC_VX: MemoryOrder = "XY "
MAPFAC_VX: description = "Map scale factor on v-grid, x direction"
MAPFAC_VX: units = ""
MAPFAC_VX: stagger = "Y"
MAPFAC_VX: coordinates = "XLONG_V XLAT_V"
float MF_VX_INV(Time,south_north_stag,west_east);
MF_VX_INV: FieldType = 104
MF_VX_INV: MemoryOrder = "XY "
MF_VX_INV: description = "Inverse map scale factor on v-grid, x direction"
MF_VX_INV: units = ""
MF_VX_INV: stagger = "Y"
MF_VX_INV: coordinates = "XLONG_V XLAT_V"
float MAPFAC_VY(Time,south_north_stag,west_east);
MAPFAC_VY: FieldType = 104
MAPFAC_VY: MemoryOrder = "XY "
MAPFAC_VY: description = "Map scale factor on v-grid, y direction"
MAPFAC_VY: units = ""
MAPFAC_VY: stagger = "Y"
MAPFAC_VY: coordinates = "XLONG_V XLAT_V"
float F(Time,south_north,west_east);
F: FieldType = 104
F: MemoryOrder = "XY "
F: description = "Coriolis sine latitude term"
F: units = "s-1"
F: stagger = ""
F: coordinates = "XLONG XLAT"
float E(Time,south_north,west_east);
E: FieldType = 104
E: MemoryOrder = "XY "
E: description = "Coriolis cosine latitude term"
E: units = "s-1"
E: stagger = ""
E: coordinates = "XLONG XLAT"
float SINALPHA(Time,south_north,west_east);
SINALPHA: FieldType = 104
SINALPHA: MemoryOrder = "XY "
SINALPHA: description = "Local sine of map rotation"
SINALPHA: units = ""
SINALPHA: stagger = ""
SINALPHA: coordinates = "XLONG XLAT"
float COSALPHA(Time,south_north,west_east);
COSALPHA: FieldType = 104
COSALPHA: MemoryOrder = "XY "
COSALPHA: description = "Local cosine of map rotation"
COSALPHA: units = ""
COSALPHA: stagger = ""
COSALPHA: coordinates = "XLONG XLAT"
float HGT(Time,south_north,west_east);
HGT: FieldType = 104
HGT: MemoryOrder = "XY "
HGT: description = "Terrain Height"
HGT: units = "m"
HGT: stagger = ""
HGT: coordinates = "XLONG XLAT"
float TSK(Time,south_north,west_east);
TSK: FieldType = 104
TSK: MemoryOrder = "XY "
TSK: description = "SURFACE SKIN TEMPERATURE"
TSK: units = "K"
TSK: stagger = ""
TSK: coordinates = "XLONG XLAT"
float P_TOP(Time);
P_TOP: FieldType = 104
P_TOP: MemoryOrder = "0 "
P_TOP: description = "PRESSURE TOP OF THE MODEL"
P_TOP: units = "Pa"
P_TOP: stagger = ""
float T00(Time);
T00: FieldType = 104
T00: MemoryOrder = "0 "
T00: description = "BASE STATE TEMPERATURE"
T00: units = "K"
T00: stagger = ""
float P00(Time);
P00: FieldType = 104
P00: MemoryOrder = "0 "
P00: description = "BASE STATE PRESURE"
P00: units = "Pa"
P00: stagger = ""
float TLP(Time);
TLP: FieldType = 104
TLP: MemoryOrder = "0 "
TLP: description = "BASE STATE LAPSE RATE"
TLP: units = ""
TLP: stagger = ""
float TISO(Time);
TISO: FieldType = 104
TISO: MemoryOrder = "0 "
TISO: description = "TEMP AT WHICH THE BASE T TURNS CONST"
TISO: units = "K"
TISO: stagger = ""
float MAX_MSTFX(Time);
MAX_MSTFX: FieldType = 104
MAX_MSTFX: MemoryOrder = "0 "
MAX_MSTFX: description = "Max map factor in domain"
MAX_MSTFX: units = ""
MAX_MSTFX: stagger = ""
float MAX_MSTFY(Time);
MAX_MSTFY: FieldType = 104
MAX_MSTFY: MemoryOrder = "0 "
MAX_MSTFY: description = "Max map factor in domain"
MAX_MSTFY: units = ""
MAX_MSTFY: stagger = ""
float RAINC(Time,south_north,west_east);
RAINC: FieldType = 104
RAINC: MemoryOrder = "XY "
RAINC: description = "ACCUMULATED TOTAL CUMULUS PRECIPITATION"
RAINC: units = "mm"
RAINC: stagger = ""
RAINC: coordinates = "XLONG XLAT"
float RAINSH(Time,south_north,west_east);
RAINSH: FieldType = 104
RAINSH: MemoryOrder = "XY "
RAINSH: description = "ACCUMULATED SHALLOW CUMULUS PRECIPITATION"
RAINSH: units = "mm"
RAINSH: stagger = ""
RAINSH: coordinates = "XLONG XLAT"
float RAINNC(Time,south_north,west_east);
RAINNC: FieldType = 104
RAINNC: MemoryOrder = "XY "
RAINNC: description = "ACCUMULATED TOTAL GRID SCALE PRECIPITATION"
RAINNC: units = "mm"
RAINNC: stagger = ""
RAINNC: coordinates = "XLONG XLAT"
float SNOWNC(Time,south_north,west_east);
SNOWNC: FieldType = 104
SNOWNC: MemoryOrder = "XY "
SNOWNC: description = "ACCUMULATED TOTAL GRID SCALE SNOW AND ICE"
SNOWNC: units = "mm"
SNOWNC: stagger = ""
SNOWNC: coordinates = "XLONG XLAT"
float GRAUPELNC(Time,south_north,west_east);
GRAUPELNC: FieldType = 104
GRAUPELNC: MemoryOrder = "XY "
GRAUPELNC: description = "ACCUMULATED TOTAL GRID SCALE GRAUPEL"
GRAUPELNC: units = "mm"
GRAUPELNC: stagger = ""
GRAUPELNC: coordinates = "XLONG XLAT"
float HAILNC(Time,south_north,west_east);
HAILNC: FieldType = 104
HAILNC: MemoryOrder = "XY "
HAILNC: description = "ACCUMULATED TOTAL GRID SCALE HAIL"
HAILNC: units = "mm"
HAILNC: stagger = ""
HAILNC: coordinates = "XLONG XLAT"
float XMB_SHALLOW(Time,south_north,west_east);
XMB_SHALLOW: FieldType = 104
XMB_SHALLOW: MemoryOrder = "XY "
XMB_SHALLOW: description = "MASSFLUX FROM SHALLOW CONVECTION (G3 only)"
XMB_SHALLOW: units = ""
XMB_SHALLOW: stagger = ""
XMB_SHALLOW: coordinates = "XLONG XLAT"
int K22_SHALLOW(Time,south_north,west_east);
K22_SHALLOW: FieldType = 106
K22_SHALLOW: MemoryOrder = "XY "
K22_SHALLOW: description = "K22 LEVEL FROM SHALLOW CONVECTION (G3 only)"
K22_SHALLOW: units = ""
K22_SHALLOW: stagger = ""
K22_SHALLOW: coordinates = "XLONG XLAT"
int KBCON_SHALLOW(Time,south_north,west_east);
KBCON_SHALLOW: FieldType = 106
KBCON_SHALLOW: MemoryOrder = "XY "
KBCON_SHALLOW: description = "KBCON LEVEL FROM SHALLOW CONVECTION (G3 only)"
KBCON_SHALLOW: units = ""
KBCON_SHALLOW: stagger = ""
KBCON_SHALLOW: coordinates = "XLONG XLAT"
int KTOP_SHALLOW(Time,south_north,west_east);
KTOP_SHALLOW: FieldType = 106
KTOP_SHALLOW: MemoryOrder = "XY "
KTOP_SHALLOW: description = "KTOP LEVEL FROM SHALLOW CONVECTION (G3 only)"
KTOP_SHALLOW: units = ""
KTOP_SHALLOW: stagger = ""
KTOP_SHALLOW: coordinates = "XLONG XLAT"
int STEPAVE_COUNT(Time);
STEPAVE_COUNT: FieldType = 106
STEPAVE_COUNT: MemoryOrder = "0 "
STEPAVE_COUNT: description = "time steps contained in averages for convective transport"
STEPAVE_COUNT: units = ""
STEPAVE_COUNT: stagger = ""
float CLDFRA(Time,bottom_top,south_north,west_east);
CLDFRA: FieldType = 104
CLDFRA: MemoryOrder = "XYZ"
CLDFRA: description = "CLOUD FRACTION"
CLDFRA: units = ""
CLDFRA: stagger = ""
CLDFRA: coordinates = "XLONG XLAT"
float SWDOWN(Time,south_north,west_east);
SWDOWN: FieldType = 104
SWDOWN: MemoryOrder = "XY "
SWDOWN: description = "DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE"
SWDOWN: units = "W m-2"
SWDOWN: stagger = ""
SWDOWN: coordinates = "XLONG XLAT"
float GLW(Time,south_north,west_east);
GLW: FieldType = 104
GLW: MemoryOrder = "XY "
GLW: description = "DOWNWARD LONG WAVE FLUX AT GROUND SURFACE"
GLW: units = "W m-2"
GLW: stagger = ""
GLW: coordinates = "XLONG XLAT"
float SWNORM(Time,south_north,west_east);
SWNORM: FieldType = 104
SWNORM: MemoryOrder = "XY "
SWNORM: description = "NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT)"
SWNORM: units = "W m-2"
SWNORM: stagger = ""
SWNORM: coordinates = "XLONG XLAT"
float ACSWUPT(Time,south_north,west_east);
ACSWUPT: FieldType = 104
ACSWUPT: MemoryOrder = "XY "
ACSWUPT: description = "ACCUMULATED UPWELLING SHORTWAVE FLUX AT TOP"
ACSWUPT: units = "J m-2"
ACSWUPT: stagger = ""
ACSWUPT: coordinates = "XLONG XLAT"
float ACSWUPTC(Time,south_north,west_east);
ACSWUPTC: FieldType = 104
ACSWUPTC: MemoryOrder = "XY "
ACSWUPTC: description = "ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP"
ACSWUPTC: units = "J m-2"
ACSWUPTC: stagger = ""
ACSWUPTC: coordinates = "XLONG XLAT"
float ACSWDNT(Time,south_north,west_east);
ACSWDNT: FieldType = 104
ACSWDNT: MemoryOrder = "XY "
ACSWDNT: description = "ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT TOP"
ACSWDNT: units = "J m-2"
ACSWDNT: stagger = ""
ACSWDNT: coordinates = "XLONG XLAT"
float ACSWDNTC(Time,south_north,west_east);
ACSWDNTC: FieldType = 104
ACSWDNTC: MemoryOrder = "XY "
ACSWDNTC: description = "ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP"
ACSWDNTC: units = "J m-2"
ACSWDNTC: stagger = ""
ACSWDNTC: coordinates = "XLONG XLAT"
float ACSWUPB(Time,south_north,west_east);
ACSWUPB: FieldType = 104
ACSWUPB: MemoryOrder = "XY "
ACSWUPB: description = "ACCUMULATED UPWELLING SHORTWAVE FLUX AT BOTTOM"
ACSWUPB: units = "J m-2"
ACSWUPB: stagger = ""
ACSWUPB: coordinates = "XLONG XLAT"
float ACSWUPBC(Time,south_north,west_east);
ACSWUPBC: FieldType = 104
ACSWUPBC: MemoryOrder = "XY "
ACSWUPBC: description = "ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM"
ACSWUPBC: units = "J m-2"
ACSWUPBC: stagger = ""
ACSWUPBC: coordinates = "XLONG XLAT"
float ACSWDNB(Time,south_north,west_east);
ACSWDNB: FieldType = 104
ACSWDNB: MemoryOrder = "XY "
ACSWDNB: description = "ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT BOTTOM"
ACSWDNB: units = "J m-2"
ACSWDNB: stagger = ""
ACSWDNB: coordinates = "XLONG XLAT"
float ACSWDNBC(Time,south_north,west_east);
ACSWDNBC: FieldType = 104
ACSWDNBC: MemoryOrder = "XY "
ACSWDNBC: description = "ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM"
ACSWDNBC: units = "J m-2"
ACSWDNBC: stagger = ""
ACSWDNBC: coordinates = "XLONG XLAT"
float ACLWUPT(Time,south_north,west_east);
ACLWUPT: FieldType = 104
ACLWUPT: MemoryOrder = "XY "
ACLWUPT: description = "ACCUMULATED UPWELLING LONGWAVE FLUX AT TOP"
ACLWUPT: units = "J m-2"
ACLWUPT: stagger = ""
ACLWUPT: coordinates = "XLONG XLAT"
float ACLWUPTC(Time,south_north,west_east);
ACLWUPTC: FieldType = 104
ACLWUPTC: MemoryOrder = "XY "
ACLWUPTC: description = "ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP"
ACLWUPTC: units = "J m-2"
ACLWUPTC: stagger = ""
ACLWUPTC: coordinates = "XLONG XLAT"
float ACLWDNT(Time,south_north,west_east);
ACLWDNT: FieldType = 104
ACLWDNT: MemoryOrder = "XY "
ACLWDNT: description = "ACCUMULATED DOWNWELLING LONGWAVE FLUX AT TOP"
ACLWDNT: units = "J m-2"
ACLWDNT: stagger = ""
ACLWDNT: coordinates = "XLONG XLAT"
float ACLWDNTC(Time,south_north,west_east);
ACLWDNTC: FieldType = 104
ACLWDNTC: MemoryOrder = "XY "
ACLWDNTC: description = "ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP"
ACLWDNTC: units = "J m-2"
ACLWDNTC: stagger = ""
ACLWDNTC: coordinates = "XLONG XLAT"
float ACLWUPB(Time,south_north,west_east);
ACLWUPB: FieldType = 104
ACLWUPB: MemoryOrder = "XY "
ACLWUPB: description = "ACCUMULATED UPWELLING LONGWAVE FLUX AT BOTTOM"
ACLWUPB: units = "J m-2"
ACLWUPB: stagger = ""
ACLWUPB: coordinates = "XLONG XLAT"
float ACLWUPBC(Time,south_north,west_east);
ACLWUPBC: FieldType = 104
ACLWUPBC: MemoryOrder = "XY "
ACLWUPBC: description = "ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM"
ACLWUPBC: units = "J m-2"
ACLWUPBC: stagger = ""
ACLWUPBC: coordinates = "XLONG XLAT"
float ACLWDNB(Time,south_north,west_east);
ACLWDNB: FieldType = 104
ACLWDNB: MemoryOrder = "XY "
ACLWDNB: description = "ACCUMULATED DOWNWELLING LONGWAVE FLUX AT BOTTOM"
ACLWDNB: units = "J m-2"
ACLWDNB: stagger = ""
ACLWDNB: coordinates = "XLONG XLAT"
float ACLWDNBC(Time,south_north,west_east);
ACLWDNBC: FieldType = 104
ACLWDNBC: MemoryOrder = "XY "
ACLWDNBC: description = "ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM"
ACLWDNBC: units = "J m-2"
ACLWDNBC: stagger = ""
ACLWDNBC: coordinates = "XLONG XLAT"
int I_ACSWUPT(Time,south_north,west_east);
I_ACSWUPT: FieldType = 106
I_ACSWUPT: MemoryOrder = "XY "
I_ACSWUPT: description = "BUCKET FOR UPWELLING SHORTWAVE FLUX AT TOP"
I_ACSWUPT: units = "J m-2"
I_ACSWUPT: stagger = ""
I_ACSWUPT: coordinates = "XLONG XLAT"
int I_ACSWUPTC(Time,south_north,west_east);
I_ACSWUPTC: FieldType = 106
I_ACSWUPTC: MemoryOrder = "XY "
I_ACSWUPTC: description = "BUCKET FOR UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP"
I_ACSWUPTC: units = "J m-2"
I_ACSWUPTC: stagger = ""
I_ACSWUPTC: coordinates = "XLONG XLAT"
int I_ACSWDNT(Time,south_north,west_east);
I_ACSWDNT: FieldType = 106
I_ACSWDNT: MemoryOrder = "XY "
I_ACSWDNT: description = "BUCKET FOR DOWNWELLING SHORTWAVE FLUX AT TOP"
I_ACSWDNT: units = "J m-2"
I_ACSWDNT: stagger = ""
I_ACSWDNT: coordinates = "XLONG XLAT"
int I_ACSWDNTC(Time,south_north,west_east);
I_ACSWDNTC: FieldType = 106
I_ACSWDNTC: MemoryOrder = "XY "
I_ACSWDNTC: description = "BUCKET FOR DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP"
I_ACSWDNTC: units = "J m-2"
I_ACSWDNTC: stagger = ""
I_ACSWDNTC: coordinates = "XLONG XLAT"
int I_ACSWUPB(Time,south_north,west_east);
I_ACSWUPB: FieldType = 106
I_ACSWUPB: MemoryOrder = "XY "
I_ACSWUPB: description = "BUCKET FOR UPWELLING SHORTWAVE FLUX AT BOTTOM"
I_ACSWUPB: units = "J m-2"
I_ACSWUPB: stagger = ""
I_ACSWUPB: coordinates = "XLONG XLAT"
int I_ACSWUPBC(Time,south_north,west_east);
I_ACSWUPBC: FieldType = 106
I_ACSWUPBC: MemoryOrder = "XY "
I_ACSWUPBC: description = "BUCKET FOR UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM"
I_ACSWUPBC: units = "J m-2"
I_ACSWUPBC: stagger = ""
I_ACSWUPBC: coordinates = "XLONG XLAT"
int I_ACSWDNB(Time,south_north,west_east);
I_ACSWDNB: FieldType = 106
I_ACSWDNB: MemoryOrder = "XY "
I_ACSWDNB: description = "BUCKET FOR DOWNWELLING SHORTWAVE FLUX AT BOTTOM"
I_ACSWDNB: units = "J m-2"
I_ACSWDNB: stagger = ""
I_ACSWDNB: coordinates = "XLONG XLAT"
int I_ACSWDNBC(Time,south_north,west_east);
I_ACSWDNBC: FieldType = 106
I_ACSWDNBC: MemoryOrder = "XY "
I_ACSWDNBC: description = "BUCKET FOR DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM"
I_ACSWDNBC: units = "J m-2"
I_ACSWDNBC: stagger = ""
I_ACSWDNBC: coordinates = "XLONG XLAT"
int I_ACLWUPT(Time,south_north,west_east);
I_ACLWUPT: FieldType = 106
I_ACLWUPT: MemoryOrder = "XY "
I_ACLWUPT: description = "BUCKET FOR UPWELLING LONGWAVE FLUX AT TOP"
I_ACLWUPT: units = "J m-2"
I_ACLWUPT: stagger = ""
I_ACLWUPT: coordinates = "XLONG XLAT"
int I_ACLWUPTC(Time,south_north,west_east);
I_ACLWUPTC: FieldType = 106
I_ACLWUPTC: MemoryOrder = "XY "
I_ACLWUPTC: description = "BUCKET FOR UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP"
I_ACLWUPTC: units = "J m-2"
I_ACLWUPTC: stagger = ""
I_ACLWUPTC: coordinates = "XLONG XLAT"
int I_ACLWDNT(Time,south_north,west_east);
I_ACLWDNT: FieldType = 106
I_ACLWDNT: MemoryOrder = "XY "
I_ACLWDNT: description = "BUCKET FOR DOWNWELLING LONGWAVE FLUX AT TOP"
I_ACLWDNT: units = "J m-2"
I_ACLWDNT: stagger = ""
I_ACLWDNT: coordinates = "XLONG XLAT"
int I_ACLWDNTC(Time,south_north,west_east);
I_ACLWDNTC: FieldType = 106
I_ACLWDNTC: MemoryOrder = "XY "
I_ACLWDNTC: description = "BUCKET FOR DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP"
I_ACLWDNTC: units = "J m-2"
I_ACLWDNTC: stagger = ""
I_ACLWDNTC: coordinates = "XLONG XLAT"
int I_ACLWUPB(Time,south_north,west_east);
I_ACLWUPB: FieldType = 106
I_ACLWUPB: MemoryOrder = "XY "
I_ACLWUPB: description = "BUCKET FOR UPWELLING LONGWAVE FLUX AT BOTTOM"
I_ACLWUPB: units = "J m-2"
I_ACLWUPB: stagger = ""
I_ACLWUPB: coordinates = "XLONG XLAT"
int I_ACLWUPBC(Time,south_north,west_east);
I_ACLWUPBC: FieldType = 106
I_ACLWUPBC: MemoryOrder = "XY "
I_ACLWUPBC: description = "BUCKET FOR UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM"
I_ACLWUPBC: units = "J m-2"
I_ACLWUPBC: stagger = ""
I_ACLWUPBC: coordinates = "XLONG XLAT"
int I_ACLWDNB(Time,south_north,west_east);
I_ACLWDNB: FieldType = 106
I_ACLWDNB: MemoryOrder = "XY "
I_ACLWDNB: description = "BUCKET FOR DOWNWELLING LONGWAVE FLUX AT BOTTOM"
I_ACLWDNB: units = "J m-2"
I_ACLWDNB: stagger = ""
I_ACLWDNB: coordinates = "XLONG XLAT"
int I_ACLWDNBC(Time,south_north,west_east);
I_ACLWDNBC: FieldType = 106
I_ACLWDNBC: MemoryOrder = "XY "
I_ACLWDNBC: description = "BUCKET FOR DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM"
I_ACLWDNBC: units = "J m-2"
I_ACLWDNBC: stagger = ""
I_ACLWDNBC: coordinates = "XLONG XLAT"
float SWUPT(Time,south_north,west_east);
SWUPT: FieldType = 104
SWUPT: MemoryOrder = "XY "
SWUPT: description = "INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT TOP"
SWUPT: units = "W m-2"
SWUPT: stagger = ""
SWUPT: coordinates = "XLONG XLAT"
float SWUPTC(Time,south_north,west_east);
SWUPTC: FieldType = 104
SWUPTC: MemoryOrder = "XY "
SWUPTC: description = "INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP"
SWUPTC: units = "W m-2"
SWUPTC: stagger = ""
SWUPTC: coordinates = "XLONG XLAT"
float SWDNT(Time,south_north,west_east);
SWDNT: FieldType = 104
SWDNT: MemoryOrder = "XY "
SWDNT: description = "INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT TOP"
SWDNT: units = "W m-2"
SWDNT: stagger = ""
SWDNT: coordinates = "XLONG XLAT"
float SWDNTC(Time,south_north,west_east);
SWDNTC: FieldType = 104
SWDNTC: MemoryOrder = "XY "
SWDNTC: description = "INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP"
SWDNTC: units = "W m-2"
SWDNTC: stagger = ""
SWDNTC: coordinates = "XLONG XLAT"
float SWUPB(Time,south_north,west_east);
SWUPB: FieldType = 104
SWUPB: MemoryOrder = "XY "
SWUPB: description = "INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT BOTTOM"
SWUPB: units = "W m-2"
SWUPB: stagger = ""
SWUPB: coordinates = "XLONG XLAT"
float SWUPBC(Time,south_north,west_east);
SWUPBC: FieldType = 104
SWUPBC: MemoryOrder = "XY "
SWUPBC: description = "INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM"
SWUPBC: units = "W m-2"
SWUPBC: stagger = ""
SWUPBC: coordinates = "XLONG XLAT"
float SWDNB(Time,south_north,west_east);
SWDNB: FieldType = 104
SWDNB: MemoryOrder = "XY "
SWDNB: description = "INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT BOTTOM"
SWDNB: units = "W m-2"
SWDNB: stagger = ""
SWDNB: coordinates = "XLONG XLAT"
float SWDNBC(Time,south_north,west_east);
SWDNBC: FieldType = 104
SWDNBC: MemoryOrder = "XY "
SWDNBC: description = "INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM"
SWDNBC: units = "W m-2"
SWDNBC: stagger = ""
SWDNBC: coordinates = "XLONG XLAT"
float LWUPT(Time,south_north,west_east);
LWUPT: FieldType = 104
LWUPT: MemoryOrder = "XY "
LWUPT: description = "INSTANTANEOUS UPWELLING LONGWAVE FLUX AT TOP"
LWUPT: units = "W m-2"
LWUPT: stagger = ""
LWUPT: coordinates = "XLONG XLAT"
float LWUPTC(Time,south_north,west_east);
LWUPTC: FieldType = 104
LWUPTC: MemoryOrder = "XY "
LWUPTC: description = "INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP"
LWUPTC: units = "W m-2"
LWUPTC: stagger = ""
LWUPTC: coordinates = "XLONG XLAT"
float LWDNT(Time,south_north,west_east);
LWDNT: FieldType = 104
LWDNT: MemoryOrder = "XY "
LWDNT: description = "INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT TOP"
LWDNT: units = "W m-2"
LWDNT: stagger = ""
LWDNT: coordinates = "XLONG XLAT"
float LWDNTC(Time,south_north,west_east);
LWDNTC: FieldType = 104
LWDNTC: MemoryOrder = "XY "
LWDNTC: description = "INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP"
LWDNTC: units = "W m-2"
LWDNTC: stagger = ""
LWDNTC: coordinates = "XLONG XLAT"
float LWUPB(Time,south_north,west_east);
LWUPB: FieldType = 104
LWUPB: MemoryOrder = "XY "
LWUPB: description = "INSTANTANEOUS UPWELLING LONGWAVE FLUX AT BOTTOM"
LWUPB: units = "W m-2"
LWUPB: stagger = ""
LWUPB: coordinates = "XLONG XLAT"
float LWUPBC(Time,south_north,west_east);
LWUPBC: FieldType = 104
LWUPBC: MemoryOrder = "XY "
LWUPBC: description = "INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM"
LWUPBC: units = "W m-2"
LWUPBC: stagger = ""
LWUPBC: coordinates = "XLONG XLAT"
float LWDNB(Time,south_north,west_east);
LWDNB: FieldType = 104
LWDNB: MemoryOrder = "XY "
LWDNB: description = "INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT BOTTOM"
LWDNB: units = "W m-2"
LWDNB: stagger = ""
LWDNB: coordinates = "XLONG XLAT"
float LWDNBC(Time,south_north,west_east);
LWDNBC: FieldType = 104
LWDNBC: MemoryOrder = "XY "
LWDNBC: description = "INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM"
LWDNBC: units = "W m-2"
LWDNBC: stagger = ""
LWDNBC: coordinates = "XLONG XLAT"
float OLR(Time,south_north,west_east);
OLR: FieldType = 104
OLR: MemoryOrder = "XY "
OLR: description = "TOA OUTGOING LONG WAVE"
OLR: units = "W m-2"
OLR: stagger = ""
OLR: coordinates = "XLONG XLAT"
float XLAT(Time,south_north,west_east);
XLAT: FieldType = 104
XLAT: MemoryOrder = "XY "
XLAT: description = "LATITUDE, SOUTH IS NEGATIVE"
XLAT: units = "degree_north"
XLAT: stagger = ""
float XLONG(Time,south_north,west_east);
XLONG: FieldType = 104
XLONG: MemoryOrder = "XY "
XLONG: description = "LONGITUDE, WEST IS NEGATIVE"
XLONG: units = "degree_east"
XLONG: stagger = ""
float XLAT_U(Time,south_north,west_east_stag);
XLAT_U: FieldType = 104
XLAT_U: MemoryOrder = "XY "
XLAT_U: description = "LATITUDE, SOUTH IS NEGATIVE"
XLAT_U: units = "degree_north"
XLAT_U: stagger = "X"
XLAT_U: coordinates = "XLONG_U XLAT_U"
float XLONG_U(Time,south_north,west_east_stag);
XLONG_U: FieldType = 104
XLONG_U: MemoryOrder = "XY "
XLONG_U: description = "LONGITUDE, WEST IS NEGATIVE"
XLONG_U: units = "degree_east"
XLONG_U: stagger = "X"
XLONG_U: coordinates = "XLONG_U XLAT_U"
float XLAT_V(Time,south_north_stag,west_east);
XLAT_V: FieldType = 104
XLAT_V: MemoryOrder = "XY "
XLAT_V: description = "LATITUDE, SOUTH IS NEGATIVE"
XLAT_V: units = "degree_north"
XLAT_V: stagger = "Y"
XLAT_V: coordinates = "XLONG_V XLAT_V"
float XLONG_V(Time,south_north_stag,west_east);
XLONG_V: FieldType = 104
XLONG_V: MemoryOrder = "XY "
XLONG_V: description = "LONGITUDE, WEST IS NEGATIVE"
XLONG_V: units = "degree_east"
XLONG_V: stagger = "Y"
XLONG_V: coordinates = "XLONG_V XLAT_V"
float ALBEDO(Time,south_north,west_east);
ALBEDO: FieldType = 104
ALBEDO: MemoryOrder = "XY "
ALBEDO: description = "ALBEDO"
ALBEDO: units = "-"
ALBEDO: stagger = ""
ALBEDO: coordinates = "XLONG XLAT"
float ALBBCK(Time,south_north,west_east);
ALBBCK: FieldType = 104
ALBBCK: MemoryOrder = "XY "
ALBBCK: description = "BACKGROUND ALBEDO"
ALBBCK: units = ""
ALBBCK: stagger = ""
ALBBCK: coordinates = "XLONG XLAT"
float EMISS(Time,south_north,west_east);
EMISS: FieldType = 104
EMISS: MemoryOrder = "XY "
EMISS: description = "SURFACE EMISSIVITY"
EMISS: units = ""
EMISS: stagger = ""
EMISS: coordinates = "XLONG XLAT"
float NOAHRES(Time,south_north,west_east);
NOAHRES: FieldType = 104
NOAHRES: MemoryOrder = "XY "
NOAHRES: description = "RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET"
NOAHRES: units = "W m{-2}"
NOAHRES: stagger = ""
NOAHRES: coordinates = "XLONG XLAT"
float TMN(Time,south_north,west_east);
TMN: FieldType = 104
TMN: MemoryOrder = "XY "
TMN: description = "SOIL TEMPERATURE AT LOWER BOUNDARY"
TMN: units = "K"
TMN: stagger = ""
TMN: coordinates = "XLONG XLAT"
float XLAND(Time,south_north,west_east);
XLAND: FieldType = 104
XLAND: MemoryOrder = "XY "
XLAND: description = "LAND MASK (1 FOR LAND, 2 FOR WATER)"
XLAND: units = ""
XLAND: stagger = ""
XLAND: coordinates = "XLONG XLAT"
float UST(Time,south_north,west_east);
UST: FieldType = 104
UST: MemoryOrder = "XY "
UST: description = "U* IN SIMILARITY THEORY"
UST: units = "m s-1"
UST: stagger = ""
UST: coordinates = "XLONG XLAT"
float PBLH(Time,south_north,west_east);
PBLH: FieldType = 104
PBLH: MemoryOrder = "XY "
PBLH: description = "PBL HEIGHT"
PBLH: units = "m"
PBLH: stagger = ""
PBLH: coordinates = "XLONG XLAT"
float HFX(Time,south_north,west_east);
HFX: FieldType = 104
HFX: MemoryOrder = "XY "
HFX: description = "UPWARD HEAT FLUX AT THE SURFACE"
HFX: units = "W m-2"
HFX: stagger = ""
HFX: coordinates = "XLONG XLAT"
float QFX(Time,south_north,west_east);
QFX: FieldType = 104
QFX: MemoryOrder = "XY "
QFX: description = "UPWARD MOISTURE FLUX AT THE SURFACE"
QFX: units = "kg m-2 s-1"
QFX: stagger = ""
QFX: coordinates = "XLONG XLAT"
float LH(Time,south_north,west_east);
LH: FieldType = 104
LH: MemoryOrder = "XY "
LH: description = "LATENT HEAT FLUX AT THE SURFACE"
LH: units = "W m-2"
LH: stagger = ""
LH: coordinates = "XLONG XLAT"
float ACHFX(Time,south_north,west_east);
ACHFX: FieldType = 104
ACHFX: MemoryOrder = "XY "
ACHFX: description = "ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE"
ACHFX: units = "J m-2"
ACHFX: stagger = ""
ACHFX: coordinates = "XLONG XLAT"
float ACLHF(Time,south_north,west_east);
ACLHF: FieldType = 104
ACLHF: MemoryOrder = "XY "
ACLHF: description = "ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE"
ACLHF: units = "J m-2"
ACLHF: stagger = ""
ACLHF: coordinates = "XLONG XLAT"
float SNOWC(Time,south_north,west_east);
SNOWC: FieldType = 104
SNOWC: MemoryOrder = "XY "
SNOWC: description = "FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER)"
SNOWC: units = ""
SNOWC: stagger = ""
SNOWC: coordinates = "XLONG XLAT"
float SR(Time,south_north,west_east);
SR: FieldType = 104
SR: MemoryOrder = "XY "
SR: description = "fraction of frozen precipitation"
SR: units = "-"
SR: stagger = ""
SR: coordinates = "XLONG XLAT"
float POTEVP(Time,south_north,west_east);
POTEVP: FieldType = 104
POTEVP: MemoryOrder = "XY "
POTEVP: description = "accumulated potential evaporation"
POTEVP: units = "W m-2"
POTEVP: stagger = ""
POTEVP: coordinates = "XLONG XLAT"
float SNOPCX(Time,south_north,west_east);
SNOPCX: FieldType = 104
SNOPCX: MemoryOrder = "XY "
SNOPCX: description = "snow phase change heat flux"
SNOPCX: units = "W m-2"
SNOPCX: stagger = ""
SNOPCX: coordinates = "XLONG XLAT"
float SOILTB(Time,south_north,west_east);
SOILTB: FieldType = 104
SOILTB: MemoryOrder = "XY "
SOILTB: description = "bottom soil temperature"
SOILTB: units = "K"
SOILTB: stagger = ""
SOILTB: coordinates = "XLONG XLAT"
int SAVE_TOPO_FROM_REAL(Time);
SAVE_TOPO_FROM_REAL: FieldType = 106
SAVE_TOPO_FROM_REAL: MemoryOrder = "0 "
SAVE_TOPO_FROM_REAL: description = "1=original topo from real/0=topo modified by WRF"
SAVE_TOPO_FROM_REAL: units = "flag"
SAVE_TOPO_FROM_REAL: stagger = ""
float AVG_FUEL_FRAC(Time,south_north,west_east);
AVG_FUEL_FRAC: FieldType = 104
AVG_FUEL_FRAC: MemoryOrder = "XY "
AVG_FUEL_FRAC: description = "fuel remaining averaged to atmospheric grid"
AVG_FUEL_FRAC: units = "1"
AVG_FUEL_FRAC: stagger = "Z"
AVG_FUEL_FRAC: coordinates = "XLONG XLAT"
float UAH(Time,south_north,west_east_stag);
UAH: FieldType = 104
UAH: MemoryOrder = "XY "
UAH: description = "wind at fire_wind_height"
UAH: units = "m/s"
UAH: stagger = "X"
UAH: coordinates = "XLONG_U XLAT_U"
float VAH(Time,south_north_stag,west_east);
VAH: FieldType = 104
VAH: MemoryOrder = "XY "
VAH: description = "wind at fire_wind_height"
VAH: units = "m/s"
VAH: stagger = "Y"
VAH: coordinates = "XLONG_V XLAT_V"
int SEED1(Time);
SEED1: FieldType = 106
SEED1: MemoryOrder = "0 "
SEED1: description = "RANDOM SEED NUMBER 1"
SEED1: units = ""
SEED1: stagger = ""
int SEED2(Time);
SEED2: FieldType = 106
SEED2: MemoryOrder = "0 "
SEED2: description = "RANDOM SEED NUMBER 2"
SEED2: units = ""
SEED2: stagger = ""
Unlimited dimension: Time

复制代码

NC数据信息如上。

chenxl · 发表于 2015-6-27 10:57:38

老师，我在运行该脚本时出现了下面的错误，我该怎样修改呢？

chenxl · 发表于 2015-6-27 12:59:55

老师，会出现这个问题是因为版本的原因吗我下载的是1.2.6 我在官网上只看到了这个的下载。

MeteoInfo · 发表于 2015-6-27 16:53:58

chenxl 发表于 2015-6-27 12:59
老师，会出现这个问题是因为版本的原因吗我下载的是1.2.6 我在官网上只看到了这个的下载。

可以在我的百度空间下载最新版本：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

MeteoInfo · 发表于 2015-6-27 17:50:31

良辰发表于 2015-6-27 09:50
王老师，刚才我查看了一下原始的nc数据，里面数据没有缺失，但是垂直高度上，U\V一致，W比UV多。如果我 ...

可以考虑直接从nc文件转arl文件

良辰 · 发表于 2015-6-27 20:22:58

MeteoInfo 发表于 2015-6-27 17:50
可以考虑直接从nc文件转arl文件

王老师，nc数据为地形追随坐标，levels如何设置？另外aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 2)中的2修改为1就可以了吗？谢谢王老师不厌其烦的指导。

MeteoInfo · 发表于 2015-6-27 21:04:00

良辰发表于 2015-6-27 20:22
王老师，nc数据为地形追随坐标，levels如何设置？另外aDH = arlDI.getDataHead(mydata.getProjectionInfo ...

应该是，你试试。下面是ARL中垂直坐标的设置：
1-sigma (fraction)
2-pressure (mb)
3-terrain (fraction)
4-hybrid (mb: offset.fraction)

chenxl · 发表于 2015-6-28 15:41:57

MeteoInfo 发表于 2015-6-27 16:53
可以在我的百度空间下载最新版本：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

老师，我已经分别用MeteoInfo和hysplit将NC数据转换为了ARL，但是为什么做不了轨迹呢?能不能帮我看看。谢谢！这些文件都传到百度云盘了。http://pan.baidu.com/s/1pJoSHrP

MeteoInfo · 发表于 2015-6-29 09:58:30

chenxl 发表于 2015-6-28 15:41
老师，我已经分别用MeteoInfo和hysplit将NC数据转换为了ARL，但是为什么做不了轨迹呢?能不能帮我看看。谢 ...

修改了一下转换脚本代码，可以生成一个ARL格式的数据。有几个难点：

1. WRF输出的是地形追随坐标，aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 1)语句中最后的1代表sigma坐标。
2. WRF输出的U, V, W等很多变量坐标是错列的（staggered），需要进行插值处理（http://www.openwfm.org/wiki/How_to_interpret_WRF_variables），MeteoInfo中的GridData类有unStagger_X()和unStagger_Y()方法来处理。U和V比较好处理，W就比较困难了，我现在还不太清楚sigma坐标中W要求的单位是什么（在pressure坐标中W的单位是hPa / m3），目前WRF输出的W的单位是m / m3。还有就是W是垂直方向上的staggered，需要在垂直层上插值。
3. 你给的WRF输出数据中不知道温度是哪个变量（T 是perturbation potential temperature (theta-t0)）、位势高度、相对湿度对应的变量也没找到。

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2015-6-29
# Purpose: Convert WRF out data to ARL data
# Note: Sample
#-----------------------------------------------------------
from org.meteoinfo.data.meteodata import MeteoDataInfo
from org.meteoinfo.data.meteodata.arl import ARLDataInfo
from org.meteoinfo.data.meteodata.arl import DataLabel
import os
#---- Set directories
dataDir = 'D:/Temp/test/hysplit-arl'
#---- Set output data file
outFile = os.path.join(dataDir, 'test7.arl')
#---- Read a GrADS data file
mydata = MeteoDataInfo()
infile = os.path.join(dataDir, 'wrfout_d01_2006-07-13_00')
print infile
mydata.openNetCDFData(infile)
print 'NetCDF file has been opened...'
#---- Set output ARL data info
arlDI = ARLDataInfo()
#---- Set variable and level list
gvars = ['PSFC','T2','U10','V10','T','U','V','W']
avars = ['PRSS','T02M','U10M','V10M','TEMP','UWND','VWND','WWND']
levels = [1, 0.9965000152587891, 0.9879999756813049, 0.9765000343322754, 0.9620000123977661, \
0.9440000057220459, 0.9214999675750732, 0.8945000171661377, 0.8547878265380859, 0.8043636083602905, \
0.7539393305778503, 0.7035150527954102, 0.6350234150886536, 0.5527187585830688, 0.478573739528656, \
0.41191452741622925, 0.35211411118507385, 0.2985895872116089, 0.2507995367050171, 0.20824146270751953, \
0.1704493761062622, 0.13699156045913696, 0.10746840387582779, 0.08151029050350189, 0.05877562612295151, \
0.03894893825054169, 0.021738944575190544, 0.0068765804171562195]
for l in levels:
arlDI.levels.add(l)
if l == 1:
arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M'])
else:
arlDI.LevelVarList.add(['TEMP','UWND','VWND','WWND'])
#---- Write ARL data file
dataInfo = mydata.getDataInfo()
arlDI.createDataFile(outFile)
arlDI.X = dataInfo.getXDimension().getValues()
arlDI.Y = dataInfo.getYDimension().getValues()
variables = dataInfo.getVariables()
tNum = dataInfo.getTimeNum()
for t in range(0, 1):
mydata.setTimeIndex(t)
atime = dataInfo.getTimes().get(t)
aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 1)
arlDI.writeIndexRecord(atime, aDH)
lidx = 0
for l in arlDI.levels:
print l
for v in arlDI.LevelVarList[lidx]:
vName = gvars[avars.index(v)]
print vName
gData = None
if lidx == 0:
mydata.setLevelIndex(lidx)
gData = mydata.getGridData(vName)
else:
if v == 'WWND':
variable = dataInfo.getVariable('U')
nlidx = variable.getZDimension().getDimValue().indexOf(l)
mydata.setLevelIndex(nlidx)
gData = mydata.getGridData(vName)
mydata.setLevelIndex(nlidx + 1)
gData1 = mydata.getGridData(vName)
gData = gData.add(gData1).mul(0.5)
else:
variable = dataInfo.getVariable(vName)
nlidx = variable.getZDimension().getDimValue().indexOf(l)
mydata.setLevelIndex(nlidx)
gData = mydata.getGridData(vName)
if v == 'PRSS':
gData = gData.div(100)
elif v == 'UWND':
gData = gData.unStagger_X()
elif v == 'VWND':
gData = gData.unStagger_Y()
aDL = DataLabel(atime)
aDL.setLevel(lidx)
aDL.setVarName(v)
aDL.setGrid(99)
aDL.setForecast(0)
arlDI.writeGridData(aDL, gData)
lidx += 1
arlDI.closeDataFile()
print 'Finished!'

#--------------------------------------------------------         
# Author: Yaqiang Wang                                           
# Date: 2015-6-29                                                
# Purpose: Convert WRF out data to ARL data   
# Note: Sample                                                    
#----------------------------------------------------------- 
from org.meteoinfo.data.meteodata import MeteoDataInfo 
from org.meteoinfo.data.meteodata.arl import ARLDataInfo 
from org.meteoinfo.data.meteodata.arl import DataLabel 
import os 
#---- Set directories 
dataDir = 'D:/Temp/test/hysplit-arl' 
#---- Set output data file 
outFile = os.path.join(dataDir, 'test7.arl') 
#---- Read a GrADS data file 
mydata = MeteoDataInfo() 
infile = os.path.join(dataDir, 'wrfout_d01_2006-07-13_00') 
print infile 
mydata.openNetCDFData(infile) 
print 'NetCDF file has been opened...' 
#---- Set output ARL data info 
arlDI = ARLDataInfo() 
#---- Set variable and level list 
gvars = ['PSFC','T2','U10','V10','T','U','V','W'] 
avars = ['PRSS','T02M','U10M','V10M','TEMP','UWND','VWND','WWND'] 
levels = [1, 0.9965000152587891, 0.9879999756813049, 0.9765000343322754, 0.9620000123977661, \ 
    0.9440000057220459, 0.9214999675750732, 0.8945000171661377, 0.8547878265380859, 0.8043636083602905, \ 
    0.7539393305778503, 0.7035150527954102, 0.6350234150886536, 0.5527187585830688, 0.478573739528656, \ 
    0.41191452741622925, 0.35211411118507385, 0.2985895872116089, 0.2507995367050171, 0.20824146270751953, \ 
    0.1704493761062622, 0.13699156045913696, 0.10746840387582779, 0.08151029050350189, 0.05877562612295151, \ 
    0.03894893825054169, 0.021738944575190544, 0.0068765804171562195] 
for l in levels: 
    arlDI.levels.add(l) 
    if l == 1: 
        arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M']) 
    else: 
        arlDI.LevelVarList.add(['TEMP','UWND','VWND','WWND']) 
#---- Write ARL data file 
dataInfo = mydata.getDataInfo() 
arlDI.createDataFile(outFile) 
arlDI.X = dataInfo.getXDimension().getValues() 
arlDI.Y = dataInfo.getYDimension().getValues() 
variables = dataInfo.getVariables() 
tNum = dataInfo.getTimeNum() 
for t in range(0, 1): 
    mydata.setTimeIndex(t) 
    atime = dataInfo.getTimes().get(t) 
    aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'GRPS', 1) 
    arlDI.writeIndexRecord(atime, aDH) 
    lidx = 0 
    for l in arlDI.levels: 
        print l 
        for v in arlDI.LevelVarList[lidx]: 
            vName = gvars[avars.index(v)] 
            print vName 
            gData = None 
            if lidx == 0: 
                mydata.setLevelIndex(lidx) 
                gData = mydata.getGridData(vName) 
            else: 
                if v == 'WWND': 
                    variable = dataInfo.getVariable('U') 
                    nlidx = variable.getZDimension().getDimValue().indexOf(l) 
                    mydata.setLevelIndex(nlidx) 
                    gData = mydata.getGridData(vName) 
                    mydata.setLevelIndex(nlidx + 1) 
                    gData1 = mydata.getGridData(vName) 
                    gData = gData.add(gData1).mul(0.5) 
                else: 
                    variable = dataInfo.getVariable(vName) 
                    nlidx = variable.getZDimension().getDimValue().indexOf(l) 
                    mydata.setLevelIndex(nlidx) 
                    gData = mydata.getGridData(vName) 
                
            if v == 'PRSS': 
                gData = gData.div(100) 
            elif v == 'UWND': 
                gData = gData.unStagger_X() 
            elif v == 'VWND': 
                gData = gData.unStagger_Y()            
            aDL = DataLabel(atime) 
            aDL.setLevel(lidx) 
            aDL.setVarName(v) 
            aDL.setGrid(99) 
            aDL.setForecast(0) 
            arlDI.writeGridData(aDL, gData) 
        lidx += 1 
arlDI.closeDataFile() 
print 'Finished!'

chenxl · 发表于 2015-6-30 14:57:04

老师，我想请教3个关于wrf-out数据错列的问题。1 wrf输出数据出现错列是由于什么导致的呢?可不可以通过wrf本身进行解决呢?
2 我另外做了一个wrf—out,这个也有错列的问题吗? 谢谢。http://pan.baidu.com/s/1bnAWd6B3 我通过NC数据转换为ARL数据做出了轨迹，但是与hysplit本身ARL数据做出的轨迹相差太大，是否是由于错列所致呢？

		自动登录	找回密码
密码			立即注册