Changes between Version 6 and Version 7 of ticket/0842_TRA_TRP

Timestamp:

2011-08-02T12:37:45+02:00 (13 years ago)

Author:

cetlod

Comment:

--

ticket/0842_TRA_TRP

@@ -42 +42 @@
 
 In the code, we ensure that the 2 coef have the same sign ( < 0 for bilaplacian operator ; > 0 for laplacian operator )
+
+{{{
+#if defined key_traldf_c3d
+!   'key_traldf_c3d' :                 aht: 3D coefficient
+#       define   fsahtt(i,j,k)   rldf * ahtt(i,j,k)
+#       define   fsahtu(i,j,k)   rldf * ahtu(i,j,k)
+#       define   fsahtv(i,j,k)   rldf * ahtv(i,j,k)
+#       define   fsahtw(i,j,k)   rldf * ahtw(i,j,k)
+#elif defined key_traldf_c2d
+!   'key_traldf_c2d' :                 aht: 2D coefficient
+#       define   fsahtt(i,j,k)   rldf * ahtt(i,j)
+#       define   fsahtu(i,j,k)   rldf * ahtu(i,j)
+#       define   fsahtv(i,j,k)   rldf * ahtv(i,j)
+#       define   fsahtw(i,j,k)   rldf * ahtw(i,j)
+#elif defined key_traldf_c1d
+!   'key_traldf_c1d' :                aht: 1D coefficient
+#       define   fsahtt(i,j,k)   rldf * ahtt(k)
+#       define   fsahtu(i,j,k)   rldf * ahtu(k)
+#       define   fsahtv(i,j,k)   rldf * ahtv(k)
+#       define   fsahtw(i,j,k)   rldf * ahtw(k)
+#else
+!   Default option :             aht: Constant coefficient
+#      define   fsahtt(i,j,k)   rldf * aht0
+#      define   fsahtu(i,j,k)   rldf * aht0
+#      define   fsahtv(i,j,k)   rldf * aht0
+#      define   fsahtw(i,j,k)   rldf * aht0
+#endif
+}}}
 
 '''(3)''' replace 3D T & S arrays by 4D array TS throughout the code [[BR]]