Hello everyone,
I am currently running VPRM on the basis of WRF-Chem V3.9.1.At present, I have successfully used pyVPRM to output the variables required for VPRM such as EVI, LSWI, etc., and write the relevant variables to vprm_input_d01. The variable information and dimensions are as follows :
dimensions:
Time = UNLIMITED ; // (32 currently)
vprm_vgcls = 8 ;
south_north = 99 ;
west_east = 112 ;
variables:
double EVI(Time, vprm_vgcls, south_north, west_east) ;
EVI:_FillValue = -9999. ;
EVI:FieldType = 104L ;
EVI:MemoryOrder = "XY" ;
EVI:coordinates = "XLONG XLAT" ;
EVI:stagger = "M" ;
double EVI_MAX(vprm_vgcls, south_north, west_east) ;
EVI_MAX:_FillValue = -9999. ;
EVI_MAX:FieldType = 104L ;
EVI_MAX:MemoryOrder = "XY" ;
EVI_MAX:coordinates = "XLONG XLAT" ;
EVI_MAX:stagger = "M" ;
double EVI_MIN(vprm_vgcls, south_north, west_east) ;
EVI_MIN:_FillValue = -9999. ;
EVI_MIN:FieldType = 104L ;
EVI_MIN:MemoryOrder = "XY" ;
EVI_MIN:coordinates = "XLONG XLAT" ;
EVI_MIN:stagger = "M" ;
double LSWI(Time, vprm_vgcls, south_north, west_east) ;
LSWI:_FillValue = -9999. ;
LSWI:FieldType = 104L ;
LSWI:MemoryOrder = "XY" ;
LSWI:coordinates = "XLONG XLAT" ;
LSWI:stagger = "M" ;
double LSWI_MAX(vprm_vgcls, south_north, west_east) ;
LSWI_MAX:_FillValue = -9999. ;
LSWI_MAX:FieldType = 104L ;
LSWI_MAX:MemoryOrder = "XY" ;
LSWI_MAX:coordinates = "XLONG XLAT" ;
LSWI_MAX:stagger = "M" ;
double LSWI_MIN(vprm_vgcls, south_north, west_east) ;
LSWI_MIN:_FillValue = -9999. ;
LSWI_MIN:FieldType = 104L ;
LSWI_MIN:MemoryOrder = "XY" ;
LSWI_MIN:coordinates = "XLONG XLAT" ;
LSWI_MIN:stagger = "M" ;
double VEGFRA_VPRM(vprm_vgcls, south_north, west_east) ;
VEGFRA_VPRM:_FillValue = -9999. ;
VEGFRA_VPRM:regrid_method = "conservative" ;
VEGFRA_VPRM:FieldType = 104L ;
VEGFRA_VPRM:MemoryOrder = "XY" ;
VEGFRA_VPRM:coordinates = "XLONG XLAT" ;
VEGFRA_VPRM:stagger = "M" ;
// global attributes:
:_NCProperties = "version=2,netcdf=4.8.1,hdf5=1.12.1" ;
:Conventions = "CF-1.8" ;
:Created_by = "Enhanced VPRM Processor" ;
:Attribute_Note = "All variables contain FieldType=104, MemoryOrder=XY" ;
In addition, I try to input EBIO_CO2OCE through wrfoce_d01, and the variable information and dimensions are as follows :
dimensions:
Time = UNLIMITED ; // (248 currently)
DateStrLen = 19 ;
south_north = 99 ;
west_east = 112 ;
variables:
double EBIO_CO2OCE(Time, south_north, west_east) ;
EBIO_CO2OCE:_FillValue = -1.e+34 ;
EBIO_CO2OCE:units = "mol km^-2 hr^-1" ;
EBIO_CO2OCE:description = "Ocean CO2 flux interpolated from CT2022 data" ;
char Times(Time, DateStrLen) ;
I specify three input files ( wrfchemi, wrfoce and vprm_input ) by namelist.input as follows ( the complete namelist.input is in the attachment ) :
io_form_auxinput5= 2
auxinput6_inname= "wrfoce_d01"
io_form_auxinput6= 2
auxinput15_inname= "vprm_input_d01"
io_form_auxinput15= 2
However, my WRF-Chem (VPRM) seems to only read wrfchemi (anthropogenic emissions), which makes CO2_ANT change, while wrfoce_d01 and vprm_input_d01 are not read, and CO2_OCE and CO2_BIO do not change accordingly. How can I solve this problem ? Is there a problem with the setting of my namelist or is the format of my input file incorrect ?
Hope to get everyone 's help and guidance, thanks!
I am currently running VPRM on the basis of WRF-Chem V3.9.1.At present, I have successfully used pyVPRM to output the variables required for VPRM such as EVI, LSWI, etc., and write the relevant variables to vprm_input_d01. The variable information and dimensions are as follows :
dimensions:
Time = UNLIMITED ; // (32 currently)
vprm_vgcls = 8 ;
south_north = 99 ;
west_east = 112 ;
variables:
double EVI(Time, vprm_vgcls, south_north, west_east) ;
EVI:_FillValue = -9999. ;
EVI:FieldType = 104L ;
EVI:MemoryOrder = "XY" ;
EVI:coordinates = "XLONG XLAT" ;
EVI:stagger = "M" ;
double EVI_MAX(vprm_vgcls, south_north, west_east) ;
EVI_MAX:_FillValue = -9999. ;
EVI_MAX:FieldType = 104L ;
EVI_MAX:MemoryOrder = "XY" ;
EVI_MAX:coordinates = "XLONG XLAT" ;
EVI_MAX:stagger = "M" ;
double EVI_MIN(vprm_vgcls, south_north, west_east) ;
EVI_MIN:_FillValue = -9999. ;
EVI_MIN:FieldType = 104L ;
EVI_MIN:MemoryOrder = "XY" ;
EVI_MIN:coordinates = "XLONG XLAT" ;
EVI_MIN:stagger = "M" ;
double LSWI(Time, vprm_vgcls, south_north, west_east) ;
LSWI:_FillValue = -9999. ;
LSWI:FieldType = 104L ;
LSWI:MemoryOrder = "XY" ;
LSWI:coordinates = "XLONG XLAT" ;
LSWI:stagger = "M" ;
double LSWI_MAX(vprm_vgcls, south_north, west_east) ;
LSWI_MAX:_FillValue = -9999. ;
LSWI_MAX:FieldType = 104L ;
LSWI_MAX:MemoryOrder = "XY" ;
LSWI_MAX:coordinates = "XLONG XLAT" ;
LSWI_MAX:stagger = "M" ;
double LSWI_MIN(vprm_vgcls, south_north, west_east) ;
LSWI_MIN:_FillValue = -9999. ;
LSWI_MIN:FieldType = 104L ;
LSWI_MIN:MemoryOrder = "XY" ;
LSWI_MIN:coordinates = "XLONG XLAT" ;
LSWI_MIN:stagger = "M" ;
double VEGFRA_VPRM(vprm_vgcls, south_north, west_east) ;
VEGFRA_VPRM:_FillValue = -9999. ;
VEGFRA_VPRM:regrid_method = "conservative" ;
VEGFRA_VPRM:FieldType = 104L ;
VEGFRA_VPRM:MemoryOrder = "XY" ;
VEGFRA_VPRM:coordinates = "XLONG XLAT" ;
VEGFRA_VPRM:stagger = "M" ;
// global attributes:
:_NCProperties = "version=2,netcdf=4.8.1,hdf5=1.12.1" ;
:Conventions = "CF-1.8" ;
:Created_by = "Enhanced VPRM Processor" ;
:Attribute_Note = "All variables contain FieldType=104, MemoryOrder=XY" ;
In addition, I try to input EBIO_CO2OCE through wrfoce_d01, and the variable information and dimensions are as follows :
dimensions:
Time = UNLIMITED ; // (248 currently)
DateStrLen = 19 ;
south_north = 99 ;
west_east = 112 ;
variables:
double EBIO_CO2OCE(Time, south_north, west_east) ;
EBIO_CO2OCE:_FillValue = -1.e+34 ;
EBIO_CO2OCE:units = "mol km^-2 hr^-1" ;
EBIO_CO2OCE:description = "Ocean CO2 flux interpolated from CT2022 data" ;
char Times(Time, DateStrLen) ;
I specify three input files ( wrfchemi, wrfoce and vprm_input ) by namelist.input as follows ( the complete namelist.input is in the attachment ) :
io_form_auxinput5= 2
auxinput6_inname= "wrfoce_d01"
io_form_auxinput6= 2
auxinput15_inname= "vprm_input_d01"
io_form_auxinput15= 2
However, my WRF-Chem (VPRM) seems to only read wrfchemi (anthropogenic emissions), which makes CO2_ANT change, while wrfoce_d01 and vprm_input_d01 are not read, and CO2_OCE and CO2_BIO do not change accordingly. How can I solve this problem ? Is there a problem with the setting of my namelist or is the format of my input file incorrect ?
Hope to get everyone 's help and guidance, thanks!