Changeset 1000 for LMDZ4/trunk/libf/dyn3dpar/disvert.F

Timestamp:

Oct 6, 2008, 10:43:22 AM (16 years ago)

Author:

Laurent Fairhead

Message:

• Modifs sur le parallelisme: masquage dans la physique
• Inclusion strato
• mise en coherence etat0
• le mode offline fonctionne maintenant en parallele,
• les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
• la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
• L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

File:

• LMDZ4/trunk/libf/dyn3dpar/disvert.F (modified) (2 diffs)

LMDZ4/trunk/libf/dyn3dpar/disvert.F

@@ -11 +11 @@
 #include "paramet.h"
 #include "iniprint.h"
+#include "logic.h"
 c
 c=======================================================================
@@ -99 +100 @@
       DO l = 1, llm
          x = 2.*asin(1.) * (FLOAT(l)-0.5) / float(llm+1)
-         dsig(l) = 1.0 + 7.0 * SIN(x)**2
+
+         IF (ok_strato) THEN
+           dsig(l) =(1.0 + 7.0 * SIN(x)**2)
+     &            *(0.5*(1.-tanh(1.*(x-asin(1.))/asin(1.))))**2        
+         ELSE
+           dsig(l) = 1.0 + 7.0 * SIN(x)**2
+         ENDIF
+
          snorm = snorm + dsig(l)
       ENDDO