Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

There are some input problems about WRF-Chem(VPRM,chem_opt=16).

Steven

New member
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" ;
......(others)
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" ;

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!
 

Attachments

  • namelist.input
    7.4 KB · Views: 0
Top