Changeset 2315 for LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/iniphysiq.F90

Timestamp:

Jun 25, 2015, 5:25:50 PM (9 years ago)

Author:

Ehouarn Millour

Message:

More on physics/dynamics separation: make a vertical_layers_mod module to contain information on the vertical discretization. This module should be used from within the physics (instead of including comvert.h from dynamics).
EM

File:

• LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/iniphysiq.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/iniphysiq.F90

@@ -3 +3 @@
 
 
-SUBROUTINE iniphysiq(iim,jjm,nlayer,punjours, pdayref,ptimestep,         &
+SUBROUTINE iniphysiq(ii,jj,nlayer,punjours, pdayref,ptimestep,           &
                      rlatu,rlonv,aire,cu,cv,                             &
                      prad,pg,pr,pcpp,iflag_phys)
@@ -13 +13 @@
                                 klon_omp_end, & ! end index of local omp subgrid
                                 klon_mpi_begin ! start indes of columns (on local mpi grid)
+  USE vertical_layers_mod, ONLY : init_vertical_layers
   USE comgeomphy, ONLY: initcomgeomphy, &
                         airephy, & ! physics grid area (m2)
@@ -28 +29 @@
   ! =======================================================================
 
+  include "dimensions.h"
+  include "comvert.h"
   include "iniprint.h"
 
@@ -36 +39 @@
   REAL, INTENT (IN) :: punjours ! length (in s) of a standard day
   INTEGER, INTENT (IN) :: nlayer ! number of atmospheric layers
-  INTEGER, INTENT (IN) :: iim ! number of atmospheric columns along longitudes
-  INTEGER, INTENT (IN) :: jjm ! number of atompsheric columns along latitudes
-  REAL, INTENT (IN) :: rlatu(jjm+1) ! latitudes of the physics grid
-  REAL, INTENT (IN) :: rlonv(iim+1) ! longitudes of the physics grid
-  REAL, INTENT (IN) :: aire(iim+1,jjm+1) ! area of the dynamics grid (m2)
-  REAL, INTENT (IN) :: cu((iim+1)*(jjm+1)) ! cu coeff. (u_covariant = cu * u)
-  REAL, INTENT (IN) :: cv((iim+1)*jjm) ! cv coeff. (v_covariant = cv * v)
+  INTEGER, INTENT (IN) :: ii ! number of atmospheric columns along longitudes
+  INTEGER, INTENT (IN) :: jj ! number of atompsheric columns along latitudes
+  REAL, INTENT (IN) :: rlatu(jj+1) ! latitudes of the physics grid
+  REAL, INTENT (IN) :: rlonv(ii+1) ! longitudes of the physics grid
+  REAL, INTENT (IN) :: aire(ii+1,jj+1) ! area of the dynamics grid (m2)
+  REAL, INTENT (IN) :: cu((ii+1)*(jj+1)) ! cu coeff. (u_covariant = cu * u)
+  REAL, INTENT (IN) :: cv((ii+1)*jj) ! cv coeff. (v_covariant = cv * v)
   INTEGER, INTENT (IN) :: pdayref ! reference day of for the simulation
   REAL, INTENT (IN) :: ptimestep !physics time step (s)
@@ -70 +73 @@
   END IF
 
-  !call init_phys_lmdz(iim,jjm+1,llm,1,(/(jjm-1)*iim+2/))
+  !call init_phys_lmdz(ii,jj+1,llm,1,(/(jj-1)*ii+2/))
   
   ! Generate global arrays on full physics grid
@@ -82 +85 @@
     cufi(1) = cu(1)
     cvfi(1) = cv(1)
-    DO j=2,jjm
-      DO i=1,iim
-        latfi((j-2)*iim+1+i)= rlatu(j)
-        lonfi((j-2)*iim+1+i)= rlonv(i)
-        cufi((j-2)*iim+1+i) = cu((j-1)*(iim+1)+i)
-        cvfi((j-2)*iim+1+i) = cv((j-1)*(iim+1)+i)
+    DO j=2,jj
+      DO i=1,ii
+        latfi((j-2)*ii+1+i)= rlatu(j)
+        lonfi((j-2)*ii+1+i)= rlonv(i)
+        cufi((j-2)*ii+1+i) = cu((j-1)*(ii+1)+i)
+        cvfi((j-2)*ii+1+i) = cv((j-1)*(ii+1)+i)
       ENDDO
     ENDDO
     ! South pole
-    latfi(klon_glo)= rlatu(jjm+1)
+    latfi(klon_glo)= rlatu(jj+1)
     lonfi(klon_glo)= 0.
-    cufi(klon_glo) = cu((iim+1)*jjm+1)
-    cvfi(klon_glo) = cv((iim+1)*jjm-iim)
+    cufi(klon_glo) = cu((ii+1)*jj+1)
+    cvfi(klon_glo) = cv((ii+1)*jj-ii)
 
     ! build airefi(), mesh area on physics grid
-    CALL gr_dyn_fi(1,iim+1,jjm+1,klon_glo,aire,airefi)
+    CALL gr_dyn_fi(1,ii+1,jj+1,klon_glo,aire,airefi)
     ! Poles are single points on physics grid
-    airefi(1)=sum(aire(1:iim,1))
-    airefi(klon_glo)=sum(aire(1:iim,jjm+1))
+    airefi(1)=sum(aire(1:ii,1))
+    airefi(klon_glo)=sum(aire(1:ii,jj+1))
 
     ! Sanity check: do total planet area match between physics and dynamics?
-    total_area_dyn=sum(aire(1:iim,1:jjm+1))
+    total_area_dyn=sum(aire(1:ii,1:jj+1))
     total_area_phy=sum(airefi(1:klon_glo))
     IF (total_area_dyn/=total_area_phy) THEN
@@ -128 +131 @@
   CALL inifis(punjours,prad,pg,pr,pcpp)
   
+  ! copy over preff , ap(), bp(), etc 
+  CALL init_vertical_layers(nlayer,preff,scaleheight, &
+                            ap,bp,presnivs,pseudoalt)
+
   ! Now generate local lon/lat/cu/cv/area arrays
   CALL initcomgeomphy