Changeset 2604 for LMDZ5/trunk/libf/dynphy_lonlat/calfis.F

Timestamp:

Jul 26, 2016, 5:37:18 PM (8 years ago)

Author:

Ehouarn Millour

Message:

Add Exner function to the call_physiq arguments (not used by the Earth physics) to harmonize physics/dynamics interface.
EM

File:

• LMDZ5/trunk/libf/dynphy_lonlat/calfis.F (modified) (4 diffs)

LMDZ5/trunk/libf/dynphy_lonlat/calfis.F

@@ -141 +141 @@
 c
       REAL zrot(iip1,jjm,llm) ! AdlC May 2014
-      REAL zufi(ngridmx,llm), zvfi(ngridmx,llm), zrfi(ngridmx,llm)
+      REAL zufi(ngridmx,llm), zvfi(ngridmx,llm)
+      REAL zrfi(ngridmx,llm) ! relative wind vorticity
       REAL ztfi(ngridmx,llm),zqfi(ngridmx,llm,nqtot)
+      REAL zpk(ngridmx,llm)
 c
       REAL pcvgu(ngridmx,llm), pcvgv(ngridmx,llm)
@@ -211 +213 @@
 
 
-c   42. pression intercouches :
+c   42. pression intercouches et fonction d'Exner:
 c
 c   -----------------------------------------------------------------
@@ -222 +224 @@
        unskap   = 1./ kappa
 c
-      DO l = 1, llmp1
+      DO l = 1, llm
+        zpk(   1,l ) = ppk(1,1,l)
         zplev( 1,l ) = pp(1,1,l)
         ig0 = 2
           DO j = 2, jjm
              DO i =1, iim
+              zpk(   ig0,l ) = ppk(i,j,l)
               zplev( ig0,l ) = pp(i,j,l)
               ig0 = ig0 +1
              ENDDO
           ENDDO
+        zpk(   ngridmx,l ) = ppk(1,jjp1,l)
         zplev( ngridmx,l ) = pp(1,jjp1,l)
       ENDDO
+        zplev( 1,llmp1 ) = pp(1,1,llmp1)
+        ig0 = 2
+          DO j = 2, jjm
+             DO i =1, iim
+              zplev( ig0,llmp1 ) = pp(i,j,llmp1)
+              ig0 = ig0 +1
+             ENDDO
+          ENDDO
+        zplev( ngridmx,llmp1 ) = pp(1,jjp1,llmp1)
 c
 c
@@ -471 +485 @@
      &                   jD_cur,jH_cur_split,zdt_split,
      &                   zplev,zplay,
-     &                   zphi,zphis,
+     &                   zpk,zphi,zphis,
      &                   presnivs,
      &                   zufi,zvfi,zrfi,ztfi,zqfi,