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 | #if defined key_s_coord |
---|
8 | !! s-coord: substitution fsdep.(,,) ==> hbat.(,) * gsig.() |
---|
9 | !! fse3.(,,) ==> hbat.(,) * esig.() |
---|
10 | # define fsdept(i,j,k) hbatt(i,j)*gsigt(k) |
---|
11 | # define fsdepu(i,j,k) hbatu(i,j)*gsigt(k) |
---|
12 | # define fsdepv(i,j,k) hbatv(i,j)*gsigt(k) |
---|
13 | # define fsdepf(i,j,k) hbatf(i,j)*gsigt(k) |
---|
14 | # define fsdepw(i,j,k) hbatt(i,j)*gsigw(k) |
---|
15 | # define fsdepuw(i,j,k) hbatu(i,j)*gsi3w(k) |
---|
16 | # define fsdepvw(i,j,k) hbatv(i,j)*gsi3w(k) |
---|
17 | # define fsde3w(i,j,k) hbatt(i,j)*gsi3w(k) |
---|
18 | # define fse3t(i,j,k) hbatt(i,j)*esigt(k) |
---|
19 | # define fse3u(i,j,k) hbatu(i,j)*esigt(k) |
---|
20 | # define fse3v(i,j,k) hbatv(i,j)*esigt(k) |
---|
21 | # define fse3f(i,j,k) hbatf(i,j)*esigt(k) |
---|
22 | # define fse3w(i,j,k) hbatt(i,j)*esigw(k) |
---|
23 | # define fse3uw(i,j,k) hbatu(i,j)*esigw(k) |
---|
24 | # define fse3vw(i,j,k) hbatv(i,j)*esigw(k) |
---|
25 | #elif defined key_partial_steps |
---|
26 | !! z-coord: substitution fsdep.(,,) ==> gdep._ps() |
---|
27 | !! fse3.(,,) ==> e3._ps() |
---|
28 | # define fsdept(i,j,k) gdept_ps(i,j,k) |
---|
29 | # define fsdepw(i,j,k) gdepw_ps(i,j,k) |
---|
30 | # define fsde3w(i,j,k) gdep3w(i,j,k) |
---|
31 | # define fsdepuw(i,j,k) gdepw(k) |
---|
32 | # define fsdepvw(i,j,k) gdepw(k) |
---|
33 | # define fse3t(i,j,k) e3t_ps(i,j,k) |
---|
34 | # define fse3u(i,j,k) e3u_ps(i,j,k) |
---|
35 | # define fse3v(i,j,k) e3v_ps(i,j,k) |
---|
36 | # define fse3f(i,j,k) e3f_ps(i,j,k) |
---|
37 | # define fse3w(i,j,k) e3w_ps(i,j,k) |
---|
38 | # define fse3uw(i,j,k) e3uw_ps(i,j,k) |
---|
39 | # define fse3vw(i,j,k) e3vw_ps(i,j,k) |
---|
40 | #else |
---|
41 | !! z-coord: substitution fsdep.(,,) ==> gdep() |
---|
42 | !! fse3.(,,) ==> e3.() |
---|
43 | # define fsdept(i,j,k) gdept(k) |
---|
44 | # define fsdepu(i,j,k) gdept(k) |
---|
45 | # define fsdepv(i,j,k) gdept(k) |
---|
46 | # define fsdepf(i,j,k) gdept(k) |
---|
47 | # define fsdepw(i,j,k) gdepw(k) |
---|
48 | # define fsdepuw(i,j,k) gdepw(k) |
---|
49 | # define fsdepvw(i,j,k) gdepw(k) |
---|
50 | # define fse3t(i,j,k) e3t(k) |
---|
51 | # define fse3u(i,j,k) e3t(k) |
---|
52 | # define fse3v(i,j,k) e3t(k) |
---|
53 | # define fse3f(i,j,k) e3t(k) |
---|
54 | # define fse3w(i,j,k) e3w(k) |
---|
55 | # define fse3uw(i,j,k) e3w(k) |
---|
56 | # define fse3vw(i,j,k) e3w(k) |
---|
57 | #endif |
---|