mozbc chemical boundary conditions do not seem to be applied

Ideas, problems or suggestions regarding the initialization of WRF Chemistry

mozbc chemical boundary conditions do not seem to be applied

Postby jwsmit12 » Tue Jul 29, 2014 5:11 pm

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.
jwsmit12
 
Posts: 8
Joined: Mon Jul 28, 2014 6:18 pm

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

Postby Schuch » Mon May 04, 2015 1:10 pm

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
Schuch
 
Posts: 15
Joined: Mon Mar 30, 2015 7:19 pm


Return to Initialization

Who is online

Users browsing this forum: No registered users and 1 guest

cron