Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_alp.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (8 weeks 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/phylmd/lmdz_thermcell_alp.F90 (modified) (14 diffs)

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

@@ -103 +103 @@
       real, parameter :: coef_m=1.            ! On considère un rendement pour alp_bl_fluct_m
       real, parameter :: coef_tke=1.          ! On considère un rendement pour alp_bl_fluct_tke
-      real :: zdp
+      REAL :: zdp
       real, dimension(ngrid) :: alp_int,dp_int
       real, dimension(ngrid) :: fm_tot
@@ -143 +143 @@
     do l=1,nlay-1
       do ig=1,ngrid
-        if (ok_lcl(ig)) then 
+        if (ok_lcl(ig)) THEN
 !ATTENTION,zw2 calcule en pplev
-!          if ((pplay(ig,l) .ge. pcon(ig)) .and. (pplay(ig,l+1) .le. pcon(ig))) then
+!          if ((pplay(ig,l) .ge. pcon(ig)) .and. (pplay(ig,l+1) .le. pcon(ig))) THEN
 !          klcl(ig)=l
 !          interp(ig)=(pcon(ig)-pplay(ig,klcl(ig)))/(pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig)))
 !          endif
-          if ((pplev(ig,l) >= pcon(ig)) .and. (pplev(ig,l+1) <= pcon(ig))) then
+          if ((pplev(ig,l) >= pcon(ig)) .and. (pplev(ig,l+1) <= pcon(ig))) THEN
           klcl(ig)=l
           interp(ig)=(pcon(ig)-pplev(ig,klcl(ig)))/(pplev(ig,klcl(ig)+1)-pplev(ig,klcl(ig)))
@@ -159 +159 @@
     do ig =1,ngrid
 !CR:REHABILITATION ZMAX CONTINU
-     if (ok_lcl(ig)) then 
+     if (ok_lcl(ig)) THEN
       rhobarz0(ig)=rhobarz(ig,klcl(ig))+(rhobarz(ig,klcl(ig)+1) &
                  -rhobarz(ig,klcl(ig)))*interp(ig)
@@ -214 +214 @@
 
   do ig=1,ngrid
-   if (ok_lcl(ig)) then
+   if (ok_lcl(ig)) THEN
      fraca0(ig)=fraca(ig,klcl(ig))+(fraca(ig,klcl(ig)+1) &
                -fraca(ig,klcl(ig)))*interp(ig)
@@ -258 +258 @@
      depth(ig)=zmax_moy(ig)-zlcl(ig)
      hmin(ig)=hmincoef*zlcl(ig)
-     if (depth(ig)>=10.) then
+     if (depth(ig)>=10.) THEN
        s2(ig)=(hcoef*depth(ig)+hmin(ig))**2
        n2(ig)=(1.-eps1)*fraca0(ig)*airephy(ig)/s2(ig)
@@ -275 +275 @@
 
 !CR: calcul de strig
-   if (iflag_strig==0) then
+   if (iflag_strig==0) THEN
       strig(:)=s_trig
-   else if (iflag_strig==1) then
+   else if (iflag_strig==1) THEN
       do ig=1,ngrid
 !         zcong_moy(ig)=zlcl(ig)+zmax_moy_coef*(zcong(ig)-zlcl(ig))
@@ -283 +283 @@
          strig(ig)=(zcong(ig)-zlcl(ig))**2  
       enddo
-   else if (iflag_strig==2) then
+   else if (iflag_strig==2) THEN
       do ig=1,ngrid 
-         if (h_trig>zlcl(ig)) then
+         if (h_trig>zlcl(ig)) THEN
          strig(ig)=(h_trig-zlcl(ig))**2
          else
@@ -296 +296 @@
    reuler=exp(1.)
    do ig=1,ngrid
-     if ( (depth(ig)>=10.) .and. (s_max(ig)>susqr2pi*reuler) ) then
+     if ( (depth(ig)>=10.) .and. (s_max(ig)>susqr2pi*reuler) ) THEN
       w_max(ig)=w0(ig)*(1.+sqrt(2.*log(s_max(ig)/susqr2pi)-log(2.*log(s_max(ig)/susqr2pi))))
       ale_bl_stat(ig)=0.5*w_max(ig)**2
@@ -319 +319 @@
   alp_bl_fluct_tke(ig)=3.*coef_m*rhobarz0(ig)*w0(ig)*fraca0(ig)*(therm_tke_max0(ig)-env_tke_max0(ig)) &
                       +3.*rhobarz0(ig)*w_conv(ig)*pbl_tke_max0(ig)
-    if (iflag_clos_bl>=2) then
+    if (iflag_clos_bl>=2) THEN
     alp_bl_conv(ig)=1.5*coef_m*rhobarz0(ig)*fraca0(ig)*(fraca0(ig)/(1.-fraca0(ig)))*w_conv(ig)* &
                      (w0(ig)**2)
@@ -371 +371 @@
       do k=1,nlay
          do ig=1,ngrid
-            if (k<=lalim_conv(ig).and.alim_star(ig,k)>1.e-10) then
+            if (k<=lalim_conv(ig).and.alim_star(ig,k)>1.e-10) THEN
                wght_th(ig,k)=alim_star(ig,k)
             endif
@@ -380 +380 @@
 !test pour prolonger la convection
       do ig=1,ngrid
-!v1d  if ((alim_star(ig,1).lt.1.e-10).and.(therm)) then
-      if ((alim_star(ig,1)<1.e-10)) then
+!v1d  if ((alim_star(ig,1).lt.1.e-10).and.(therm)) THEN
+      if ((alim_star(ig,1)<1.e-10)) THEN
       lalim_conv(ig)=1
       wght_th(ig,1)=1.
@@ -399 +399 @@
       do l=2,nlay
         do ig=1,ngrid
-           if(l<=lmax(ig)) THEN
+           IF(l<=lmax(ig)) THEN
            zdp=pplay(ig,l-1)-pplay(ig,l)
            alp_int(ig)=alp_int(ig)+0.5*rhobarz(ig,l)*wth3(ig,l)*zdp
@@ -407 +407 @@
       enddo
 
-      if (iflag_coupl>=3 .and. iflag_coupl<=5) then
+      if (iflag_coupl>=3 .and. iflag_coupl<=5) THEN
       do ig=1,ngrid
 !valeur integree de alp_bl * 0.5:
-        if (dp_int(ig)>0.) then
+        if (dp_int(ig)>0.) THEN
         alp_bl(ig)=alp_int(ig)/dp_int(ig)
         endif
@@ -424 +424 @@
 
 
-      return
+      RETURN
       end
 END MODULE lmdz_thermcell_alp