Zero'ing surface fluxes below tabular icebergs#388
Conversation
…the ocean below large icebergs (with an area greater than half of the grid cell). The SW, LW, evap, lprec and sensible fluxes are set to zero. The frazil is set to zero, and the appropiate amount of freshwater flux is added to evap. This was done to limit sea ice growth below tabular icebergs.
|
Will this option break conservation of heat for the climate model? When you zero out the fluxes to the ocean, are they instead passed to the icebergs? |
MJHarrison-GFDL
left a comment
There was a problem hiding this comment.
@StephenGriffies @sternalon Yes. Thislresults in a coupled model heat sink. If we want to modify the fluxes here, then we need to make the corresponding change to the flux on the ATM side of the coupler.
| "The viscosity of the icebergs", units="m2 s-1",default=1.0e10) | ||
| call get_param(param_file, mod, "DENSITY_ICEBERGS", OS%density_iceberg, & | ||
| "A typical density of icebergs.", units="kg m-3", default=917.0) | ||
| call get_param(param_file, mod, "Lat_fusion", OS%Lat_fusion, & |
There was a problem hiding this comment.
Run-time parameters are upper case by convention.
There was a problem hiding this comment.
Turns out we already have this parameter called LATENT_HEAT_FUSION so I'm using that...
|
|
||
| !Zero'ing out other fluxes under the tabular icebergs | ||
| do j=jsd,jed ; do i=isd,ied | ||
| if (fluxes%frac_shelf_h(i,j) > 0.5) then !Only applying for ice shelf covering most of cell |
There was a problem hiding this comment.
Will this affect any other experiments? If we make this "0.5" a run-time parameter we could at least turn off the zero'ing.
Removed the coord_slight module and all calls to it, and obsoleted all run-time parameters that are exclusively related to it. This code was an attempt from 2015 to define an appropriate hybrid vertical coordinate for global climate modeling, but it never worked very well (usually falling apart in the second year), and it has not been used in any publication or active model for many years. The test case that exercised this coordinate in the MOM6-examples test suite is also being removed via MOM6-examples PR mom-ocean#388. The coord_SLight code is being eliminated altogether now to simplify the MOM6 code base and reduce the volume of untested and unused code. All answers in all known MOM6 configurations in active use are bitwise identical, although there is a remote chance that someone somewhere might be using the SLIGHT coordinate.
Removed the coord_slight module and all calls to it, and obsoleted all run-time parameters that are exclusively related to it. This code was an attempt from 2015 to define an appropriate hybrid vertical coordinate for global climate modeling, but it never worked very well (usually falling apart in the second year), and it has not been used in any publication or active model for many years. The test case that exercised this coordinate in the MOM6-examples test suite is also being removed via MOM6-examples PR #388. The coord_SLight code is being eliminated altogether now to simplify the MOM6 code base and reduce the volume of untested and unused code. All answers in all known MOM6 configurations in active use are bitwise identical, although there is a remote chance that someone somewhere might be using the SLIGHT coordinate.
A piece of code has been added which zeros out the surface fluxes to the ocean below large icebergs (with an area greater than half of the grid cell).
The SW, LW, evap, lprec and sensible fluxes are set to zero.
The frazil is set to zero, and the appropiate amount of freshwater flux is added to evap.
This was done to limit sea ice growth below tabular icebergs.