Changeset 5098 for LMDZ6/branches/Amaury_dev/libf/dyn3d

Timestamp:

Jul 22, 2024, 6:53:44 PM (6 months ago)

Author:

abarral

Message:

Remove CRAY key (obsolete calls to functions that don't exist anymore, bugs in some implementations, irrelevant now)
Replace usage of CPP_XIOS key by using_xios logical
Remove always unused testcpu bits
Replace most uses of CPP_StratAer by the corresponding logical defined in lmdz_cppkeys_wrapper.F90 [this breaks iso compilation because phyiso doesn't define all aerosols - to be fixed later]
Replaces uses of include "yomcst.h" by the lmdz_yomcst.f90 module in .[fF]90 files

Location:

LMDZ6/branches/Amaury_dev/libf/dyn3d

Files:

• vlsplt.F (modified) (9 diffs)
• vlspltqs.F (modified) (9 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.F

@@ -210 +210 @@
 
             DO ij=iip2+1,ip1jm
-#ifdef CRAY
-               dxq(ij,l)=
-     ,         cvmgp(dxqu(ij-1)+dxqu(ij),0.,dxqu(ij-1)*dxqu(ij))
-#else
                IF(dxqu(ij-1)*dxqu(ij)>0) THEN
                   dxq(ij,l)=dxqu(ij-1)+dxqu(ij)
@@ -220 +216 @@
                   dxq(ij,l)=0.
                ENDIF
-#endif
                dxq(ij,l)=0.5*dxq(ij,l)
                dxq(ij,l)=
@@ -274 +269 @@
 c   calcul des flux a gauche et a droite
 
-#ifdef CRAY
-
-      DO l=1,llm
-       DO ij=iip2,ip1jm-1
-          zdum(ij,l)=cvmgp(1.-u_m(ij,l)/masse(ij,l,iq),
-     ,                     1.+u_m(ij,l)/masse(ij+1,l,iq),
-     ,                     u_m(ij,l))
-          zdum(ij,l)=0.5*zdum(ij,l)
-          u_mq(ij,l)=cvmgp(
-     ,                q(ij,l,iq)+zdum(ij,l)*dxq(ij,l),
-     ,                q(ij+1,l,iq)-zdum(ij,l)*dxq(ij+1,l),
-     ,                u_m(ij,l))
-          u_mq(ij,l)=u_m(ij,l)*u_mq(ij,l)
-       ENDDO
-      ENDDO
-#else
 c   on cumule le flux correspondant a toutes les mailles dont la masse
 c   au travers de la paroi pENDant le pas de temps.
@@ -307 +286 @@
        ENDDO
       ENDDO
-#endif
 
 c       go to 9999
@@ -883 +861 @@
       INTEGER ifils,iq2 ! CRisi
 
-      LOGICAL testcpu
-      SAVE testcpu
-
 #ifdef BIDON
       REAL temps0,temps1,second
       SAVE temps0,temps1
 
-      DATA testcpu/.false./
       DATA temps0,temps1/0.,0./
 #endif
@@ -899 +873 @@
       !write(*,*) 'vlz 923: entree'
 
-#ifdef BIDON
-      IF(testcpu) THEN
-         temps0=second(0.)
-      ENDIF
-#endif
       DO l=2,llm
          DO ij=1,ip1jmp1
@@ -913 +882 @@
       DO l=2,llm-1
          DO ij=1,ip1jmp1
-#ifdef CRAY
-            dzq(ij,l)=0.5*
-     ,      cvmgp(dzqw(ij,l)+dzqw(ij,l+1),0.,dzqw(ij,l)*dzqw(ij,l+1))
-#else
             IF(dzqw(ij,l)*dzqw(ij,l+1)>0.) THEN
                 dzq(ij,l)=0.5*(dzqw(ij,l)+dzqw(ij,l+1))
@@ -922 +887 @@
                 dzq(ij,l)=0.
             ENDIF
-#endif
             dzqmax=pente_max*min(adzqw(ij,l),adzqw(ij,l+1))
             dzq(ij,l)=sign(min(abs(dzq(ij,l)),dzqmax),dzq(ij,l))
@@ -934 +898 @@
       ENDDO
 
-#ifdef BIDON
-      IF(testcpu) THEN
-         temps1=temps1+second(0.)-temps0
-      ENDIF
-#endif
 c ---------------------------------------------------------------
 c   .... calcul des termes d'advection verticale  .......

LMDZ6/branches/Amaury_dev/libf/dyn3d/vlspltqs.F

@@ -55 +55 @@
       REAL temps1,temps2,temps3
       REAL zzpbar, zzw
-      LOGICAL testcpu
-      SAVE testcpu
       SAVE temps1,temps2,temps3
 
       REAL qmin,qmax
       DATA qmin,qmax/0.,1.e33/
-      DATA testcpu/.false./
       DATA temps1,temps2,temps3/0.,0.,0./
 
@@ -220 +217 @@
       INTEGER ifils,iq2 ! CRisi
 
-      Logical first,testcpu
-      SAVE first,testcpu
+      Logical first
+      SAVE first
 
       REAL      SSUM
@@ -228 +225 @@
 
 
-      DATA first,testcpu/.true.,.false./
+      DATA first/.true./
 
       IF(first) THEN
@@ -278 +275 @@
 
             DO ij=iip2+1,ip1jm
-#ifdef CRAY
-               dxq(ij,l)=
-     ,         cvmgp(dxqu(ij-1)+dxqu(ij),0.,dxqu(ij-1)*dxqu(ij))
-#else
                IF(dxqu(ij-1)*dxqu(ij)>0) THEN
                   dxq(ij,l)=dxqu(ij-1)+dxqu(ij)
@@ -288 +281 @@
                   dxq(ij,l)=0.
                ENDIF
-#endif
                dxq(ij,l)=0.5*dxq(ij,l)
                dxq(ij,l)=
@@ -341 +333 @@
 c   calcul des flux a gauche et a droite
 
-#ifdef CRAY
-c--pas encore modification sur Qsat
-      DO l=1,llm
-       DO ij=iip2,ip1jm-1
-          zdum(ij,l)=cvmgp(1.-u_m(ij,l)/masse(ij,l,iq),
-     ,                     1.+u_m(ij,l)/masse(ij+1,l,iq),
-     ,                     u_m(ij,l))
-          zdum(ij,l)=0.5*zdum(ij,l)
-          u_mq(ij,l)=cvmgp(
-     ,                q(ij,l,iq)+zdum(ij,l)*dxq(ij,l),
-     ,                q(ij+1,l,iq)-zdum(ij,l)*dxq(ij+1,l),
-     ,                u_m(ij,l))
-          u_mq(ij,l)=u_m(ij,l)*u_mq(ij,l)
-       ENDDO
-      ENDDO
-#else
 c   on cumule le flux correspondant a toutes les mailles dont la masse
 c   au travers de la paroi pENDant le pas de temps.
@@ -373 +349 @@
        ENDDO
       ENDDO
-#endif
 
 
@@ -581 +556 @@
       REAL qpns,qpsn,dyn1,dys1,dyn2,dys2,newmasse,fn,fs
 c     REAL newq,oldmasse
-      Logical first,testcpu
+      Logical first
       REAL temps0,temps1,temps2,temps3,temps4,temps5
       SAVE temps0,temps1,temps2,temps3,temps4,temps5
-      SAVE first,testcpu
+      SAVE first
 
       REAL convpn,convps,convmpn,convmps
@@ -598 +573 @@
       REAL      SSUM
 
-      DATA first,testcpu/.true.,.false./
+      DATA first/.true./
       DATA temps0,temps1,temps2,temps3,temps4,temps5/0.,0.,0.,0.,0.,0./