Fix TEMPO aerosol initialization #69
Merged
Conversation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
…AP or the monthly climatology data. * Adds a config variable, config_tempo_rap, to the core_init_atmosphere Registry that must be set to true in the init_atmosphere namelist to use RAP aerosols. * Adds a logical function that check for config_tempo_rap = true and non-zero values of nwfa and nifa from RAP. This function will control whether to use RAP aerosols or the monthly climatology.
information from the first guess field (e.g., RAP) instead of the monthly climatology. This aerosol information, nwfa and nifa, is needed for the aerosol-aware TEMPO microphysics.
Use_rap_aerosols will now check for existing aerosols in the initial conditions if lbc_state is not present. Otherwise, it will check for existing aerosols in lbc_state (lateral boundary condition data type)
can be added to IC/LBC files.
lbc_state. * For qv, the lbc index is changed from index_qv to index_lbc_qv and state is changed to lbc_state * This is done for all LBC variables * Fixes a bug that caused an LBC variable (nwfa) to populate the incorrect variable name (lbc_nccn) since index_nwfa is NOT equal to index_lbc_nwfa
|
Thanks @AndersJensen-NOAA for adding options for how to initialize water- and ice-friendly aerosols, and for debugging and testing. I just ran tests with the CONUS 3-km MPAS domain, and results are good. Aerosol initialization is working for all of the following situations: (1) Initial conditions produced from a parent model that provides aerosols (config_tempo_rap = true) |
AndersJensen-NOAA
changed the title
| @@ -771,8 +771,8 @@ subroutine init_atm_thompson_aerosols_lbc(timestamp,timestart,block,mesh,diag,st | |||
| call mpas_pool_get_dimension(mesh,'nVertLevels' ,nVertLevels ) | |||
| call mpas_pool_get_dimension(mesh,'nMonths' ,nMonths ) | |||
|
|
|||
| call mpas_pool_get_dimension(state,'index_nifa',index_nifa) | |||
There was a problem hiding this comment.
This is a change to NCAR's code. I emailed Laura Fowler about this change.
There was a problem hiding this comment.
Are you hoping to hear back before this is merged?
There was a problem hiding this comment.
This has now been merged into NCAR's 8.2.3:
MPAS-Dev#1249
jderrico-noaa
requested review from
barlage,
joeolson42,
tanyasmirnova and
mdtoyNOAA
mdtoyNOAA
approved these changes
Nov 15, 2024
joeolson42
approved these changes
Nov 16, 2024
ghost
pushed a commit
that referenced
this pull request
Feb 11, 2025
* Adds code to determine whether to use aerosols (nwfa and nifa) from RAP or the monthly climatology data. * Adds a config variable, config_tempo_rap, to the core_init_atmosphere Registry that must be set to true in the init_atmosphere namelist to use RAP aerosols. * Adds a logical function that check for config_tempo_rap = true and non-zero values of nwfa and nifa from RAP. This function will control whether to use RAP aerosols or the monthly climatology. * Modifies logic for creation of ICs/LBCs so that users can use aerosol information from the first guess field (e.g., RAP) instead of the monthly climatology. This aerosol information, nwfa and nifa, is needed for the aerosol-aware TEMPO microphysics. * Adds optional argument lbc_state to the use_rap_aerosols function Use_rap_aerosols will now check for existing aerosols in the initial conditions if lbc_state is not present. Otherwise, it will check for existing aerosols in lbc_state (lateral boundary condition data type) * Fix missing parentheses in if statement * Initialize tempo_aerosolaware package so that nwfa and nifa can be added to IC/LBC files. * Changes lbc code to use lbc-specific indexes when adding variables to lbc_state. * For qv, the lbc index is changed from index_qv to index_lbc_qv and state is changed to lbc_state * This is done for all LBC variables * Fixes a bug that caused an LBC variable (nwfa) to populate the incorrect variable name (lbc_nccn) since index_nwfa is NOT equal to index_lbc_nwfa * Fix typos
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR fixes issue where initial water- and ice-friendly aerosol values from the RAP model are being overwritten by the monthly climatology. This PR adds logic to allow users to use either initial aerosols from a first guess model field if they exist or values from a monthly climatology.
To use RAP aerosol initial conditions, set
config_tempo_rap = .true.in preproc_stages. If first guess values do not exist (global max values of zero), then the monthly climatology will be used as a backup.