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 tags/nemo_v3_2/nemo_v3_2/NEMO/OFF_SRC/DOM – NEMO

source: tags/nemo_v3_2/nemo_v3_2/NEMO/OFF_SRC/DOM/domzgr_substitute.h90 @ 1878

Last change on this file since 1878 was 1878, checked in by flavoni, 14 years ago

initial test for nemogcm
File size: 6.5 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  fsdept_b(i,j,k)  (fsdept_0(i,j,k)*(1.+sshb(i,j)*mut(i,j,k)))
49#   define  fsdepw_b(i,j,k)  (fsdepw_0(i,j,k)*(1.+sshb(i,j)*mut(i,j,k)))
50#   define  fsde3w_b(i,j,k)  (fsde3w_0(i,j,k)*(1.+sshb(i,j)*mut(i,j,k))-sshb(i,j))
51#   define  fse3t_b(i,j,k)   (fse3t_0(i,j,k)*(1.+sshb(i,j)*mut(i,j,k)))
52#   define  fse3u_b(i,j,k)   (fse3u_0(i,j,k)*(1.+sshu_b(i,j)*muu(i,j,k)))
53#   define  fse3v_b(i,j,k)   (fse3v_0(i,j,k)*(1.+sshv_b(i,j)*muv(i,j,k)))
54#   define  fse3f_b(i,j,k)   (fse3f_0(i,j,k)*(1.+sshf_b(i,j)*muf(i,j,k)))
55#   define  fse3w_b(i,j,k)   (fse3w_0(i,j,k)*(1.+sshb(i,j)*mut(i,j,k)))
56#   define  fse3uw_b(i,j,k)  (fse3uw_0(i,j,k)*(1.+sshu_b(i,j)*muu(i,j,k)))
57#   define  fse3vw_b(i,j,k)  (fse3vw_0(i,j,k)*(1.+sshv_b(i,j)*muv(i,j,k)))
58
59#   define  fsdept_n(i,j,k)  (fsdept_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
60#   define  fsdepw_n(i,j,k)  (fsdepw_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
61#   define  fsde3w_n(i,j,k)  (fsde3w_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k))-sshn(i,j))
62#   define  fse3t_n(i,j,k)   (fse3t_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
63#   define  fse3u_n(i,j,k)   (fse3u_0(i,j,k)*(1.+sshu_n(i,j)*muu(i,j,k)))
64#   define  fse3v_n(i,j,k)   (fse3v_0(i,j,k)*(1.+sshv_n(i,j)*muv(i,j,k)))
65#   define  fse3f_n(i,j,k)   (fse3f_0(i,j,k)*(1.+sshf_n(i,j)*muf(i,j,k)))
66#   define  fse3w_n(i,j,k)   (fse3w_0(i,j,k)*(1.+sshn(i,j)*mut(i,j,k)))
67#   define  fse3uw_n(i,j,k)  (fse3uw_0(i,j,k)*(1.+sshu_n(i,j)*muu(i,j,k)))
68#   define  fse3vw_n(i,j,k)  (fse3vw_0(i,j,k)*(1.+sshv_n(i,j)*muv(i,j,k)))
69
70#   define  fse3t_m(i,j,k)   (fse3t_0(i,j,k)*(1.+ssh_m(i,j)*mut(i,j,k)))
71
72#   define  fsdept_a(i,j,k)  (fsdept_0(i,j,k)*(1.+ssha(i,j)*mut(i,j,k)))
73#   define  fsdepw_a(i,j,k)  (fsdepw_0(i,j,k)*(1.+ssha(i,j)*mut(i,j,k)))
74#   define  fsde3w_a(i,j,k)  (fsde3w_0(i,j,k)*(1.+ssha(i,j)*mut(i,j,k))-ssha(i,j))
75#   define  fse3t_a(i,j,k)   (fse3t_0(i,j,k)*(1.+ssha(i,j)*mut(i,j,k)))
76#   define  fse3u_a(i,j,k)   (fse3u_0(i,j,k)*(1.+sshu_a(i,j)*muu(i,j,k)))
77#   define  fse3v_a(i,j,k)   (fse3v_0(i,j,k)*(1.+sshv_a(i,j)*muv(i,j,k)))
78#   define  fse3f_a(i,j,k)   (fse3f_0(i,j,k)*(1.+sshf_a(i,j)*muf(i,j,k)))
79#   define  fse3w_a(i,j,k)   (fse3w_0(i,j,k)*(1.+ssha(i,j)*mut(i,j,k)))
80#   define  fse3uw_a(i,j,k)  (fse3uw_0(i,j,k)*(1.+sshu_a(i,j)*muu(i,j,k)))
81#   define  fse3vw_a(i,j,k)  (fse3vw_0(i,j,k)*(1.+sshv_a(i,j)*muv(i,j,k)))
82
83#else
84! z- or s-coordinate (1D or 3D + no time dependency) use reference in all cases
85#   define  fsdept(i,j,k)  fsdept_0(i,j,k)
86#   define  fsdepw(i,j,k)  fsdepw_0(i,j,k)
87#   define  fsde3w(i,j,k)  fsde3w_0(i,j,k)
88#   define  fse3t(i,j,k)   fse3t_0(i,j,k)
89#   define  fse3u(i,j,k)   fse3u_0(i,j,k)
90#   define  fse3v(i,j,k)   fse3v_0(i,j,k)
91#   define  fse3f(i,j,k)   fse3f_0(i,j,k)
92#   define  fse3w(i,j,k)   fse3w_0(i,j,k)
93#   define  fse3uw(i,j,k)  fse3uw_0(i,j,k)
94#   define  fse3vw(i,j,k)  fse3vw_0(i,j,k)
95
96#   define  fsdept_b(i,j,k)  fsdept_0(i,j,k)
97#   define  fsdepw_b(i,j,k)  fsdepw_0(i,j,k)
98#   define  fsde3w_b(i,j,k)  fsde3w_0(i,j,k)
99#   define  fse3t_b(i,j,k)   fse3t_0(i,j,k)
100#   define  fse3u_b(i,j,k)   fse3u_0(i,j,k)
101#   define  fse3v_b(i,j,k)   fse3v_0(i,j,k)
102#   define  fse3f_b(i,j,k)   fse3f_0(i,j,k)
103#   define  fse3w_b(i,j,k)   fse3w_0(i,j,k)
104#   define  fse3uw_b(i,j,k)  fse3uw_0(i,j,k)
105#   define  fse3vw_b(i,j,k)  fse3vw_0(i,j,k)
106
107#   define  fsdept_n(i,j,k)  fsdept_0(i,j,k)
108#   define  fsdepw_n(i,j,k)  fsdepw_0(i,j,k)
109#   define  fsde3w_n(i,j,k)  fsde3w_0(i,j,k)
110#   define  fse3t_n(i,j,k)   fse3t_0(i,j,k)
111#   define  fse3u_n(i,j,k)   fse3u_0(i,j,k)
112#   define  fse3v_n(i,j,k)   fse3v_0(i,j,k)
113#   define  fse3f_n(i,j,k)   fse3f_0(i,j,k)
114#   define  fse3w_n(i,j,k)   fse3w_0(i,j,k)
115#   define  fse3uw_n(i,j,k)  fse3uw_0(i,j,k)
116#   define  fse3vw_n(i,j,k)  fse3vw_0(i,j,k)
117
118#   define  fse3t_m(i,j,k)   fse3t_0(i,j,k)
119
120#   define  fsdept_a(i,j,k)  fsdept_0(i,j,k)
121#   define  fsdepw_a(i,j,k)  fsdepw_0(i,j,k)
122#   define  fsde3w_a(i,j,k)  fsde3w_0(i,j,k)
123#   define  fse3t_a(i,j,k)   fse3t_0(i,j,k)
124#   define  fse3u_a(i,j,k)   fse3u_0(i,j,k)
125#   define  fse3v_a(i,j,k)   fse3v_0(i,j,k)
126#   define  fse3f_a(i,j,k)   fse3f_0(i,j,k)
127#   define  fse3w_a(i,j,k)   fse3w_0(i,j,k)
128#   define  fse3uw_a(i,j,k)  fse3uw_0(i,j,k)
129#   define  fse3vw_a(i,j,k)  fse3vw_0(i,j,k)
130#endif
131   !!----------------------------------------------------------------------
132   !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009)
133   !! $Id: domzgr_substitute.h90 1640 2009-10-07 13:37:08Z cetlod $
134   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
135   !!----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.