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).
Provide a clear and concise overview of the new feature requested.
Currently dry deposition velocities calculated in HEMCO for seaflux species like ACET aren’t output with the DryDep diagnostic. The values over the ocean are all zero. Velocities can be output with HEMCO diagnostics, but the values being output are odd (orders of magnitude from what's expected), suggesting a unit conversion error perhaps.
To address absence of seaflux species dry dep velocities over the ocean, I’ve simply moved where DryDepVel is archived from drydep_mod.F90 to hco_interface_gc_mod.F90. This places it in the same location as Archive_DryDep and Archive_DryDepMix. I’ve also moved the SatDiagnDryDepVel diagnostic. I haven’t moved the ConcAboveSurf drydepvel diagnostic, as I can’t quite make sense of what values it’s collecting in drydep_mod.F90 and I’m also not familiar with this diagnostic.
I’ve tested that the updated code works by comparing velocities for acetone (ACET) from the DryDep and SatDiagn collection (images below). In the original code, values are 0.1 cm/s over land and zero over the ocean, whereas in the updated code, the values calculated in HEMCO extension hcox_seaflux_mod.F90 are also output. I’ve also compared NO2 (images below) to confirm that the updated code doesn’t affect velocities of other species with no seaflux.
I’ve uploaded a patch that includes 2 code commits. The first updates the DryDep collection and the second updates the SatDiagn collection.
The patch also includes a very minor code update that changes indices in hco_interface_gc_mod.F90 from (:,:) to (I,J), as these occur within a lat, long loop.
Your name
Eloise Marais
Your affiliation
UCL, UK
Provide a clear and concise overview of the new feature requested.
Currently dry deposition velocities calculated in HEMCO for seaflux species like ACET aren’t output with the DryDep diagnostic. The values over the ocean are all zero. Velocities can be output with HEMCO diagnostics, but the values being output are odd (orders of magnitude from what's expected), suggesting a unit conversion error perhaps.
To address absence of seaflux species dry dep velocities over the ocean, I’ve simply moved where DryDepVel is archived from drydep_mod.F90 to hco_interface_gc_mod.F90. This places it in the same location as Archive_DryDep and Archive_DryDepMix. I’ve also moved the SatDiagnDryDepVel diagnostic. I haven’t moved the ConcAboveSurf drydepvel diagnostic, as I can’t quite make sense of what values it’s collecting in drydep_mod.F90 and I’m also not familiar with this diagnostic.
I’ve tested that the updated code works by comparing velocities for acetone (ACET) from the DryDep and SatDiagn collection (images below). In the original code, values are 0.1 cm/s over land and zero over the ocean, whereas in the updated code, the values calculated in HEMCO extension hcox_seaflux_mod.F90 are also output. I’ve also compared NO2 (images below) to confirm that the updated code doesn’t affect velocities of other species with no seaflux.
I’ve uploaded a patch that includes 2 code commits. The first updates the DryDep collection and the second updates the SatDiagn collection.
The patch also includes a very minor code update that changes indices in hco_interface_gc_mod.F90 from (:,:) to (I,J), as these occur within a lat, long loop.
Reference(s)
No response
Will you be implementing this feature yourself?
Yes
Additional information
Patch to apply changes that move drydepvel output from drydep_mod to hco_interface_gc_mod: move-drydepvel-drydep-satdiagn.patch.zip