From e56542f636e046681a7cea3ef3ff1c9c799bd029 Mon Sep 17 00:00:00 2001 From: Ravan Ahmadov Date: Thu, 2 Mar 2017 18:27:26 -0700 Subject: [PATCH] TYPE: bug fixes SOURCE: internal, Doug Lowe (UK), Paolo Tuccella (Italy) DESCRIPTION OF CHANGES: 1) registry.chem is changed to reduce the memory footprint and remove unused 3D arrays. Some commented lines have been removed from registry too. 2) Several modules are changed so chemistry option RACM_SOA_VBS_AQCHEM_KPP (chem_opt=109) works properly. These fixes are mostly related to CASE statements for the RACM_SOA_VBS_AQCHEM_KPP option. LIST OF MODIFIED FILES: modified: ../Registry/registry.chem modified: ../chem/aerosol_driver.F modified: ../chem/chemics_init.F modified: ../chem/cloudchem_driver.F modified: ../chem/module_bioemi_megan2.F TESTS CONDUCTED: Code has been compiled and run successfully on NOAA HPC. The regression test is in progress. --- Registry/registry.chem | 38 +++++++++++++------------------------ chem/aerosol_driver.F | 2 +- chem/chemics_init.F | 10 +++++----- chem/cloudchem_driver.F | 2 +- chem/module_bioemi_megan2.F | 4 ++-- 5 files changed, 22 insertions(+), 34 deletions(-) diff --git a/Registry/registry.chem b/Registry/registry.chem index 5b56a49dc7..bfa540e458 100644 --- a/Registry/registry.chem +++ b/Registry/registry.chem @@ -3660,23 +3660,13 @@ state real c4h9ooh ikjftb chem 1 - i0{12}rhusdf=(bdy_ state real meo2 ikjftb chem 1 - i0{12}rhusdf=(bdy_interp:dt) "meo2" "" "" state real meoh ikjftb chem 1 - i0{12}rhusdf=(bdy_interp:dt) "meoh" "" "" state real meo2no2 ikjftb chem 1 - i0{12}rhusdf=(bdy_interp:dt) "meo2no2" "" "" - #End MATCH-MPI MAINZ global chemistry -# time averaged stuff -#state integer STEPAVE_COUNT - misc 1 - r "STEPAVE_COUNT" "time steps contained in averages for convective transport" "" -#state real RAINCV_A ij misc 1 - r "RAINCV_A" "taveragd TIME-STEP CUMULUS PRECIPITATION" "mm" -#state real RAINCV_B ij misc 1 - r "RAINCV_B" "taveragd TIME-STEP CUMULUS PRECIPITATION" "mm" -#state real GD_CLOUD_A ikj misc 1 - r "GD_CLOUD_A" "taveragd CLOUD WATER MIXING RAIO IN GD CLOUD" "kg kg-1" -#state real GD_CLOUD2_A ikj misc 1 - r "GD_CLOUD2_A" "taveragd cloud ice mix ratio in GD" "kg kg-1" -#state real GD_CLOUD_B ikj misc 1 - r "GD_CLOUD_B" "taveragd CLOUD WATER MIXING RAIO IN GD CLOUD" "kg kg-1" -#state real GD_CLOUD2_B ikj misc 1 - r "GD_CLOUD2_B" "taveragd cloud ice mix ratio in GD" "kg kg-1" state integer STEPBIOE - misc 1 - r "STEPBIOE" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN BIOGENIC EMIS CALLS" "NA" state integer STEPPHOT - misc 1 - r "STEPPHOT" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN PHOTOLYSIS CALLS" "NA" state integer STEPCHEM - misc 1 - r "STEPCHEM" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN CHEM MECH CALLS" "NA" state integer STEPFIREPL - misc 1 - r "STEPFIREPL" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN CALLSTO PLUMERISE FIRE ROUTINE" "NA" - # rconfig character emi_inname namelist,chem 1 "wrfchemi_d_" - "name of chem emissions infile" "" "" rconfig character fireemi_inname namelist,chem 1 "wrffirechemi_d_" - "name of chem fire emissions infile" "" "" @@ -3753,7 +3743,7 @@ rconfig integer gas_drydep_opt namelist,chem max_domains rconfig integer aer_drydep_opt namelist,chem max_domains 1 rh "aer_drydep_opt" "" "" rconfig integer diagnostic_chem namelist,chem max_domains 1 rh "diagnostic_chem" "" "" rconfig integer aer_aerodynres_opt namelist,chem max_domains 1 rh "aer_aerodynres_opt" "" "" -rconfig integer emiss_opt namelist,chem max_domains 4 rh "emiss_opt" "" "" +rconfig integer emiss_opt namelist,chem max_domains 0 rh "emiss_opt" "" "" rconfig integer emiss_opt_vol namelist,chem max_domains 0 rh "emiss_opt_vol" "" "" rconfig integer dust_opt namelist,chem 1 0 rh "dust_opt" "" "" rconfig integer dust_schme namelist,chem 1 2 rh "dust_schme" "" "" @@ -3774,6 +3764,7 @@ rconfig integer aer_ra_feedback namelist,chem max_domains rconfig integer aer_op_opt namelist,chem max_domains 1 rh "aer_op_opt" "" "" rconfig integer opt_pars_out namelist,chem 1 0 h "opt_pars_out" "" "" rconfig integer diagnostic_dep namelist,chem max_domains 0 rh "diagnostic_dep" "" "" +rconfig integer diagnostic_soa namelist,chem max_domains 0 rh "diagnostic_soa" "" "" # aircraft emissions rconfig integer aircraft_emiss_opt namelist,chem max_domains 0 rh "aircraft_emiss_opt" "" "" @@ -3910,6 +3901,9 @@ package saprc99_kpp chem_opt==195 - chem:o3,h2o2,no,no2,no3,n2o5,hono package saprc99_mosaic_4bin_vbs2_kpp chem_opt==198 - chem:o3,h2o2,no,no2,no3,n2o5,hono,hno3,hno4,so2,h2so4,co,hcho,ccho,rcho,acet,mek,hcooh,meoh,etoh,cco_oh,rco_oh,gly,mgly,bacl,cres,bald,isoprod,methacro,mvk,prod2,dcb1,dcb2,dcb3,ethene,isoprene,c2h6,c3h8,c2h2,c3h6,alk3,alk4,alk5,aro1,aro2,ole1,ole2,terp,rno3,nphe,phen,pan,pan2,pbzn,ma_pan,co2,cco_ooh,rco_o2,rco_ooh,xn,xc,ho,ho2,c_o2,cooh,rooh,ro2_r,r2o2,ro2_n,cco_o2,bzco_o2,ma_rco3,sesq,pcg1_b_c,pcg2_b_c,pcg1_b_o,pcg2_b_o,opcg1_b_c,opcg1_b_o,pcg1_f_c,pcg2_f_c,pcg1_f_o,pcg2_f_o,opcg1_f_c,opcg1_f_o,psd1,psd2,nh3,hcl,nume,den,ant1_c,ant1_o,biog1_c,biog1_o,ch4,so4_a01,no3_a01,cl_a01,nh4_a01,na_a01,oin_a01,oc_a01,bc_a01,hysw_a01,water_a01,pcg1_b_c_a01,pcg2_b_c_a01,pcg1_b_o_a01,pcg2_b_o_a01,opcg1_b_c_a01,opcg1_b_o_a01,pcg1_f_c_a01,pcg2_f_c_a01,pcg1_f_o_a01,pcg2_f_o_a01,opcg1_f_c_a01,opcg1_f_o_a01,ant1_c_a01,ant1_o_a01,biog1_c_a01,biog1_o_a01,num_a01,so4_a02,no3_a02,cl_a02,nh4_a02,na_a02,oin_a02,oc_a02,bc_a02,hysw_a02,water_a02,pcg1_b_c_a02,pcg2_b_c_a02,pcg1_b_o_a02,pcg2_b_o_a02,opcg1_b_c_a02,opcg1_b_o_a02,pcg1_f_c_a02,pcg2_f_c_a02,pcg1_f_o_a02,pcg2_f_o_a02,opcg1_f_c_a02,opcg1_f_o_a02,ant1_c_a02,ant1_o_a02,biog1_c_a02,biog1_o_a02,num_a02,so4_a03,no3_a03,cl_a03,nh4_a03,na_a03,oin_a03,oc_a03,bc_a03,hysw_a03,water_a03,pcg1_b_c_a03,pcg2_b_c_a03,pcg1_b_o_a03,pcg2_b_o_a03,opcg1_b_c_a03,opcg1_b_o_a03,pcg1_f_c_a03,pcg2_f_c_a03,pcg1_f_o_a03,pcg2_f_o_a03,opcg1_f_c_a03,opcg1_f_o_a03,ant1_c_a03,ant1_o_a03,biog1_c_a03,biog1_o_a03,num_a03,so4_a04,no3_a04,cl_a04,nh4_a04,na_a04,oin_a04,oc_a04,bc_a04,hysw_a04,water_a04,pcg1_b_c_a04,pcg2_b_c_a04,pcg1_b_o_a04,pcg2_b_o_a04,opcg1_b_c_a04,opcg1_b_o_a04,pcg1_f_c_a04,pcg2_f_c_a04,pcg1_f_o_a04,pcg2_f_o_a04,opcg1_f_c_a04,opcg1_f_o_a04,ant1_c_a04,ant1_o_a04,biog1_c_a04,biog1_o_a04,num_a04,ca_a01,ca_a02,ca_a03,ca_a04,co3_a01,co3_a02,co3_a03,co3_a04 +# KPP mechanism from MATCH-MPI Mainz used for global chemistry +package nmhc9_kpp chem_opt==200 - chem:o3,h2o2,ch4,op1,hcho,ch3oh,co,hno3,no3,n2o5,hno4,no,no2,isopr,mvk,iso2,isooh,mvko2,mvkooh,ison,aca,acol,hcooh,mpan,naca,pan,pa,paa,mglo,c2h6,eTooh,ald,c3h8,pRooh,acet,acooh,eTo2,pRo2,aco2,c3h6,c3h6ooh,c2h4,c4h10,buooh,mek,mekooh,mEcoco,c3h6o2,c4h9o2,meko2,onit,pRono2,ch3o2,acetol,acetp,aceto2,ch3cooh,c4h9ooh,mEoh,ho,ho2,mEo2,mEo2no2 + # using an empirical soa scheme for anthropogenic and biomass burning sources #package mozart_mosaic_4bin_vbs0_kpp chem_opt==201 - chem:o3,o1d_cb4,o,no,no2,no3,n2o5,hno3,hno4,so2,ho,ho2,h2o2,sulf,co,hcho,ch3ooh,ch3o2,ch4,h2,eo2,ch3cooh,c2h4,n2o,ch3oh,aco3,acet,mgly,paa,gly,c3h6ooh,pan,mpan,macr,mvk,c3h6,etooh,prooh,acetp,xooh,onitr,isooh,acetol,glyald,mek,eto2,open,alkooh,mekooh,tolooh,terpooh,ald,mco3,c2h5oh,eo,c2h6,c3h8,pro2,po2,aceto2,bigene,bigalk,eneo2,alko2,isopr,iso2,mvko2,mvkooh,hydrald,xo2,c10h16,terpo2,tol,cres,to2,xoh,onit,isopn,dms,nh3,nume,den,voca,vocbb,smpa,smpbb,biog1_c,biog1_o,meko2,so4_a01,no3_a01,cl_a01,nh4_a01,na_a01,oin_a01,oc_a01,bc_a01,hysw_a01,water_a01,smpa_a01,smpbb_a01,biog1_c_a01,biog1_o_a01,num_a01,so4_a02,no3_a02,cl_a02,nh4_a02,na_a02,oin_a02,oc_a02,bc_a02,hysw_a02,water_a02,smpa_a02,smpbb_a02,biog1_c_a02,biog1_o_a02,num_a02,so4_a03,no3_a03,cl_a03,nh4_a03,na_a03,oin_a03,oc_a03,bc_a03,hysw_a03,water_a03,smpa_a03,smpbb_a03,biog1_c_a03,biog1_o_a03,num_a03,so4_a04,no3_a04,cl_a04,nh4_a04,na_a04,oin_a04,oc_a04,bc_a04,hysw_a04,water_a04,smpa_a04,smpbb_a04,biog1_c_a04,biog1_o_a04,num_a04,ca_a01,ca_a02,ca_a03,ca_a04,co3_a01,co3_a02,co3_a03,co3_a04 package mozart_mosaic_4bin_kpp chem_opt==201 - chem:o3,o1d_cb4,o,no,no2,no3,n2o5,hno3,hno4,so2,ho,ho2,h2o2,sulf,co,hcho,hcooh,c2h2,hoch2oo,ch3ooh,ch3o2,ch4,h2,eo2,ch3cooh,c2h4,n2o,ch3oh,aco3,acet,mgly,paa,gly,c3h6ooh,pan,mpan,macr,mvk,c3h6,etooh,prooh,acetp,xooh,onitr,isooh,acetol,glyald,mek,eto2,open,alkooh,mekooh,tolooh,terpooh,ald,mco3,c2h5oh,eo,c2h6,c3h8,pro2,po2,aceto2,bigene,bigalk,eneo2,alko2,isopr,iso2,mvko2,mvkooh,hydrald,xo2,apin,bpin,limon,myrc,bcary,terprod1,terprod2,terp2o2,terp2ooh,nterpo2,terpo2,tol,cres,to2,onit,isopn,dms,mbo,mboo2,hmprop,hmpropo2,mboooh,mbono3o2,nh3,nume,den,voca,vocbb,smpa,smpbb,biog1_c,biog1_o,benzene,phen,bepomuc,benzo2,pheno2,pheno,phenooh,c6h5o2,c6h5ooh,benzooh,bigald1,bigald2,bigald3,bigald4,malo2,tepomuc,bzoo,bzooh,bald,acbzo2,dicarbo2,mdialo2,xyl,xylol,xylolo2,xylolooh,xyleno2,xylenooh,pbznit,hono,meko2,so4_a01,no3_a01,smpa_a01,smpbb_a01,glysoa_sfc_a01,biog1_c_a01,biog1_o_a01,cl_a01,co3_a01,nh4_a01,na_a01,ca_a01,oin_a01,oc_a01,bc_a01,hysw_a01,water_a01,num_a01,so4_a02,no3_a02,smpa_a02,smpbb_a02,glysoa_sfc_a02,biog1_c_a02,biog1_o_a02,cl_a02,co3_a02,nh4_a02,na_a02,ca_a02,oin_a02,oc_a02,bc_a02,hysw_a02,water_a02,num_a02,so4_a03,no3_a03,smpa_a03,smpbb_a03,glysoa_sfc_a03,biog1_c_a03,biog1_o_a03,cl_a03,co3_a03,nh4_a03,na_a03,ca_a03,oin_a03,oc_a03,bc_a03,hysw_a03,water_a03,num_a03,so4_a04,no3_a04,smpa_a04,smpbb_a04,glysoa_sfc_a04,biog1_c_a04,biog1_o_a04,cl_a04,co3_a04,nh4_a04,na_a04,ca_a04,oin_a04,oc_a04,bc_a04,hysw_a04,water_a04,num_a04 @@ -3921,17 +3915,6 @@ package saprc99_mosaic_8bin_vbs2_aq_kpp chem_opt==203 - chem:o3,h2o2, package saprc99_mosaic_8bin_vbs2_kpp chem_opt==204 - chem:o3,h2o2,no,no2,no3,n2o5,hono,hno3,hno4,so2,h2so4,co,hcho,ccho,rcho,acet,mek,hcooh,meoh,etoh,cco_oh,rco_oh,gly,mgly,bacl,cres,bald,isoprod,methacro,mvk,prod2,dcb1,dcb2,dcb3,ethene,isoprene,c2h6,c3h8,c2h2,c3h6,alk3,alk4,alk5,aro1,aro2,ole1,ole2,terp,rno3,nphe,phen,pan,pan2,pbzn,ma_pan,co2,cco_ooh,rco_o2,rco_ooh,xn,xc,ho,ho2,c_o2,cooh,rooh,ro2_r,r2o2,ro2_n,cco_o2,bzco_o2,ma_rco3,sesq,pcg1_b_c,pcg2_b_c,pcg1_b_o,pcg2_b_o,opcg1_b_c,opcg1_b_o,pcg1_f_c,pcg2_f_c,pcg1_f_o,pcg2_f_o,opcg1_f_c,opcg1_f_o,psd1,psd2,nh3,hcl,nume,den,ant1_c,biog1_c,ant2_c,biog2_c,biog3_c,biog1_o,biog2_o,ant3_c,ant4_c,bgas,agas,ch4,so4_a01,no3_a01,cl_a01,nh4_a01,na_a01,oin_a01,oc_a01,bc_a01,hysw_a01,water_a01,pcg1_b_c_a01,pcg1_b_o_a01,opcg1_b_c_a01,opcg1_b_o_a01,pcg1_f_c_a01,pcg1_f_o_a01,opcg1_f_c_a01,opcg1_f_o_a01,ant1_c_a01,biog1_c_a01,ant2_c_a01,biog2_c_a01,biog3_c_a01,biog1_o_a01,biog2_o_a01,ant3_c_a01,ant4_c_a01,num_a01,so4_a02,no3_a02,cl_a02,nh4_a02,na_a02,oin_a02,oc_a02,bc_a02,hysw_a02,water_a02,pcg1_b_c_a02,pcg1_b_o_a02,opcg1_b_c_a02,opcg1_b_o_a02,pcg1_f_c_a02,pcg1_f_o_a02,opcg1_f_c_a02,opcg1_f_o_a02,ant1_c_a02,biog1_c_a02,ant2_c_a02,biog2_c_a02,biog3_c_a02,biog1_o_a02,biog2_o_a02,ant3_c_a02,ant4_c_a02,num_a02,so4_a03,no3_a03,cl_a03,nh4_a03,na_a03,oin_a03,oc_a03,bc_a03,hysw_a03,water_a03,pcg1_b_c_a03,pcg1_b_o_a03,opcg1_b_c_a03,opcg1_b_o_a03,pcg1_f_c_a03,pcg1_f_o_a03,opcg1_f_c_a03,opcg1_f_o_a03,ant1_c_a03,biog1_c_a03,ant2_c_a03,biog2_c_a03,biog3_c_a03,biog1_o_a03,biog2_o_a03,ant3_c_a03,ant4_c_a03,num_a03,so4_a04,no3_a04,cl_a04,nh4_a04,na_a04,oin_a04,oc_a04,bc_a04,hysw_a04,water_a04,pcg1_b_c_a04,pcg1_b_o_a04,opcg1_b_c_a04,opcg1_b_o_a04,pcg1_f_c_a04,pcg1_f_o_a04,opcg1_f_c_a04,opcg1_f_o_a04,ant1_c_a04,biog1_c_a04,ant2_c_a04,biog2_c_a04,biog3_c_a04,biog1_o_a04,biog2_o_a04,ant3_c_a04,ant4_c_a04,num_a04,so4_a05,no3_a05,cl_a05,nh4_a05,na_a05,oin_a05,oc_a05,bc_a05,hysw_a05,water_a05,pcg1_b_c_a05,pcg1_b_o_a05,opcg1_b_c_a05,opcg1_b_o_a05,pcg1_f_c_a05,pcg1_f_o_a05,opcg1_f_c_a05,opcg1_f_o_a05,ant1_c_a05,biog1_c_a05,ant2_c_a05,biog2_c_a05,biog3_c_a05,biog1_o_a05,biog2_o_a05,ant3_c_a05,ant4_c_a05,num_a05,so4_a06,no3_a06,cl_a06,nh4_a06,na_a06,oin_a06,oc_a06,bc_a06,hysw_a06,water_a06,pcg1_b_c_a06,pcg1_b_o_a06,opcg1_b_c_a06,opcg1_b_o_a06,pcg1_f_c_a06,pcg1_f_o_a06,opcg1_f_c_a06,opcg1_f_o_a06,ant1_c_a06,biog1_c_a06,ant2_c_a06,biog2_c_a06,biog3_c_a06,biog1_o_a06,biog2_o_a06,ant3_c_a06,ant4_c_a06,num_a06,so4_a07,no3_a07,cl_a07,nh4_a07,na_a07,oin_a07,oc_a07,bc_a07,hysw_a07,water_a07,pcg1_b_c_a07,pcg1_b_o_a07,opcg1_b_c_a07,opcg1_b_o_a07,pcg1_f_c_a07,pcg1_f_o_a07,opcg1_f_c_a07,opcg1_f_o_a07,ant1_c_a07,biog1_c_a07,ant2_c_a07,biog2_c_a07,biog3_c_a07,biog1_o_a07,biog2_o_a07,ant3_c_a07,ant4_c_a07,num_a07,so4_a08,no3_a08,cl_a08,nh4_a08,na_a08,oin_a08,oc_a08,bc_a08,hysw_a08,water_a08,pcg1_b_c_a08,pcg1_b_o_a08,opcg1_b_c_a08,opcg1_b_o_a08,pcg1_f_c_a08,pcg1_f_o_a08,opcg1_f_c_a08,opcg1_f_o_a08,ant1_c_a08,biog1_c_a08,ant2_c_a08,biog2_c_a08,biog3_c_a08,biog1_o_a08,biog2_o_a08,ant3_c_a08,ant4_c_a08,num_a08 - - - - - - - - -# KPP mechanism from MATCH-MPI Mainz used for global chemistry -package nmhc9_kpp chem_opt==200 - chem:o3,h2o2,ch4,op1,hcho,ch3oh,co,hno3,no3,n2o5,hno4,no,no2,isopr,mvk,iso2,isooh,mvko2,mvkooh,ison,aca,acol,hcooh,mpan,naca,pan,pa,paa,mglo,c2h6,eTooh,ald,c3h8,pRooh,acet,acooh,eTo2,pRo2,aco2,c3h6,c3h6ooh,c2h4,c4h10,buooh,mek,mekooh,mEcoco,c3h6o2,c4h9o2,meko2,onit,pRono2,ch3o2,acetol,acetp,aceto2,ch3cooh,c4h9ooh,mEoh,ho,ho2,mEo2,mEo2no2 - # GOCART Packages (no Ozone chemistry in this choice) package gocart_simple chem_opt==300 - chem:so2,sulf,dms,msa,p25,bc1,bc2,oc1,oc2,dust_1,dust_2,dust_3,dust_4,dust_5,seas_1,seas_2,seas_3,seas_4,p10 @@ -3971,6 +3954,7 @@ package cri_mosaic_4bin_aq_kpp chem_opt==611 - chem:dms, #emission package definitions # +package noantemis emiss_opt==0 - - package eradm emiss_opt==2 - emis_ant:e_iso,e_so2,e_no,e_co,e_eth,e_hc3,e_hc5,e_hc8,e_xyl,e_ol2,e_olt,e_oli,e_tol,e_csl,e_hcho,e_ald,e_ket,e_ora2,e_nh3 package eradmsorg emiss_opt==3 - emis_ant:e_iso,e_so2,e_no,e_no2,e_co,e_eth,e_hc3,e_hc5,e_hc8,e_xyl,e_ol2,e_olt,e_oli,e_tol,e_csl,e_hcho,e_ald,e_ket,e_ora2,e_nh3,e_pm25i,e_pm25j,e_pm_10,e_eci,e_ecj,e_orgi,e_orgj,e_so4i,e_so4j,e_no3i,e_no3j,e_naaj,e_naai,e_orgi_a,e_orgj_a,e_orgi_bb,e_orgj_bb package ecbmz_mosaic emiss_opt==4 - emis_ant:e_iso,e_so2,e_no,e_co,e_eth,e_hc3,e_hc5,e_hc8,e_xyl,e_ol2,e_olt,e_oli,e_tol,e_csl,e_hcho,e_ald,e_ket,e_ora2,e_nh3,e_no2,e_ch3oh,e_c2h5oh,e_pm25i,e_pm25j,e_eci,e_ecj,e_orgi,e_orgj,e_so4i,e_so4j,e_no3i,e_no3j,e_so4c,e_no3c,e_orgc,e_ecc @@ -3985,7 +3969,7 @@ package esaprcnov emiss_opt==13 - emis_ant: # For CB05 mechanism based on CBMZ speciation, used with emiss_inpt_opt = 102 package ecb05_opt1 emiss_opt==14 - emis_ant:e_no2,e_xyl,e_tol,e_terp,e_so2,e_ora2,e_olt,e_oli,e_ol2,e_no,e_nh3,e_iso,e_hcl,e_hcho,e_eth,e_csl,e_co,e_ch3oh,e_c2h5oh,e_ald,e_aldx,e_hc3,e_hc5,e_hc8,e_ket,e_pm25i,e_pm25j,e_eci,e_ecj,e_orgi,e_orgj,e_so4i,e_so4j,e_no3i,e_no3j,e_so4c,e_no3c,e_orgc,e_ecc,e_pm10 # For CB05 emissions inventory base on CB05 speciation, used with emiss_inpt_opt = 101 -package ecb05_opt2 emiss_opt==15 - emis_ant:e_acet,e_par,e_alk1,e_alk2,e_alk3,e_alk4,e_alk5,e_tol,e_xyl,e_bald,e_ald2,e_ccooh,e_co,e_cres,e_eth,e_etha,e_gly,e_form,e_hcooh,e_iprod,e_isop,e_macr,e_mek,e_meoh,e_meo2,e_etoh,e_mgly,e_nh3,e_hcl,e_no,e_no2,e_iole,e_ole,e_phen,e_prod2,e_aldx,e_rcooh,e_so2,e_psulf,e_terp,e_pm25i,e_pm25j,e_eci,e_ecj,e_orgi,e_orgj,e_so4i,e_so4j,e_no3i,e_no3j,e_so4c,e_no3c,e_orgc,e_ecc,e_pm10 +package ecb05_opt2 emiss_opt==15 - emis_ant:e_acet,e_par,e_alk3,e_alk4,e_alk5,e_tol,e_xyl,e_bald,e_ald2,e_ccooh,e_co,e_cres,e_eth,e_etha,e_gly,e_form,e_hcooh,e_iprod,e_isop,e_macr,e_mek,e_meoh,e_meo2,e_etoh,e_mgly,e_nh3,e_hcl,e_no,e_no2,e_iole,e_ole,e_phen,e_prod2,e_aldx,e_so2,e_psulf,e_terp,e_pm25i,e_pm25j,e_eci,e_ecj,e_orgi,e_orgj,e_so4i,e_so4j,e_no3i,e_no3j,e_so4c,e_no3c,e_orgc,e_ecc,e_pm10 # Anthropogenic CO2, CO and CH4 emissions: package eco2 emiss_opt==16 - emis_ant:e_co2,e_co2tst,e_co @@ -4011,14 +3995,18 @@ package wesely gas_drydep_opt==1 - - # diagnostic packages, first for deposition velocities (original package for Mozart) # package depvel1 diagnostic_chem==1 - dvel:dvel_o3,dvel_no,dvel_no2,dvel_nh3,dvel_hno3,dvel_hno4,dvel_h2o2,dvel_co,dvel_ch3ooh,dvel_hcho,dvel_ch3oh,dvel_eo2,dvel_ald,dvel_ch3cooh,dvel_acet,dvel_mgly,dvel_gly,dvel_paa,dvel_pooh,dvel_pan,dvel_mpan,dvel_mco3,dvel_mvkooh,dvel_c2h5oh,dvel_etooh,dvel_prooh,dvel_acetp,dvel_onit,dvel_onitr,dvel_isooh,dvel_acetol,dvel_glyald,dvel_hydrald,dvel_alkooh,dvel_mekooh,dvel_tolooh,dvel_xooh,dvel_so2,dvel_so4,dvel_terpooh,dvel_cvasoaX,dvel_cvasoa1,dvel_cvasoa2,dvel_cvasoa3,dvel_cvasoa4,dvel_cvbsoaX,dvel_cvbsoa1,dvel_cvbsoa2,dvel_cvbsoa3,dvel_cvbsoa4,ddmass_o3,ddmass_no,ddmass_no2,ddmass_nh3,ddmass_hno3,ddmass_hno4,ddmass_h2o2,ddmass_co,ddmass_ch3ooh,ddmass_hcho,ddmass_ch3oh,ddmass_eo2,ddmass_ald,ddmass_ch3cooh,ddmass_acet,ddmass_mgly,ddmass_gly,ddmass_paa,ddmass_pooh,ddmass_pan,ddmass_mpan,ddmass_mco3,ddmass_mvkooh,ddmass_c2h5oh,ddmass_etooh,ddmass_prooh,ddmass_acetp,ddmass_onit,ddmass_onitr,ddmass_isooh,ddmass_acetol,ddmass_glyald,ddmass_hydrald,ddmass_alkooh,ddmass_mekooh,ddmass_tolooh,ddmass_xooh,ddmass_so2,ddmass_so4,ddmass_terpooh,ddmass_cvasoaX,ddmass_cvasoa1,ddmass_cvasoa2,ddmass_cvasoa3,ddmass_cvasoa4,ddmass_cvbsoaX,ddmass_cvbsoa1,ddmass_cvbsoa2,ddmass_cvbsoa3,ddmass_cvbsoa4,ddmass_so4_a01,ddmass_no3_a01,ddmass_cl_a01,ddmass_nh4_a01,ddmass_na_a01,ddmass_oin_a01,ddmass_oc_a01,ddmass_bc_a01,ddmass_smpa_a01,ddmass_smpbb_a01,ddmass_glysoa_a01,ddmass_biog1_c_a01,ddmass_biog1_o_a01,ddmass_asoaX_a01,ddmass_asoa1_a01,ddmass_asoa2_a01,ddmass_asoa3_a01,ddmass_asoa4_a01,ddmass_bsoaX_a01,ddmass_bsoa1_a01,ddmass_bsoa2_a01,ddmass_bsoa3_a01,ddmass_bsoa4_a01,ddmass_so4_a02,ddmass_no3_a02,ddmass_cl_a02,ddmass_nh4_a02,ddmass_na_a02,ddmass_oin_a02,ddmass_oc_a02,ddmass_bc_a02,ddmass_smpa_a02,ddmass_smpbb_a02,ddmass_glysoa_a02,ddmass_biog1_c_a02,ddmass_biog1_o_a02,ddmass_asoaX_a02,ddmass_asoa1_a02,ddmass_asoa2_a02,ddmass_asoa3_a02,ddmass_asoa4_a02,ddmass_bsoaX_a02,ddmass_bsoa1_a02,ddmass_bsoa2_a02,ddmass_bsoa3_a02,ddmass_bsoa4_a02,ddmass_so4_a03,ddmass_no3_a03,ddmass_cl_a03,ddmass_nh4_a03,ddmass_na_a03,ddmass_oin_a03,ddmass_oc_a03,ddmass_bc_a03,ddmass_smpa_a03,ddmass_smpbb_a03,ddmass_glysoa_a03,ddmass_biog1_c_a03,ddmass_biog1_o_a03,ddmass_asoaX_a03,ddmass_asoa1_a03,ddmass_asoa2_a03,ddmass_asoa3_a03,ddmass_asoa4_a03,ddmass_bsoaX_a03,ddmass_bsoa1_a03,ddmass_bsoa2_a03,ddmass_bsoa3_a03,ddmass_bsoa4_a03,ddmass_so4_a04,ddmass_no3_a04,ddmass_cl_a04,ddmass_nh4_a04,ddmass_na_a04,ddmass_oin_a04,ddmass_oc_a04,ddmass_bc_a04,ddmass_smpa_a04,ddmass_smpbb_a04,ddmass_glysoa_a04,ddmass_biog1_c_a04,ddmass_biog1_o_a04,ddmass_asoaX_a04,ddmass_asoa1_a04,ddmass_asoa2_a04,ddmass_asoa3_a04,ddmass_asoa4_a04,ddmass_bsoaX_a04,ddmass_bsoa1_a04,ddmass_bsoa2_a04,ddmass_bsoa3_a04,ddmass_bsoa4_a04,ddmass_ca_a01,ddmass_ca_a02,ddmass_ca_a03,ddmass_ca_a04,ddmass_co3_a01,ddmass_co3_a02,ddmass_co3_a03,ddmass_co3_a04,ddmass_so4_cw01,ddmass_no3_cw01,ddmass_cl_cw01,ddmass_nh4_cw01,ddmass_na_cw01,ddmass_oin_cw01,ddmass_oc_cw01,ddmass_bc_cw01,ddmass_smpa_cw01,ddmass_smpbb_cw01,ddmass_glysoa_cw01,ddmass_biog1_c_cw01,ddmass_biog1_o_cw01,ddmass_asoaX_cw01,ddmass_asoa1_cw01,ddmass_asoa2_cw01,ddmass_asoa3_cw01,ddmass_asoa4_cw01,ddmass_bsoaX_cw01,ddmass_bsoa1_cw01,ddmass_bsoa2_cw01,ddmass_bsoa3_cw01,ddmass_bsoa4_cw01,ddmass_so4_cw02,ddmass_no3_cw02,ddmass_cl_cw02,ddmass_nh4_cw02,ddmass_na_cw02,ddmass_oin_cw02,ddmass_oc_cw02,ddmass_bc_cw02,ddmass_smpa_cw02,ddmass_smpbb_cw02,ddmass_glysoa_cw02,ddmass_biog1_c_cw02,ddmass_biog1_o_cw02,ddmass_asoaX_cw02,ddmass_asoa1_cw02,ddmass_asoa2_cw02,ddmass_asoa3_cw02,ddmass_asoa4_cw02,ddmass_bsoaX_cw02,ddmass_bsoa1_cw02,ddmass_bsoa2_cw02,ddmass_bsoa3_cw02,ddmass_bsoa4_cw02,ddmass_so4_cw03,ddmass_no3_cw03,ddmass_cl_cw03,ddmass_nh4_cw03,ddmass_na_cw03,ddmass_oin_cw03,ddmass_oc_cw03,ddmass_bc_cw03,ddmass_smpa_cw03,ddmass_smpbb_cw03,ddmass_glysoa_cw03,ddmass_biog1_c_cw03,ddmass_biog1_o_cw03,ddmass_asoaX_cw03,ddmass_asoa1_cw03,ddmass_asoa2_cw03,ddmass_asoa3_cw03,ddmass_asoa4_cw03,ddmass_bsoaX_cw03,ddmass_bsoa1_cw03,ddmass_bsoa2_cw03,ddmass_bsoa3_cw03,ddmass_bsoa4_cw03,ddmass_so4_cw04,ddmass_no3_cw04,ddmass_cl_cw04,ddmass_nh4_cw04,ddmass_na_cw04,ddmass_oin_cw04,ddmass_oc_cw04,ddmass_bc_cw04,ddmass_smpa_cw04,ddmass_smpbb_cw04,ddmass_glysoa_cw04,ddmass_biog1_c_cw04,ddmass_biog1_o_cw04,ddmass_asoaX_cw04,ddmass_asoa1_cw04,ddmass_asoa2_cw04,ddmass_asoa3_cw04,ddmass_asoa4_cw04,ddmass_bsoaX_cw04,ddmass_bsoa1_cw04,ddmass_bsoa2_cw04,ddmass_bsoa3_cw04,ddmass_bsoa4_cw04,ddmass_ca_cw01,ddmass_ca_cw02,ddmass_ca_cw03,ddmass_ca_cw04,ddmass_co3_cw01,ddmass_co3_cw02,ddmass_co3_cw03,ddmass_co3_cw04 + +# Diagnostic packages for the SOA schemes +package diag_vbs diagnostic_soa==108 - state:br_rto + # package gunther1 bio_emiss_opt==1 - - package beis314 bio_emiss_opt==2 - - package megan2 bio_emiss_opt==3 - state:mebio_isop,mebio_apin,mebio_bcar,mebio_acet,mebio_mbo,mebio_no,msebio_isop,mlai,pftp_bt,pftp_nt,pftp_sb,pftp_hb,mtsa,mswdown,EFmegan # Biospheric CO2 and CH4 emissions -package ebioco2 bio_emiss_opt==16 - vprm_in:vegfra_vprm,evi,evi_min,evi_max,lswi,lswi_max,lswi_min;eghg_bio:ebio_gee,ebio_res,ebio_co2oce -package ebioghg bio_emiss_opt==17 - vprm_in:vegfra_vprm,evi,evi_min,evi_max,lswi,lswi_max,lswi_min;wet_in:cpool,wetmap,t_ann;eghg_bio:ebio_gee,ebio_res,ebio_co2oce,ebio_ch4wet,ebio_ch4soil,ebio_ch4term +package ebioco2 bio_emiss_opt==16 - state:rad_vprm,lambda_vprm,alpha_vprm,resp_vprm;vprm_in:vegfra_vprm,evi,evi_min,evi_max,lswi,lswi_max,lswi_min;eghg_bio:ebio_gee,ebio_res,ebio_co2oce +package ebioghg bio_emiss_opt==17 - state:rad_vprm,lambda_vprm,alpha_vprm,resp_vprm;vprm_in:vegfra_vprm,evi,evi_min,evi_max,lswi,lswi_max,lswi_min;wet_in:cpool,wetmap,t_ann;eghg_bio:ebio_gee,ebio_res,ebio_co2oce,ebio_ch4wet,ebio_ch4soil,ebio_ch4term # the following arrays only needed for burn options package biomassb biomass_burn_opt==1 - state:mean_fct_agtf,mean_fct_agef,mean_fct_agsv,mean_fct_aggr,firesize_agtf,firesize_agef,firesize_agsv,firesize_aggr;ebu:ebu_no,ebu_no2,ebu_co,ebu_co2,ebu_eth,ebu_hc3,ebu_hc5,ebu_hc8,ebu_ete,ebu_olt,ebu_oli,ebu_pm25,ebu_pm10,ebu_dien,ebu_iso,ebu_api,ebu_lim,ebu_tol,ebu_csl,ebu_hcho,ebu_ald,ebu_ket,ebu_macr,ebu_ora1,ebu_ora2,ebu_so2,ebu_nh3,ebu_oc,ebu_bc,ebu_sulf,ebu_dms;ebu_in:ebu_in_no,ebu_in_no2,ebu_in_co,ebu_in_co2,ebu_in_eth,ebu_in_hc3,ebu_in_hc5,ebu_in_hc8,ebu_in_ete,ebu_in_olt,ebu_in_oli,ebu_in_pm25,ebu_in_pm10,ebu_in_dien,ebu_in_iso,ebu_in_api,ebu_in_lim,ebu_in_tol,ebu_in_csl,ebu_in_hcho,ebu_in_ald,ebu_in_xyl,ebu_in_ket,ebu_in_macr,ebu_in_ora1,ebu_in_ora2,ebu_in_so2,ebu_in_nh3,ebu_in_oc,ebu_in_bc,ebu_in_sulf,ebu_in_dms,ebu_in_ash diff --git a/chem/aerosol_driver.F b/chem/aerosol_driver.F index c16c330090..a804ef01f3 100755 --- a/chem/aerosol_driver.F +++ b/chem/aerosol_driver.F @@ -492,7 +492,7 @@ SUBROUTINE sum_pm_driver ( config_flags, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) - CASE (RACM_SOA_VBS_KPP) + CASE (RACM_SOA_VBS_KPP,RACM_SOA_VBS_AQCHEM_KPP) CALL wrf_debug(15,'sum_pm_driver: calling sum_pm_soa_vbs') CALL sum_pm_soa_vbs ( & alt, chem, h2oaj, h2oai, & diff --git a/chem/chemics_init.F b/chem/chemics_init.F index d3f3ae62ff..aec477c9a2 100755 --- a/chem/chemics_init.F +++ b/chem/chemics_init.F @@ -143,12 +143,12 @@ subroutine chem_init (id,chem,emis_ant,scalar,dt,bioemdt,photdt,chemdt,stepbioe, #ifdef CHEM_DBG_I call print_chem_species_index( config_flags%chem_opt ) #endif - program_name = "* PROGRAM:WRF/CHEM " // TRIM(release_version) // " MODEL" + program_name = "* PROGRAM:WRF-CHEM " // TRIM(release_version) // " MODEL" call wrf_message("*********************************************************************") call wrf_message(program_name) call wrf_message("* *") -call wrf_message("* PLEASE REPORT ANY BUGS TO WRF/CHEM HELP at *") +call wrf_message("* PLEASE REPORT ANY BUGS TO WRF-CHEM HELP at *") call wrf_message("* *") call wrf_message("* wrfchemhelp.gsd@noaa.gov *") call wrf_message("* *") @@ -1746,7 +1746,7 @@ subroutine chem_init (id,chem,emis_ant,scalar,dt,bioemdt,photdt,chemdt,stepbioe, endif chem(its:ite,kts:min(kte,kde-1),jts:jte,:)=max(chem(its:ite,kts:min(kte,kde-1),jts:jte,:),epsilc) ! - CASE (RACM_SOA_VBS_AQCHEM_KPP,CB05_SORG_VBS_AQ_KPP) + CASE (CB05_SORG_VBS_AQ_KPP) CALL wrf_debug(15,'call MADE/SOA_VBS aerosols initialization') call aerosols_sorgam_vbs_init(chem,convfac,z_at_w, & @@ -1773,7 +1773,7 @@ subroutine chem_init (id,chem,emis_ant,scalar,dt,bioemdt,photdt,chemdt,stepbioe, endif chem(its:ite,kts:min(kte,kde-1),jts:jte,:)=max(chem(its:ite,kts:min(kte,kde-1),jts:jte,:),epsilc) - CASE (RACM_SOA_VBS_KPP) + CASE (RACM_SOA_VBS_KPP,RACM_SOA_VBS_AQCHEM_KPP) CALL wrf_debug(15,'call MADE/SOA_VBS aerosols initialization') call aerosols_soa_vbs_init(chem,convfac,z_at_w, & @@ -3342,7 +3342,7 @@ subroutine print_chem_species_index( chem_opt ) print*,p_ac0,"ac0" print*,p_corn,"corn" !!! TUCCELLA -case (RACM_SOA_VBS_AQCHEM_KPP) + case (RACM_SOA_VBS_AQCHEM_KPP) print*,p_so4aj,"so4aj" print*,p_so4ai,"so4ai" print*,p_nh4aj,"nh4aj" diff --git a/chem/cloudchem_driver.F b/chem/cloudchem_driver.F index 0b216818c3..d3116bf249 100644 --- a/chem/cloudchem_driver.F +++ b/chem/cloudchem_driver.F @@ -249,7 +249,7 @@ SUBROUTINE cloudchem_driver( & CASE ( RADM2SORG_AQCHEM, RACMSORG_AQCHEM_KPP, RACM_ESRLSORG_AQCHEM_KPP, & - CB05_SORG_AQ_KPP ) + CB05_SORG_AQ_KPP,RACM_SOA_VBS_AQCHEM_KPP ) call wrf_debug(15, & 'cloudchem_driver calling sorgam_aqchem_driver') diff --git a/chem/module_bioemi_megan2.F b/chem/module_bioemi_megan2.F index 50d7ef275d..2fd74a54d5 100644 --- a/chem/module_bioemi_megan2.F +++ b/chem/module_bioemi_megan2.F @@ -1137,7 +1137,7 @@ SUBROUTINE bio_emissions_megan2(id,config_flags,ktau,dtstep, & END DO - CASE (RACM_SOA_VBS_AQCHEM_KPP,CB05_SORG_AQ_KPP) + CASE (CB05_SORG_AQ_KPP) DO icount = 1, n_megan2cb05 IF ( p_of_cb05 (icount) .NE. non_react ) THEN @@ -1400,7 +1400,7 @@ SUBROUTINE bio_emissions_megan2(id,config_flags,ktau,dtstep, & END DO - CASE (RACM_SOA_VBS_KPP) + CASE (RACM_SOA_VBS_KPP,RACM_SOA_VBS_AQCHEM_KPP) DO icount = 1, n_megan2racmSOA