Opened 3 years ago

Last modified 2 years ago

#1759 assigned Bug

diatmb.F90 : small error in z- and zps coordinate, wrong in s-coordinate

Reported by: gm Owned by: acc
Priority: low Milestone:
Component: OCE Version: trunk
Severity: Keywords: diatmb
Cc: acc



It is only a diagnostic issue which concern both v3.6_stable and trunk.
The top, middle and bottom output of tracer and horizontal velocity done diatmb.F90 module has several issues that can be rather easily solved (see below).
It should also be generalized to all available vertical coordinate and take into account ice-shelf cavities.

Analysis & Fix

1- Better to use mkbt which is directly the level needed (with 1 on land) instead of mbathy array.

2- small  BUG for u- and v-points mid and bottom calculation: the use of mbathy is wrong !
==⇒>>> mbku and mbkv should be used instead

3- mid depth diag. case: here only the z-coordinate is OK.

For zps there will be a small error at the boundary.
For s-coordinate, the diag is just wrong. Either add "IF(  ln_sco )   CALL ctl_stop(….)" in dia_tmb_init or adapt the diag to s-coordinate.

4- surface diag. case : in case of ice shelf cavities, the top of the ocean is not the level  k=1, but the bottom go the ice shelves. use mikt, miku or mikv (for t-, u-, or v-point)  

Commit History (0)

(No commits)

Change History (2)

comment:1 Changed 2 years ago by clevy

  • Cc acc added
  • Owner changed from nemo to acc
  • Status changed from new to assigned

comment:2 Changed 2 years ago by nicolasmartin

  • Version changed from nemo_v3_6_STABLE to trunk
Note: See TracTickets for help on using tickets.