Opened 12 years ago

Closed 5 years ago

Last modified 3 years ago

#337 closed Enhancement (fixed)

partial cells and vvl

Reported by: jchanut Owned by: nemo
Priority: low Milestone:
Component: OCE Version: release-3.0
Severity: Keywords: VVL release-3.0
Cc: Review:
MP ready?:
Progress:

Description

Here are a couple of changes I have made to made both key_vvl and partial cells compatible.
The idea is still to perform a vertical interpolation of tracers at the last level as in the constant volume case. My feeling was to perform this interpolation not on horizontal surfaces but on z* surfaces (z*(i,j,k)=e3t_0(k)+sshn(i,j), this is my feeling so it is questionable). One still has to add the Jacobian correction (depending on the sla gradient) everywhere including at the bottom.
There are certainly other ways to do this.

I created a new routine in dynhpg.F90: hpg_zps_vvl dedicated to that special case (key_vvl, ln_zps=.true.). It should however work in the more general case where ln_zco or ln_zps=.true., with or without key_vvl.
One more thing: this only works with key_sigma_vvl and bugs on the vertical scale factors definition corrected in domvvl (see Gurvan's tickets).

Jerome

Commit History (0)

(No commits)

Attachments (1)

PC-VVL.tar (100.0 KB) - added by jchanut 12 years ago.
Patch VVL-Partial cells

Download all attachments as: .zip

Change History (8)

Changed 12 years ago by jchanut

Patch VVL-Partial cells

comment:1 Changed 12 years ago by jchanut

The definition of z* surfaces I gave is wrong. This is rather:
z*(i,j,k)=gdept_0(k)*(1.+sshn(i,j)/H(i,j)), where H(i,j)=1./mut(i,j,:) the total depth with key_sigma_vvl.

comment:2 Changed 5 years ago by clevy

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

comment:3 Changed 4 years ago by nicolasmartin

  • Keywords VVL added; vvl removed

comment:4 Changed 4 years ago by nicolasmartin

  • Keywords nemo_v3* added

comment:5 Changed 4 years ago by nicolasmartin

  • Keywords cells removed

comment:6 Changed 4 years ago by nicolasmartin

  • Keywords partial removed

comment:7 Changed 3 years ago by nemo

  • Keywords release-3.0 added; nemo_v3* removed
Note: See TracTickets for help on using tickets.