Page 1 of 1

mozbc chemical boundary conditions do not seem to be applied

PostPosted: Tue Jul 29, 2014 5:11 pm
by jwsmit12
Hello:

The mozbc utility developed by NCAR implements MOZART GTCM chemistry and aerosol boundary and initial conditions. My model domain is over 40 W-25E and 20N-20S. I have attached a picture of WRF-Chem carbon monoxide output which should unrealistically low mixing ratios of CO on the extremities of the model domain.

I have emailed Stacy Walters at NCAR but wondering if anyone was familiar with this issue and how to resolve it.

I am using the chem_opt=301 which is GOCART aerosols and RACM chemistry.
In my WRF-Chem namelist have_bcs_chem = 1. I am thinking that this indicates bc's are being applied. Are there any other chem namelist parameters I need to pay attention to.

My mozbc input for GOCART is the following:
&control

do_bc = .true.
do_ic = .true.
domain = 1
dir_wrf = '/scratch1/jwsmith/get_int_wps_testjws/'
dir_moz = '/data/jwsmith/MOZART/'
fn_moz = 'mozart4geos5-20140707131233409239.nc'

spc_map = 'BC1 -> .4143*CB1;1.e9', 'BC2 -> .4143*CB2;1.e9',
'OC1 -> .4143*OC1;1.e9', 'OC2 -> .4143*OC2;1.e9',
'DUST_1 -> 1.1738*[DUST1];1.e9', 'DUST_2 -> .939*[DUST2];1.e9',
'DUST_3 -> .2348*[DUST2]+.939*[DUST3];1.e9',
'DUST_4 -> .2348*[DUST3]+.5869*[DUST4];1.e9', 'DUST_5 -> .5869*[DUST4];1.e9',
'SEAS_1 -> 2.*SA1;1.e9', 'SEAS_2 -> 2.*SA2;1.e9',
'SEAS_3 -> 2.*SA3;1.e9', 'SEAS_4 -> 2.*SA4;1.e9'


For RACM:

&control

do_bc = .true.
do_ic = .true.
domain = 1

dir_wrf = '/scratch1/jwsmith/get_int_wps_testjws/'
dir_moz = '/data/jwsmith/MOZART/'
fn_moz = 'mozart4geos5-20140707131233409239.nc'

spc_map = 'co -> CO',
'o3 -> O3',
'so2 -> SO2',
'no -> NO',
'no2 -> NO2',
'h2o2 -> H2O2',
'no3 -> NO3',
'n2o5 -> N2O5',
'hno3 -> HNO3',
'nh3 -> NH3',
'hcho -> CH2O',
'pan -> PAN',
'tpan -> MPAN',
'tol -> .75*TOLUENE',
'xyl -> .25*TOLUENE',
'hno4 -> HO2NO2',
'eth -> C2H6',
'hc3 -> C3H8',
'hc5 -> BIGALK',
'iso -> ISOP',
'macr -> MACR',
'ald -> GLYALD+CH3CHO',
'ket -> CH3COCH3',
'mgly -> CH3COCHO',
'hket -> HYAC',
'csl -> CRESOL',
'ete -> C2H4',
'olt -> C3H6+MVK+0.5*BIGENE',
'oli -> .5*BIGENE',
'paa -> CH3COOOH',
'onit -> ONIT+ONITR',
'ora2 -> CH3COOH',
'op1 -> CH3OOH',
'ch4 -> CH4',
'sulf -> SO4;0.30173e6'
/


Please let me know if you all have found solutions.

Re: mozbc chemical boundary conditions do not seem to be app

PostPosted: Mon May 04, 2015 1:10 pm
by Schuch
Do you run './mozbc < GOCART.inp' and then './mozbc < RACM.inp'?

Try this:
1. save a RACM_GOCART.inp in /data/jwsmith/MOZART with

Code: Select all
&control

do_bc = .true.
do_ic = .true.
domain = 1

dir_wrf = '/scratch1/jwsmith/get_int_wps_testjws/'
dir_moz = '/data/jwsmith/MOZART/'
fn_moz = 'mozart4geos5-20140707131233409239.nc'

spc_map = 'co -> CO',
'o3 -> O3',
'so2 -> SO2',
'no -> NO',
'no2 -> NO2',
'h2o2 -> H2O2',
'no3 -> NO3',
'n2o5 -> N2O5',
'hno3 -> HNO3',
'nh3 -> NH3',
'hcho -> CH2O',
'pan -> PAN',
'tpan -> MPAN',
'tol -> .75*TOLUENE',
'xyl -> .25*TOLUENE',
'hno4 -> HO2NO2',
'eth -> C2H6',
'hc3 -> C3H8',
'hc5 -> BIGALK',
'iso -> ISOP',
'macr -> MACR',
'ald -> GLYALD+CH3CHO',
'ket -> CH3COCH3',
'mgly -> CH3COCHO',
'hket -> HYAC',
'csl -> CRESOL',
'ete -> C2H4',
'olt -> C3H6+MVK+0.5*BIGENE',
'oli -> .5*BIGENE',
'paa -> CH3COOOH',
'onit -> ONIT+ONITR',
'ora2 -> CH3COOH',
'op1 -> CH3OOH',
'ch4 -> CH4',
'sulf -> SO4;0.30173e6',
'BC1 -> .4143*CB1;1.e9',
'BC2 -> .4143*CB2;1.e9',
'OC1 -> .4143*OC1;1.e9',
'OC2 -> .4143*OC2;1.e9',
'DUST_1 -> 1.1738*[DUST1];1.e9',
'DUST_2 -> .939*[DUST2];1.e9',
'DUST_3 -> .2348*[DUST2]+.939*[DUST3];1.e9',
'DUST_4 -> .2348*[DUST3]+.5869*[DUST4];1.e9',
'DUST_5 -> .5869*[DUST4];1.e9',
'SEAS_1 -> 2.*SA1;1.e9',
'SEAS_2 -> 2.*SA2;1.e9',
'SEAS_3 -> 2.*SA3;1.e9',
'SEAS_4 -> 2.*SA4;1.e9'
/


2. run './mozbc < RACM_GOCART.inp > log.txt', in case of error this log.txt can be usefull