Changeset 1176 for LMDZ4/branches/LMDZ4-dev/libf/dyn3d/inidissip.F

Timestamp:

Jun 11, 2009, 10:54:10 AM (15 years ago)

Author:

Laurent Fairhead

Message:

Modif pour la compilation avec fcm sur Vargas SD
Corrections relatives au controle de histins et de l'utilisation des
options INST(X) pour l'appel au routines de IOIPSL FH
Mise en place de diagnostics sur les critères CFL pour l'advection
de traceurs. FH
Controle dans les .def de la dépendance verticale de l'efficacite
de la diffusion. Actif pour le moment uniquement avec ok_strato=y et llm=39. FH
LF

File:

• LMDZ4/branches/LMDZ4-dev/libf/dyn3d/inidissip.F (modified) (4 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/inidissip.F

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
       SUBROUTINE inidissip ( lstardis,nitergdiv,nitergrot,niterh  ,
@@ -18 +18 @@
 #include "comvert.h"
 #include "control.h"
+#include "logic.h"
 
       LOGICAL lstardis
@@ -29 +30 @@
       INTEGER l,ij,idum,ii
       REAL tetamin
+      REAL pseudoz
 
       REAL ran1
@@ -174 +176 @@
 c   --------------------------------------------------
 
-      DO l=1,llm
-         zvert(l)=1.
-      ENDDO
-
-      fact=2.
-c
-      DO l = 1, llm
-         zz      = 1. - preff/presnivs(l)
-         zvert(l)= fact -( fact-1.)/( 1.+zz*zz )
-      ENDDO
+      if (ok_strato .and. llm==39) then
+         do l=1,llm
+            pseudoz=8.*log(preff/presnivs(l))
+            zvert(l)=1+
+     s      (tanh((pseudoz-dissip_zref)/dissip_deltaz)+1.)/2.
+     s      *(dissip_factz-1.)
+         enddo 
+      else
+         DO l=1,llm
+            zvert(l)=1.
+         ENDDO
+         fact=2.
+         DO l = 1, llm
+            zz      = 1. - preff/presnivs(l)
+            zvert(l)= fact -( fact-1.)/( 1.+zz*zz )
+         ENDDO
+      endif