Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  • Copy seq_drydep_mod.F90 from  $CCSMROOT/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. The Henry's law constant table consists of the following 6 columns, respectively: KH298 (M atm-1), dH/R (K), K1298, dH1/R (K), K2298, dH2/R (K). These values are used to calculate the Henry's law constant (KH) and for acids and bases the effective Henry's law constant (Heff) using the following formulas:
    • 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.
  • Add new species into gas_wetdep_list in user_nl_cam. Note: if the default gas_wetdep_list is not already in your user_nl_cam file, you will need to build your case, copy the default gas_wetdep_list from $CASEROOT/CaseDocs/atm_in, and then add your new species to the end of this list.

...

  • Copy seq_drydep_mod.F90 from  $CCSMROOT/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.

...