Dear @cjcoats, thank you for your reply.
this is the ncdump -h output:
ncdump -h GRIDBDY2D_d01.nc
netcdf GRIDBDY2D_d01 {
dimensions:
TSTEP = 1 ;
DATE-TIME = 2 ;
LAY = 1 ;
VAR = 28 ;
PERIM = 104 ;
variables:
int TFLAG(TSTEP, VAR, DATE-TIME) ;
TFLAG:units = "<YYYYDDD,HHMMSS>" ;
TFLAG:long_name = "TFLAG " ;
TFLAG:var_desc = "Timestep-valid flags: (1) YYYYDDD or (2) HHMMSS " ;
float LAT(TSTEP, LAY, PERIM) ;
LAT:long_name = "LAT " ;
LAT:units = "DEGREES " ;
LAT:var_desc = "latitude (south negative) " ;
float LON(TSTEP, LAY, PERIM) ;
LON:long_name = "LON " ;
LON:units = "DEGREES " ;
LON:var_desc = "longitude (west negative) " ;
float MSFX2(TSTEP, LAY, PERIM) ;
MSFX2:long_name = "MSFX2 " ;
MSFX2:units = "(M/M)**2 " ;
MSFX2:var_desc = "squared map-scale factor (CROSS) " ;
float HT(TSTEP, LAY, PERIM) ;
HT:long_name = "HT " ;
HT:units = "M " ;
HT:var_desc = "terrain elevation " ;
float DLUSE(TSTEP, LAY, PERIM) ;
DLUSE:long_name = "DLUSE " ;
DLUSE:units = "CATEGORY " ;
DLUSE:var_desc = "dominant land use category from classification MODIS NOAH " ;
float LWMASK(TSTEP, LAY, PERIM) ;
LWMASK:long_name = "LWMASK " ;
LWMASK:units = "CATEGORY " ;
LWMASK:var_desc = "land-water mask (1=land, 0=water) " ;
float PURB(TSTEP, LAY, PERIM) ;
PURB:long_name = "PURB " ;
PURB:units = "PERCENT " ;
PURB:var_desc = "urban percent of cell based on land " ;
float LUFRAC_01(TSTEP, LAY, PERIM) ;
LUFRAC_01:long_name = "LUFRAC_01 " ;
LUFRAC_01:units = "FRACTION " ;
LUFRAC_01:var_desc = "MODIS NOAH: Evergreen Needleleaf Forest " ;
float LUFRAC_02(TSTEP, LAY, PERIM) ;
LUFRAC_02:long_name = "LUFRAC_02 " ;
LUFRAC_02:units = "FRACTION " ;
LUFRAC_02:var_desc = "MODIS NOAH: Evergreen Broadleaf Forest " ;
float LUFRAC_03(TSTEP, LAY, PERIM) ;
LUFRAC_03:long_name = "LUFRAC_03 " ;
LUFRAC_03:units = "FRACTION " ;
LUFRAC_03:var_desc = "MODIS NOAH: Deciduous Needleleaf Forest " ;
float LUFRAC_04(TSTEP, LAY, PERIM) ;
LUFRAC_04:long_name = "LUFRAC_04 " ;
LUFRAC_04:units = "FRACTION " ;
LUFRAC_04:var_desc = "MODIS NOAH: Deciduous Broadleaf Forest " ;
float LUFRAC_05(TSTEP, LAY, PERIM) ;
LUFRAC_05:long_name = "LUFRAC_05 " ;
LUFRAC_05:units = "FRACTION " ;
LUFRAC_05:var_desc = "MODIS NOAH: Mixed Forests " ;
float LUFRAC_06(TSTEP, LAY, PERIM) ;
LUFRAC_06:long_name = "LUFRAC_06 " ;
LUFRAC_06:units = "FRACTION " ;
LUFRAC_06:var_desc = "MODIS NOAH: Closed Shrublands " ;
float LUFRAC_07(TSTEP, LAY, PERIM) ;
LUFRAC_07:long_name = "LUFRAC_07 " ;
LUFRAC_07:units = "FRACTION " ;
LUFRAC_07:var_desc = "MODIS NOAH: Open Shrublands " ;
float LUFRAC_08(TSTEP, LAY, PERIM) ;
LUFRAC_08:long_name = "LUFRAC_08 " ;
LUFRAC_08:units = "FRACTION " ;
LUFRAC_08:var_desc = "MODIS NOAH: Woody Savannas " ;
float LUFRAC_09(TSTEP, LAY, PERIM) ;
LUFRAC_09:long_name = "LUFRAC_09 " ;
LUFRAC_09:units = "FRACTION " ;
LUFRAC_09:var_desc = "MODIS NOAH: Savannas " ;
float LUFRAC_10(TSTEP, LAY, PERIM) ;
LUFRAC_10:long_name = "LUFRAC_10 " ;
LUFRAC_10:units = "FRACTION " ;
LUFRAC_10:var_desc = "MODIS NOAH: Grasslands " ;
float LUFRAC_11(TSTEP, LAY, PERIM) ;
LUFRAC_11:long_name = "LUFRAC_11 " ;
LUFRAC_11:units = "FRACTION " ;
LUFRAC_11:var_desc = "MODIS NOAH: Permanent Wetlands " ;
float LUFRAC_12(TSTEP, LAY, PERIM) ;
LUFRAC_12:long_name = "LUFRAC_12 " ;
LUFRAC_12:units = "FRACTION " ;
LUFRAC_12:var_desc = "MODIS NOAH: Croplands " ;
float LUFRAC_13(TSTEP, LAY, PERIM) ;
LUFRAC_13:long_name = "LUFRAC_13 " ;
LUFRAC_13:units = "FRACTION " ;
LUFRAC_13:var_desc = "MODIS NOAH: Urban and Built-Up " ;
float LUFRAC_14(TSTEP, LAY, PERIM) ;
LUFRAC_14:long_name = "LUFRAC_14 " ;
LUFRAC_14:units = "FRACTION " ;
LUFRAC_14:var_desc = "MODIS NOAH: Cropland-Natural Vegetation Mosaic " ;
float LUFRAC_15(TSTEP, LAY, PERIM) ;
LUFRAC_15:long_name = "LUFRAC_15 " ;
LUFRAC_15:units = "FRACTION " ;
LUFRAC_15:var_desc = "MODIS NOAH: Snow and Ice " ;
float LUFRAC_16(TSTEP, LAY, PERIM) ;
LUFRAC_16:long_name = "LUFRAC_16 " ;
LUFRAC_16:units = "FRACTION " ;
LUFRAC_16:var_desc = "MODIS NOAH: Barren or Sparsely Vegetated " ;
float LUFRAC_17(TSTEP, LAY, PERIM) ;
LUFRAC_17:long_name = "LUFRAC_17 " ;
LUFRAC_17:units = "FRACTION " ;
LUFRAC_17:var_desc = "MODIS NOAH: Water " ;
float LUFRAC_18(TSTEP, LAY, PERIM) ;
LUFRAC_18:long_name = "LUFRAC_18 " ;
LUFRAC_18:units = "FRACTION " ;
LUFRAC_18:var_desc = "MODIS NOAH: Wooded Tundra " ;
float LUFRAC_19(TSTEP, LAY, PERIM) ;
LUFRAC_19:long_name = "LUFRAC_19 " ;
LUFRAC_19:units = "FRACTION " ;
LUFRAC_19:var_desc = "MODIS NOAH: Mixed Tundra " ;
float LUFRAC_20(TSTEP, LAY, PERIM) ;
LUFRAC_20:long_name = "LUFRAC_20 " ;
LUFRAC_20:units = "FRACTION " ;
LUFRAC_20:var_desc = "MODIS NOAH: Barren Tundra " ;
float LUFRAC_21(TSTEP, LAY, PERIM) ;
LUFRAC_21:long_name = "LUFRAC_21 " ;
LUFRAC_21:units = "FRACTION " ;
LUFRAC_21:var_desc = "MODIS NOAH: Lake (if from WRFv3.8 or later) " ;
// global attributes:
:IOAPI_VERSION = "ioapi-3.2: $Id: init3.F90 98 2018-04-05 14:35:07Z coats $ " ;
:EXEC_ID = "mcip " ;
:FTYPE = 2 ;
:CDATE = 2019182 ;
:CTIME = 175644 ;
:WDATE = 2019182 ;
:WTIME = 175644 ;
:SDATE = 2015266 ;
:STIME = 0 ;
:TSTEP = 0 ;
:NTHIK = 1 ;
:NCOLS = 25 ;
:NROWS = 25 ;
:NLAYS = 1 ;
:NVARS = 28 ;
:GDTYP = 2 ;
:P_ALP = -26.5 ;
:P_BET = -13.5 ;
:P_GAM = -45.9000015258789 ;
:XCENT = -45.9000015258789 ;
:YCENT = -23.550500869751 ;
:XORIG = -198000. ;
:YORIG = -132500. ;
:XCELL = 9000. ;
:YCELL = 9000. ;
:VGTYP = 7 ;
:VGTOP = 5000.f ;
:VGLVLS = 1.f, 0.9983f ;
:GDNAM = "GRIDOUT_BR_9km_M" ;
:UPNAM = "GRIDOUT " ;
:VAR-LIST = "LAT LON MSFX2 HT DLUSE LWMASK PURB LUFRAC_01 LUFRAC_02 LUFRAC_03 LUFRAC_04 LUFRAC_05 LUFRAC_06 LUFRAC_07 LUFRAC_08 LUFRAC_09 LUFRAC_10 LUFRAC_11 LUFRAC_12 LUFRAC_13 LUFRAC_14 LUFRAC_15 LUFRAC_16 LUFRAC_17 LUFRAC_18 LUFRAC_19 LUFRAC_20 LUFRAC_21 " ;
:FILEDESC = "US EPA COMMUNITY MULTISCALE AIR QUALITY MODEL METEOROLOGY-CHEMISTRY INTERFACE PROCESSOR MCIP V4.5 FROZEN 03/16/2018 INPUT METEOROLOGY DATA FROM WRF ARW V4.1 INPUT RUN INITIALIZED: 2015-07-22-00:00:00.0000 CUMULUS PARAMETERIZATION: Multi-Scale Kain-Fritsch SHALLOW CONVECTION: No shallow convection MICROPHYSICS: Morrison 2-moment LONGWAVE RADIATION: RRTMg SHORTWAVE RADIATION: Goddard PBL SCHEME: YSU SURFACE LAYER SCHEME: Revised MM5 (Jimenez) LAND-SURFACE SCHEME: NOAH Land-Surface Model URBAN MODEL: No urban physics LAND USE CLASSIFICATION: MODIFIED_IGBP_MODIS 3D ANALYSIS NUDGING: GRID WIND COEFF: 6.000E-04 s-1 TEMP COEFF: 3.000E-04 s-1 MOIS COEFF: 1.000E-05 s-1 GEOP COEFF: not applicable SFC ANALYSIS NUDGING: OFF WIND COEFF: not applicable TEMP COEFF: not applicable MOIS COEFF: not applicable OBS NUDGING: OFF WIND COEFF: not applicable TEMP COEFF: not applicable MOIS COEFF: not applicable EARTH RADIUS ASSUMED IN MCIP: 6370000.000 m " ;
:HISTORY = "" ;
data:
TFLAG =
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0,
0, 0 ;
}
Also this is the ncdump -h of my WRF file:
ncdump -h wrfout_d01_2015-09-01_September_W014
netcdf wrfout_d01_2015-09-01_September_W014 {
dimensions:
Time = UNLIMITED ; // (726 currently)
DateStrLen = 19 ;
west_east = 190 ;
south_north = 160 ;
bottom_top = 34 ;
bottom_top_stag = 35 ;
soil_layers_stag = 4 ;
west_east_stag = 191 ;
south_north_stag = 161 ;
seed_dim_stag = 2 ;
variables:
// global attributes:
:TITLE = " OUTPUT FROM WRF V4.1 MODEL" ;
:START_DATE = "2015-09-01_00:00:00" ;
:SIMULATION_START_DATE = "2015-07-22_00:00:00" ;
:WEST-EAST_GRID_DIMENSION = 191 ;
:SOUTH-NORTH_GRID_DIMENSION = 161 ;
:BOTTOM-TOP_GRID_DIMENSION = 35 ;
:DX = 9000.f ;
:DY = 9000.f ;
:AERCU_OPT = 0 ;
:AERCU_FCT = 1.f ;
:IDEAL_CASE = 0 ;
:DIFF_6TH_SLOPEOPT = 0 ;
:AUTO_LEVELS_OPT = 2 ;
:DIFF_6TH_THRESH = 0.1f ;
:DZBOT = 50.f ;
:DZSTRETCH_S = 1.3f ;
:DZSTRETCH_U = 1.1f ;
:SKEBS_ON = 0 ;
:SPEC_BDY_FINAL_MU = 1 ;
:USE_Q_DIABATIC = 0 ;
:GRIDTYPE = "C" ;
:DIFF_OPT = 1 ;
:KM_OPT = 4 ;
:DAMP_OPT = 0 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 10 ;
:RA_LW_PHYSICS = 4 ;
:RA_SW_PHYSICS = 2 ;
:SF_SFCLAY_PHYSICS = 1 ;
:SF_SURFACE_PHYSICS = 2 ;
:BL_PBL_PHYSICS = 1 ;
:CU_PHYSICS = 11 ;
:SF_LAKE_PHYSICS = 0 ;
:SURFACE_INPUT_SOURCE = 3 ;
:SST_UPDATE = 0 ;
:GRID_FDDA = 1 ;
:GFDDA_INTERVAL_M = 180 ;
:GFDDA_END_H = 99999999 ;
:GRID_SFDDA = 2 ;
:SGFDDA_INTERVAL_M = 180 ;
:SGFDDA_END_H = 99999999 ;
:HYPSOMETRIC_OPT = 2 ;
:USE_THETA_M = 1 ;
:GWD_OPT = 0 ;
:SF_URBAN_PHYSICS = 0 ;
:SF_SURFACE_MOSAIC = 0 ;
:SF_OCEAN_PHYSICS = 0 ;
:SHCU_PHYSICS = 0 ;
:MFSHCONV = 0 ;
:FEEDBACK = 1 ;
:SMOOTH_OPTION = 0 ;
:SWRAD_SCAT = 1.f ;
:W_DAMPING = 0 ;
:DT = 25.f ;
:RADT = 5.f ;
:BLDT = 0.f ;
:CUDT = 0.f ;
:AER_OPT = 0 ;
:SWINT_OPT = 0 ;
:AER_TYPE = 1 ;
:AER_AOD550_OPT = 1 ;
:AER_ANGEXP_OPT = 1 ;
:AER_SSA_OPT = 1 ;
:AER_ASY_OPT = 1 ;
:AER_AOD550_VAL = 0.12f ;
:AER_ANGEXP_VAL = 1.3f ;
:AER_SSA_VAL = 0.85f ;
:AER_ASY_VAL = 0.9f ;
:MOIST_ADV_OPT = 1 ;
:SCALAR_ADV_OPT = 1 ;
:TKE_ADV_OPT = 1 ;
:DIFF_6TH_OPT = 0 ;
:DIFF_6TH_FACTOR = 0.12f ;
:FGDT = 0.f ;
:GUV = 0.0006f ;
:GT = 0.0003f ;
:GQ = 1.e-05f ;
:IF_RAMPING = 1 ;
:DTRAMP_MIN = 60.f ;
:GUV_SFC = 0.0006f ;
:GT_SFC = 0.0003f ;
:GQ_SFC = 1.e-05f ;
:RINBLW = 5000.f ;
:OBS_NUDGE_OPT = 0 ;
:BUCKET_MM = -1.f ;
:BUCKET_J = -1.f ;
:PREC_ACC_DT = 0.f ;
:ISFTCFLX = 0 ;
:ISHALLOW = 0 ;
:ISFFLX = 1 ;
:ICLOUD = 1 ;
:ICLOUD_CU = 2 ;
:TRACER_PBLMIX = 1 ;
:SCALAR_PBLMIX = 0 ;
:YSU_TOPDOWN_PBLMIX = 0 ;
:GRAV_SETTLING = 0 ;
:DFI_OPT = 0 ;
:SIMULATION_INITIALIZATION_TYPE = "REAL-DATA CASE" ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 190 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 191 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 160 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 161 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 34 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 35 ;
:GRID_ID = 1 ;
:PARENT_ID = 1 ;
:I_PARENT_START = 1 ;
:J_PARENT_START = 1 ;
:PARENT_GRID_RATIO = 1 ;
:CEN_LAT = -19.70001f ;
:CEN_LON = -45.90039f ;
:TRUELAT1 = -13.5f ;
:TRUELAT2 = -26.5f ;
:MOAD_CEN_LAT = -19.70001f ;
:STAND_LON = -45.9f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 0.f ;
:JULYR = 2015 ;
:JULDAY = 203 ;
:MAP_PROJ = 1 ;
:MAP_PROJ_CHAR = "Lambert Conformal" ;
:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
:NUM_LAND_CAT = 21 ;
:ISWATER = 17 ;
:ISLAKE = 21 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:HYBRID_OPT = 0 ;
:ETAC = 0.f ;
}
In my past MCIP cases, I have never seen the SOI_CRO and LUFRAC_CRO being written in a separate file. Is this normal?
According to the "ckfile3.f " code, my VGTYP3 is being set to -999, but in the header of WRF and MCIP files, the VGTYP3 =7, which is supported by IOAPI.
IF ( VGTYP3(FID) .EQ. IMISS3 ) THEN ! "other" -- legal but unusual
WRITE( MESG, 94010 )
& 'WARNING: Vertical grid/coordinate type:',
& VGTYP3(FID),
& '"MISSING" in file "' //
& TRIM( FLIST3( FID ) ) // '"'
CALL M3WARN( 'CKFILE3', 0, 0, MESG )
ELSE IF ( ( VGTYP3(FID) .NE. VGSGPN3 ) .AND.
& ( VGTYP3(FID) .NE. VGSGPH3 ) .AND.
& ( VGTYP3(FID) .NE. VGSIGZ3 ) .AND.
& ( VGTYP3(FID) .NE. VGPRES3 ) .AND.
& ( VGTYP3(FID) .NE. VGZVAL3 ) .AND.
& ( VGTYP3(FID) .NE. VGHVAL3 ) .AND.
& ( VGTYP3(FID) .NE. VGWRFEM ) .AND.
& ( VGTYP3(FID) .NE. VGWRFNM ) ) THEN
What do you think is causing this problem?
Thank you so much for your time and attention.
Rizzieri