Changeset 1654 for LMDZ5/trunk/libf/dyn3d/calfis.F

Timestamp:

Sep 24, 2012, 5:07:18 PM (12 years ago)

Author:

aslmd

Message:

---

• Changements dans calfis pour appeler la physique generique planetaire NB: dans le futur il faudrait un appel commun mais actuellement differences dans les arguments...
• Changements dans dynetat0 pour s'accommoder des tableaux "controle" terrestre et planeto
• Passage de la variable pseudoalt dans comvert.h (compatibilite planeto)
• arch : ajout des fichiers pour faire fonctionner sur ciclad avec ifort

verifie : compatibilite bit-par-bit des resultats et initialisations, ainsi que 1+1=2 avec unpun.sh

---

• Changes in calfis so that generic physics for planets could be interfaced NB: ideally a common call should be implemented but for now there are differences in arguments...
• Changes in dynetat0 to cope with difference in "controle" between earth and planeto
• Transfer of variable pseudoalt into comvert.h (this is only used by planeto and done for compatibility)
• arch: added files to allow use of ifort compiler on ciclad cluster

checked : bit-by-bit compatibility of results and starts, as well as 1+1=2 with unpun.sh

---

Author : AS equipe planeto

File:

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

LMDZ5/trunk/libf/dyn3d/calfis.F

@@ -179 +179 @@
 !      REAL rdayvrai
       REAL, intent(in):: jD_cur, jH_cur
+
+      LOGICAL tracerdyn
+
 c
 c-----------------------------------------------------------------------
@@ -459 +462 @@
       zdqfic(:,:,:)=0.
 
-      if (planet_type=="earth") then
 #ifdef CPP_PHYS
 
@@ -467 +469 @@
          debut_split=debut.and.isplit==1
          lafin_split=lafin.and.isplit==nsplit_phys
+
+      if (planet_type=="earth") then
 
          CALL physiq (ngridmx,
@@ -495 +499 @@
      .             PVteta)
 
+      else if ( planet_type=="generic" ) then
+
+         CALL physiq (ngridmx,     !! ngrid
+     .             llm,            !! nlayer
+     .             nqtot,          !! nq
+     .             tname,          !! tracer names from dynamical core (given in infotrac)
+     .             debut_split,    !! firstcall 
+     .             lafin_split,    !! lastcall
+     .             float(day_ini), !! pday <-- day_ini (dans temps.h)
+     .             jH_cur_split,   !! ptime "fraction of day"
+     .             zdt_split,      !! ptimestep
+     .             zplev,          !! pplev
+     .             zplay,          !! pplay
+     .             zphi,           !! pphi
+     .             zufi,           !! pu
+     .             zvfi,           !! pv
+     .             ztfi,           !! pt
+     .             zqfi,           !! pq
+     .             flxwfi,         !! pw !! or 0. anyway this is for diagnostic. not used in physiq.
+     .             zdufi,          !! pdu
+     .             zdvfi,          !! pdv
+     .             zdtfi,          !! pdt
+     .             zdqfi,          !! pdq
+     .             zdpsrf,         !! pdpsrf
+     .             tracerdyn)      !! tracerdyn <-- utilite ???
+
+      endif ! of if (planet_type=="earth")
+
          zufi(:,:)=zufi(:,:)+zdufi(:,:)*zdt_split
          zvfi(:,:)=zvfi(:,:)+zdvfi(:,:)*zdt_split
@@ -509 +541 @@
 #endif
 ! of #ifdef CPP_PHYS
-      endif ! of if (planet_type=="earth")
 
       zdufi(:,:)=zdufic(:,:)/nsplit_phys