Add melt/freeze potential via MCT cap#77
Merged
alperaltuntas merged 13 commits intoAug 8, 2018
Conversation
sfc_state%melt_potential is now accumulated over time_in_thermo_cycle.
Conflicts: src/core/MOM.F90 src/core/MOM_variables.F90
Member
|
There are answer changes in the following 6 fields. I think these changes are expected, but double checking just to make sure: o2x_Fioo_q |
Collaborator
Author
|
Thanks @alperaltuntas. Did you run both G and C compsets? I would expect these changes in the G compset, but not in the C compset. |
Member
|
Yes, the differences appears in C compset as well. It look like, in C compset, all six files used to be zero, but now are non-zero. |
alperaltuntas
pushed a commit
that referenced
this pull request
Apr 2, 2026
Updates for candidate 200918
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
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
2 participants
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 implements melt/freeze potential via MCT cap.
Both melt and freeze potential (defined as negative and positive, respectively) are passed to the coupler via variable o2x_Fioo_q, see type cpl_indices_type in MOM6/config_src/mct_driver/ocn_cpl_indices.F90. Freeze is computed using the frazil formation algorithm that was already implemented in MOM6 (see subroutine
make_frazilinMOM6/src/parameterizations/vertical/MOM_diabatic_aux.F90).Frazilis computed during the beginning and end of diabatic processes (seeMOM6/src/parameterizations/vertical/MOM_diabatic_driver.F90). However, unlike frazil, melt potential is computed insideextract_surface_state, which is a subroutine inMOM.F90. This subroutine is call at the end of subroutinestep_MOM, also inMOM.F90. Melt potential is defined as follows:Melt Potential = rho * c_p * Hmix * (T_m(1)-T_f(1)) / Delta_time,where Delta_time is the ocean coupled time step, T_f(1) is the potential temperature at the surface, rho is the mean ocean density, c_p is the heat capacity of sea water, and Hmix is the thickness over which melt potential is comuted (user control this via HMIX_SFC_PROP) . We have tested this using HMIX_SFC_PROP = dz(1) = 10m. Melt Potential is accumulated (cumulative sum) over CS%time_in_thermo_cycle, which is the ocean coupling time step. It is passed to the cap via the surface type in J/m^2. In subroutine ocn_export (MOM6/config_src/mct_driver/ocn_cap_methods.F90), melt potential is converted to W/m^2.
TODO: