Changeset 1279 for LMDZ4/trunk/libf/phylmd/hines_gwd.F

Timestamp:

Dec 10, 2009, 10:02:56 AM (15 years ago)

Author:

Laurent Fairhead

Message:

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/hines_gwd.F (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/hines_gwd.F

@@ +1 @@
+!
+! $Id$
+!
       SUBROUTINE HINES_GWD(NLON,NLEV,DTIME,paphm1x, papm1x,
      I      rlat,tx,ux,vx,
@@ -1666 +1669 @@
 C  the variances.
 C
-      DO 80 N = 1,NAZ
-        DO 70 I = IL1,IL2
-          IF (I_ALPHA(I,N).LT.0.)  THEN
-            WRITE (6,*) 
-            WRITE (6,*) '******************************'
-            WRITE (6,*) 'Hines integral I_ALPHA < 0 '
-            WRITE (6,*) '  longitude I=',I
-            WRITE (6,*) '  azimuth   N=',N
-            WRITE (6,*) '  level   LEV=',LEV
-            WRITE (6,*) '  I_ALPHA =',I_ALPHA(I,N)
-            WRITE (6,*) '  V_ALPHA =',V_ALPHA(I,LEV,N)
-            WRITE (6,*) '  M_ALPHA =',M_ALPHA(I,LEV,N)
-            WRITE (6,*) '  Q_ALPHA =',V_ALPHA(I,LEV,N) / BVFB(I)
-            WRITE (6,*) '  QM      =',V_ALPHA(I,LEV,N) / BVFB(I) 
-     ^                                * M_ALPHA(I,LEV,N)
-            WRITE (6,*) '******************************'
-          END IF
- 70     CONTINUE
- 80   CONTINUE
+c      DO 80 N = 1,NAZ
+c        DO 70 I = IL1,IL2
+c          IF (I_ALPHA(I,N).LT.0.)  THEN
+c            WRITE (6,*) 
+c            WRITE (6,*) '******************************'
+c            WRITE (6,*) 'Hines integral I_ALPHA < 0 '
+c            WRITE (6,*) '  longitude I=',I
+c            WRITE (6,*) '  azimuth   N=',N
+c            WRITE (6,*) '  level   LEV=',LEV
+c            WRITE (6,*) '  I_ALPHA =',I_ALPHA(I,N)
+c            WRITE (6,*) '  V_ALPHA =',V_ALPHA(I,LEV,N)
+c            WRITE (6,*) '  M_ALPHA =',M_ALPHA(I,LEV,N)
+c            WRITE (6,*) '  Q_ALPHA =',V_ALPHA(I,LEV,N) / BVFB(I)
+c            WRITE (6,*) '  QM      =',V_ALPHA(I,LEV,N) / BVFB(I) 
+c     ^                                * M_ALPHA(I,LEV,N)
+c            WRITE (6,*) '******************************'
+c          END IF
+c 70     CONTINUE
+c 80   CONTINUE
 C
       RETURN