Changeset 7311 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_substitute.h90

Timestamp:

2016-11-23T07:30:18+01:00 (8 years ago)

Author:

cbricaud

Message:

changes in CRS branch for traldf to allow to run CRS only for output of physics ( no key_top and no use of TOP_SRC

File:

• branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_substitute.h90 (modified) (1 diff)

branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_substitute.h90

@@ -35 +35 @@
 #      define   fsahtw(i,j,k)   rldf * aht0
 #endif
+
+   !!----------------------------------------------------------------------
+
+#if defined key_crs
+#if defined key_traldf_c3d
+!   'key_traldf_c3d' :                 aht: 3D coefficient
+#       define   fsahtt_crs(i,j,k)   rldf * ahtt_crs(i,j,k)
+#       define   fsahtu_crs(i,j,k)   rldf * ahtu_crs(i,j,k) * r_fact_lap_crs(i,j,k)
+#       define   fsahtv_crs(i,j,k)   rldf * ahtv_crs(i,j,k)
+#       define   fsahtw_crs(i,j,k)   rldf * ahtw_crs(i,j,k)
+#elif defined key_traldf_c2d
+!   'key_traldf_c2d' :                 aht: 2D coefficient
+#       define   fsahtt_crs(i,j,k)   rldf * ahtt_crs(i,j)
+#       define   fsahtu_crs(i,j,k)   rldf * ahtu_crs(i,j) * r_fact_lap_crs(i,j,k)
+#       define   fsahtv_crs(i,j,k)   rldf * ahtv_crs(i,j)
+#       define   fsahtw_crs(i,j,k)   rldf * ahtw_crs(i,j)
+#elif defined key_traldf_c1d
+!   'key_traldf_c1d' :                aht: 1D coefficient
+#       define   fsahtt_crs(i,j,k)   rldf * ahtt_crs(k)
+#       define   fsahtu_crs(i,j,k)   rldf * ahtu_crs(k) * r_fact_lap_crs(i,j,k)
+#       define   fsahtv_crs(i,j,k)   rldf * ahtv_crs(k)
+#       define   fsahtw_crs(i,j,k)   rldf * ahtw_crs(k)
+#else
+!   Default option :             aht: Constant coefficient
+#      define   fsahtt_crs(i,j,k)   rldf * aht0
+#      define   fsahtu_crs(i,j,k)   rldf * aht0 * r_fact_lap_crs(i,j,k)
+#      define   fsahtv_crs(i,j,k)   rldf * aht0
+#      define   fsahtw(_crsi,j,k)   rldf * aht0
+#endif
+#endif