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

Timestamp:

Apr 17, 2016, 10:23:21 PM (9 years ago)

Author:

emillour

Message:

Mars GCM:

• fix for 1D in writediagfi to enable writing at "ecritphy" rate.
• removed iniprint.h from phymars/dyn1d since it is in "misc"
• Some code cleanup in anticipation of future updates:
  • changed variable names in comgeomphy.F90: give them more explicit names: rlond => longitude , rlatd => latitude, airephy => cell_area, cuphy => dx , cvphy => dy
  • removed long(), lati() and area() from comgeomfi_h.F90, use longitude(), latitude() and cell_are() from comgeomphy.F90 instead

EM

File:

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

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

@@ -1 +1 @@
 subroutine getslopes(ngrid,geopot)
     
-use comgeomfi_h, only: long, lati
+use comgeomphy, only: longitude, latitude ! in radians
 use slope_mod, only: theta_sl, psi_sl
 use comcstfi_h, only: g, rad, pi
@@ -40 +40 @@
    do i=1,nbp_lon
       topogrid(i,j)=geopot(ig0+i)/g
-      latigrid(i,j)=lati(ig0+i)
-      longgrid(i,j)=long(ig0+i)
+      latigrid(i,j)=latitude(ig0+i)
+      longgrid(i,j)=longitude(ig0+i)
    enddo
 enddo
 !poles :
 topogrid(:,1) = geopot(1)/g
-latigrid(:,1) = lati(1)
-longgrid(:,1) = long(1)
+latigrid(:,1) = latitude(1)
+longgrid(:,1) = longitude(1)
 topogrid(:,nbp_lat) = geopot(ngrid)/g
-latigrid(:,nbp_lat) = lati(ngrid)
-longgrid(:,nbp_lat) = long(ngrid)
+latigrid(:,nbp_lat) = latitude(ngrid)
+longgrid(:,nbp_lat) = longitude(ngrid)