From d07387b269436388ece09ffd7698c9f924139312 Mon Sep 17 00:00:00 2001 From: Andrew Kiss <31054815+aekiss@users.noreply.github.com> Date: Wed, 7 Aug 2024 15:34:06 +1000 Subject: [PATCH] use DT_THERM=6*DT and DIABATIC_FIRST = False - https://github.com/COSIMA/ACCESS-OM3/issues/138 Revert to 1mo run --- MOM_input | 101 +++++++++++++++++++----------------------------- nuopc.runconfig | 2 +- nuopc.runseq | 2 +- 3 files changed, 42 insertions(+), 63 deletions(-) diff --git a/MOM_input b/MOM_input index a0f67aa..c140eb6 100644 --- a/MOM_input +++ b/MOM_input @@ -1,66 +1,40 @@ ! This input file provides the adjustable run-time parameters for version 6 of the Modular Ocean Model (MOM6). ! === module MOM === -DT = 1800.0 - ! "[s] - ! The (baroclinic) dynamics time step. The time-step that - ! is actually used will be an integer fraction of the - ! forcing time-step (DT_FORCING in ocean-only mode or the - ! coupling timestep in coupled mode.)" - -DT_THERM = 3600.0 - ! "[s] default = 3600.0 - ! The thermodynamic and tracer advection time step. - ! Ideally DT_THERM should be an integer multiple of DT - ! and less than the forcing or coupling time-step. - ! By default DT_THERM is set to DT. - ! The minimum amount of time in seconds between - ! calculations of depth-space diagnostics. Making this - ! larger than DT_THERM reduces the performance penalty - ! of regridding to depth online." - -HFREEZE = 10.0 - ! "[m] default = -1.0 - ! If HFREEZE > 0, melt potential will be computed. The actual depth - ! over which melt potential is computed will be min(HFREEZE, OBLD), - ! where OBLD is the boundary layer depth. If HFREEZE <= 0 (default), - ! melt potential will not be computed." - -C_P = 3992.0 - ! "[J kg-1 K-1] default = 3991.86795711963 - ! The heat capacity of sea water, approximated as a - ! constant. This is only used if ENABLE_THERMODYNAMICS is - ! true. The default value is from the TEOS-10 definition - ! of conservative temperature." - -USE_REGRIDDING = True - ! "[Boolean] default = False - ! If True, use the ALE algorithm (regridding/remapping). - ! If False, use the layered isopycnal algorithm." - -THICKNESSDIFFUSE = True - ! "[Boolean] default = False - ! If true, interfaces or isopycnal surfaces are diffused, - ! depending on the value of FULL_THICKNESSDIFFUSE." - -THICKNESSDIFFUSE_FIRST = True - ! "[Boolean] default = False - ! If true, do thickness diffusion before dynamics. - ! This is only used if THICKNESSDIFFUSE is true." - -DTBT_RESET_PERIOD = 0.0 - ! "[s] default = 7200.0 - ! The period between recalculations of DTBT (if DTBT <= 0). - ! If DTBT_RESET_PERIOD is negative, DTBT is set based - ! only on information available at initialization. If 0, - ! DTBT will be set every dynamics time step. The default - ! is set by DT_THERM. This is only used if SPLIT is true." - -FRAZIL = True - ! "[Boolean] default = False - ! If true, water freezes if it gets too cold, and the - ! the accumulated heat deficit is returned in the - ! surface state. FRAZIL is only used if +DIABATIC_FIRST = False ! [Boolean] default = False + ! If true, apply diabatic and thermodynamic processes, including buoyancy + ! forcing and mass gain or loss, before stepping the dynamics forward. +USE_REGRIDDING = True ! [Boolean] default = False + ! If True, use the ALE algorithm (regridding/remapping). If False, use the + ! layered isopycnal algorithm. +THICKNESSDIFFUSE = True ! [Boolean] default = False + ! If true, isopycnal surfaces are diffused with a Laplacian coefficient of KHTH. +THICKNESSDIFFUSE_FIRST = True ! [Boolean] default = False + ! If true, do thickness diffusion or interface height smoothing before dynamics. + ! This is only used if THICKNESSDIFFUSE or APPLY_INTERFACE_FILTER is true. +DT = 1800.0 ! [s] + ! The (baroclinic) dynamics time step. The time-step that is actually used will + ! be an integer fraction of the forcing time-step (DT_FORCING in ocean-only mode + ! or the coupling timestep in coupled mode.) +DT_THERM = 10800.0 ! [s] default = 1800.0 + ! The thermodynamic and tracer advection time step. Ideally DT_THERM should be + ! an integer multiple of DT and less than the forcing or coupling time-step, + ! unless THERMO_SPANS_COUPLING is true, in which case DT_THERM can be an integer + ! multiple of the coupling timestep. By default DT_THERM is set to DT. +THERMO_SPANS_COUPLING = True +HFREEZE = 10.0 ! [m] default = -1.0 + ! If HFREEZE > 0, melt potential will be computed. The actual depth over which + ! melt potential is computed will be min(HFREEZE, OBLD), where OBLD is the + ! boundary layer depth. If HFREEZE <= 0 (default), melt potential will not be + ! computed. +DTBT_RESET_PERIOD = 10800.0 ! [s] default = 3600.0 + ! The period between recalculations of DTBT (if DTBT <= 0). If DTBT_RESET_PERIOD + ! is negative, DTBT is set based only on information available at + ! initialization. If 0, DTBT will be set every dynamics time step. The default + ! is set by DT_THERM. This is only used if SPLIT is true. +FRAZIL = True ! [Boolean] default = False + ! If true, water freezes if it gets too cold, and the accumulated heat deficit + ! is returned in the surface state. FRAZIL is only used if ! ENABLE_THERMODYNAMICS is true. ! If true, apply geothermal heating." @@ -74,7 +48,12 @@ BOUND_SALINITY = True ! If true, limit salinity to being positive. (The sea-ice ! model may ask for more salt than is available and ! drive the salinity negative otherwise.)" - +C_P = 3992.0 + ! "[J kg-1 K-1] default = 3991.86795711963 + ! The heat capacity of sea water, approximated as a + ! constant. This is only used if ENABLE_THERMODYNAMICS is + ! true. The default value is from the TEOS-10 definition + ! of conservative temperature." USE_PSURF_IN_EOS = True ! "[Boolean] default = True ! If true, always include the surface pressure contributions in equation of diff --git a/nuopc.runconfig b/nuopc.runconfig index 795f40f..c88d5dd 100644 --- a/nuopc.runconfig +++ b/nuopc.runconfig @@ -287,7 +287,7 @@ CLOCK_attributes:: history_ymd = -999 ice_cpl_dt = 99999 #not used lnd_cpl_dt = 99999 #not used - ocn_cpl_dt = 3600 + ocn_cpl_dt = 1800 #ignored (coupling timestep set by nuopc.runseq) unless stop_option is nsteps restart_n = 1 restart_option = ndays restart_ymd = -999 diff --git a/nuopc.runseq b/nuopc.runseq index b3138da..abcbf4b 100644 --- a/nuopc.runseq +++ b/nuopc.runseq @@ -1,5 +1,5 @@ runSeq:: -@3600 +@1800 MED med_phases_aofluxes_run MED med_phases_prep_ocn_accum MED med_phases_ocnalb_run