Changeset 2399 for LMDZ5/trunk/libf/phylmd/hgardfou.F90

Timestamp:

Nov 20, 2015, 5:23:28 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Follow-up from commit 2395: get rid of rlon and rlat, longitude_deg and latitude_deg (from module geometry_mod) should be used instead. Longitudes and latitudes are no longer loaded from startphy.nc but inherited from dynamics (and compatibility with values in startphy.nc is checked). This will change bench results because of roundoffs differences between the two.
EM

File:

• LMDZ5/trunk/libf/phylmd/hgardfou.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/hgardfou.F90

@@ -2 +2 @@
 ! $Id$
 SUBROUTINE hgardfou(t, tsol, text,abortphy)
-  USE dimphy
-  USE phys_state_var_mod
-  USE indice_sol_mod
+  USE dimphy, ONLY: klon, klev
+  USE phys_state_var_mod, ONLY: pctsrf
+  USE geometry_mod, ONLY: longitude_deg, latitude_deg
+  USE indice_sol_mod, ONLY: nbsrf
   USE print_control_mod, ONLY: lunout
   IMPLICIT NONE
@@ -55 +56 @@
       DO i = 1, jbad
         WRITE (lunout, *) 'i,k,temperature,lon,lat,pourc ter,lic,oce,sic =', &
-          jadrs(i), k, zt(jadrs(i)), rlon(jadrs(i)), rlat(jadrs(i)), &
-          (pctsrf(jadrs(i),nsrf), nsrf=1, nbsrf)
+          jadrs(i), k, zt(jadrs(i)), longitude_deg(jadrs(i)), &
+          latitude_deg(jadrs(i)),(pctsrf(jadrs(i),nsrf), nsrf=1, nbsrf)
       END DO
     END IF
@@ -75 +76 @@
       DO i = 1, jbad
         WRITE (lunout, *) 'i,k,temperature,lon,lat,pourc ter,lic,oce,sic =', &
-          jadrs(i), k, zt(jadrs(i)), rlon(jadrs(i)), rlat(jadrs(i)), &
-          (pctsrf(jadrs(i),nsrf), nsrf=1, nbsrf)
+          jadrs(i), k, zt(jadrs(i)), longitude_deg(jadrs(i)), &
+          latitude_deg(jadrs(i)), (pctsrf(jadrs(i),nsrf), nsrf=1, nbsrf)
       END DO
     END IF
@@ -101 +102 @@
         WRITE (lunout, *) &
           'i,nsrf,temperature,lon,lat,pourc ter,lic,oce,sic =', jadrs(i), &
-          nsrf, zt(jadrs(i)), rlon(jadrs(i)), rlat(jadrs(i)), &
-          pctsrf(jadrs(i), nsrf)
+          nsrf, zt(jadrs(i)), longitude_deg(jadrs(i)), &
+          latitude_deg(jadrs(i)), pctsrf(jadrs(i), nsrf)
       END DO
     END IF
@@ -122 +123 @@
         WRITE (lunout, *) &
           'i,nsrf,temperature,lon,lat,pourc ter,lic,oce,sic =', jadrs(i), &
-          nsrf, zt(jadrs(i)), rlon(jadrs(i)), rlat(jadrs(i)), &
-          pctsrf(jadrs(i), nsrf)
+          nsrf, zt(jadrs(i)), longitude_deg(jadrs(i)), &
+          latitude_deg(jadrs(i)), pctsrf(jadrs(i), nsrf)
       END DO
     END IF