Changeset 2399 for LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/init_ssrf_m.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/dynlonlat_phylonlat/phylmd/init_ssrf_m.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/init_ssrf_m.F90

@@ -3 +3 @@
 !*******************************************************************************
 
-  USE indice_sol_mod
+  USE indice_sol_mod, ONLY: is_ter, is_oce, is_oce, is_lic, epsfra
   USE dimphy,             ONLY: klon, zmasq
-  USE phys_state_var_mod, ONLY: pctsrf, rlon, rlat
+  USE phys_state_var_mod, ONLY: pctsrf
+  USE geometry_mod, ONLY : longitude_deg, latitude_deg
   USE grid_atob_m,        ONLY: grille_m
   USE ioipsl,             ONLY: flininfo, flinopen, flinget, flinclo
@@ -45 +46 @@
 
 !--- Physical grid points coordinates
-  DO j=2,jjm; rlat((j-2)*iim+2:(j-1)*iim+1)=rlatu(j);    END DO
-  DO j=2,jjm; rlon((j-2)*iim+2:(j-1)*iim+1)=rlonv(1:im); END DO
-  rlat(1) = pi/2.; rlat(klon) = - pi/2.; rlat(:)=rlat(:)/deg2rad
-  rlon(1) = 0.0;   rlon(klon) = 0.0;     rlon(:)=rlon(:)/deg2rad
+  DO j=2,jjm; latitude_deg((j-2)*iim+2:(j-1)*iim+1)=rlatu(j);    END DO
+  DO j=2,jjm; longitude_deg((j-2)*iim+2:(j-1)*iim+1)=rlonv(1:im); END DO
+  latitude_deg(1) = pi/2.; latitude_deg(klon) = - pi/2.
+  latitude_deg(:)=latitude_deg(:)/deg2rad
+  longitude_deg(1) = 0.0;   longitude_deg(klon) = 0.0;
+  longitude_deg(:)=longitude_deg(:)/deg2rad
 
 ! Compute ground geopotential, sub-cells quantities and possibly the mask.