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

LMDZ5/trunk/libf/phylmd/add_phys_tend.F90

@@ -14 +14 @@
 !======================================================================
 
-use dimphy
-use phys_local_var_mod
-use phys_state_var_mod
-use print_control_mod, only: prt_level
+USE dimphy, ONLY: klon, klev
+USE phys_local_var_mod, ONLY: u_seri, v_seri, ql_seri, qs_seri, q_seri, &
+                              t_seri
+USE phys_state_var_mod, ONLY: ftsol, zdu
+USE geometry_mod, ONLY: longitude_deg, latitude_deg
+USE print_control_mod, ONLY: prt_level
 IMPLICIT none
   include "YOMCST.h"
@@ -106 +108 @@
          i=jadrs(j)
          if(prt_level.ge.debug_level) THEN
-          print*,'PLANTAGE POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
+          print*,'PLANTAGE POUR LE POINT i lon lat =',&
+                 i,longitude_deg(i),latitude_deg(i),text
           print*,'l    T     dT       Q     dQ    '
           DO k = 1, klev
@@ -124 +127 @@
         i=jqadrs(j)
           if(prt_level.ge.debug_level) THEN
-           print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
+           print*,'WARNING  : EAU POUR LE POINT i lon lat =',&
+                  i,longitude_deg(i),latitude_deg(i),text
            print*,'l    T     dT       Q     dQ    '
            DO k = 1, klev
@@ -205 +209 @@
          k=kadrs(j)
          if(prt_level.ge.debug_level) THEN
-          print*,'PLANTAGE2 POUR LE POINT i itap rlon rlat txt jbad zdt t',i,itap,rlon(i),rlat(i),text,jbad, &
+          print*,'PLANTAGE2 POUR LE POINT i itap lon lat txt jbad zdt t',&
+                 i,itap,longitude_deg(i),latitude_deg(i),text,jbad, &
        &        zdt(i,k),t_seri(i,k)-zdt(i,k)
 !!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN 
@@ -222 +227 @@
          k=kqadrs(j)
          if(prt_level.ge.debug_level) THEN
-          print*,'WARNING  : EAU2 POUR LE POINT i itap rlon rlat txt jqbad zdq q zdql ql',i,itap,rlon(i),rlat(i),text,jqbad,&
+          print*,'WARNING  : EAU2 POUR LE POINT i itap lon lat txt jqbad zdq q zdql ql',&
+                 i,itap,longitude_deg(i),latitude_deg(i),text,jqbad,&
        &        zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)
 !!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN