Changeset 960 for LMDZ4/trunk/libf/dyn3d/calfis.F

Timestamp:

May 27, 2008, 6:46:50 PM (17 years ago)

Author:

lsce

Message:

• Ajoute du parametre config_inca dans conf_gcm.F config_inca='none'(sans INCA, par defaut) config_inca='chem'(avec INCA config chemie) config_inca='aero'(avec INCA config aerosol)
• Menage parmis les cles CPP INCA
• Enleve le calcul d'omega dans calfis.F et active le calcul correspondant dans advtrac.F(avant uniquement pour INCA).

JG

File:

• LMDZ4/trunk/libf/dyn3d/calfis.F (modified) (7 diffs)

LMDZ4/trunk/libf/dyn3d/calfis.F

@@ -22 +22 @@
      $                  pdteta,
      $                  pdq,
-     $                  pw,
-#ifdef INCA
      $                  flxw,
-#endif
      $                  clesphy0,
      $                  pdufi,
@@ -121 +118 @@
       REAL pdq(iip1,jjp1,llm,nqmx)
 c
-      REAL pw(iip1,jjp1,llm)
-
       REAL pps(iip1,jjp1)
       REAL pp(iip1,jjp1,llmp1)
@@ -151 +146 @@
       REAL pcvgu(ngridmx,llm), pcvgv(ngridmx,llm)
       REAL pcvgt(ngridmx,llm), pcvgq(ngridmx,llm,2)
-c
-      REAL pvervel(ngridmx,llm)
 c
       REAL zdufi(ngridmx,llm),zdvfi(ngridmx,llm)
@@ -169 +162 @@
       REAL PVteta(ngridmx,ntetaSTD)
 c
-#ifdef INCA
-      REAL flxw(iip1,jjp1,llm)
-      REAL flxwfi(ngridmx,llm)
-#endif
+      REAL flxw(iip1,jjp1,llm)  ! Flux de masse verticale sur la grille dynamique
+      REAL flxwfi(ngridmx,llm)  ! Flux de masse verticale sur la grille physiq
 c
       
@@ -328 +319 @@
 
 c   ....  Calcul de la vitesse  verticale  ( en Pa*m*s  ou Kg/s )  ....
-c
-      DO l=1,llm
-        pvervel(1,l)=pw(1,1,l) * g /apoln
-        ig0=2
-       DO j=2,jjm
-           DO i = 1, iim
-              pvervel(ig0,l) = pw(i,j,l) * g * unsaire(i,j)
-              ig0 = ig0 + 1
-           ENDDO
-       ENDDO
-        pvervel(ig0,l)=pw(1,jjp1,l) * g /apols
-      ENDDO
+c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux 
+c    de masse est calclue dans advtrac.F  
+c      DO l=1,llm
+c        pvervel(1,l)=pw(1,1,l) * g /apoln
+c        ig0=2
+c       DO j=2,jjm
+c           DO i = 1, iim
+c              pvervel(ig0,l) = pw(i,j,l) * g * unsaire(i,j)
+c              ig0 = ig0 + 1
+c           ENDDO
+c       ENDDO
+c        pvervel(ig0,l)=pw(1,jjp1,l) * g /apols
+c      ENDDO
 
 c
@@ -442 +434 @@
      $           ntetaSTD,rtetaSTD,PVteta)
 c
-#ifdef INCA
+c On change de grille, dynamique vers physiq, pour le flux de masse verticale
       CALL gr_dyn_fi(llm,iip1,jjp1,ngridmx,flxw,flxwfi)
-#endif
-
 
 c-----------------------------------------------------------------------
@@ -470 +460 @@
      .             ztfi,
      .             zqfi,
-     .             pvervel,
-#ifdef INCA
      .             flxwfi,
-#endif
      .             zdufi,
      .             zdvfi,