Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/phylmd/vdif_kcay.F90

Timestamp:

Jul 23, 2024, 3:29:36 PM (8 weeks ago)

Author:

abarral

Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/vdif_kcay.F90 (modified) (11 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/vdif_kcay.F90

@@ -338 +338 @@
   END IF
 
-  ! (abd 5 2)      print*,'LONG0=',long0
+  ! (abd 5 2)      PRINT*,'LONG0=',long0
 
   ! -----------------------------------------------------------------------
@@ -385 +385 @@
       ! abd
       ! if(ilev.le.57.and.ilev.ge.37) then
-      ! print*,'L=',ilev,'   GN=',gn,'  SM=',sm(igrid,ilev)
+      ! PRINT*,'L=',ilev,'   GN=',gn,'  SM=',sm(igrid,ilev)
       ! endif
       ! --->
@@ -408 +408 @@
         snstable = snstable*snstable
 
-        ! abde       print*,'SN ',ilev,sn(1,ilev),snstable
+        ! abde       PRINT*,'SN ',ilev,sn(1,ilev),snstable
         IF (sn(igrid,ilev)<snstable) THEN
           sn(igrid, ilev) = snstable
@@ -505 +505 @@
       mcstat = sqrt(m2cstat)
 
-      ! abde      print*,'M2 L=',ilev,mpre(igrid,ilev),mcstat
+      ! abde      PRINT*,'M2 L=',ilev,mpre(igrid,ilev),mcstat
 
       ! -----{puis on ecrit la valeur de q qui annule l'equation de m
@@ -534 +534 @@
       m(igrid, ilev) = mcstat
       ! abd       if(ilev.le.57.and.ilev.ge.37) then
-      ! print*,'L=',ilev,'   M2=',m2(igrid,ilev),m2cstat,
+      ! PRINT*,'L=',ilev,'   M2=',m2(igrid,ilev),m2cstat,
       ! s     'N2=',n2(igrid,ilev)
       ! abd       endif
@@ -562 +562 @@
       ! abd
       ! if(ilev.le.57.and.ilev.ge.37) then
-      ! print*,'L=',ilev,'   GN=',gn,'  SM=',sm(igrid,ilev)
+      ! PRINT*,'L=',ilev,'   GN=',gn,'  SM=',sm(igrid,ilev)
       ! endif
 
@@ -595 +595 @@
       sssq = sssq + (plev(1,ilev-1)-plev(1,ilev+1))*q2(1, ilev)
     END DO
-    ! print*,'Q2moy avant',sssq/sss
-    ! print*,'Q2q20 ',(q2(1,ilev),ilev=1,10)
-    ! print*,'Q2km0 ',(km(1,ilev),ilev=1,10)
+    ! PRINT*,'Q2moy avant',sssq/sss
+    ! PRINT*,'Q2q20 ',(q2(1,ilev),ilev=1,10)
+    ! PRINT*,'Q2km0 ',(km(1,ilev),ilev=1,10)
     ! ! C'est quoi ca qu'etait dans l'original???
     ! do igrid=1,ngrid
@@ -604 +604 @@
     ! q2s=q2
     ! do iii=1,10
-    ! call vdif_q2(dt,g,rconst,plev,temp,km,q2)
+    ! CALL vdif_q2(dt,g,rconst,plev,temp,km,q2)
     ! do ilev=1,klev+1
     ! write(iii+49,*) q2(1,ilev),zlev(1,ilev)
@@ -611 +611 @@
     ! stop
     ! do ilev=1,klev
-    ! print*,zlev(1,ilev),q2s(1,ilev),q2(1,ilev)
+    ! PRINT*,zlev(1,ilev),q2s(1,ilev),q2(1,ilev)
     ! enddo
     ! q2s=q2-q2s
     ! do ilev=1,klev
-    ! print*,q2s(1,ilev),zlev(1,ilev)
+    ! PRINT*,q2s(1,ilev),zlev(1,ilev)
     ! enddo
     DO ilev = 2, klev - 1
@@ -648 +648 @@
           snstable = snstable*snstable
 
-          ! abde      print*,'SN ',ilev,sn(1,ilev),snstable
+          ! abde      PRINT*,'SN ',ilev,sn(1,ilev),snstable
           IF (sn(igrid,ilev)<snstable) THEN
             sn(igrid, ilev) = snstable
@@ -667 +667 @@
       END DO
     END DO
-    ! print*,'Q2km1 ',(km(1,ilev),ilev=1,10)
+    ! PRINT*,'Q2km1 ',(km(1,ilev),ilev=1,10)
 
   END IF