source: trunk/NEMO/TOP_SRC/passivetrc_substitute.h90 @ 719

Last change on this file since 719 was 719, checked in by ctlod, 13 years ago

get back to the nemo_v2_3 version for trunk

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.3 KB
Line 
1   !!----------------------------------------------------------------------
2   !!                    ***  passivetrc_substitute.h90   ***
3   !!----------------------------------------------------------------------
4   !! ** purpose : Statement function file: to be include in all routines
5   !!              concerning passive tracer model
6   !!----------------------------------------------------------------------
7   !!  TOP 1.0,  LOCEAN-IPSL (2005)
8   !! $Header$
9   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
10   !!----------------------------------------------------------------------
11
12    ! Lateral eddy diffusivity coefficient for passive tracer:
13    ! ========================================================
14#if defined key_traldf_c3d || defined key_off_degrad
15!   'key_traldf_c3d' :                 aht: 3D coefficient
16#       define fsahtrt(i,j,k)  trcrat * ahtt(i,j,k)
17#       define fsahtru(i,j,k)  trcrat * ahtu(i,j,k)
18#       define fsahtrv(i,j,k)  trcrat * ahtv(i,j,k)
19#       define fsahtrw(i,j,k)  trcrat * ahtw(i,j,k)
20#elif defined key_traldf_c2d
21!   'key_traldf_c2d' :                 aht: 2D coefficient
22#       define fsahtrt(i,j,k)  trcrat * ahtt(i,j)
23#       define fsahtru(i,j,k)  trcrat * ahtu(i,j)
24#       define fsahtrv(i,j,k)  trcrat * ahtv(i,j)
25#       define fsahtrw(i,j,k)  trcrat * ahtw(i,j)
26#elif defined key_traldf_c1d
27!   'key_traldf_c1d' :                 aht: 1D coefficient
28#       define fsahtrt(i,j,k)  trcrat * ahtt(k)
29#       define fsahtru(i,j,k)  trcrat * ahtu(k)
30#       define fsahtrv(i,j,k)  trcrat * ahtv(k)
31#       define fsahtrw(i,j,k)  trcrat * ahtw(k)
32#else
33!   Default option :             aht: Constant coefficient
34#       define fsahtrt(i,j,k)  ahtrc0
35#       define fsahtru(i,j,k)  ahtrc0
36#       define fsahtrv(i,j,k)  ahtrc0
37#       define fsahtrw(i,j,k)  ahtrc0
38#endif
39
40    ! Eddy induced velocity  coefficient for passive tracer:
41    ! ========================================================
42#if defined key_traldf_c3d || defined key_off_degrad
43!   'key_traldf_c3d' :                 eiv: 3D coefficient
44#       define fsaeitru(i,j,k)  trcrat * aeiu(i,j,k)
45#       define fsaeitrv(i,j,k)  trcrat * aeiv(i,j,k)
46#       define fsaeitrw(i,j,k)  trcrat * aeiw(i,j,k)
47#elif defined key_traldf_c2d
48!   'key_traldf_c2d' :                 eiv: 2D coefficient
49#       define fsaeitru(i,j,k)  trcrat * aeiu(i,j)
50#       define fsaeitrv(i,j,k)  trcrat * aeiv(i,j)
51#       define fsaeitrw(i,j,k)  trcrat * aeiw(i,j)
52#elif defined key_traldf_c1d
53!   'key_traldf_c1d' :                 eiv: 1D coefficient
54#       define fsaeitru(i,j,k)  trcrat * aeiu(k)
55#       define fsaeitrv(i,j,k)  trcrat * aeiv(k)
56#       define fsaeitrw(i,j,k)  trcrat * aeiw(k)
57#else
58!   Default option :             eiv: Constant coefficient
59#       define fsaeitru(i,j,k)  aeiv0
60#       define fsaeitrv(i,j,k)  aeiv0
61#       define fsaeitrw(i,j,k)  aeiv0
62#endif
63    ! mixing for passive tracer:
64    ! ========================================================
65#if defined key_trc_zdfddm
66!   'key_trc_zdfddm' :    avs
67#       define fstravs(i,j,k)   avs(i,j,k)
68#else
69!   Defautl option :   avs = avt
70#       define fstravs(i,j,k)   avt(i,j,k)
71#endif
72
73! ========================================================
74#include "domzgr_substitute.h90"
75#include "ldfeiv_substitute.h90"
76#include "ldftra_substitute.h90"
77#include "vectopt_loop_substitute.h90"
Note: See TracBrowser for help on using the repository browser.