New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#1400 (No isopycnal mixing with vvl ?) – NEMO

Opened 8 years ago

Closed 8 years ago

Last modified 8 years ago

#1400 closed Bug (invalid)

No isopycnal mixing with vvl ?

Reported by: jchanut Owned by: rfurner??
Priority: low Milestone:
Component: OCE Version: trunk
Severity: Keywords:
Cc:

Description

In ldfslp.F90, revision 4488 added the following condition

IF ( ln_traldf_iso .OR. ln_dynldf_iso ) THEN
...
compute the isopycnal slopes

ELSEIF ( lk_vvl ) THEN 
! set the slope of diffusion to the slope of s-surfaces 
... set isopycnal slopes to horizontal slopes

ENDIF

It means that with key_vvl isopycnal mixing reverts to horizontal mixing. Any reason for this ?

Commit History (0)

(No commits)

Change History (3)

comment:1 Changed 8 years ago by clevy

  • Owner changed from NEMO team to rfurner??

comment:2 Changed 8 years ago by timgraham

  • Resolution set to invalid
  • Status changed from new to closed

I've discussed this with Rachel and we don't think there is a problem.

The first part of the statement (compute the isopycnal slopes) will always be done if

IF ( ln_traldf_iso .OR. ln_dynldf_iso ) THEN

is true regardless of whether lk_vvl is set.

The part after the ELSEIF is only used if both ln_traldf_iso and ln_dynldf_iso are false and lk_vvl is true. This is used in the case of S-coordinates and geopotential mixing with VVL (without VVL these are set as part of the init stage so they don't need to be amended.

Also note that if using a mix of geopotential and isoneutral (one for tracer and one for momentum) then the isoneutral slopes are used for the first mixing, then the slopes are over written with the geopotential at the start of the module for the next mixing, and then set back to isoneutral with this routine at the start of the next time step (you don’t need to be able to do the vice versa, geopotential then isoneutral, as Gurvan doesn’t think you should need to and so it calls stop if you try!)

comment:3 Changed 8 years ago by jchanut

You are totally right. My Mistake, sorry.
I think I got mistaken by the statement ' Horizontal mixing in s-coordinate: slope = slope of s-surfaces' in the initialization part of my ocean.output

Note: See TracTickets for help on using tickets.