Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add temperature and salinity to ocean extrapolation #124

Open
wants to merge 8 commits into
base: develop
Choose a base branch
from
Open

Add temperature and salinity to ocean extrapolation #124

Show file tree
Hide file tree
Changes from 6 commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
0d7e8e0
Modified ocean extrapolation routine to take in a 3d ocean cavity mask
cshafer Aug 26, 2024
93c7ad7
Modified which ocean layers are picked for creating TFocean and added…
cshafer Aug 26, 2024
932736e
Revert "Modified which ocean layers are picked for creating TFocean a…
matthewhoffman Aug 28, 2024
5aa744c
Handle seafloor and lakes in extrapolation module
matthewhoffman Aug 29, 2024
75c1940
Replace verbose per-iteration log message with a summary
matthewhoffman Aug 29, 2024
bba3a42
Add salt + temp to ocean extrapolation
alexolinhager Sep 5, 2024
735add8
Add thermal forcing calculation routine
alexolinhager Sep 16, 2024
8edf178
2d from 3d fields moved to mpas_li_iceshelf_melt.F
alexolinhager Nov 4, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
32 changes: 18 additions & 14 deletions components/mpas-albany-landice/src/Registry.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -454,16 +454,16 @@
possible_values="any non-negative value"
/>
<nml_option name="config_ocean_data_extrapolation" type="logical" default_value=".false." units="unitless"
description="If true, extrapolate ocean data (temperature, salinity, thermal forcing) from external source into underneath the ice draft."
description="If true, extrapolate ocean data (temperature and salinity, or thermal forcing) from external source into underneath the ice draft."
possible_values=".true. or .false."
/>
<nml_option name="config_ocean_data_extrap_ncells_extra" type="integer" default_value="10" units="unitless"
/>
<nml_option name="config_recalculate_thermal_forcing" type="logical" defaualt_value=".false." units="unitless"
description="If true, ocean temperature and salinity are extrapolated under ice drafts and to glacier termini, where thermal forcing is then calculated according to options inconfig_thermal_forcing_param_type. If false, MALI will look for an external thermal forcing field to extrapolate instead, and the subsequent 2d thermal forcing field will be defined at the grounding line. Option is only valid if config_ocean_data_extrapolation is true."
possible_values=".true. or .false."
/>
<nml_option name="config_ocean_data_extrap_ncells_extra" type="integer" default_value="10" units="unitless"
description="number of extra cells for over-extrapolation into grounded ice" possible_values="any non-negative value"
/>
<nml_option name="config_invalid_value_TF" type="real" default_value="1.0e36" units="unitless"
description="value assigned to indicate invalid thermal forcing value when config_ocean_data_extrapolation is set to true"
possible_values="Any real value"
/>
<nml_option name="config_weight_value_cell" type="real" default_value="0.9" units="unitless"
description="weight used to smooth horizontal extrapolation of ocean data field. Value close to 1 implies more weight of the current cell value versus the averaged value of the neighbouring cells around the cell"
possible_values="any real value between 0 and 1"
Expand Down Expand Up @@ -810,6 +810,7 @@
<var name="ismip6Runoff" packages="ismip6GroundedFaceMelt" />
<var name="ismip6_2dThermalForcing" packages="ismip6GroundedFaceMelt" />
<var name="origOceanMaskHoriz" packages="extrapOceanData" />
<var name="orig3dOceanCavityMask" packages="extrapOceanData" /> <!-- CAS 8/16/2024 -->
</stream>

<!-- An alternate way to allow the HO variables to exist in a separate file.
Expand Down Expand Up @@ -915,15 +916,14 @@
<var name="ismip6_2dThermalForcingCurrent" packages="ismip6GroundedFaceMelt" />
<var name="forcingTimeStamp" packages="ismip6GroundedFaceMelt" />
<!-- these variables are just for the ocean thermal forcing re-extrapolation scheme -->
<var name="origOceanMaskHoriz" packages="extrapOceanData" />
<var name="origOceanMaskHoriz" packages="extrapOceanData" />
<var name="orig3dOceanCavityMask" packages="extrapOceanData" /> <!-- CAS 8/16/2024 -->
<var name="validOceanMask" packages="extrapOceanData" />
<var name="validOceanMaskOrig" packages="extrapOceanData" />
<var name="availOceanMask" packages="extrapOceanData" />
<var name="growOceanMaskHoriz" packages="extrapOceanData" />
<var name="seedOceanMaskHoriz" packages="extrapOceanData" />
<var name="seedOceanMaskHorizInit" packages="extrapOceanData" />
<var name="TFoceanOld" packages="extrapOceanData" />
<var name="TFocean" packages="extrapOceanData" />
</stream>


Expand Down Expand Up @@ -1689,6 +1689,10 @@ is the value of that variable from the *previous* time level!
<var name="origOceanMaskHoriz" type="integer" dimensions="nCells" units="none"
description="2D mask for original valid ocean data (e.g., ISMIP6 original TF data)"
/>
<!-- CAS 8/16/2024 -->
<var name="orig3dOceanCavityMask" type="integer" dimensions="nISMIP6OceanLayers nCells Time" units="none"
description="3D mask for original valid ocean data that includes cavities (e.g., SORRM original TF data)"
/>
<var name="validOceanMask" type="integer" dimensions="nISMIP6OceanLayers nCells Time" units="none"
description="3D mask for indicating where ocean data is present (i.e., seed mask)"
/>
Expand All @@ -1707,11 +1711,11 @@ is the value of that variable from the *previous* time level!
<var name="seedOceanMaskHorizInit" type="integer" dimensions="nCells Time" units="none"
description="initial 2D seedMask before the over-extrapolation, for benchmark/debug purpose"
/>
<var name="TFoceanOld" type="real" dimensions="nISMIP6OceanLayers nCells Time" units="deg. C"
default_value="1.0"
<var name="MPAS_3dOceanTemperature" type="real" dimensions="nISMIP6OceanLayers nCells Time" units="deg. C"
description="MPAS-Ocean potential temperature to be extrapolated to all grounding lines if config_recalcuate_thermal_forcing is true."
/>
<var name="TFocean" type="real" dimensions="nISMIP6OceanLayers nCells Time" units="deg. C"
default_value="0.0"
<var name="MPAS_3dOceanSalinity" type="real" dimensions="nISMIP6OceanLayers nCells Time" units="g/kg"
description="MPAS-Ocean salinity to be extrapolated to all grounding lines if config_recalcuate_thermal_forcing is true."
/>
</var_struct>
<!-- ================ -->
Expand Down
2 changes: 1 addition & 1 deletion components/mpas-albany-landice/src/mode_forward/Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -82,7 +82,7 @@ mpas_li_velocity_external.o: Interface_velocity_solver.o

mpas_li_bedtopo.o: mpas_li_advection.o

mpas_li_ocean_extrap.o:
mpas_li_ocean_extrap.o: mpas_li_calving.o

Interface_velocity_solver.o:

Expand Down
Loading