Re: Use of mozbc utility with WACCM forecasts
845 views
Skip to first unread message
Gabriele Pfister
Jun 7, 2018, 11:34:42 AM6/7/18
to Yiannis Proestos, bar...@ucar.edu, Louisa Emmons, wrf-chem-mozbc
Hello Yiannis,
yes, the mozbc tool will work as is also with WACCM output. The main changes in the longer-lived gas phase speciation is that the lumped toluene has been split up into toluene, xylene and benzene. You will have to change this in the RADM2 or RACM mapping accordingly.
WACCM also uses a different aerosol scheme compared to MOZART-4 and the mapping MADE-SORGAM would also need to be changed. Unfortunately we cannot yet provide you with such a mapping. The only aerosol mapping we can provide at the moment is to the MOSAIC-8bin scheme. I am attaching an example mozbc namelist file for MOZART-MOSAIC-8bin, which hopefully gets you on the right path.
Gabriele
&control
do_bc = .true.
do_ic = .true.
domain = 1
dir_wrf = '/glade/scratch/pfister/WRFreal_WACCM/'
dir_moz = './'
fn_moz = 'test_WACMM_0001.nc'
moz_var_suffix = ''
def_missing_var = .true.
spc_map ='o3 -> O3', 'n2o -> N2O', 'no -> NO',
'no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',
'n2o5 -> N2O5', 'h2o2 -> H2O2',
'ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH',
'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',
'ald -> CH3CHO', 'acet -> CH3COCH3', 'mgly -> CH3COCHO',
'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',
'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8',
'c2h5oh -> C2H5OH', 'c10h16 -> MTERP',
'isopr -> ISOP',
'acetol -> HYAC',
'mek -> MEK',
'bigene -> BIGENE', 'bigalk -> BIGALK',
'tol -> TOLUENE',
'benzene -> BENZENE','xylenes ->XYLENES',
'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2',
'oc_a01->0.0093*pom_a1+0.7510*soa1_a2+0.0093*soa1_a1+0.7510*soa2_a2+0.0093*soa2_a1+0.7510*soa3_a2+0.0093*soa3_a1+0.7510*soa4_a2+0.0093*soa4_a1+0.7510*soa5_a2+0.0093*soa5_a1;1.e9',
'oc_a02->0.1123*pom_a1+0.2376*soa1_a2+0.1123*soa1_a1+0.2376*soa2_a2+0.1123*soa2_a1+0.2376*soa3_a2+0.1123*soa3_a1+0.2376*soa4_a2+0.1123*soa4_a1+0.2376*soa5_a2+0.1123*soa5_a1;1.e9',
'oc_a03->0.3835*pom_a1+0.0113*soa1_a2+0.3835*soa1_a1+0.0133*soa2_a2+0.3835*soa2_a1+0.0113*soa3_a2+0.3838*soa3_a1+0.0113*soa4_a2+0.3838*soa4_a1+0.0113*soa5_a2+0.3838*soa5_a1;1.e9',
'oc_a04->0.3783*pom_a1+0.0001*soa1_a2+0.3783*soa1_a1+0.0001*soa2_a2+0.3783*soa2_a1+0.0001*soa3_a2+0.3783*soa3_a1+0.0001*soa4_a2+0.3783*soa4_a1+0.0001*soa5_a2+0.3783*soa5_a1;1.e9',
'oc_a05->0.1077*pom_a1+0.0000*soa1_a2+0.1077*soa1_a1+0.0000*soa2_a2+0.1077*soa2_a1+0.0000*soa3_a2+0.1077*soa3_a1+0.0000*soa4_a2+0.1077*soa4_a1+0.0000*soa5_a2+0.1077*soa5_a1;1.e9',
'oc_a06->0.0087*pom_a1+0.0000*soa1_a2+0.0087*soa1_a1+0.0000*soa2_a2+0.0087*soa2_a1+0.0000*soa3_a2+0.0087*soa3_a1+0.0000*soa4_a2+0.0087*soa4_a1+0.0000*soa5_a2+0.0087*soa5_a1;1.e9',
'oc_a07->0.0002*pom_a1+0.0000*soa1_a2+0.0002*soa1_a1+0.0000*soa2_a2+0.0002*soa2_a1+0.0000*soa3_a2+0.0002*soa3_a1+0.0000*soa4_a2+0.0002*soa4_a1+0.0000*soa5_a2+0.0002*soa5_a1;1.e9',
'oc_a08->0.0000*pom_a1+0.0000*soa1_a2+0.0000*soa1_a1+0.0000*soa2_a2+0.0000*soa2_a1+0.0000*soa3_a2+0.0000*soa3_a1+0.0000*soa4_a2+0.0000*soa4_a1+0.0000*soa5_a2+0.0000*soa5_a1;1.e9',
'bc_a01->0.0093*bc_a1+0.0093*bc_a4;1.e9',
'bc_a02->0.1123*bc_a1+0.1123*bc_a4;1.e9',
'bc_a03->0.3835*bc_a1+0.3835*bc_a4;1.e9',
'bc_a04->0.3783*bc_a1+0.3783*bc_a4;1.e9',
'bc_a05->0.1077*bc_a1+0.1077*bc_a4;1.e9',
'bc_a06->0.0087*bc_a1+0.0087*bc_a4;1.e9',
'bc_a07->0.0002*bc_a1+0.0002*bc_a4;1.e9',
'bc_a08->0.0000*bc_a1+0.0000*bc_a4;1.e9',
'so4_a01->0.7510*so4_a2+0.0093*so4_a1+0.0000*so4_a3;1.e9',
'so4_a02->0.2376*so4_a2+0.1123*so4_a1+0.0000*so4_a3;1.e9',
'so4_a03->0.0113*so4_a2+0.3835*so4_a1+0.0000*so4_a3;1.e9',
'so4_a04->0.0001*so4_a2+0.3783*so4_a1+0.0002*so4_a3;1.e9',
'so4_a05->0.0000*so4_a2+0.1077*so4_a1+0.0061*so4_a3;1.e9',
'so4_a06->0.0000*so4_a2+0.0087*so4_a1+0.0934*so4_a3;1.e9',
'so4_a07->0.0000*so4_a2+0.0002*so4_a1+0.4020*so4_a3;1.e9',
'so4_a08->0.0000*so4_a2+0.0000*so4_a1+0.4983*so4_a3;1.e9',
'nh4_a01->0.1410*so4_a2+0.0033*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a02->0.0446*so4_a2+0.0017*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a03->0.0021*so4_a2+0.0210*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a04->0.0000*so4_a2+0.0720*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a05->0.0000*so4_a2+0.0202*so4_a1+0.0011*so4_a3;1.e9',
'nh4_a06->0.0000*so4_a2+0.0001*so4_a1+0.0175*so4_a3;1.e9',
'nh4_a07->0.0000*so4_a2+0.0000*so4_a1+0.0755*so4_a3;1.e9',
'nh4_a08->0.0000*so4_a2+0.0000*so4_a1+0.0935*so4_a3;1.e9',
'no3_a01->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a02->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a03->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a04->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a05->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a06->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a07->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'no3_a08->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
'na_a01->0.2954*ncl_a2+0.0037*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a02->0.0935*ncl_a2+0.0442*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a03->0.0045*ncl_a2+0.1509*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a04->0.0000*ncl_a2+0.1488*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a05->0.0000*ncl_a2+0.0424*ncl_a1+0.0024*ncl_a3;1.e9',
'na_a06->0.0000*ncl_a2+0.0034*ncl_a1+0.0367*ncl_a3;1.e9',
'na_a07->0.0000*ncl_a2+0.0000*ncl_a1+0.1582*ncl_a3;1.e9',
'na_a08->0.0000*ncl_a2+0.0000*ncl_a1+0.1960*ncl_a3;1.e9',
'cl_a01->0.4555*ncl_a2+0.0056*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a02->0.1441*ncl_a2+0.0681*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a03->0.0068*ncl_a2+0.2326*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a04->0.0000*ncl_a2+0.2295*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a05->0.0000*ncl_a2+0.0654*ncl_a1+0.0037*ncl_a3;1.e9',
'cl_a06->0.0000*ncl_a2+0.0055*ncl_a1+0.0567*ncl_a3;1.e9',
'cl_a07->0.0000*ncl_a2+0.0001*ncl_a1+0.2439*ncl_a3;1.e9',
'cl_a08->0.0000*ncl_a2+0.0000*ncl_a1+0.3023*ncl_a3;1.e9',
'oin_a01->0.7510*dst_a2+0.0093*dst_a1+0.0000*dst_a3;1.e9',
'oin_a02->0.2376*dst_a2+0.1123*dst_a1+0.0000*dst_a3;1.e9',
'oin_a03->0.0113*dst_a2+0.3835*dst_a1+0.0000*dst_a3;1.e9',
'oin_a04->0.0001*dst_a2+0.3783*dst_a1+0.0002*dst_a3;1.e9',
'oin_a05->0.0000*dst_a2+0.1077*dst_a1+0.0061*dst_a3;1.e9',
'oin_a06->0.0000*dst_a2+0.0087*dst_a1+0.0934*dst_a3;1.e9',
'oin_a07->0.0000*dst_a2+0.0002*dst_a1+0.4020*dst_a3;1.e9',
'oin_a08->0.0000*dst_a2+0.0000*dst_a1+0.4983*dst_a3;1.e9',
'num_a01->0.9502*num_a2+0.2509*num_a1+0.0000*num_a3;1.0',
'num_a02->0.0494*num_a2+0.4626*num_a1+0.0000*num_a3;1.0',
'num_a03->0.0004*num_a2+0.2470*num_a1+0.0007*num_a3;1.0',
'num_a04->0.0000*num_a2+0.0377*num_a1+0.0232*num_a3;1.0',
'num_a05->0.0000*num_a2+0.0016*num_a1+0.1886*num_a3;1.0',
'num_a06->0.0000*num_a2+0.0000*num_a1+0.4372*num_a3;1.0',
'num_a07->0.0000*num_a2+0.0000*num_a1+0.2935*num_a3;1.0',
'num_a08->0.0000*num_a2+0.0000*num_a1+0.0566*num_a3;1.0'
/
On Thu, Jun 7, 2018 at 4:15 AM, Yiannis Proestos <y.pro...@cyi.ac.cy> wrote:
Dear Dr. Barth and Dr. Pfister,
I have been using the mozbc utility to process chemical initial/boundary conditions for WRF-Chem simulations using the MOZART-4 forecasts. However, since the MOZART-4 is not operational any more and has been replaced by the WACCM system, I have a couple of questions regarding the processing of ic/bc from WACCM model output using mozbc.
1) Is the mozbc utility still applicable for processing the WACCM forecast output? At least, from a quick look of the mozbc source code I didn't notice any hard coded values that will limit its usage on the WACCM output.
2) Comparing the chemical species between MOZART-4 and WACC, there are some differences, so I was wondering whether there is any available mapping between lets say RADM2/MADE-SORGAM (or RACM/MADE-SORGAM) similar to the conversion examples you have kindly provided for MOZART-4 species (e.g., http://www.acom.ucar.edu/wrf-chem/ConversionMozartWRFChem.pdf).
Thanks in advance for any comments regarding my inquiries and many thanks for providing to the modelling community useful tools such as mozbc.
Best regards,
Yiannis
----------------------------------------------------------------------------------------------------------------
Yiannis Proestos, Ph.D. | Computational Scientist-Climate Modelling and Climate Change Impacts
The Cyprus Institute | Energy, Environment and Water Research Center (EEWRC)
20 Konstantinou Kavafi Street | 2121 Aglantzia, Nicosia, Cyprus
Phone: +357-22208634 | Skype: yproatcyi
y.pro...@cyi.ac.cy | https://www.cyi.ac.cy
----------------------------------------------------------------------------------------------------------------
================================
Gabriele Pfister, Deputy Director
Atmospheric Chemistry Observations & ModelingGabriele Pfister, Deputy Director
Liz Coleman
Sep 3, 2018, 6:07:11 AM9/3/18
to wrf-chem-mozbc, y.pro...@cyi.ac.cy, bar...@ucar.edu, emm...@ucar.edu
Hi all,
I'm just wondering if anyone has had any success running mozbc for WACCM forecast?
The WACCM data that I downloaded didn't have a complete chemical dataset (eg missing N2O for a start) and I wish to create a BC file for use with CRIMECH.
Any information would be most welcome,
Kind Regards,
Liz
Gabriele Pfister
Sep 4, 2018, 11:22:51 AM9/4/18
to Liz Coleman, wrf-chem-mozbc, Yiannis Proestos, Mary Barth, Louisa Emmons
Dear Liz,
sorry, N2O had not been part of the output but we will add it starting with tomorrow's forecast. Are there other species you are missing?
Best regards,
Gabriele
Gabriele Pfister
Sep 4, 2018, 11:44:12 AM9/4/18
to Liz Coleman, wrf-chem-mozbc, Yiannis Proestos, Mary Barth, Louisa Emmons
Liz,
Louisa just pointed out to me that while N2O is not in the 3-hour output, it is included in the 6-hour output that is provided at:
https://www.acom.ucar.edu/waccm/DATA/
ncdump -h f.e20.FWSD.f09_f09_mg17.beta08_misc05_cam5_4_170.forecast.001.cam.h3.2018-09-01-00000.nc |grep N2O
float N2O(time, lev, lat, lon) ;
N2O:mdims = 1 ;
N2O:units = "mol/mol" ;
N2O:long_name = "N2O concentration" ;
N2O:cell_methods = "time: mean" ;
Gabriele
Reply all
Reply to author
Forward
0 new messages