geoschem / GCHP

The "superproject" wrapper repository for GCHP, the high-performance instance of the GEOS-Chem chemical-transport model.
https://gchp.readthedocs.io
Other
23 stars 27 forks source link

Transporting species initialized in a gridded component outside GeosChem #424

Closed ShreyaSharma2023 closed 4 months ago

ShreyaSharma2023 commented 4 months ago

Your name

Shreya Sharma

Your affiliation

MIT LAE

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

I am trying to advect species which have been initialized outside geoschem in another gridded component for MAM. For more information, refer to https://github.com/ShreyaSharma2023/GCHP.git (branch: one-way)

The bundles ([species]A[mode], for eg, SU_A_AIT) have been made 'Friendly to dynamics'. I have also added the species in geoschem_config.yaml and to species_database.yml in the run directory. However, I don't seem to see any change in the results before and after doing these steps, which makes me think the species are not getting transported.

Is there another way to have these species transported in GCHP?

lizziel commented 4 months ago

Hi Shreya, I think I found the missing step for connecting declared friendlies to advection. We have this code in GCHP_GridCompMod.F90:

    ! AdvCore Tracers                                                                                                     
    !----------------                                                                                                     
    call ESMF_StateGet( GIM(ADV), 'TRADV', BUNDLE, RC=STATUS )
    _VERIFY(STATUS)

    call MAPL_GridCompGetFriendlies(GCS(CHEM), "DYNAMICS", BUNDLE, RC=STATUS )
    _VERIFY(STATUS)

It looks like you need to add one of these calls for each of the gridded components you want to connect to advection.

ShreyaSharma2023 commented 4 months ago

Hi Lizzie Thanks! Adding this code seems to work. I do get differences in the output now, and the log file also lists all the species in the bundle as being advected. I assume I do not need to add the individual species in geoschem_config.yaml and to species_database.yml now? Because I added only one species to them, but I see changes in the values of all species in the bundle.

lizziel commented 4 months ago

Correct. These species are part of a different gridded component so should not need to be in any of the GEOS-Chem configuration files. I'm glad it works!