Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/climb_hq_mod.F90

Timestamp:

Jul 24, 2024, 4:23:34 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmdiso/climb_hq_mod.F90 (modified) (12 diffs)

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

@@ -119 +119 @@
 
 #ifdef ISO
-    real, DIMENSION(klon,2:klev) :: gamaxt_tmp
-    real, DIMENSION(klon,klev) :: xt_tmp
-    real, DIMENSION(klon,klev) ::  Ccoef_XT_tmp,Dcoef_XT_tmp
-    real, DIMENSION(klon) ::  Acoef_XT_tmp,Bcoef_XT_tmp
-    integer ixt
+    REAL, DIMENSION(klon,2:klev) :: gamaxt_tmp
+    REAL, DIMENSION(klon,klev) :: xt_tmp
+    REAL, DIMENSION(klon,klev) ::  Ccoef_XT_tmp,Dcoef_XT_tmp
+    REAL, DIMENSION(klon) ::  Acoef_XT_tmp,Bcoef_XT_tmp
+    INTEGER ixt
 #endif 
 ! Include
@@ -132 +132 @@
 #ifdef ISO
 #ifdef ISOVERIF
-   if (iso_eau.gt.0) THEN
+   IF (iso_eau.gt.0) THEN
        do k = 1, klev
           DO i = 1, knon
@@ -315 +315 @@
        enddo ! do ixt=1,ntraciso
 #ifdef ISOVERIF
-        if (iso_eau.gt.0) THEN
+        IF (iso_eau.gt.0) THEN
          do k = 1, klev
           DO i = 1, knon
@@ -375 +375 @@
 !****************************************************************************************
 !!! jyg le 07/02/2012
-!!jyg       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .EQ.1) THEN
        IF (mod(iflag_pbl_split,10) >=1) THEN
 !!! nrlmd le 02/05/2011
@@ -548 +548 @@
     REAL, DIMENSION(ntraciso,klon,klev)               :: xt_new
     REAL, DIMENSION(ntraciso,klon), INTENT(IN)        :: flx_xt1
-    integer ixt
+    INTEGER ixt
 !#ifdef ISOVERIF
 !    integer iso_verif_noNaN_nostop
@@ -571 +571 @@
        DO k = 1, klev
         DO i = 1, knon  
-         if (iso_eau.gt.0) THEN
+         IF (iso_eau.gt.0) THEN
            CALL iso_verif_egalite(xt_old(iso_eau,i,k), &
                q_old(i,k),'climb_hq_mod 421')
@@ -598 +598 @@
 
 !!! jyg le 07/02/2012
-!!jyg       IF (mod(iflag_pbl_split,2) .eq.1) THEN
+!!jyg       IF (mod(iflag_pbl_split,2) .EQ.1) THEN
        IF (mod(iflag_pbl_split,10) >=1) THEN
 !!! nrlmd le 02/05/2011
@@ -674 +674 @@
        DO i = 1, knon     
         do ixt=1,ntraciso
-         if (iso_verif_noNaN_nostop(xt_new(ixt,i,k),'climb_hq 507').eq.1) THEN
+         IF (iso_verif_noNaN_nostop(xt_new(ixt,i,k),'climb_hq 507').EQ.1) THEN
              WRITE(*,*) 'Acoef_XT(ixt,i)=',Acoef_XT(ixt,i)
              WRITE(*,*) 'Bcoef_XT(ixt,i)=',Bcoef_XT(ixt,i)
              WRITE(*,*) 'flx_xt1(ixt,i)=',flx_xt1(ixt,i)
-             if (k.ge.2) THEN
+             IF (k.ge.2) THEN
                WRITE(*,*) 'Ccoef_XT(ixt,i,k)=',Ccoef_XT(ixt,i,k)
                WRITE(*,*) 'Dcoef_XT(ixt,i,k)=',Dcoef_XT(ixt,i,k)
@@ -689 +689 @@
 #endif
 #ifdef ISOVERIF
-     if (iso_eau.gt.0) THEN
+     IF (iso_eau.gt.0) THEN
         CALL iso_verif_egalite_vect2D( &
                  xt_new,q_new, &
                  'climb_hq_mod 504',ntraciso,klon,klev)
      endif !if (iso_eau.gt.0) THEN
-     if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) THEN
+     IF ((iso_HDO.gt.0).AND.(iso_O18.gt.0)) THEN
         do i=1,klon 
          do k=1,klev 
-            if (q_new(i,k).gt.ridicule) THEN
-               if (iso_verif_o18_aberrant_nostop( &
+            IF (q_new(i,k).gt.ridicule) THEN
+               IF (iso_verif_o18_aberrant_nostop( &
                     xt_new(iso_HDO,i,k)/q_new(i,k), &
                     xt_new(iso_O18,i,k)/q_new(i,k), &
-                    'climb_hq_mod 690').eq.1) THEN
+                    'climb_hq_mod 690').EQ.1) THEN
                   WRITE(*,*) 'i,k,q_new(i,k)=',i,k,q_new(i,k)
                   stop
@@ -708 +708 @@
           enddo !k=1,klev
          enddo  !i=1,klon
-        endif !if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) THEN
+        endif !if ((iso_HDO.gt.0).AND.(iso_O18.gt.0)) THEN
 #endif
 #endif   
@@ -770 +770 @@
 #endif
 #ifdef ISOVERIF
-        if (iso_eau.gt.0) THEN
+        IF (iso_eau.gt.0) THEN
            CALL iso_verif_egalite(d_xt(iso_eau,i,k), &
                d_q(i,k),'climb_hq_mod 503')
@@ -787 +787 @@
 #ifdef ISOVERIF
 !        WRITE(*,*) 'climb_hq_mod 758: d_xt,d_q=',d_xt(iso_eau,1,1),d_q(1,1)
-        if (iso_eau.gt.0) THEN
+        IF (iso_eau.gt.0) THEN
          CALL iso_verif_egalite_vect2D( &
                 d_xt,d_q, &