Changeset 130 for trunk/libf

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

Files:

• dyn3d/gcm.F (modified) (2 diffs)
• dyn3d/leapfrog.F (modified) (3 diffs)
• dyn3d/sortvarc.F (modified) (3 diffs)
• dyn3dpar/gcm.F (modified) (2 diffs)
• dyn3dpar/leapfrog_p.F (modified) (3 diffs)
• dyn3dpar/sortvarc.F (modified) (3 diffs)
• phytitan/optci.F (modified) (1 diff)
• phytitan/optci_1pt.F (modified) (2 diffs)
• phytitan/pg3.F (modified) (1 diff)

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
 

trunk/libf/dyn3dpar/gcm.F ¶

@@ -82 +82 @@
 
       REAL zdtvr
-c      INTEGER nbetatmoy, nbetatdem,nbetat
-      INTEGER nbetatmoy, nbetatdem
 
 c   variables dynamiques
@@ -416 +414 @@
 #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/dyn3dpar/leapfrog_p.F ¶

@@ -71 +71 @@
       
       real zqmin,zqmax
-      INTEGER nbetatmoy, nbetatdem,nbetat
 
 c   variables dynamiques
@@ -1519 +1518 @@
              if (leapf.or.(.not.leapf.and.(.not.forward))) then
 c$OMP BARRIER
-c$OMP MASTER
-              nbetat = nbetatdem
-c$OMP END MASTER
-c$OMP BARRIER
 
 ! ADAPTATION GCM POUR CP(T)
@@ -1747 +1742 @@
 c              IF(MOD(itau,iecri*day_step).EQ.0) THEN
 c$OMP BARRIER
-c$OMP MASTER
-                nbetat = nbetatdem
-c$OMP END MASTER
-c$OMP BARRIER
 
 ! ADAPTATION GCM POUR CP(T)

trunk/libf/dyn3dpar/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
 

trunk/libf/phytitan/optci.F ¶

@@ -74 +74 @@
 
       integer iopti,iwarning     ! iopti: premier appel, une seule boucle sur les l.d'o.
-      integer ig_,seulmtunpt
+      integer ig,seulmtunpt
       save iopti,iwarning,seulmtunpt
       data iopti,iwarning,seulmtunpt/0,0,0/

trunk/libf/phytitan/optci_1pt.F ¶

@@ -318 +318 @@
 C     ??FLAG? HERE MUST BE WATCHED CAREFULLY
                      U=COLDEN(J)*6.02204E23/BMU
-          if(ig.eq.1.and.k.eq.nspecv/2) print*,'Avant GAS2'
+c         if(ig.eq.1.and.k.eq.nspecv/2) print*,'Avant GAS2'
                      if((ylellouch).or.(.not.hcnrad)) then
                        CALL GAS2_NOHCN(J, KGAS,TBAR,PBAR,U,TAU2)
@@ -324 +324 @@
                        CALL GAS2(J, KGAS,TBAR,PBAR,U,TAU2)
                      endif
-          if(ig.eq.1.and.k.eq.nspecv/2) print*,'Apres GAS2'
+c         if(ig.eq.1.and.k.eq.nspecv/2) print*,'Apres GAS2'
                      TAUGAS=TAUGAS+TAU2
        ENDIF

trunk/libf/phytitan/pg3.F ¶

@@ -1204 +1204 @@
          
 
-c        p0=.3
+c        p0=0.3
          p0=1.
-         do i=1,nz
+         do i=1,nz-1
           if (pb(i).lt.p0.and.pb(i+1).gt.p0) zalt0=(z(i)+z(i+1))/2.
          enddo