Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume.F90

Timestamp:

Jul 20, 2024, 12:00:23 PM (2 months ago)

Author:

abarral

Message:

remove fixed-form ^{\s+& remaining in .f90,.F90}

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_plume.F90 (modified) (8 diffs)

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

@@ -6 +6 @@
 
       SUBROUTINE thermcell_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,  &
-     &           zlev,pplev,pphi,zpspsk,alim_star,alim_star_tot,  &
-     &           lalim,f0,detr_star,entr_star,f_star,csc,ztva,  &
-     &           ztla,zqla,zqta,zha,zw2,w_est,ztva_est,zqsatth,lmix,lmix_bis,linter &
-    &           ,lev_out,lunout1,igout)
+             zlev,pplev,pphi,zpspsk,alim_star,alim_star_tot,  &
+             lalim,f0,detr_star,entr_star,f_star,csc,ztva,  &
+             ztla,zqla,zqta,zha,zw2,w_est,ztva_est,zqsatth,lmix,lmix_bis,linter &
+             ,lev_out,lunout1,igout)
 !     &           ,lev_out,lunout1,igout,zbuoy,zbuoyjam)
 !--------------------------------------------------------------------------
@@ -227 +227 @@
         ztva_est(ig,l) = ztva_est(ig,l)/zpspsk(ig,l)
         ztva_est(ig,l) = ztva_est(ig,l)*(1.+RETV*(zqta(ig,l-1)  &
-     &      -zqla_est(ig,l))-zqla_est(ig,l))
+        -zqla_est(ig,l))-zqla_est(ig,l))
  
 
@@ -265 +265 @@
         coefzlmel=Min(1.,(zlmelup-zltdwn)/zdz)
         zbuoyjam(ig,l)=1.*RG*(coefzlmel*(ztva_est(ig,l)- & 
-    &   ztv(ig,lt))/ztv(ig,lt)+(1.-coefzlmel)*(ztva_est(ig,l)- &
-    &   ztv(ig,lt-1))/ztv(ig,lt-1))+0.*zbuoy(ig,l)
+     ztv(ig,lt))/ztv(ig,lt)+(1.-coefzlmel)*(ztva_est(ig,l)- &
+     ztv(ig,lt-1))/ztv(ig,lt-1))+0.*zbuoy(ig,l)
 
 !------------------------------------------------
@@ -303 +303 @@
 
           detr_star(ig,l)=f_star(ig,l)*zdz             &
-    &     *( mix0 * 0.1 / (zalpha+0.001)               &
-    &     + MAX(detr_min, -afact*zbetalpha*zbuoyjam(ig,l)/zw2m   &
-    &     + detr_q_coef*(zdqt(ig,l)/zw2m)**detr_q_power))
+       *( mix0 * 0.1 / (zalpha+0.001)               &
+       + MAX(detr_min, -afact*zbetalpha*zbuoyjam(ig,l)/zw2m   &
+       + detr_q_coef*(zdqt(ig,l)/zw2m)**detr_q_power))
 
           if ( iflag_thermals_ed == 20 ) then
              entr_star(ig,l)=f_star(ig,l)*zdz* (         &
-    &          mix0 * 0.1 / (zalpha+0.001)               &
-    &        + zbetalpha*MAX(entr_min,                   &
-    &        afact*zbuoyjam(ig,l)/zw2m - fact_epsilon))
+            mix0 * 0.1 / (zalpha+0.001)               &
+          + zbetalpha*MAX(entr_min,                   &
+          afact*zbuoyjam(ig,l)/zw2m - fact_epsilon))
           else
              entr_star(ig,l)=f_star(ig,l)*zdz* (         &
-    &          mix0 * 0.1 / (zalpha+0.001)               &
-    &        + zbetalpha*MAX(entr_min,                   &
-    &        afact*zbuoy(ig,l)/zw2m - fact_epsilon))
+            mix0 * 0.1 / (zalpha+0.001)               &
+          + zbetalpha*MAX(entr_min,                   &
+          afact*zbuoy(ig,l)/zw2m - fact_epsilon))
           endif
           
@@ -326 +326 @@
         endif
         f_star(ig,l+1)=f_star(ig,l)+alim_star(ig,l)+entr_star(ig,l)  &
-     &              -detr_star(ig,l)
+                -detr_star(ig,l)
 
       endif
@@ -341 +341 @@
            Zsat=.false.
            ztla(ig,l)=(f_star(ig,l)*ztla(ig,l-1)+  &
-     &            (alim_star(ig,l)+entr_star(ig,l))*zthl(ig,l))  &
-     &            /(f_star(ig,l+1)+detr_star(ig,l))
+              (alim_star(ig,l)+entr_star(ig,l))*zthl(ig,l))  &
+              /(f_star(ig,l+1)+detr_star(ig,l))
            zqta(ig,l)=(f_star(ig,l)*zqta(ig,l-1)+  &
-     &            (alim_star(ig,l)+entr_star(ig,l))*po(ig,l))  &
-     &            /(f_star(ig,l+1)+detr_star(ig,l))
+              (alim_star(ig,l)+entr_star(ig,l))*po(ig,l))  &
+              /(f_star(ig,l+1)+detr_star(ig,l))
 
         endif
@@ -362 +362 @@
            zha(ig,l) = ztva(ig,l)
            ztva(ig,l) = ztva(ig,l)*(1.+RETV*(zqta(ig,l)  &
-     &              -zqla(ig,l))-zqla(ig,l))
+                -zqla(ig,l))-zqla(ig,l))
            zbuoy(ig,l)=RG*(ztva(ig,l)-ztv(ig,l))/ztv(ig,l)
            zdz=zlev(ig,l+1)-zlev(ig,l)
@@ -393 +393 @@
         if (zw2(ig,l+1)<0.) then
            linter(ig)=(l*(zw2(ig,l+1)-zw2(ig,l))  &
-     &               -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l))
+                 -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l))
            zw2(ig,l+1)=0.
 !+CR:04/05/12:correction calcul linter pour calcul de zmax continu
         elseif (f_star(ig,l+1)<0.) then
            linter(ig)=(l*(f_star(ig,l+1)-f_star(ig,l))  &
-     &               -f_star(ig,l))/(f_star(ig,l+1)-f_star(ig,l))
+                 -f_star(ig,l))/(f_star(ig,l+1)-f_star(ig,l))
            zw2(ig,l+1)=0.
 !fin CR:04/05/12