Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

How to run ozone chemistry with MAM3 physics in CESM1.2

C

chen63@illinois_edu

Member
Which compset should I use if I want to run ozone chemistry with MAM3 physiscs for aerosols?Another question is how I can run prescribed aerosols with MAM3 physics in CESM1.2? I saw that this function is now available.Thanks!
 
hannay

hannay

Cecile Hannay
AMWG Liaison
Staff member
About MAM prescribed, see:
https://bb.cgd.ucar.edu/node/1002386#comment-1006283
 
T

tilmes

Member
there are two compsets in CESM1.2 that allow you to run trop-mozart chemistry: B_2000_MOZMAM_CN, or F_2000_MOZMAM_CN, depending on what you are interested in. There are also other compsets if you want to run with both tropospheric and stratospheric chemistry.Simone
 
C

chen63@illinois_edu

Member
I have several new questions regarding the compset F_2000_MOZMAM_CN for ozone chemistry. Could someone please help me on these questions? Thanks a lot!(1) I just notice that co2vmr in the compset F_2000_MOZMAM_CN is set to 0.000001e-6  by default. Could anyone please give some hint what it is set this way? Should I change it? (2) Another question is about the grid settings of compset F_2000_MOZMAM_CN. Can it only be run under 1.9*2.5, or can it be run for any other finer resolutions? Do I need to prepare my own emission files for other resolutions as I noticed that the current emission files in the input data directory for srf_emis_specifier   are only for 1.9*2.5?(3) One more question is about the gas species in this compset. According to Emmons et al., 2010, there are 85 gaseous species included in MOZART 4. However, when I looked at the atm_in, the gaseous species included are specified in srf_emis_specifier. Are these all the gaseous species included? How about the reaction envolving other gaseous species?srf_emis_specifier             = 'BIGALK   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.BIGALK.surface.1.9x2.5_c110426.nc',         'BIGENE   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.BIGENE.surface.1.9x2.5_c110426.nc',         'C10H16   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C10H16.surface.1.9x2.5_c110426.nc',         'C2H2     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C2H2.surface.1.9x2.5_c110426.nc',         'C2H4     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C2H4.surface.1.9x2.5_c110426.nc',         'C2H5OH   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C2H5OH.surface.1.9x2.5_c110426.nc',         'C2H6     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C2H6.surface.1.9x2.5_c110426.nc',         'C3H6     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C3H6.surface.1.9x2.5_c110426.nc',         'C3H8     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.C3H8.surface.1.9x2.5_c110426.nc',         'CH2O     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CH2O.surface.1.9x2.5_c110426.nc',         'CH3CHO   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CH3CHO.surface.1.9x2.5_c110426.nc',         'CH3CN    -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CH3CN.surface.1.9x2.5_c110426.nc',         'CH3COCH3 -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CH3COCH3.surface.1.9x2.5_c110426.nc',         'CH3COOH  -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CH3COOH.surface.1.9x2.5_c110426.nc',         'CH3OH    -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CH3OH.surface.1.9x2.5_c110426.nc',         'CO       -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.CO.surface.1.9x2.5_c110426.nc',         'DMS      -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/aerocom_mam3_dms_surf_2000_c120315.nc',         'HCN      -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.HCN.surface.1.9x2.5_c110426.nc',         'HCOOH    -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.HCOOH.surface.1.9x2.5_c110426.nc',         'ISOP     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.ISOP.surface.1.9x2.5_c110426.nc',         'MEK      -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.MEK.surface.1.9x2.5_c110426.nc',         'NH3      -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/emis_NH3_2000_c111014.nc',         'NO       -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.NO.surface.1.9x2.5_c110426.nc',         'SO2      -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so2_surf_2000_c120315.nc',         'SOAG     -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_soag_1.5_surf_2000_c120315.nc',         'TOLUENE  -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/emis/1992-2010/emissions.TOLUENE.surface.1.9x2.5_c110426.nc',         'bc_a1    -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_bc_surf_2000_c120315.nc',         'num_a1   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a1_surf_2000_c120315.nc',         'num_a2   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_num_a2_surf_2000_c120315.nc',         'pom_a1   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_oc_surf_2000_c120315.nc',         'so4_a1   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a1_surf_2000_c120315.nc',         'so4_a2   -> /project/projectdirs/ccsm1/inputdata/atm/cam/chem/trop_mozart_aero/emis/ar5_mam3_so4_a2_surf_2000_c120315.nc' But there are more species listed in  gas_wetdep_list                = 'ALKOOH', 'C2H5OH', 'C2H5OOH', 'C3H7OOH', 'CH2O', 'CH3CHO', 'CH3CN', 'CH3COCHO', 'CH3COOH', 'CH3COOOH', 'CH3OH',         'CH3OOH', 'EOOH', 'GLYALD', 'H2O2', 'H2SO4', 'HCN', 'HCOOH', 'HNO3', 'HO2NO2', 'HYAC', 'HYDRALD',         'ISOPNO3', 'ISOPOOH', 'MACR', 'MACROOH', 'MEKOOH', 'MVK', 'NH3', 'ONIT', 'ONITR', 'POOH', 'Pb',         'ROOH', 'SO2', 'TERPOOH', 'TOLOOH', 'XOOH' 
 
T

tilmes

Member
1) co2vmr is set to that value to inidate that CO2 is used from the lower boundary file. This is very confusing and does not really make much sense but that's how it is set in this model version. Therefore, don't change it.2) This version was tested for the 1.9x2.5 resolutions and indeed all the emission data sets are for that resolution. You can run a higer resolution and the input datasets will be interpolated to the other resolution. Depending on what your research question is, you may want to prodcude your own emissions. 3) You find the full chemistry mechanism in CaseDocs/chem_mech.in . In surf_emis_specifier, all the componets are listed that have emissions, but of course this is not the case for all species, same as for the deposition.
 
C

chen63@illinois_edu

Member
Thanks for your replies during the holiday. Could you please explain to me more on the follow-up questions?(1) Which is the lower boundary file? How is it specified in the namelist?(2) Do you mean that the model will automatically interpolate the input datasets to the other resolution when I create the case for another resolution (e.g., 0.23*0.31 or 0.47*0.63)? If the answer is yea, if I need to reproduce some of the emission files, such as CO, can I just produce the emission file with 1.9*2.5 resolution, but create the case with finer resolutions by letting the model interpolate the data?Thanks a lot!
 
T

tilmes

Member
1) here is an example of the lower boundary condition file, as defined in the namelist:flbc_file              = '/glade/p/cesm/cseg//inputdata/atm/waccm/lb/LBC_1765-2500_1.9x2.5_CMIP5_RCP45_za_c120127.nc'2) yes the model will interpolate, but of course the information are still for a 2x resolution, so if you have higher resolution emissions that would be better. You may have to also apply some tuning if you change the resolution. Po-Lun has done high resolution MAM3 runs, so you can ask him about this as well.
 
You must log in or register to reply here.
Top