ldfdyn_substitute.h90 1.8 KB

12345678910111213141516171819202122232425262728293031323334353637383940
  1. !!----------------------------------------------------------------------
  2. !! *** ldfdyn_substitute.h90 ***
  3. !!----------------------------------------------------------------------
  4. !! ** purpose : substitute fsahm., the lateral eddy viscosity 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: ldfdyn_substitute.h90 2528 2010-12-27 17:33:53Z rblod $
  10. !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
  11. !!----------------------------------------------------------------------
  12. !!
  13. !! fsahmt, fsahmf - used for laplaian operator only
  14. !! fsahmu, fsahmv - used for bilaplacian operator only
  15. !!
  16. #if defined key_dynldf_c3d
  17. ! ' key_dynldf_c3d' : 3D coefficient
  18. # define fsahmt(i,j,k) ahm1(i,j,k)
  19. # define fsahmf(i,j,k) ahm2(i,j,k)
  20. # define fsahmu(i,j,k) ahm3(i,j,k)
  21. # define fsahmv(i,j,k) ahm4(i,j,k)
  22. #elif defined key_dynldf_c2d
  23. ! ' key_dynldf_c2d' : 2D coefficient
  24. # define fsahmt(i,j,k) ahm1(i,j)
  25. # define fsahmf(i,j,k) ahm2(i,j)
  26. # define fsahmu(i,j,k) ahm3(i,j)
  27. # define fsahmv(i,j,k) ahm4(i,j)
  28. #elif defined key_dynldf_c1d
  29. ! 'key_dynldf_c1d' : 1D coefficient
  30. # define fsahmt(i,j,k) ahm1(k)
  31. # define fsahmf(i,j,k) ahm1(k)
  32. # define fsahmu(i,j,k) ahm2(k)
  33. # define fsahmv(i,j,k) ahm2(k)
  34. #else
  35. ! default option : Constant coefficient
  36. # define fsahmt(i,j,k) ahm0
  37. # define fsahmf(i,j,k) ahm0
  38. # define fsahmu(i,j,k) ahm0_blp
  39. # define fsahmv(i,j,k) ahm0_blp
  40. #endif