Changeset 2346 for LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd

Timestamp:

Aug 21, 2015, 5:13:46 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Physics/dynamics separation:

• remove all references to dimensions.h from physics. nbp_lon (==iim) , nbp_lat (==jjm+1) and nbp_lev (==llm) from mod_grid_phy_lmdz should be used instead.
• added module regular_lonlat_mod in phy_common to store information about the global (lon-lat) grid cell boundaries and centers.

EM

File:

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

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

@@ -3 +3 @@
 
 
-SUBROUTINE iniphysiq(ii,jj,nlayer,punjours, pdayref,ptimestep,           &
-                     rlatu,rlonv,aire,cu,cv,                             &
+SUBROUTINE iniphysiq(ii,jj,nlayer,punjours, pdayref,ptimestep, &
+                     rlatu,rlatv,rlonu,rlonv,aire,cu,cv,       &
                      prad,pg,pr,pcpp,iflag_phys)
   USE dimphy, ONLY: klev ! number of atmospheric levels
@@ -35 +35 @@
   USE phystokenc_mod, ONLY: init_phystokenc
   USE phyaqua_mod, ONLY: iniaqua
+  USE regular_lonlat_mod, ONLY : init_regular_lonlat, &
+                                 east, west, north, south, &
+                                 north_east, north_west, &
+                                 south_west, south_east
   IMPLICIT NONE
 
@@ -44 +48 @@
   include "dimensions.h"
   include "comvert.h"
+  include "comconst.h"
   include "iniprint.h"
   include "temps.h"
@@ -57 +62 @@
   INTEGER, INTENT (IN) :: jj ! number of atompsheric columns along latitudes
   REAL, INTENT (IN) :: rlatu(jj+1) ! latitudes of the physics grid
+  REAL, INTENT (IN) :: rlatv(jj) ! latitude boundaries of the physics grid
   REAL, INTENT (IN) :: rlonv(ii+1) ! longitudes of the physics grid
+  REAL, INTENT (IN) :: rlonu(ii+1) ! longitude boundaries 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)
@@ -71 +78 @@
   REAL :: total_area_phy, total_area_dyn
 
+  ! boundaries, on global grid
+  REAL,ALLOCATABLE :: boundslon_reg(:,:)
+  REAL,ALLOCATABLE :: boundslat_reg(:,:)
 
   ! global array, on full physics grid:
@@ -87 +97 @@
   !call init_phys_lmdz(ii,jj+1,llm,1,(/(jj-1)*ii+2/))
   
+  ! init regular global longitude-latitude grid points and boundaries
+  ALLOCATE(boundslon_reg(ii,2))
+  ALLOCATE(boundslat_reg(jj+1,2))
+  
+  DO i=1,ii
+   boundslon_reg(i,east)=rlonu(i) 
+   boundslon_reg(i,west)=rlonu(i+1) 
+  ENDDO
+
+  boundslat_reg(1,north)= PI/2 
+  boundslat_reg(1,south)= rlatv(1)
+  DO j=2,jj
+   boundslat_reg(j,north)=rlatv(j-1) 
+   boundslat_reg(j,south)=rlatv(j) 
+  ENDDO
+  boundslat_reg(jj+1,north)= rlatv(jj) 
+  boundslat_reg(jj+1,south)= -PI/2
+
+  ! Write values in module regular_lonlat_mod
+  CALL init_regular_lonlat(ii,jj+1, rlonv(1:ii), rlatu, &
+                           boundslon_reg, boundslat_reg)
+
   ! Generate global arrays on full physics grid
   ALLOCATE(latfi(klon_glo),lonfi(klon_glo),cufi(klon_glo),cvfi(klon_glo))