Add tracers for freshwater sources to MPAS-Ocean #7087
Conversation
cbegeman
added
mpas-ocean
BFB
TestingThis branch has been tested again master 2481382 on chrys. The |
|
Further discussion is located at E3SM-Ocean-Discussion#108 |
|
@cbegeman, as we discussed, this will need a stealth test. Here are the existing ones: |
| @@ -282,6 +282,7 @@ | |||
| "SMS_D_Ld1.T62_oQU240wLI.GMPAS-IAF-PISMF.mpaso-impl_top_drag", | |||
| "SMS_D_Ld1.T62_oQU240.GMPAS-IAF.mpaso-harmonic_mean_drag", | |||
| "SMS_D_Ld1.T62_oQU240.GMPAS-IAF.mpaso-upwind_advection", | |||
| "SMS_D_Ld1.T62_oQU240.GMPAS-IAF.mpaso-freshwater_tracers", | |||
There was a problem hiding this comment.
It might be worth considering:
| "SMS_D_Ld1.T62_oQU240.GMPAS-IAF.mpaso-freshwater_tracers", | |
| "SMS_D_Ld1.T62_oQU240wLI.GMPAS-IAF-PISMF.mpaso-freshwater_tracers", |
to also test the ice-shelf-melt freshwater tracer. What do you think?
| else | ||
| call mpas_log_write('indexSnow is not associated') |
There was a problem hiding this comment.
I think this output (here and elsewhere in this file) would happen at every time step or at least every coupling time step and would be confusing to users. I think we probably want to remove it.
| else | ||
| call mpas_log_write('indexLandIceFreshwater or landIceFreshwaterFlux is not associated') | ||
| endif | ||
| call mpas_log_write('end fw assign sflux') |
There was a problem hiding this comment.
I think this should probably also be removed.
| !call mpas_pool_print_summary(tracersSurfaceFluxPool, MPAS_POOL_FIELD, .true.) | ||
|
|
|
|
||
| if ( groupName == 'freshwaterTracers' ) then | ||
| if (config_use_tracerGroup_surface_restoring) & | ||
| call mpas_log_write('surface restoring not supported for freshwaterTracers group') |
There was a problem hiding this comment.
I think these (here and below) should be critical errors, not just messages in the log file:
| call mpas_log_write('surface restoring not supported for freshwaterTracers group') | |
| call mpas_log_write('surface restoring not supported for freshwaterTracers group', messageType=MPAS_LOG_CRIT) |
This PR creates a new group of freshwater tracers, each associated with a source of fresh water mass to the ocean (including sea ice meltwater, which has a salinity > 0 PSU). All freshwater sources to the ocean are included. The value assigned to each tracers is its volumetric concentration.
All tracer tendency terms are included except the nonlocal flux term because it would require a new
nonLocalSurfaceTracerFluxstruct array to separate the salinity flux into its constituent mass source terms. That term is also not included for other passive tracers.These freshwater tracers do not have any sinks and thus their concentrations do not equilibrate. A possible extension of this work would be adding sinks proportional to freshwater tracer concentrations times negative mass fluxes at the ocean surface. This feature is not generally included in the literature so we have chosen not to include it here.
This feature is off by default. To activate freshwater tracers, set
config_use_freshwaterTracers=.true.. All freshwater tracers are activated together, but the user may choose to output only a subset of these tracers.[BFB][stealth]