Changeset 4064 for LMDZ6/trunk/libf/dyn3d/vlsplt.F

Timestamp:

Jan 25, 2022, 7:14:39 PM (2 years ago)

Author:

dcugnet

Message:

• minor fixes (unused variables suppressed, comas after a WRITE() statement, etc.)
• parser routines taken from version 7 of ​https://svn.lmd.jussieu.fr/tracers-parser
• few changes in infotrac, and few fixes of (at least) the sequential version:
  • uadv and vadv were deallocated twice (fix was lost by mistake just before last commit)
  • in [( dum(im), im=1, nm)] implicit loops, ifort evaluates "dum(im)" even if nm==0,

resulting in a crash, "im" being unitialized.

File:

• LMDZ6/trunk/libf/dyn3d/vlsplt.F (modified) (13 diffs)

LMDZ6/trunk/libf/dyn3d/vlsplt.F

@@ -38 +38 @@
 c   ---------
 c
-      INTEGER i,ij,l,j,ii
-      INTEGER ijlqmin,iqmin,jqmin,lqmin
-c
-      REAL zm(ip1jmp1,llm,nqtot),newmasse
+      INTEGER ij,l
+c
+      REAL zm(ip1jmp1,llm,nqtot)
       REAL mu(ip1jmp1,llm)
       REAL mv(ip1jm,llm)
       REAL mw(ip1jmp1,llm+1)
-      REAL zq(ip1jmp1,llm,nqtot),zz
-      REAL dqx(ip1jmp1,llm),dqy(ip1jmp1,llm),dqz(ip1jmp1,llm)
-      REAL second,temps0,temps1,temps2,temps3
-      REAL ztemps1,ztemps2,ztemps3
+      REAL zq(ip1jmp1,llm,nqtot)
       REAL zzpbar, zzw
-      LOGICAL testcpu
-      SAVE testcpu
-      SAVE temps1,temps2,temps3
-      INTEGER iminn,imaxx
       INTEGER ifils,iq2 ! CRisi
 
       REAL qmin,qmax
       DATA qmin,qmax/0.,1.e33/
-      DATA testcpu/.false./
-      DATA temps1,temps2,temps3/0.,0.,0./
-
 
         zzpbar = 0.5 * pdt
@@ -157 +146 @@
 c   ----------
       REAL masse(ip1jmp1,llm,nqtot),pente_max
-      REAL u_m( ip1jmp1,llm ),pbarv( iip1,jjm,llm)
+      REAL u_m( ip1jmp1,llm )
       REAL q(ip1jmp1,llm,nqtot)
-      REAL w(ip1jmp1,llm)
       INTEGER iq ! CRisi
 c
@@ -169 +157 @@
 c
       REAL new_m,zu_m,zdum(ip1jmp1,llm)
-      REAL sigu(ip1jmp1),dxq(ip1jmp1,llm),dxqu(ip1jmp1)
+c      REAL sigu(ip1jmp1)
+      REAL dxq(ip1jmp1,llm),dxqu(ip1jmp1)
       REAL zz(ip1jmp1)
       REAL adxqu(ip1jmp1),dxqmax(ip1jmp1,llm)
@@ -178 +167 @@
       INTEGER ifils,iq2 ! CRisi
 
-      Logical extremum,first,testcpu
-      SAVE first,testcpu
-
-      REAL      SSUM
-      REAL temps0,temps1,temps2,temps3,temps4,temps5,second
-      SAVE temps0,temps1,temps2,temps3,temps4,temps5
-
-      REAL z1,z2,z3
-
-      DATA first,testcpu/.true.,.false./
-
-      IF(first) THEN
-         temps1=0.
-         temps2=0.
-         temps3=0.
-         temps4=0.
-         temps5=0.
-         first=.false.
-      ENDIF
+      Logical first
+      SAVE first
+      DATA first/.true./
 
 c   calcul de la pente a droite et a gauche de la maille
@@ -432 +405 @@
          ENDDO
       ENDIF  ! n0.gt.0 
-9999    continue
+c9999    continue
 
 
@@ -536 +509 @@
       REAL masse(ip1jmp1,llm,nqtot),pente_max
       REAL masse_adv_v( ip1jm,llm)
-      REAL q(ip1jmp1,llm,nqtot), dq( ip1jmp1,llm)
+      REAL q(ip1jmp1,llm,nqtot)
       INTEGER iq ! CRisi
 c
@@ -545 +518 @@
 c
       REAL airej2,airejjm,airescb(iim),airesch(iim)
-      REAL dyq(ip1jmp1,llm),dyqv(ip1jm),zdvm(ip1jmp1,llm)
+      REAL dyq(ip1jmp1,llm),dyqv(ip1jm)
       REAL adyqv(ip1jm),dyqmax(ip1jmp1)
       REAL qbyv(ip1jm,llm)
 
-      REAL qpns,qpsn,appn,apps,dyn1,dys1,dyn2,dys2,newmasse,fn,fs
+      REAL qpns,qpsn,dyn1,dys1,dyn2,dys2,newmasse,fn,fs
+c     REAL appn apps
 c     REAL newq,oldmasse
-      Logical extremum,first,testcpu
-      REAL temps0,temps1,temps2,temps3,temps4,temps5,second
-      SAVE temps0,temps1,temps2,temps3,temps4,temps5
-      SAVE first,testcpu
+      LOGICAL first
+      SAVE first
 
       REAL convpn,convps,convmpn,convmps
@@ -570 +542 @@
       REAL      SSUM
 
-      DATA first,testcpu/.true.,.false./
-      DATA temps0,temps1,temps2,temps3,temps4,temps5/0.,0.,0.,0.,0.,0./
+      DATA first/.true./
 
       !write(*,*) 'vly 578: entree, iq=',iq
@@ -905 +876 @@
 c   ---------
 c
-      INTEGER i,ij,l,j,ii
+      INTEGER ij,l
 c
       REAL wq(ip1jmp1,llm+1),newmasse
@@ -918 +889 @@
       SAVE testcpu
 
-      REAL temps0,temps1,temps2,temps3,temps4,temps5,second
-      SAVE temps0,temps1,temps2,temps3,temps4,temps5
-      REAL      SSUM
+#ifdef BIDON
+      REAL temps0,temps1,second
+      SAVE temps0,temps1
 
       DATA testcpu/.false./
-      DATA temps0,temps1,temps2,temps3,temps4,temps5/0.,0.,0.,0.,0.,0./
+      DATA temps0,temps1/0.,0./
+#endif
 
 c    On oriente tout dans le sens de la pression c'est a dire dans le
@@ -1083 +1055 @@
       real zzq(iip1,jjp1,llm)
 
+#ifdef isminmax
       integer imin,jmin,lmin,ijlmin
       integer imax,jmax,lmax,ijlmax
@@ -1088 +1061 @@
       integer ismin,ismax
 
-#ifdef isminismax
       call scopy (ip1jmp1*llm,zq,1,zzq,1)
 
@@ -1116 +1088 @@
 #endif
       return
-9999  format(a20,'  q(',i3,',',i2,',',i2,')=',e12.5,e12.5)
+c9999  format(a20,'  q(',i3,',',i2,',',i2,')=',e12.5,e12.5)
       end