Changeset 130 for trunk/libf/dyn3d

Timestamp:

May 25, 2011, 3:45:43 PM (14 years ago)

Author:

slebonnois

Message:

Sebastien:

• correction de bugs dans phytitan suite a compil avec -debug sur gnome
• elimination de nbetat* dans gcm.F et leapfrog.F (et dans le parallele)
• correction rday -> rjour dans sortvarc.F
• ajustement des arch-GNOME*

Location:

trunk/libf/dyn3d

Files:

• gcm.F (modified) (2 diffs)
• leapfrog.F (modified) (3 diffs)
• sortvarc.F (modified) (3 diffs)

trunk/libf/dyn3d/gcm.F ¶

@@ -85 +85 @@
 
       REAL zdtvr
-      INTEGER nbetatmoy, nbetatdem,nbetat
 
 c   variables dynamiques
@@ -397 +396 @@
 #endif
 
-c  nombre d'etats dans les fichiers demarrage et histoire
-      nbetatdem = nday / iecri
-      nbetatmoy = nday / periodav + 1
-
       if (iflag_phys.eq.1) then
       ! these initialisations have already been done (via iniacademic)

trunk/libf/dyn3d/leapfrog.F ¶

@@ -68 +68 @@
 
       real zqmin,zqmax
-      INTEGER nbetatmoy, nbetatdem,nbetat
 
 c   variables dynamiques
@@ -694 +693 @@
              ! Ehouarn: output only during LF or Backward Matsuno
              if (leapf.or.(.not.leapf.and.(.not.forward))) then
-              nbetat = nbetatdem
 ! ADAPTATION GCM POUR CP(T)
               call tpot2t(ijp1llm,teta,temp,pk)
@@ -821 +819 @@
               IF(MOD(itau,iecri         ).EQ.0) THEN
 c              IF(MOD(itau,iecri*day_step).EQ.0) THEN
-                nbetat = nbetatdem
 ! ADAPTATION GCM POUR CP(T)
                 call tpot2t(ijp1llm,teta,temp,pk)

trunk/libf/dyn3d/sortvarc.F ¶

@@ -53 +53 @@
       REAL cosphi(ip1jm),omegcosp(ip1jm)
       REAL dtvrs1j,rjour,heure,radsg,radomeg
-      REAL rday, massebxy(ip1jm,llm)
+      REAL massebxy(ip1jm,llm)
       INTEGER  l, ij, imjmp1
 
@@ -143 +143 @@
 
       IF (firstcal.and.resetvarc) then
-       rday = REAL(INT(time-jD_ref-jH_ref))
-         PRINT 3500, itau, rday, heure,time
+         PRINT 3500, itau, rjour, heure,time
          PRINT*,'WARNING!!! On recalcule les valeurs initiales de :'
          PRINT*,'ptot,rmsdpdt,etot,ztot,stot,rmsv,ang'
@@ -164 +163 @@
       firstcal = .false.
 
-      PRINT 3500, itau, rday, heure, time
+      PRINT 3500, itau, rjour, heure, time
       PRINT 4000, ptot,rmsdpdt,etot,ztot,stot,rmsv,ang