corrrect pool name in init_atm_thompson_aerosols_lbc #1249
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
mgduda
requested changes
Dec 13, 2024
There was a problem hiding this comment.
With the changes in this PR, I'm seeing the following error (with the ifx compiler, after changing -check all to -check bounds,pointers in the DEBUG flags):
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
libpthread-2.31.s 0000149FF33BA8C0 Unknown Unknown Unknown
init_atmosphere_m 00000000005E9DD6 init_atm_thompson 783 mpas_init_atm_thompson_aerosols.F
init_atmosphere_m 0000000000475A29 init_atm_setup_ca 334 mpas_init_atm_cases.F
init_atmosphere_m 0000000000471F4A init_atm_core_run 92 mpas_init_atm_core.F
init_atmosphere_m 00000000004110F3 mpas_run 416 mpas_subdriver.F
init_atmosphere_m 000000000040E2A7 mpas 20 mpas.F
init_atmosphere_m 000000000040E24D Unknown Unknown Unknown
libc-2.31.so 0000149FEF27329D __libc_start_main Unknown Unknown
init_atmosphere_m 000000000040E17A Unknown Unknown Unknown
If we switch to the lbc_state pool, perhaps index_nifa should be index_lbc_nifa, and similarly for nwfa.
…the dimension of index_nifa and index_nwfa from the pool lbc_state instead of state. Corrected the call to subroutine init_atm_thompson_aerosols_lbc accordingly. * In ./src/core_init_atmosphere/mpas_init_atm_thompson_aerosols.F, replaced the index name of nifa and nwfa from index_nifa (index_nwfa) to index_lbc_nifa (index_lbc_nwfa) from the pool lbc_state.
ldfowler58
force-pushed
the
NSF-MPAS/correct_thompson_aerosols_lbc
branch
from
907ec80 to
1720638
Compare
|
I first verified that I was going a similar error message when computing the lateral boundary conditions for nwfa and nifa. After switching to using the index index_lbc_nifa and index_lbc_nwfa, I was able to buid the lateral boundary condition files without further issues. I updated the sourcecode accordingly. |
mgduda
approved these changes
Dec 18, 2024
mgduda
added a commit
that referenced
this pull request
May 23, 2025
This merge addresses several issues in the MPAS-Atmosphere model and in the MPAS infrastructure. Specific changes include: * Correction of the pool from which lbc_scalar constituent indices are obtained in the init_atm_thompson_aerosols_lbc routine. Rather than obtaining index_nifa and index_nwfa from the state pool, the indices of lbc_nifa and lbc_nwfa should be obtained from the lbc_state pool. (PR #1249) * Correction to the computation of the soil temperature (TSLB) in the Noah-MP land surface scheme through the addition of initialization of the soil liquid water (SH2O) in the noahmp_init subroutine in module mpas_atmphys_lsm_noahmpinit.F prior to calling NoahmpInitMain. (PR #1244) * Correction of the units of the fields 'greenfrac', 'shdmin', 'shdmax', 'vegfra', and 'albedo12m' from "unitless" to "percent" in the init_atmosphere and atmosphere core Registry.xml files. Also, a correction to the spelling of 'greenness' in several places. (PR #1248) * Removal of a duplicate allocation of indexToEdgeID % array in the mpas_io_setup_edge_block_fields routine that was the source of a memory leak. (PR #1258) * Fix for a memory leak in mpas_block_creator_build_cell_halos by deallocating the cellLimitField field before the routine returns. (PR #1264) * Fix for a bug in the logic for determining when decompositions can be reused by the SMIOL library. In almost any practical situation, however, this bug created no issues. (PR #1288) * Changes in the init_atmosphere core to provide more reliable error messages in case config_nfglevels is not set to a value that is at least as large as the number of vertical levels in the first-guess intermediate file. (PR #1291) * Correction of the loop for Noah-MP snow initialization, capping snow water equivalent maximum at 2000 mm. (PR #1300) * Fix for a bug in the horizontal 2nd-order filter for the CAM upper absorbing layer, where the wrong level in the kdiff field was being used when enforcing a lower-bound on kdiff. This absorbing layer is active only when config_mpas_cam_coef > 0.0. (PR #1302) * Fix in the mountain wave idealized test case initialization when multiple MPI tasks are used. The 'xc' variable, which represents the center-point location of the mountain, was previously computed based on the maximum xCell values local to an MPI task, leading to inconsistent values on each MPI rank. By finding the maximum of xCell over all MPI ranks and ensuring that all MPI ranks use this global maximum, the terrain field is computed consistently between serial and parallel runs of the init_atmosphere_model program for the mountain wave test case (config_init_case = 6). (PR #1312) * Correction to the calculation of the 2-meter diagnostics (T2M, TH2M, and Q2) when using the Noah-MP land surface scheme. While the computation of 2-meter diagnostics is the same for Noah and Noah-MP over oceans, it is different between the two land surface schemes over land. In Noah-MP, the 2-meter diagnostics are weighted as functions of their respective diagnostics over bare soil and over vegetation. The updated diagnostics for Noah and Noah-MP are now computed in the new file mpas_atmphys_sfc_diagnostics.F. (PR #1242) * Fix to provide consistency in the ringing behavior of recurring alarms after their reference time has been adjusted with a call to mpas_adjust_alarm_to_reference_time. Now, adjusting the reference time for an alarm will always leave that alarm in a state such that it is considered by the mpas_is_alarm_ringing routine to be ringing at the current time. With this fix, limited-area simulations can be restarted at times between LBC updates, provided the reference_time attribute for the 'lbc_in' stream is set to the simulation initial time in the streams.atmosphere file. (PR #1290). * Correction of an indexing error for rvcuten in code blocks specific to the Grell-Freitas scheme in the convection driver. Specifically, in the convection_from_MPAS and convection_to_MPAS routines, rvcuten used (k,k) as indexing in a loop, where (k,i) is needed. Since the Grell-Freitas scheme does not provide momentum tendencies, the changes in this merge have no impact on results. (PR #1283)
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.
In ./src/core_init_atmosphere/mpas_init_atm_thompson_aerosols.F, corrected the pool name to get the dimension for index_nifa and index_nwfa in subroutine init_atm_thompson_aerosols_lbc. In its original release, index_nifa and index_nwfa
was obtained from the pool "state' instead of the pool "lbc_state".
This PR replaces the pool name "state" with the pool name "lbc_state", and modified the list of arguments needed in subroutine init_atm_thompson_aerosols_lbc.