Changeset 127 for trunk/libf/dyn3dpar/leapfrog_p.F

Timestamp:

May 24, 2011, 1:26:29 PM (14 years ago)

Author:

emillour

Message:

Ehouarn: suite de l'implementation de la discretisation verticale,
avec quelques mises a jour pour concorder avec la version terrestre.
-> Finalement, on met un flag "disvert_type" pour fixer la discretisation

disvert_type==1 (defaut si planet_type=="earth") pour cas terrestre
disvert_type==2 (defaut si planet_type!="earth") pour cas planeto (z2sig.def)

-> au passage, pour rester en phase avec modele terrestre on renomme

disvert_terre en disvert (le disvert "alternatif" demeure 'disvert_noterre')

File:

• trunk/libf/dyn3dpar/leapfrog_p.F (modified) (3 diffs)

trunk/libf/dyn3dpar/leapfrog_p.F

@@ -271 +271 @@
 
       CALL pression ( ip1jmp1, ap, bp, ps, p       )
-      if (planet_type.eq."earth") then
+      if (disvert_type==1) then
         CALL exner_hyb( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
-      else
+      else ! we assume that we are in the disvert_type==2 case
         CALL exner_milieu( ip1jmp1, ps, p, beta, pks, pk, pkf )
       endif
@@ -751 +751 @@
 
 c$OMP BARRIER
-         if (planet_type.eq."earth") then
-          CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
-         else
-          CALL exner_milieu_p( ip1jmp1, ps, p, beta, pks, pk, pkf )
+         if (disvert_type==1) then
+           CALL exner_hyb_p(  ip1jmp1, ps, p,alpha,beta,pks, pk, pkf )
+         else ! we assume that we are in the disvert_type==2 case
+           CALL exner_milieu_p( ip1jmp1, ps, p, beta, pks, pk, pkf )
          endif
 c$OMP BARRIER
@@ -1105 +1105 @@
         CALL pression_p ( ip1jmp1, ap, bp, ps, p                  )
 c$OMP BARRIER
-        if (planet_type.eq."earth") then
+        if (disvert_type==1) then
           CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
-        else
+        else ! we assume that we are in the disvert_type==2 case
           CALL exner_milieu_p( ip1jmp1, ps, p, beta, pks, pk, pkf )
         endif