Skip to content

*Fixes to get clamped dumbbell cases to run more correctly.#738

Merged
adcroft merged 3 commits into
mom-ocean:dev/gfdlfrom
ESMG:dev/esmg
Mar 30, 2018
Merged

*Fixes to get clamped dumbbell cases to run more correctly.#738
adcroft merged 3 commits into
mom-ocean:dev/gfdlfrom
ESMG:dev/esmg

Conversation

@kshedstrom
Copy link
Copy Markdown
Collaborator

  • normal_trans needs to include appropriate dx or dy.
  • Remove the masking on specified OBC segments (required when reading in
    array values from files).
  • Remove unused dx_Cv_obc and dy_CU_obc variables.
  • Clean up continuity for cases with more than one kind of OBC.
  • Changes answers on "simple" OBC cases.

- normal_trans needs to include appropriate dx or dy.
- Remove the masking on specified OBC segments (required when reading in
  array values from files).
- Remove unused dx_Cv_obc and dy_CU_obc variables.
- Clean up continuity for cases with more than one kind of OBC.
- Changes answers on "simple" OBC cases.
@MJHarrison-GFDL
Copy link
Copy Markdown
Contributor

Conflict with SIS2

src/SIS_transcribe_grid.F90:84.36:

SG%dy_Cu_obc(I,j) = dG%dy_Cu_obc(I+ido,j+jdo)
                                1

Error: 'dy_cu_obc' at (1) is not a member of the 'dyn_horgrid_type' structure

@kshedstrom
Copy link
Copy Markdown
Collaborator Author

Well, SIS2 isn't using those either. I'll do a pull request there as soon as my current CCS run finishes.

@adcroft
Copy link
Copy Markdown
Collaborator

adcroft commented Mar 30, 2018

This has to be handled manually because of the implied SIS2 update. Needless to say the automated tests failed on SIS2: https://gitlab.gfdl.noaa.gov/ogrp/MOM6/pipelines/3646

@adcroft adcroft merged commit f9b30b3 into mom-ocean:dev/gfdl Mar 30, 2018
OlgaSergienko pushed a commit to OlgaSergienko/MOM6 that referenced this pull request Nov 18, 2024
Added SQG vertical structure in Varmix to provide vertical profile for diffusivities
- added function calc_sqg_struct in MOM_lateral_mixing_coeffs to compute sqg_struct
- added sqg_expo to set the exponent of sqg_struct
- to use sqg_struct for the backscatter, set BS_use_sqg=true, sqg_expo>0., and BS_EBT_power=0.
- if SQG_USE_MEKE=True, use the eddy length scale from MEKE to compute sqg_struct
- added eddy length scale Le in MEKE if SQG_USE_MEKE=True
- added MEKE%Le into restart file if SQG_USE_MEKE=True
- added MEKE in Varmix
- registered N2_u and N2_v diagnostics when SQG_EXPO>0

(cherry picked from commit 6d3df0541c33d6f6d1f9fcb695f1a1eb961ec1b3)

* Compute vertical structures for khth, khtr, backscatter, and kdgl90 all in VarMix

- Vertical structures including khth_struct, khtr_struct, BS_struct, and kdgl90_struct are now computed in VarMix
- Each diffusivity/viscosity have two vertical structure options, equivalent barotropic (EBT) and surface quasigeostrophic (SQG) mode structures
- KHTH_USE_EBT_STRUCT, KHTR_USE_EBT_STRUCT, KDGL90_USE_EBT_STRUCT and BS_EBT_POWER parameters, which already existed, still control whether to use the EBT structure for khth, khtr, kdgl90, and backscatter, respectively
- Added KHTH_USE_SQG_STRUCT, KHTR_USE_SQG_STRUCT, KDGL90_USE_SQG_STRUCT and BS_USE_SQG parameters to control whether to use the SQG structure for khth, khtr, kdgl90, and backscatter, respectively
- If neither EBT nor SQG is called, no vertical structure will be used for that diffusivity/viscosity
- An error will be called if both EBT and SQG structures are called for the same diffusivity/viscosity
- Added dt as an input of calc_resoln_function. dt is needed for calc_sqg_struct called in calc_resoln_function

Bug fixes
- Changed ()**0.5 to sqrt in sqg_struct, which solves the dimension error
- Added if statement for using BS_struct in MOM_hor_visc
- Added if statement for CS%sqg_struct<=0.
- Changed the subroundoff of Coriolis parameter to 1e-40 in calc_sqg_struct
- Changed parameter BS_USE_SQG to BS_USE_SQG_STRUCT

---------

Co-authored-by: Alistair Adcroft <adcroft@users.noreply.github.com>
alex-huth pushed a commit to alex-huth/MOM6 that referenced this pull request Apr 14, 2025
Added SQG vertical structure in Varmix to provide vertical profile for diffusivities
- added function calc_sqg_struct in MOM_lateral_mixing_coeffs to compute sqg_struct
- added sqg_expo to set the exponent of sqg_struct
- to use sqg_struct for the backscatter, set BS_use_sqg=true, sqg_expo>0., and BS_EBT_power=0.
- if SQG_USE_MEKE=True, use the eddy length scale from MEKE to compute sqg_struct
- added eddy length scale Le in MEKE if SQG_USE_MEKE=True
- added MEKE%Le into restart file if SQG_USE_MEKE=True
- added MEKE in Varmix
- registered N2_u and N2_v diagnostics when SQG_EXPO>0

(cherry picked from commit 6d3df0541c33d6f6d1f9fcb695f1a1eb961ec1b3)

* Compute vertical structures for khth, khtr, backscatter, and kdgl90 all in VarMix

- Vertical structures including khth_struct, khtr_struct, BS_struct, and kdgl90_struct are now computed in VarMix
- Each diffusivity/viscosity have two vertical structure options, equivalent barotropic (EBT) and surface quasigeostrophic (SQG) mode structures
- KHTH_USE_EBT_STRUCT, KHTR_USE_EBT_STRUCT, KDGL90_USE_EBT_STRUCT and BS_EBT_POWER parameters, which already existed, still control whether to use the EBT structure for khth, khtr, kdgl90, and backscatter, respectively
- Added KHTH_USE_SQG_STRUCT, KHTR_USE_SQG_STRUCT, KDGL90_USE_SQG_STRUCT and BS_USE_SQG parameters to control whether to use the SQG structure for khth, khtr, kdgl90, and backscatter, respectively
- If neither EBT nor SQG is called, no vertical structure will be used for that diffusivity/viscosity
- An error will be called if both EBT and SQG structures are called for the same diffusivity/viscosity
- Added dt as an input of calc_resoln_function. dt is needed for calc_sqg_struct called in calc_resoln_function

Bug fixes
- Changed ()**0.5 to sqrt in sqg_struct, which solves the dimension error
- Added if statement for using BS_struct in MOM_hor_visc
- Added if statement for CS%sqg_struct<=0.
- Changed the subroundoff of Coriolis parameter to 1e-40 in calc_sqg_struct
- Changed parameter BS_USE_SQG to BS_USE_SQG_STRUCT

---------

Co-authored-by: Alistair Adcroft <adcroft@users.noreply.github.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants