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 branches/DEV_r1837_MLF/NEMO/OPA_SRC/DOM – NEMO

source: branches/DEV_r1837_MLF/NEMO/OPA_SRC/DOM/domzgr_substitute.h90 @ 2005

Last change on this file since 2005 was 2005, checked in by mlelod, 14 years ago

ticket: #663 MLF: second part (local compatibility essentially)

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 5.0 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 :  1.0  !  2005-10  (A. Beckmann, G. Madec) generalisation to all coord.
8   !!            3.1  !  2009-02  (G. Madec, M. Leclair)  pure z* coordinate
9   !!----------------------------------------------------------------------
10#if defined key_zco
11! reference for pure z-coordinate (1D - no i,j and time dependency)
12#   define  fsdept_0(i,j,k)  gdept_0(k)
13#   define  fsdepw_0(i,j,k)  gdepw_0(k)
14#   define  fsde3w_0(i,j,k)  gdepw_0(k)
15#   define  fse3t_0(i,j,k)   e3t_0(k)
16#   define  fse3u_0(i,j,k)   e3t_0(k)
17#   define  fse3v_0(i,j,k)   e3t_0(k)
18#   define  fse3f_0(i,j,k)   e3t_0(k)
19#   define  fse3w_0(i,j,k)   e3w_0(k)
20#   define  fse3uw_0(i,j,k)  e3w_0(k)
21#   define  fse3vw_0(i,j,k)  e3w_0(k)
22#else
23! reference for s- or zps-coordinate (3D no time dependency)
24#   define  fsdept_0(i,j,k)  gdept(i,j,k)
25#   define  fsdepw_0(i,j,k)  gdepw(i,j,k)
26#   define  fsde3w_0(i,j,k)  gdep3w(i,j,k)
27#   define  fse3t_0(i,j,k)   e3t(i,j,k)
28#   define  fse3u_0(i,j,k)   e3u(i,j,k)
29#   define  fse3v_0(i,j,k)   e3v(i,j,k)
30#   define  fse3f_0(i,j,k)   e3f(i,j,k)
31#   define  fse3w_0(i,j,k)   e3w(i,j,k)
32#   define  fse3uw_0(i,j,k)  e3uw(i,j,k)
33#   define  fse3vw_0(i,j,k)  e3vw(i,j,k)
34#endif
35#if defined key_vvl
36! s* or z*-coordinate (3D + time dependency) + use of additional now arrays (..._1)
37#   define  fsdept(i,j,k)  gdept_1(i,j,k)
38#   define  fsdepw(i,j,k)  gdepw_1(i,j,k)
39#   define  fsde3w(i,j,k)  gdep3w_1(i,j,k)
40#   define  fse3t(i,j,k)   e3t_1(i,j,k)
41#   define  fse3u(i,j,k)   e3u_1(i,j,k)
42#   define  fse3v(i,j,k)   e3v_1(i,j,k)
43#   define  fse3f(i,j,k)   e3f_1(i,j,k)
44#   define  fse3w(i,j,k)   e3w_1(i,j,k)
45#   define  fse3uw(i,j,k)  e3uw_1(i,j,k)
46#   define  fse3vw(i,j,k)  e3vw_1(i,j,k)
47
48#   define  fse3t_b(i,j,k)   e3t_b(i,j,k)
49#   define  fse3u_b(i,j,k)   e3u_b(i,j,k)
50#   define  fse3v_b(i,j,k)   e3v_b(i,j,k)
51#   define  fse3uw_b(i,j,k)  (fse3uw_0(i,j,k)*(1.+sshu_b(i,j)*muu(i,j,k)))
52#   define  fse3vw_b(i,j,k)  (fse3vw_0(i,j,k)*(1.+sshv_b(i,j)*muv(i,j,k)))
53
54#   define  fsdept_n(i,j,k)  (fsdept_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
55#   define  fsdepw_n(i,j,k)  (fsdepw_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
56#   define  fsde3w_n(i,j,k)  (fsde3w_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k))-sshn(i,j))
57#   define  fse3t_n(i,j,k)   (fse3t_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
58#   define  fse3u_n(i,j,k)   (fse3u_0(i,j,k)*(1.+sshu_n(i,j)*muu(i,j,k)))
59#   define  fse3v_n(i,j,k)   (fse3v_0(i,j,k)*(1.+sshv_n(i,j)*muv(i,j,k)))
60#   define  fse3f_n(i,j,k)   (fse3f_0(i,j,k)*(1.+sshf_n(i,j)*muf(i,j,k)))
61#   define  fse3w_n(i,j,k)   (fse3w_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
62#   define  fse3uw_n(i,j,k)  (fse3uw_0(i,j,k)*(1.+sshu_n(i,j)*muu(i,j,k)))
63#   define  fse3vw_n(i,j,k)  (fse3vw_0(i,j,k)*(1.+sshv_n(i,j)*muv(i,j,k)))
64
65#   define  fse3t_m(i,j,k)   e3t_m(i,j,k)
66
67#   define  fse3t_a(i,j,k)   (fse3t_0(i,j,k)*(1.+ssha(i,j)*mut(i,j,k)))
68#   define  fse3u_a(i,j,k)   (fse3u_0(i,j,k)*(1.+sshu_a(i,j)*muu(i,j,k)))
69#   define  fse3v_a(i,j,k)   (fse3v_0(i,j,k)*(1.+sshv_a(i,j)*muv(i,j,k)))
70
71#else
72! z- or s-coordinate (1D or 3D + no time dependency) use reference in all cases
73#   define  fsdept(i,j,k)  fsdept_0(i,j,k)
74#   define  fsdepw(i,j,k)  fsdepw_0(i,j,k)
75#   define  fsde3w(i,j,k)  fsde3w_0(i,j,k)
76#   define  fse3t(i,j,k)   fse3t_0(i,j,k)
77#   define  fse3u(i,j,k)   fse3u_0(i,j,k)
78#   define  fse3v(i,j,k)   fse3v_0(i,j,k)
79#   define  fse3f(i,j,k)   fse3f_0(i,j,k)
80#   define  fse3w(i,j,k)   fse3w_0(i,j,k)
81#   define  fse3uw(i,j,k)  fse3uw_0(i,j,k)
82#   define  fse3vw(i,j,k)  fse3vw_0(i,j,k)
83
84#   define  fse3t_b(i,j,k)   fse3t_0(i,j,k)
85#   define  fse3u_b(i,j,k)   fse3u_0(i,j,k)
86#   define  fse3v_b(i,j,k)   fse3v_0(i,j,k)
87#   define  fse3uw_b(i,j,k)  fse3uw_0(i,j,k)
88#   define  fse3vw_b(i,j,k)  fse3vw_0(i,j,k)
89
90#   define  fsdept_n(i,j,k)  fsdept_0(i,j,k)
91#   define  fsdepw_n(i,j,k)  fsdepw_0(i,j,k)
92#   define  fsde3w_n(i,j,k)  fsde3w_0(i,j,k)
93#   define  fse3t_n(i,j,k)   fse3t_0(i,j,k)
94#   define  fse3u_n(i,j,k)   fse3u_0(i,j,k)
95#   define  fse3v_n(i,j,k)   fse3v_0(i,j,k)
96#   define  fse3f_n(i,j,k)   fse3f_0(i,j,k)
97#   define  fse3w_n(i,j,k)   fse3w_0(i,j,k)
98#   define  fse3uw_n(i,j,k)  fse3uw_0(i,j,k)
99#   define  fse3vw_n(i,j,k)  fse3vw_0(i,j,k)
100
101#   define  fse3t_m(i,j,k)   fse3t_0(i,j,k)
102
103#   define  fse3t_a(i,j,k)   fse3t_0(i,j,k)
104#   define  fse3u_a(i,j,k)   fse3u_0(i,j,k)
105#   define  fse3v_a(i,j,k)   fse3v_0(i,j,k)
106#endif
107   !!----------------------------------------------------------------------
108   !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009)
109   !! $Id$
110   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
111   !!----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.