source: branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_substitute.h90 @ 6498

Last change on this file since 6498 was 6498, checked in by timgraham, 5 years ago

Merge head of nemo_v3_6_STABLE into package branch

File size: 1.8 KB
Line 
1   !!----------------------------------------------------------------------
2   !!                    *** ldftra_substitute.h90  ***
3   !!----------------------------------------------------------------------
4   !! ** purpose :   substitute fsaht. the eddy diffusivity coeff.
5   !!      with a constant or 1D or 2D or 3D array, using CPP macro.
6   !!----------------------------------------------------------------------
7   !!----------------------------------------------------------------------
8   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
9   !! $Id$
10   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
11   !!----------------------------------------------------------------------
12#if defined key_traldf_c3d
13!   'key_traldf_c3d' :                 aht: 3D coefficient
14#       define   fsahtt(i,j,k)   rldf * ahtt(i,j,k)
15#       define   fsahtu(i,j,k)   rldf * ahtu(i,j,k) * r_fact_lap(i,j,k)
16#       define   fsahtv(i,j,k)   rldf * ahtv(i,j,k)
17#       define   fsahtw(i,j,k)   rldf * ahtw(i,j,k)
18#elif defined key_traldf_c2d
19!   'key_traldf_c2d' :                 aht: 2D coefficient
20#       define   fsahtt(i,j,k)   rldf * ahtt(i,j)
21#       define   fsahtu(i,j,k)   rldf * ahtu(i,j) * r_fact_lap(i,j,k)
22#       define   fsahtv(i,j,k)   rldf * ahtv(i,j)
23#       define   fsahtw(i,j,k)   rldf * ahtw(i,j)
24#elif defined key_traldf_c1d
25!   'key_traldf_c1d' :                aht: 1D coefficient
26#       define   fsahtt(i,j,k)   rldf * ahtt(k)
27#       define   fsahtu(i,j,k)   rldf * ahtu(k) * r_fact_lap(i,j,k)
28#       define   fsahtv(i,j,k)   rldf * ahtv(k)
29#       define   fsahtw(i,j,k)   rldf * ahtw(k)
30#else
31!   Default option :             aht: Constant coefficient
32#      define   fsahtt(i,j,k)   rldf * aht0
33#      define   fsahtu(i,j,k)   rldf * aht0 * r_fact_lap(i,j,k)
34#      define   fsahtv(i,j,k)   rldf * aht0
35#      define   fsahtw(i,j,k)   rldf * aht0
36#endif
Note: See TracBrowser for help on using the repository browser.