Changeset 1520 for LMDZ5/trunk/libf/dyn3d/leapfrog.F

Timestamp:

May 23, 2011, 1:37:09 PM (13 years ago)

Author:

Ehouarn Millour

Message:

Implementation of a different vertical discretization (from/for planets, but
can in principle also be used for Earth).
Choice of vertical discretization is set by flag 'disvert_type';
'disvert_type=1' is Earth standard (default; ie set to 1 if
planet_type=="earth") case.
With 'disvert_type=2', approximate altitude of layers and reference atmospheric
scale height must be given using an input file ("z2sig.def", first line
should give scale height, in km, following lines must specify the altitude,
in km above surface, of mid-layers, one per line; see disvert_noterre.F).

Checked that these changes do not impact on 'bench' results, on Vargas.

EM.

File:

• LMDZ5/trunk/libf/dyn3d/leapfrog.F (modified) (5 diffs)

LMDZ5/trunk/libf/dyn3d/leapfrog.F

@@ -172 +172 @@
 
       character*80 dynhist_file, dynhistave_file
-      character(len=20) :: modname
+      character(len=*),parameter :: modname="leapfrog"
       character*80 abort_message
 
@@ -193 +193 @@
       itaufin   = nday*day_step
       itaufinp1 = itaufin +1
-      modname="leapfrog"
       
 
@@ -211 +210 @@
       dq(:,:,:)=0.
       CALL pression ( ip1jmp1, ap, bp, ps, p       )
-      CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
+      if (disvert_type==1) then
+        CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
+      else ! we assume that we are in the disvert_type==2 case
+        CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf )
+      endif
 
 c-----------------------------------------------------------------------
@@ -359 +362 @@
 
          CALL pression (  ip1jmp1, ap, bp, ps,  p      )
-         CALL exner_hyb(  ip1jmp1, ps, p,alpha,beta,pks, pk, pkf )
+         if (disvert_type==1) then
+           CALL exner_hyb(  ip1jmp1, ps, p,alpha,beta,pks, pk, pkf )
+         else ! we assume that we are in the disvert_type==2 case
+           CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf )
+         endif
 
 !           rdaym_ini  = itau * dtvr / daysec
@@ -463 +470 @@
 
         CALL pression ( ip1jmp1, ap, bp, ps, p                  )
-        CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
+        if (disvert_type==1) then
+          CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
+        else ! we assume that we are in the disvert_type==2 case
+          CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf )
+        endif