Changeset 1529 for trunk/LMDZ.GENERIC/libf/dynlonlat_phylonlat/phystd/iniphysiq_mod.F90

Timestamp:

Apr 5, 2016, 10:51:51 AM (9 years ago)

Author:

emillour

Message:

Generic GCM: Towards a cleaner separation between physics and dynamics

• Got rid of references to "dimensions.h" from physics packages: use nbp_lon (=iim), nbp_lat (=jjp1) and nbp_lev from module mod_grid_phy_lmdz (in phy_common) instead.
• Removed module "comhdiff_mod.F90", as it is only used by module surf_heat_transp_mod.F90, moved module variables there.
• Addedin "surf_heat_transp_mod" local versions of some arrays and routines (from dyn3d) required to compute gradient, divergence, etc. on the global dynamics grid. As before, the slab ocean only works in serial.

EM

File:

• trunk/LMDZ.GENERIC/libf/dynlonlat_phylonlat/phystd/iniphysiq_mod.F90 (modified) (7 diffs)

trunk/LMDZ.GENERIC/libf/dynlonlat_phylonlat/phystd/iniphysiq_mod.F90

@@ -4 +4 @@
 
 subroutine iniphysiq(ii,jj,nlayer,punjours, pdayref,ptimestep,           &
-                     rlatu,rlatv,rlonu,rlonv,aire,cu,cv,                 &
+               rlatudyn,rlatvdyn,rlonudyn,rlonvdyn,airedyn,cudyn,cvdyn,  &
                      prad,pg,pr,pcpp,iflag_phys)
 
@@ -21 +21 @@
                        rlond, & ! longitudes
                        rlatd ! latitudes
+use surf_heat_transp_mod, only: ini_surf_heat_transp
 use infotrac, only : nqtot ! number of advected tracers
 use planete_mod, only: ini_planete_mod
@@ -29 +30 @@
                               north_east, north_west, &
                               south_west, south_east
+use ioipsl_getin_p_mod, only: getin_p
 
 implicit none
 include "dimensions.h"
+include "paramet.h"
+include "comgeom.h"
 include "iniprint.h"
 
@@ -43 +47 @@
 integer,intent(in) :: ii ! number of atmospheric coulumns 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) :: 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)
-real,intent(in) :: cv((ii+1)*jj) ! cv coeff. (v_covariant = cv * v)
+real,intent(in) :: rlatudyn(jj+1) ! latitudes of the physics grid
+real,intent(in) :: rlatvdyn(jj) ! latitude boundaries of the physics grid
+real,intent(in) :: rlonvdyn(ii+1) ! longitudes of the physics grid
+real,intent(in) :: rlonudyn(ii+1) ! longitude boundaries of the physics grid
+real,intent(in) :: airedyn(ii+1,jj+1) ! area of the dynamics grid (m2)
+real,intent(in) :: cudyn((ii+1)*(jj+1)) ! cu coeff. (u_covariant = cu * u)
+real,intent(in) :: cvdyn((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)
@@ -60 +64 @@
 real :: total_area_phy, total_area_dyn
 real :: pi
+logical :: ok_slab_ocean
 
 ! boundaries, on global grid
@@ -99 +104 @@
   
 DO i=1,ii
-   boundslon_reg(i,east)=rlonu(i) 
-   boundslon_reg(i,west)=rlonu(i+1) 
+   boundslon_reg(i,east)=rlonudyn(i) 
+   boundslon_reg(i,west)=rlonudyn(i+1) 
 ENDDO
 
 boundslat_reg(1,north)= PI/2 
-boundslat_reg(1,south)= rlatv(1)
+boundslat_reg(1,south)= rlatvdyn(1)
 DO j=2,jj
-   boundslat_reg(j,north)=rlatv(j-1) 
-   boundslat_reg(j,south)=rlatv(j) 
+   boundslat_reg(j,north)=rlatvdyn(j-1) 
+   boundslat_reg(j,south)=rlatvdyn(j) 
 ENDDO
-boundslat_reg(jj+1,north)= rlatv(jj) 
+boundslat_reg(jj+1,north)= rlatvdyn(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, &
+CALL init_regular_lonlat(ii,jj+1, rlonvdyn(1:ii), rlatudyn, &
                          boundslon_reg, boundslat_reg)
 
 ! Generate global arrays on full physics grid
 allocate(latfi(klon_glo),lonfi(klon_glo),cufi(klon_glo),cvfi(klon_glo))
-latfi(1)=rlatu(1)
+latfi(1)=rlatudyn(1)
 lonfi(1)=0.
-cufi(1) = cu(1)
-cvfi(1) = cv(1)
+cufi(1) = cudyn(1)
+cvfi(1) = cvdyn(1)
 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)
+    latfi((j-2)*ii+1+i)= rlatudyn(j)
+    lonfi((j-2)*ii+1+i)= rlonvdyn(i)
+    cufi((j-2)*ii+1+i) = cudyn((j-1)*(ii+1)+i)
+    cvfi((j-2)*ii+1+i) = cvdyn((j-1)*(ii+1)+i)
   ENDDO
 ENDDO
-latfi(klon_glo)= rlatu(jj+1)
+latfi(klon_glo)= rlatudyn(jj+1)
 lonfi(klon_glo)= 0.
-cufi(klon_glo) = cu((ii+1)*jj+1)
-cvfi(klon_glo) = cv((ii+1)*jj-ii)
+cufi(klon_glo) = cudyn((ii+1)*jj+1)
+cvfi(klon_glo) = cvdyn((ii+1)*jj-ii)
 
 ! build airefi(), mesh area on physics grid
 allocate(airefi(klon_glo))
-CALL gr_dyn_fi(1,ii+1,jj+1,klon_glo,aire,airefi)
+CALL gr_dyn_fi(1,ii+1,jj+1,klon_glo,airedyn,airefi)
 ! Poles are single points on physics grid
-airefi(1)=sum(aire(1:ii,1))
-airefi(klon_glo)=sum(aire(1:ii,jj+1))
+airefi(1)=sum(airedyn(1:ii,1))
+airefi(klon_glo)=sum(airedyn(1:ii,jj+1))
 
 ! Sanity check: do total planet area match between physics and dynamics?
-total_area_dyn=sum(aire(1:ii,1:jj+1))
+total_area_dyn=sum(airedyn(1:ii,1:jj+1))
 total_area_phy=sum(airefi(1:klon_glo))
 IF (total_area_dyn/=total_area_phy) THEN
@@ -174 +179 @@
 call ini_planete_mod(nlayer,preff,ap,bp)
 
+! for slab ocean, copy over some arrays
+ok_slab_ocean=.false. ! default value
+call getin_p("ok_slab_ocean",ok_slab_ocean)
+if (ok_slab_ocean) then
+  call ini_surf_heat_transp(ip1jm,ip1jmp1,unsairez,fext,unsaire, &
+                            cu,cuvsurcv,cv,cvusurcu,aire,apoln,apols, &
+                            aireu,airev)
+endif
+
 ! copy some fundamental parameters to physics 
 ! and do some initializations