Changeset 5160 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90

Timestamp:

Aug 3, 2024, 2:56:58 PM (7 weeks ago)

Author:

abarral

Message:

Put .h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90 (modified) (17 diffs)

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

@@ -202 +202 @@
   INTEGER il
 
-!    print *,'cv3_incrcount : noconv_stop ',noconv_stop
-!    print *,'cv3_incrcount in, sig(1,nd) ',sig(1,nd)
+!    PRINT *,'cv3_incrcount : noconv_stop ',noconv_stop
+!    PRINT *,'cv3_incrcount in, sig(1,nd) ',sig(1,nd)
     IF(ok_convstop) THEN
       DO il = 1, len
@@ -215 +215 @@
       END DO
     ENDIF  ! (ok_convstop)
-!    print *,'cv3_incrcount out, sig(1,nd) ',sig(1,nd)
+!    PRINT *,'cv3_incrcount out, sig(1,nd) ',sig(1,nd)
 
 
@@ -280 +280 @@
                  (p(i,k-1)-p(i,k))/ph(i, k)        !convect3
 
-! c        print *,' gz(',k,')',gz(i,k),' tvx',tvx,' tvy ',tvy
+! c        PRINT *,' gz(',k,')',gz(i,k),' tvx',tvx,' tvy ',tvy
 
 ! ori         gz(i,k)=gz(i,k-1)+hrd*(tv(i,k-1)+tv(i,k))
@@ -611 +611 @@
   ENDIF
   IF (prt_level >= 10) THEN
-    print *,'cv3_feed : iflag(1), pfeed(1), plcl(1), wghti(1,k) ', &
+    PRINT *,'cv3_feed : iflag(1), pfeed(1), plcl(1), wghti(1,k) ', &
                         iflag(1), pfeed(1), plcl(1), (wghti(1,k),k=1,10)
   ENDIF
@@ -1582 +1582 @@
 #endif
   IF (prt_level >= 10) THEN
-    print *,'cv3_undilute2.0. icvflag_Tpa, t(1,k), q(1,k), qs(1,k) ', &
+    PRINT *,'cv3_undilute2.0. icvflag_Tpa, t(1,k), q(1,k), qs(1,k) ', &
                         icvflag_Tpa, (k, t(1,k), q(1,k), qs(1,k), k = 1,nl)
   ENDIF
@@ -1775 +1775 @@
             ENDIF ! (tg .gt. Tx)
 
-!!     print *,' j, k, Um, U, Ux, aa, bb, discr, dd, ddelta ', j, k, Um, U, Ux, aa, bb, discr, dd, ddelta
-!!     print *,' j, k, ah0(i), ahg, tg, qg, tp(i,k), ff ', j, k, ah0(i), ahg, tg, qg, tp(i,k), ff
+!!     PRINT *,' j, k, Um, U, Ux, aa, bb, discr, dd, ddelta ', j, k, Um, U, Ux, aa, bb, discr, dd, ddelta
+!!     PRINT *,' j, k, ah0(i), ahg, tg, qg, tp(i,k), ff ', j, k, ah0(i), ahg, tg, qg, tp(i,k), ff
           END IF ! (k>=(icbs(i)+1))
         END DO ! i = 1, ncum
@@ -1839 +1839 @@
             Tp(i,k) = tg + (ah0(i) - ahg)/ &
                     (cpd + (cl-cpd)*qta(i,k-1) + (lv(i,k)+frac(i,k)*lf(i,k))*dqgdT)
-!!   print *,'undilute2 iterations k, Tp(i,k), ah0(i), ahg ', &
+!!   PRINT *,'undilute2 iterations k, Tp(i,k), ah0(i), ahg ', &
 !!                                 k, Tp(i,k), ah0(i), ahg
           END IF ! (k>=(icbs(i)+1))
@@ -2216 +2216 @@
 !   =========================================================================
   IF (prt_level >= 10) THEN
-    print *,'cv3_undilute2.1. tp(1,k), tvp(1,k) ', &
+    PRINT *,'cv3_undilute2.1. tp(1,k), tvp(1,k) ', &
                           (k, tp(1,k), tvp(1,k), k = 1,nl)
   ENDIF
@@ -3696 +3696 @@
 
 ! ------------------------------------------------------
-IF (prt_level >= 10) print *,' ->cv3_unsat, iflag(1) ', iflag(1)
+IF (prt_level >= 10) PRINT *,' ->cv3_unsat, iflag(1) ', iflag(1)
 
 smallestreal=tiny(smallestreal)
@@ -4114 +4114 @@
 !    water(il,i)=revap*revap
 !    prec(il,i)=revap*revap
-!!        print *,' i,b6,c6,revap,evap(il,i),water(il,i),wdtrain(il) ', &
+!!        PRINT *,' i,b6,c6,revap,evap(il,i),water(il,i),wdtrain(il) ', &
 !!                 i,b6,c6,revap,evap(il,i),water(il,i),wdtrain(il)
 !!---end jyg---
@@ -4139 +4139 @@
           b6 = bfac*50.*sigd(il)*(ph(il,i)-ph(il,i+1))*sigt*afac1
           c6 = prec(il, i+1) + 0.5*bfac*wdtrain(il)
-! print *,'bfac,sigd(il),sigt,afac1 ',bfac,sigd(il),sigt,afac1
-! print *,'prec(il,i+1),wdtrain(il) ',prec(il,i+1),wdtrain(il)
-! print *,'b6,c6,b6*b6+4.*c6 ',b6,c6,b6*b6+4.*c6
+! PRINT *,'bfac,sigd(il),sigt,afac1 ',bfac,sigd(il),sigt,afac1
+! PRINT *,'prec(il,i+1),wdtrain(il) ',prec(il,i+1),wdtrain(il)
+! PRINT *,'b6,c6,b6*b6+4.*c6 ',b6,c6,b6*b6+4.*c6
           IF (c6>b6*b6+1.E-20) THEN
             revap = 2.*c6/(b6+sqrt(b6*b6+4.*c6))
@@ -4163 +4163 @@
 !JYG---   Dans tous les cas, evaporation = [tt ce qui entre dans la couche i]
 ! moins [tt ce qui sort de la couche i]
-! print *, 'evap avec ice'
+! PRINT *, 'evap avec ice'
           evap(il, i) = (wdtrain(il)+sigd(il)*wt(il,i)*(prec(il,i+1)-prec(il,i))) / &
                         (sigd(il)*(ph(il,i)-ph(il,i+1))*100.)
@@ -4227 +4227 @@
             water(il, i) = 0.
           END IF
-! print *, 'evap sans ice'
+! PRINT *, 'evap sans ice'
           evap(il, i) = (wdtrain(il)+sigd(il)*wt(il,i)*(water(il,i+1)-water(il,i)))/ &
                         (sigd(il)*(ph(il,i)-ph(il,i+1))*100.)
@@ -5645 +5645 @@
 !-----------------------------------------------------------
 
-!!   print *,'yield, i, amp1, ad', i, amp1(1), ad(1)
+!!   PRINT *,'yield, i, amp1, ad', i, amp1(1), ad(1)
 
     DO il = 1, ncum
@@ -6798 +6798 @@
 
     IF (prt_level >= 5) THEN
-      print *,' CV3_YIELD : alpha_qpos ',alpha_qpos(1)
+      PRINT *,' CV3_YIELD : alpha_qpos ',alpha_qpos(1)
     ENDIF
 
@@ -7022 +7022 @@
           END IF
         END IF
-! c        print *,'cbmf',il,i,k,cbmf(il),wghti(il,k)
+! c        PRINT *,'cbmf',il,i,k,cbmf(il),wghti(il,k)
       END DO
     END DO
@@ -7035 +7035 @@
           dnwd(il, i) = dnwd(il, i) + dn1(il, k, i)
         END IF
-! c         print *,'upwd',il,i,k,inb(il),upwd(il,i),m(il,k),up1(il,k,i)
+! c         PRINT *,'upwd',il,i,k,inb(il),upwd(il,i),m(il,k),up1(il,k,i)
       END DO
     END DO