Changeset 1541 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

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/dynphy_lonlat/phymars/iniphysiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/iniphysiq_mod.F90

@@ -16 +16 @@
 
 use comgeomphy, only : initcomgeomphy, &
-                       airephy, & ! physics grid area (m2)
-                       cuphy, & ! cu coeff. (u_covariant = cu * u)
-                       cvphy, & ! cv coeff. (v_covariant = cv * v)
-                       rlond, & ! longitudes
-                       rlatd ! latitudes
+                       cell_area, & ! physics grid area (m2)
+                       dx, & ! cu coeff. (u_covariant = cu * u)
+                       dy, & ! cv coeff. (v_covariant = cv * v)
+                       longitude, & ! longitudes (rad)
+                       latitude ! latitudes (rad)
 use infotrac, only : nqtot ! number of advected tracers
 use comgeomfi_h, only: ini_fillgeom
@@ -164 +164 @@
       
 offset=klon_mpi_begin-1
-airephy(1:klon_omp)=airefi(offset+klon_omp_begin:offset+klon_omp_end)
-cuphy(1:klon_omp)=cufi(offset+klon_omp_begin:offset+klon_omp_end)
-cvphy(1:klon_omp)=cvfi(offset+klon_omp_begin:offset+klon_omp_end)
-rlond(1:klon_omp)=lonfi(offset+klon_omp_begin:offset+klon_omp_end)
-rlatd(1:klon_omp)=latfi(offset+klon_omp_begin:offset+klon_omp_end)
+cell_area(1:klon_omp)=airefi(offset+klon_omp_begin:offset+klon_omp_end)
+dx(1:klon_omp)=cufi(offset+klon_omp_begin:offset+klon_omp_end)
+dy(1:klon_omp)=cvfi(offset+klon_omp_begin:offset+klon_omp_end)
+longitude(1:klon_omp)=lonfi(offset+klon_omp_begin:offset+klon_omp_end)
+latitude(1:klon_omp)=latfi(offset+klon_omp_begin:offset+klon_omp_end)
 
 ! copy some fundamental parameters to physics 
@@ -175 +175 @@
                          day_ini,hour_ini,punjours,ptimestep, &
                          prad,pg,pr,pcpp)
-call ini_fillgeom(klon_omp,rlatd,rlond,airephy)
+call ini_fillgeom(klon_omp,latitude,longitude,cell_area)
 call conf_phys(klon_omp,nlayer,nqtot)