Changeset 1176 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar

Timestamp:

Jun 11, 2009, 10:54:10 AM (16 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

Location:

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar

Files:

• comconst.h (modified) (3 diffs)
• conf_gcm.F (modified) (6 diffs)
• inidissip.F (modified) (4 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/comconst.h

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
 !-----------------------------------------------------------------------
@@ -7 +7 @@
       COMMON/comconst/im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl,           &
      & dtvr,daysec,                                                     &
-     & pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg
+     & pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg              &
+     &                   ,dissip_factz,dissip_deltaz,dissip_zref
+
 
       INTEGER im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl
@@ -13 +15 @@
       REAL pi,dtphys,dtdiss,rad,r,cpp,kappa
       REAL cotot,unsim,g,omeg
+      REAL dissip_factz,dissip_deltaz,dissip_zref
 
 !-----------------------------------------------------------------------

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/conf_gcm.F

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
 c
@@ -44 +44 @@
 !#include "clesphys.h"
 #include "iniprint.h"
+#include "comconst.h"
 
 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
@@ -277 +278 @@
        tetatemp  = 7200.
        CALL getin('tetatemp',tetatemp )
+
+! Parametres controlant la variation sur la verticale des constantes de
+! dissipation.
+! Pour le moment actifs uniquement dans la version a 39 niveaux
+! avec ok_strato=y
+
+       dissip_factz=4.
+       dissip_deltaz=10.
+       dissip_zref=30.
+       CALL getin('dissip_factz',dissip_factz )
+       CALL getin('dissip_deltaz',dissip_deltaz )
+       CALL getin('dissip_zref',dissip_zref )
 
 !Config  Key  = coefdis
@@ -590 +603 @@
       write(lunout,*)' read_start = ', read_start
       write(lunout,*)' iflag_phys = ', iflag_phys
+      write(lunout,*)' iphysiq = ', iphysiq
       write(lunout,*)' clonn = ', clonn 
       write(lunout,*)' clatt = ', clatt
@@ -776 +790 @@
 !Config  Desc = activation de la version strato
 !Config  Def  = .FALSE.
-!Config  Help = active la version stratosph�rique de LMDZ de F. Lott
+!Config  Help = active la version stratosphérique de LMDZ de F. Lott
 
       ok_strato=.FALSE.
@@ -813 +827 @@
       write(lunout,*)' read_start = ', read_start
       write(lunout,*)' iflag_phys = ', iflag_phys
-      write(lunout,*)' clon = ', clon 
+      write(lunout,*)' iphysiq = ', iphysiq
+      write(lunout,*)' clon = ', clon
       write(lunout,*)' clat = ', clat
       write(lunout,*)' grossismx = ', grossismx

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/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