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

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/dyn3d

Files:

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

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

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
 c
@@ -30 +30 @@
 #include "ener.h"
 #include "description.h"
+#include "iniprint.h"
 
 c-------------------------------------------------------------------
@@ -74 +75 @@
       DATA dum/.true./
 
+      integer,save :: countcfl=0
+      real cflx(ip1jmp1,llm)
+      real cfly(ip1jm,llm)
+      real cflz(ip1jmp1,llm)
+      real, save :: cflxmax(llm),cflymax(llm),cflzmax(llm)
 
       IF(iadvtr.EQ.0) THEN
@@ -139 +145 @@
          ENDDO
 
+
+c-------------------------------------------------------------------
+! Calcul des criteres CFL en X, Y et Z
+c-------------------------------------------------------------------
+
+      if (countcfl == 0. ) then
+          cflxmax(:)=0.
+          cflymax(:)=0.
+          cflzmax(:)=0.
+      endif
+
+      countcfl=countcfl+iapp_tracvl
+      cflx(:,:)=0.
+      cfly(:,:)=0.
+      cflz(:,:)=0.
+      do l=1,llm
+         do ij=iip2,ip1jm-1
+            if (pbarug(ij,l)>=0.) then
+                cflx(ij,l)=pbarug(ij,l)*dtvr/masse(ij,l)
+            else
+                cflx(ij,l)=-pbarug(ij,l)*dtvr/masse(ij+1,l)
+            endif 
+         enddo
+      enddo
+      do l=1,llm
+         do ij=iip2,ip1jm-1,iip1
+            cflx(ij+iip1,l)=cflx(ij,l)
+         enddo
+      enddo
+
+      do l=1,llm
+         do ij=1,ip1jm
+            if (pbarvg(ij,l)>=0.) then
+                cfly(ij,l)=pbarvg(ij,l)*dtvr/masse(ij,l)
+            else
+                cfly(ij,l)=-pbarvg(ij,l)*dtvr/masse(ij+iip1,l)
+            endif 
+         enddo
+      enddo
+
+      do l=2,llm
+         do ij=1,ip1jm
+            if (wg(ij,l)>=0.) then
+                cflz(ij,l)=wg(ij,l)*dtvr/masse(ij,l)
+            else
+                cflz(ij,l)=-wg(ij,l)*dtvr/masse(ij,l-1)
+            endif 
+         enddo
+      enddo
+
+      do l=1,llm
+         cflxmax(l)=max(cflxmax(l),maxval(cflx(:,l)))
+         cflymax(l)=max(cflymax(l),maxval(cfly(:,l)))
+         cflzmax(l)=max(cflzmax(l),maxval(cflz(:,l)))
+      enddo
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Par defaut, on sort le diagnostic des CFL tous les jours.
+! Si on veut le sortir a chaque pas d'advection en cas de plantage 
+!     if (countcfl==iapp_tracvl) then
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+      if (countcfl==day_step) then
+         do l=1,llm
+         write(lunout,*) 'L, CFLmax '
+     s   ,l,maxval(cflx(:,l)),maxval(cfly(:,l)),maxval(cflz(:,l))
+         enddo
+         countcfl=0
+      endif
+    
 c-------------------------------------------------------------------
 c   Advection proprement dite (Modification Le Croller (07/2001)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/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/dyn3d/conf_gcm.F

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
 c
@@ -38 +38 @@
 #include "serre.h"
 #include "comdissnew.h"
+#include "comconst.h"
 
 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
@@ -267 +268 @@
        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

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