domzgr_substitute.h90 4.8 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122
  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_vvl
  11. ! s* or z*-coordinate (3D + time dependency) + use of additional now arrays (..._n)
  12. # define fse3t_b(i,j,k) e3t_b(i,j,k)
  13. # define fse3u_b(i,j,k) e3u_b(i,j,k)
  14. # define fse3v_b(i,j,k) e3v_b(i,j,k)
  15. # define fse3uw_b(i,j,k) e3uw_b(i,j,k)
  16. # define fse3vw_b(i,j,k) e3vw_b(i,j,k)
  17. # define fsdept_n(i,j,k) gdept_n(i,j,k)
  18. # define fsdepw_n(i,j,k) gdepw_n(i,j,k)
  19. # define fsde3w_n(i,j,k) gdep3w_n(i,j,k)
  20. # define fse3t_n(i,j,k) e3t_n(i,j,k)
  21. # define fse3u_n(i,j,k) e3u_n(i,j,k)
  22. # define fse3v_n(i,j,k) e3v_n(i,j,k)
  23. # define fse3f_n(i,j,k) e3f_n(i,j,k)
  24. # define fse3w_n(i,j,k) e3w_n(i,j,k)
  25. # define fse3uw_n(i,j,k) e3uw_n(i,j,k)
  26. # define fse3vw_n(i,j,k) e3vw_n(i,j,k)
  27. # define fsdept_b(i,j,k) gdept_b(i,j,k)
  28. # define fsdepw_b(i,j,k) gdepw_b(i,j,k)
  29. # define fse3w_b(i,j,k) e3w_b(i,j,k)
  30. # define fse3t_a(i,j,k) e3t_a(i,j,k)
  31. # define fse3u_a(i,j,k) e3u_a(i,j,k)
  32. # define fse3v_a(i,j,k) e3v_a(i,j,k)
  33. # define fse3t_m(i,j) e3t_m(i,j)
  34. # define hu_a(i,j) ehu_a(i,j)
  35. # define hur_a(i,j) ehur_a(i,j)
  36. # define hv_a(i,j) ehv_a(i,j)
  37. # define hvr_a(i,j) ehvr_a(i,j)
  38. # define hu_b(i,j) ehu_b(i,j)
  39. # define hur_b(i,j) ehur_b(i,j)
  40. # define hv_b(i,j) ehv_b(i,j)
  41. # define hvr_b(i,j) ehvr_b(i,j)
  42. ! This part should be removed one day ...
  43. ! ... In that case all occurence of the above statement functions
  44. ! have to be replaced in the code by xxx_n
  45. # define fsdept(i,j,k) gdept_n(i,j,k)
  46. # define fsdepw(i,j,k) gdepw_n(i,j,k)
  47. # define fsde3w(i,j,k) gdep3w_n(i,j,k)
  48. # define fse3t(i,j,k) e3t_n(i,j,k)
  49. # define fse3u(i,j,k) e3u_n(i,j,k)
  50. # define fse3v(i,j,k) e3v_n(i,j,k)
  51. # define fse3f(i,j,k) e3f_n(i,j,k)
  52. # define fse3w(i,j,k) e3w_n(i,j,k)
  53. # define fse3uw(i,j,k) e3uw_n(i,j,k)
  54. # define fse3vw(i,j,k) e3vw_n(i,j,k)
  55. #else
  56. ! z- or s-coordinate (1D or 3D + no time dependency) use reference in all cases
  57. # define fse3t_b(i,j,k) e3t_0(i,j,k)
  58. # define fse3u_b(i,j,k) e3u_0(i,j,k)
  59. # define fse3v_b(i,j,k) e3v_0(i,j,k)
  60. # define fse3uw_b(i,j,k) e3uw_0(i,j,k)
  61. # define fse3vw_b(i,j,k) e3vw_0(i,j,k)
  62. # define fsdept_n(i,j,k) gdept_0(i,j,k)
  63. # define fsdepw_n(i,j,k) gdepw_0(i,j,k)
  64. # define fsde3w_n(i,j,k) gdep3w_0(i,j,k)
  65. # define fse3t_n(i,j,k) e3t_0(i,j,k)
  66. # define fse3u_n(i,j,k) e3u_0(i,j,k)
  67. # define fse3v_n(i,j,k) e3v_0(i,j,k)
  68. # define fse3f_n(i,j,k) e3f_0(i,j,k)
  69. # define fse3w_n(i,j,k) e3w_0(i,j,k)
  70. # define fse3uw_n(i,j,k) e3uw_0(i,j,k)
  71. # define fse3vw_n(i,j,k) e3vw_0(i,j,k)
  72. # define fsdept_b(i,j,k) gdept_0(i,j,k)
  73. # define fsdepw_b(i,j,k) gdepw_0(i,j,k)
  74. # define fse3w_b(i,j,k) e3w_0(i,j,k)
  75. # define fse3t_a(i,j,k) e3t_0(i,j,k)
  76. # define fse3u_a(i,j,k) e3u_0(i,j,k)
  77. # define fse3v_a(i,j,k) e3v_0(i,j,k)
  78. # define fse3t_m(i,j) e3t_0(i,j,1)
  79. # define hu_a(i,j) hu(i,j)
  80. # define hur_a(i,j) hur(i,j)
  81. # define hv_a(i,j) hv(i,j)
  82. # define hvr_a(i,j) hvr(i,j)
  83. # define hu_b(i,j) hu(i,j)
  84. # define hur_b(i,j) hur(i,j)
  85. # define hv_b(i,j) hv(i,j)
  86. # define hvr_b(i,j) hvr(i,j)
  87. ! This part should be removed one day ...
  88. ! ... In that case all occurence of the above statement functions
  89. ! have to be replaced in the code by xxx_n
  90. # define fsdept(i,j,k) gdept_0(i,j,k)
  91. # define fsdepw(i,j,k) gdepw_0(i,j,k)
  92. # define fsde3w(i,j,k) gdep3w_0(i,j,k)
  93. # define fse3t(i,j,k) e3t_0(i,j,k)
  94. # define fse3u(i,j,k) e3u_0(i,j,k)
  95. # define fse3v(i,j,k) e3v_0(i,j,k)
  96. # define fse3f(i,j,k) e3f_0(i,j,k)
  97. # define fse3w(i,j,k) e3w_0(i,j,k)
  98. # define fse3uw(i,j,k) e3uw_0(i,j,k)
  99. # define fse3vw(i,j,k) e3vw_0(i,j,k)
  100. #endif
  101. !!----------------------------------------------------------------------
  102. !! NEMO/OPA 3.3 , NEMO Consortium (2010)
  103. !! $Id: domzgr_substitute.h90 4488 2014-02-06 10:43:09Z rfurner $
  104. !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
  105. !!----------------------------------------------------------------------