Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

horizontal diffusivity depth dependent hack #8

Open
riquitorres opened this issue Jun 8, 2022 · 0 comments
Open

horizontal diffusivity depth dependent hack #8

riquitorres opened this issue Jun 8, 2022 · 0 comments

Comments

@riquitorres
Copy link
Contributor

In mod_fabm_3D.F in the advection subroutine around lines 1660 there is a modification of the horizontal viscosity that is dependent on depth. More specifically, it is reduced on the shelf (5% of original value) and increases linearly between 200-800 to its nominal value.

!*******************************Jadon*************************************
VISCOF=(FACT*0.5_SP*(VISCOFF(IA)*NN_HVC(IA)+VISCOFF(IB)*NN_HVC(IB)) +
FM1*0.5_SP*(NN_HVC(IA)+NN_HVC(IB)))
      if (min(h(ia),h(ib))>800.0) then
      elseif (min(h(ia),h(ib))<200.0) then
        VISCOF=VISCOF*0.05
      else
        VISCOF=VISCOF*(1-(800.0-min(h(ia),h(ib)))/600.0*0.95)
      endif
!*****************************************************************************

I believe this was done to solve instabilities issues in the South China Sea implementation shelf break for nutrients but this shouldn't be there by default... and it should certainly be documented otherwise. There is no equivalent treatment in the salinity advection code

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

No branches or pull requests

1 participant