search

geoschem / geos-chem

GEOS-Chem "Science Codebase" repository. Contains GEOS-Chem science routines, run directory generation scripts, and interface code. This repository is used as a submodule within the GCClassic and GCHP wrappers, as well as in other modeling contexts (external ESMs).
http://geos-chem.org
Other
167 stars 163 forks source link

Number concentration and mass concentration #2432

Closed ChenBHXMU closed 2 months ago

ChenBHXMU commented 2 months ago

Your name

Baihua Chen

Your affiliation

Xiamen University

Please provide a clear and concise description of your question or discussion topic.

Dear @yantosca ,

I hope this message finds you well. While reading hcox_seasalt_mod.F90, I came across this piece of code. I think that SALT and SALT_N should represent mass concentration and number concentration, respectively. I have two question: 1. Will both types of concentration be saved? 2. Which type of concentration is recorded in the GEOSChem.Metrics.xxxxx_xxxxz.nc4 file in the OutputDir folder?

Could you kindly provide guidance or recommend any resources that could help me with this process?

Thank you very much for your assistance.

Best regards, Baihua Chen

图片1
yantosca commented 2 months ago

Thanks for writing @ChenBHXMU.

  1. I believe the SALT_N fields will be saved to HEMCO diagnostics if the marinePOA option is activated:

    ! Create marine density diagnostics only if marine POA enabled
IF ( HcoState%MarinePOA ) THEN

   CALL Diagn_Create ( HcoState   = HcoState,              &
                       cName      = 'SEASALT_DENS_PHOBIC', &
                       ExtNr      = Inst%ExtNrSS,         &
                       Cat        = -1,                    &
                       Hier       = -1,                    &
                       HcoID      = Inst%IDTMOPO,          &
                       SpaceDim   = 2,                     &
                       OutUnit    = 'number_dens',         &
                       AutoFill   = 0,                     &
                       Trgt2D     = Inst%NDENS_MOPO,       &
                       COL = HcoState%Diagn%HcoDiagnIDManual, &
                       RC         = RC                      )
   IF ( RC /= HCO_SUCCESS ) THEN
       CALL HCO_ERROR( 'ERROR 20', RC, THISLOC=LOC )
       RETURN
   ENDIF

   CALL Diagn_Create ( HcoState   = HcoState,              &
                       cName      = 'SEASALT_DENS_PHILIC', &
                       ExtNr      = Inst%ExtNrSS,          &
                       Cat        = -1,                    &
                       Hier       = -1,                    &
                       HcoID      = Inst%IDTMOPI,          &
                       SpaceDim   = 2,                     &
                       OutUnit    = 'number_dens',         &
                       AutoFill   = 0,                     &
                       Trgt2D     = Inst%NDENS_MOPI,       &
                       COL = HcoState%Diagn%HcoDiagnIDManual, &
                       RC         = RC                      )
   IF ( RC /= HCO_SUCCESS ) THEN
       CALL HCO_ERROR( 'ERROR 21', RC, THISLOC=LOC )
       RETURN
   ENDIF

ENDIF
  2. The Metrics collection only saves gas-phase fields needed to compute the mean OH concentration, CH4 lifetime, and methyl chloroform lifetime.
ChenBHXMU commented 2 months ago

Thank you very much for your assistance. @yantosca