...
- Copy seq_drydep_mod.F90 from $CESMROOT/cime/src/drivers/mct/shr/ to your $CASEROOT/SourceMods/src.share directory
- Edit the fortran code in seq_drydep_mod.F90. The example below maps tracer XHNO3 to HNO3. This only works if the "test" species is also simulated.
case( 'XHNO3' )
test_name = 'HNO3'
...
- integer :: new_species_ndx
- logical :: new_species_dd
...
'
...
- Add new species into drydep_list in user_nl_cam. Note: if the default drydep_list is not already in your user_nl_cam file, you will need to build your case, copy the default drydep_list from $CASEROOT/CaseDocs/drv_flds_in, and then add your new species to the end of this list.
...
- Copy seq_drydep_mod.F90 from $CESMROOT/cime/src/drivers/mct/shr/ to your $CASEROOT/SourceMods/src.share directory.
- Edit the fortran code in seq_drydep_mod.F90.
- In this code, there are several arrays containing: 1) species names, 2) reactivity factors (f0), 3) henry's law constants, and 4) molecular weights. Add your new species characteristics at the end of each of these arrays. (See updating Wet deposition for a larger description of how these henry's law constants are calculated.
Update the variable "maxspc" to be equal or greater than the total number of species you are dry depositing (i.e., the maximum number of species you will dry deposit).
- Update the variable "n_species_table" to the total number of species listed in these arrays.
- Copy mo_drydep.F90 from $CCSMROOT/components/cam/src/chemistry/mozart/ to your $CASEROOT/SourceMods/src.cam directory.
- Edit the fortran code in mo_drydep.F90
- .
- initialize your new variables
- integer :: new_species_ndx
- logical :: new_species_dd
In subroutine dvel_inti_xactive - new_species_ndx = get_spc_ndx( 'new_species' ) new_species_dd = has_drydep ( 'new_species' )
- Add new species into drydep_list in user_nl_cam. Note: if the default drydep_list is not already in your user_nl_cam file, you will need to build your case, copy the default drydep_list from $CASEROOT/CaseDocs/drv_flds_in, and then add your new species to the end of this list.
...