Changeset 72 for trunk/src/tb.h

Timestamp:

09/30/11 17:32:34 (13 years ago)

Author:

jbrlod

Message:

test du qs

File:

• trunk/src/tb.h (modified) (1 diff)

trunk/src/tb.h

@@ -41 +41 @@
       // pour generer le tourbillon
       /*  YS1=YS1+exp(-((Yi+1.-(3.*NX)/4.)*(Yi+1.-(3.*NX)/4.)+ 
-                    (Yj+1.-(3.*NY)/4.)*(Yj+1.-(3.*NY)/4.))/4.)
-                    *(1.-exp(-double(Yk+1.-NZ)/NZ))*(1.-exp(-double(Yt-TU)/24.));*/
-    }
+          (Yj+1.-(3.*NY)/4.)*(Yj+1.-(3.*NY)/4.))/4.)
+          *(1.-exp(-double(Yk+1.-NZ)/NZ))*(1.-exp(-double(Yt-TU)/24.));*/
+#ifdef EDDY
+      if (Yt<27)
+      //pour generer un tourbillon à Yi=19 et Yj=13
+      YS1=YS1+4*exp(-0.75*((Yi - 17.*(NX/32.))*(Yi - 17.*(NX/32.))+ 
+                    (Yj - 10.*(NY/22.))*(Yj - 10.*(NY/22.))))
+        *(1.-exp(-double(Yk+1.-NZ)/NZ))*(1.-exp(-double(Yt-TU)/24.));
+#endif
+        }
   }
   //
 }
-//===========================================================================
-//                         methode  backward
+    //===========================================================================
+    //                         methode  backward
 
 backward (YREAL x1,YREAL x2,YREAL x3,YREAL x4)