Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (2 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90 (modified) (15 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90

@@ -355 +355 @@
   ! PRINT*, 'wake initialisations'
 !#ifdef ISOVERIF
-!        write(*,*) 'wake 358: entree'
+!        WRITE(*,*) 'wake 358: entree'
 !#endif
 
@@ -398 +398 @@
   ! alpk = 0.05
 
- if (first) then
-
+ if (first) THEN
   igout = klon/2+1/klon
 
@@ -686 +685 @@
   DO k = 1, klev
     DO i = 1, klon
-      ! write(*,*)'wake 1',i,k,rd,te(i,k)
+      ! WRITE(*,*)'wake 1',i,k,rd,te(i,k)
       rho(i, k) = p(i, k)/(rd*te(i,k))
-      ! write(*,*)'wake 2',rho(i,k)
+      ! WRITE(*,*)'wake 2',rho(i,k)
       IF (k==1) THEN
-        ! write(*,*)'wake 3',i,k,rd,te(i,k)
+        ! WRITE(*,*)'wake 3',i,k,rd,te(i,k)
         rhoh(i, k) = ph(i, k)/(rd*te(i,k))
-        ! write(*,*)'wake 4',i,k,rd,te(i,k)
+        ! WRITE(*,*)'wake 4',i,k,rd,te(i,k)
         zhh(i, k) = 0
       ELSE
-        ! write(*,*)'wake 5',rd,(te(i,k)+te(i,k-1))
+        ! WRITE(*,*)'wake 5',rd,(te(i,k)+te(i,k-1))
         rhoh(i, k) = ph(i, k)*2./(rd*(te(i,k)+te(i,k-1)))
-        ! write(*,*)'wake 6',(-rhoh(i,k)*RG)+zhh(i,k-1)
+        ! WRITE(*,*)'wake 6',(-rhoh(i,k)*RG)+zhh(i,k-1)
         zhh(i, k) = (ph(i,k)-ph(i,k-1))/(-rhoh(i,k)*rg) + zhh(i, k-1)
       END IF
-      ! write(*,*)'wake 7',ppi(i,k)
+      ! WRITE(*,*)'wake 7',ppi(i,k)
       the(i, k) = te(i, k)/ppi(i, k)
       thu(i, k) = (te(i,k)-deltatw(i,k)*sigmaw(i))/ppi(i, k)
       tu(i, k) = te(i, k) - deltatw(i, k)*sigmaw(i)
       qu(i, k) = qe(i, k) - deltaqw(i, k)*sigmaw(i)
-      ! write(*,*)'wake 8',(rd*(te(i,k)+deltatw(i,k)))
+      ! WRITE(*,*)'wake 8',(rd*(te(i,k)+deltatw(i,k)))
       rhow(i, k) = p(i, k)/(rd*(te(i,k)+deltatw(i,k)))
       dth(i, k) = deltatw(i, k)/ppi(i, k)
@@ -713 +712 @@
         enddo !do ixt=1,ntraciso
 #ifdef ISOVERIF
-        if (iso_eau.gt.0) then
+        if (iso_eau.gt.0) THEN
               CALL iso_verif_egalite(deltaqw(i,k),deltaxtw(iso_eau,i,k),'wake 723a')
               CALL iso_verif_egalite(qu(i,k),xtu(iso_eau,i,k),'wake 723b')
         endif
-        if (iso_HDO.gt.0) then
+        if (iso_HDO.gt.0) THEN
                if (iso_verif_aberrant_enc_choix_nostop(xtu(iso_hdo,i,k),qu(i,k),ridicule,deltalim, &
-                        'wake 723c xtu').eq.1) then
+                        'wake 723c xtu').eq.1) THEN
                      stop
                endif
@@ -1456 +1455 @@
         enddo !do ixt=1,ntraciso
 #ifdef ISOVERIF
-      if (iso_eau.gt.0) then
+      if (iso_eau.gt.0) THEN
         CALL iso_verif_egalite(deltaqw(i,k-1),deltaxtw(iso_eau,i,k-1),'wake 1460a')
         CALL iso_verif_egalite(deltaqw(i,k),deltaxtw(iso_eau,i,k),'wake 1460b')
@@ -1504 +1503 @@
         enddo !do ixt=1,ntraciso
 #ifdef ISOVERIF
-      if (iso_eau.gt.0) then
+      if (iso_eau.gt.0) THEN
         CALL iso_verif_egalite(d_q1(i,k),d_xt1(iso_eau,i,k),'wake 1502a')
         CALL iso_verif_egalite(d_q2(i,k),d_xt2(iso_eau,i,k),'wake 1502b')
@@ -1632 +1631 @@
         enddo !do ixt=1,ntraciso
 #ifdef ISOVERIF
-      if (iso_eau.gt.0) then
+      if (iso_eau.gt.0) THEN
         CALL iso_verif_egalite(dqke(i,k),dxtKE(iso_eau,i,k),'wake 1621a')
         CALL iso_verif_egalite(dqdwn(i,k),dxtdwn(iso_eau,i,k),'wake 1621b')
@@ -1663 +1662 @@
           ! Jingmei
 
-          ! write(lunout,*)'wake.F ',i,k, dtimesub,d_deltat_gw(i,k),
+          ! WRITE(lunout,*)'wake.F ',i,k, dtimesub,d_deltat_gw(i,k),
           ! &  Tgw(i,k),deltatw(i,k)
           d_deltat_gw(i, k) = d_deltat_gw(i, k) - tgw(i, k)*deltatw(i, k)* &
             dtimesub
-          ! write(lunout,*)'wake.F ',i,k, dtimesub,d_deltatw(i,k)
+          ! WRITE(lunout,*)'wake.F ',i,k, dtimesub,d_deltatw(i,k)
           ff(i) = d_deltatw(i, k)/dtimesub
 
@@ -1709 +1708 @@
         enddo !do ixt=1,ntraciso
 #ifdef ISOVERIF
-      if (iso_eau.gt.0) then
+      if (iso_eau.gt.0) THEN
         CALL iso_verif_egalite(dqke(i,k),dxtKE(iso_eau,i,k),'wake 1692a')
         CALL iso_verif_egalite(deltaqw(i,k),deltaxtw(iso_eau,i,k),'wake 1692b')
@@ -1728 +1727 @@
 #ifdef ISO
 #ifdef ISOVERIF
-      if (iso_eau.gt.0) then
+      if (iso_eau.gt.0) THEN
              CALL iso_verif_egalite_vect2D(d_deltaxtw,d_deltaqw, &
                  'wake 1359',ntraciso,klon,klev)
@@ -1798 +1797 @@
 #ifdef ISO
 #ifdef ISOVERIF
-      if (iso_eau.gt.0) then
+      if (iso_eau.gt.0) THEN
         DO k= 1,klev 
           DO i = 1,klon 
@@ -1805 +1804 @@
           enddo ! DO i = 1,klon     
         enddo ! DO k= 1,klev 
-      endif !if (iso_eau.gt.0) then
+      endif !if (iso_eau.gt.0) THEN
 #endif
 #endif
@@ -1840 +1839 @@
 #ifdef ISO
 #ifdef ISOVERIF
-      write(*,*) 'wake 1859'
-      if (iso_eau.gt.0) then
+      WRITE(*,*) 'wake 1859'
+      if (iso_eau.gt.0) THEN
         DO k= 1,klev 
           DO i = 1,klon 
@@ -1848 +1847 @@
           enddo ! DO i = 1,klon     
         enddo ! DO k= 1,klev 
-      endif !if (iso_eau.gt.0) then     
-      if (iso_hdo.gt.0) then
+      endif !if (iso_eau.gt.0) THEN
+      if (iso_hdo.gt.0) THEN
       CALL iso_verif_aberrant_enc_vect2D( &
             xte,qe, &
@@ -2232 +2231 @@
 #ifdef ISO
 #ifdef ISOVERIF
-      if (iso_hdo.gt.0) then
+      if (iso_hdo.gt.0) THEN
       CALL iso_verif_aberrant_enc_vect2D( &
               xtu,qu, &