Skip to content

Static ice shelf initialize#3

Closed
MJHarrison-GFDL wants to merge 2 commits into
OlgaSergienko:ice_dynamicsfrom
MJHarrison-GFDL:static_ice_shelf_initialize
Closed

Static ice shelf initialize#3
MJHarrison-GFDL wants to merge 2 commits into
OlgaSergienko:ice_dynamicsfrom
MJHarrison-GFDL:static_ice_shelf_initialize

Conversation

@MJHarrison-GFDL
Copy link
Copy Markdown

Hydrostatic initialization under ice shelf cavities

  • Iteratively solve for the initial surface displacement under ice shelf cavities based on the applied pressure using the initial guess layered T/S profile above the level of the ice shelf displacement.
  • This option should deprecate the use of TRIM_IC_FOR_P_SURF
  • Turing off both TRIM_IC_FOR_P_SURF and DEPRESS_INITIAL_SURFACE to activate this option.
  • Existing ice shelf implementations (e.g. ISOMIP) should be reconfigured to use this option.
  • This option was tested using https://github.com/MJHarrison-GFDL/MOM6-IceShelf-examples/tree/47e13d7e9
    where the model was found to produce a more balanced initialization.

 - Iteratively solve for the initial ice shelf
   displacement in ice shelf cavities using
   the initial layer T/S values.
 - This change should make TRIM_IC_FOR_PSURF option
   unnecessay. Existing implementations (e.g. ISOMIP)
   should turn off this option in order to exercise
   this feature.
Hallberg-NOAA pushed a commit that referenced this pull request Dec 14, 2021
* reads in porous topography parameters from CHANNEL_LIST_FILE

*new module to compute curve fit for porous topography

*porous constraints used to modify continuity_PPM, CoriolisAdv, and Rayleigh bottom channel drag
Hallberg-NOAA added a commit that referenced this pull request Dec 14, 2021
(+) porous topography implementation
OlgaSergienko pushed a commit that referenced this pull request Jan 14, 2022
  Use the por_face_area[UV] in the effective thickness calculations in
zonal_face_thickness and merid_face_thickness, so that they are more consistent
with their use elsewhere in the code for the relative weights in calculating the
barotropic accelerations.  Because these por_face_area arrays are still 1 in all
test cases, the answers are unchanged in any test cases from before a few weeks
ago, but there could be answer changes in cases that are using the very recently
added capability (in PR #3) to set fractional face areas.  This change was
discussed with Sam Ditkovsky, and agreed that there is no reason to keep the
ability to recover the previous answers in any cases that use the recently added
partial face width option.

  This commit also expanded the comments describing the h_u and h_v arguments to
btcalc(), zonal_face_thickness(), and merid_face_thickness() routines, the
diag_h[uv] elements of the accel_diag_ptrs type and the h_u and h_v elements of
the BT_cont_type.

  All answers and output are bitwise identical in the MOM6-examples test suite
and TC tests, but answer changes are possible in cases using a very recently
added code option.
OlgaSergienko added a commit that referenced this pull request Feb 3, 2025
* Removed second registration of sfc_mass_flux with register_diag_field.

* Registered sfc_mass_flux to runs with static and dynamic ice sheets.
OlgaSergienko added a commit that referenced this pull request Feb 18, 2025
* Removed second registration of sfc_mass_flux with register_diag_field.

* Registered sfc_mass_flux to runs with static and dynamic ice sheets.
OlgaSergienko added a commit that referenced this pull request Mar 3, 2025
* Commit for coupling of land to ocean adot

* Remove unused ice sheet enabled flag

* add adot to Shelf restart

* Register sfc mass flux outside of initialize fluxes

* Gfdl cryo merge (#3)

* Removed second registration of sfc_mass_flux with register_diag_field.

* Registered sfc_mass_flux to runs with static and dynamic ice sheets.

* This commit removed the second registration of 'sfc_mass_flux' in 'register_diag_field()'

* Testing shelf_sfc_mass_flux restart.

* The earlier removed 'register_diag_field()' for 'mass_flux' has been added back and 'register_restart_field()' of 'fluxes_in%shelf_sfc_mass_flux' is moved before 'restore_state()'

* Corrected the line length error.

* Added back 'CS%id_mass_flux = register_diag_field()'

---------

Co-authored-by: NJSchlegel <nschlegel@gmail.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant