Changeset 1528 for trunk/LMDZ.MARS/libf/phymars/getslopes.F90

Timestamp:

Apr 2, 2016, 4:09:43 PM (9 years ago)

Author:

emillour

Message:

Mars GCM:

• 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.
• Added "ioipsl_getin_p_mod.F90" (getin_p routine) in phy_common to correctly read in parameters from *.def files in a parallel environment.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/getslopes.F90 (modified) (6 diffs)

trunk/LMDZ.MARS/libf/phymars/getslopes.F90

@@ -3 +3 @@
 use comgeomfi_h, only: long, lati
 use slope_mod, only: theta_sl, psi_sl
-USE comcstfi_h
+use comcstfi_h, only: g, rad, pi
+use mod_phys_lmdz_para, only: is_parallel
+use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
 implicit none
 
-#include "dimensions.h"
 
 ! This routine computes slope inclination and orientation for the GCM (callslope=.true. in callphys.def)
@@ -16 +17 @@
 integer,intent(in) :: ngrid ! nnumber of atmospheric columns
 real,intent(in) :: geopot(ngrid)     ! geopotential on phy grid
-real topogrid(iim,jjm+1) ! topography on lat/lon grid with poles and only one -180/180 point
-real latigrid(iim,jjm+1),longgrid(iim,jjm+1) ! meshgrid of latitude and longitude values (radians)
+real topogrid(nbp_lon,nbp_lat) ! topography on lat/lon grid with poles and only one -180/180 point
+real latigrid(nbp_lon,nbp_lat),longgrid(nbp_lon,nbp_lat) ! meshgrid of latitude and longitude values (radians)
 real theta_val ! slope inclination
 real psi_val   ! slope orientation
-real gradx(iim,jjm+1) ! x: latitude-wise topography gradient,  increasing northward
-real grady(iim,jjm+1) ! y: longitude-wise topography gradient, increasing westward
+real gradx(nbp_lon,nbp_lat) ! x: latitude-wise topography gradient,  increasing northward
+real grady(nbp_lon,nbp_lat) ! y: longitude-wise topography gradient, increasing westward
 integer i,j,ig0
 integer id2,idm1 ! a trick to compile testphys1d with debug option
 
+if (is_parallel) then
+  ! This routine only works in serial mode so stop now.
+  write(*,*) "getslopes Error: this routine is not designed to run in parallel"
+  stop
+endif
+
 id2  = 2
-idm1 = iim-1
+idm1 = nbp_lon-1
 
 ! rearrange topography on a 2d array
-do j=2,jjm
-   ig0= 1+(j-2)*iim
-   do i=1,iim
+do j=2,nbp_lat-1
+   ig0= 1+(j-2)*nbp_lon
+   do i=1,nbp_lon
       topogrid(i,j)=geopot(ig0+i)/g
       latigrid(i,j)=lati(ig0+i)
@@ -41 +48 @@
 latigrid(:,1) = lati(1)
 longgrid(:,1) = long(1)
-topogrid(:,jjm+1) = geopot(ngrid)/g
-latigrid(:,jjm+1) = lati(ngrid)
-longgrid(:,jjm+1) = long(ngrid)
+topogrid(:,nbp_lat) = geopot(ngrid)/g
+latigrid(:,nbp_lat) = lati(ngrid)
+longgrid(:,nbp_lat) = long(ngrid)
 
 
@@ -49 +56 @@
 ! compute topography gradient
 ! topogrid and rad are both in meters
-do j=2,jjm
-   do i=1,iim
+do j=2,nbp_lat-1
+   do i=1,nbp_lon
      gradx(i,j) = (topogrid(i,j+1) - topogrid(i,j-1)) / (latigrid(i,j+1)-latigrid(i,j-1))
      gradx(i,j) = gradx(i,j) / rad
    enddo
-   grady(1,j) = (topogrid(id2,j) - topogrid(iim,j)) / (2*pi+longgrid(id2,j)-longgrid(iim,j)) 
+   grady(1,j) = (topogrid(id2,j) - topogrid(nbp_lon,j)) / (2*pi+longgrid(id2,j)-longgrid(nbp_lon,j)) 
    grady(1,j) = grady(1,j) / rad
-   grady(iim,j) = (topogrid(1,j) - topogrid(idm1,j)) / (2*pi+longgrid(1,j)-longgrid(idm1,j)) 
-   grady(iim,j) = grady(iim,j) / rad
-   do i=2,iim-1
+   grady(nbp_lon,j) = (topogrid(1,j) - topogrid(idm1,j)) / (2*pi+longgrid(1,j)-longgrid(idm1,j)) 
+   grady(nbp_lon,j) = grady(nbp_lon,j) / rad
+   do i=2,nbp_lon-1
      grady(i,j) = (topogrid(i+1,j) - topogrid(i-1,j)) / (longgrid(i+1,j)-longgrid(i-1,j)) 
      grady(i,j) = grady(i,j) / rad
@@ -66 +73 @@
 gradx(:,1) = 0.
 grady(:,1) = 0.
-gradx(:,jjm+1) = 0.
-grady(:,jjm+1) = 0.
+gradx(:,nbp_lat) = 0.
+grady(:,nbp_lat) = 0.
 
 
@@ -74 +81 @@
 theta_sl(:) = 0.
 psi_sl(:)   = 0.
-do j=2,jjm
-   do i=1,iim
+do j=2,nbp_lat-1
+   do i=1,nbp_lon
    
-     ig0= 1+(j-2)*iim
+     ig0= 1+(j-2)*nbp_lon
    
      theta_val=atan(sqrt( (gradx(i,j))**2 + (grady(i,j))**2 ))