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.
domzgr_substitute.h90 in NEMO/trunk/src/SWE – NEMO

source: NEMO/trunk/src/SWE/domzgr_substitute.h90 @ 13237

Last change on this file since 13237 was 12983, checked in by techene, 4 years ago

SWE r12822 from Antoine updated with qco and rk3 options : draft 1

File size: 1.7 KB
Line 
1!!----------------------------------------------------------------------
2!!                    ***  domzgr_substitute.h90   ***
3!!----------------------------------------------------------------------
4!! ** purpose :   substitute fsdep. and fse.., the vert. depth and scale
5!!      factors depending on the vertical coord. used, using CPP macro.
6!!----------------------------------------------------------------------
7!! History :  4.2  !  2020-02  (S. Techene, G. Madec)  star coordinate
8!!----------------------------------------------------------------------
9!! NEMO/OCE 4.2 , NEMO Consortium (2020)
10!! $Id$
11!! Software governed by the CeCILL license (see ./LICENSE)
12!!----------------------------------------------------------------------
13#if defined key_qco
14#   define  e3t(i,j,k,t)   (e3t_0(i,j,k)*(1._wp+r3t(i,j,t)))
15#   define  e3u(i,j,k,t)   (e3u_0(i,j,k)*(1._wp+r3u(i,j,t)))
16#   define  e3v(i,j,k,t)   (e3v_0(i,j,k)*(1._wp+r3v(i,j,t)))
17#   define  e3f(i,j,k)     (e3f_0(i,j,k)*(1._wp+r3f(i,j)))
18#   define  e3w(i,j,k,t)   (e3w_0(i,j,k)*(1._wp+r3t(i,j,t)))
19#   define  e3uw(i,j,k,t)  (e3uw_0(i,j,k)*(1._wp+r3u(i,j,t)))
20#   define  e3vw(i,j,k,t)  (e3vw_0(i,j,k)*(1._wp+r3v(i,j,t)))
21#   define  ht(i,j)        (ht_0(i,j)+ssh(i,j,Kmm))
22#   define  hu(i,j,t)      (hu_0(i,j)*(1._wp+r3u(i,j,t)))
23#   define  hv(i,j,t)      (hv_0(i,j)*(1._wp+r3v(i,j,t)))
24#   define  r1_hu(i,j,t)   (r1_hu_0(i,j)/(1._wp+r3u(i,j,t)))
25#   define  r1_hv(i,j,t)   (r1_hv_0(i,j)/(1._wp+r3v(i,j,t)))
26#   define  gdept(i,j,k,t) (gdept_0(i,j,k)*(1._wp+r3t(i,j,t)))
27#   define  gdepw(i,j,k,t) (gdepw_0(i,j,k)*(1._wp+r3t(i,j,t)))
28#   define  gde3w(i,j,k)   (gdept_0(i,j,k)*(1._wp+r3t(i,j,Kmm))-ssh(i,j,Kmm))
29#endif
30!!----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.