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

Timestamp:

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

Author:

abarral

Message:

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

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_alp.F90 (modified) (5 diffs)

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

@@ -5 +5 @@
 
       SUBROUTINE thermcell_alp(ngrid,nlay,ptimestep  &                         ! in
-     &                  ,pplay,pplev  &                                        ! in
-     &                  ,fm0,entr0,lmax  &                                     ! in
-     &                  ,pbl_tke,pctsrf,omega,airephy &                        ! in
-     &                  ,zw2,fraca &                                           ! in
-     &                  ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &  ! in
+                    ,pplay,pplev  &                                        ! in
+                    ,fm0,entr0,lmax  &                                     ! in
+                    ,pbl_tke,pctsrf,omega,airephy &                        ! in
+                    ,zw2,fraca &                                           ! in
+                    ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &  ! in
 !
-     &                  ,zcong,ale_bl,alp_bl,lalim_conv,wght_th &                    ! out
-     &                  ,zlcl,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &   ! out
-     &                  ,n2,s2,strig,ale_bl_stat &                                   ! out
-     &                  ,therm_tke_max,env_tke_max &                           ! out
-     &                  ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &          ! out
-     &                  ,alp_bl_conv,alp_bl_stat &                             ! out
-     &)
+                    ,zcong,ale_bl,alp_bl,lalim_conv,wght_th &                    ! out
+                    ,zlcl,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &   ! out
+                    ,n2,s2,strig,ale_bl_stat &                                   ! out
+                    ,therm_tke_max,env_tke_max &                           ! out
+                    ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &          ! out
+                    ,alp_bl_conv,alp_bl_stat &                             ! out
+  )
 
       USE indice_sol_mod
@@ -161 +161 @@
      if (ok_lcl(ig)) then 
       rhobarz0(ig)=rhobarz(ig,klcl(ig))+(rhobarz(ig,klcl(ig)+1) &
- &               -rhobarz(ig,klcl(ig)))*interp(ig)
+                 -rhobarz(ig,klcl(ig)))*interp(ig)
       zlcl(ig)=(pplev(ig,1)-pcon(ig))/(rhobarz0(ig)*RG)
       zlcl(ig)=min(zlcl(ig),zmax(ig))   ! Si zlcl > zmax alors on pose zlcl = zmax
@@ -200 +200 @@
 !-----Calcul de la TKE transport�e par les thermiques : therm_tke_max
    call thermcell_tke_transport(ngrid,nlay,ptimestep,fm0,entr0,  &  ! in
-  &           rg,pplev,therm_tke_max)                               ! out
+             rg,pplev,therm_tke_max)                               ! out
 !   print *,' thermcell_tke_transport -> '   !!jyg
 
@@ -216 +216 @@
    if (ok_lcl(ig)) then
      fraca0(ig)=fraca(ig,klcl(ig))+(fraca(ig,klcl(ig)+1) &
- &             -fraca(ig,klcl(ig)))*interp(ig)
+               -fraca(ig,klcl(ig)))*interp(ig)
      w0(ig)=zw2(ig,klcl(ig))+(zw2(ig,klcl(ig)+1) &
- &         -zw2(ig,klcl(ig)))*interp(ig)
+           -zw2(ig,klcl(ig)))*interp(ig)
      w_conv(ig)=w_ls(ig,klcl(ig))+(w_ls(ig,klcl(ig)+1) &
- &             -w_ls(ig,klcl(ig)))*interp(ig)
+               -w_ls(ig,klcl(ig)))*interp(ig)
      therm_tke_max0(ig)=therm_tke_max(ig,klcl(ig)) &
- &                     +(therm_tke_max(ig,klcl(ig)+1)-therm_tke_max(ig,klcl(ig)))*interp(ig)
+                       +(therm_tke_max(ig,klcl(ig)+1)-therm_tke_max(ig,klcl(ig)))*interp(ig)
      env_tke_max0(ig)=env_tke_max(ig,klcl(ig))+(env_tke_max(ig,klcl(ig)+1) &
- &                   -env_tke_max(ig,klcl(ig)))*interp(ig)
+                     -env_tke_max(ig,klcl(ig)))*interp(ig)
      pbl_tke_max0(ig)=pbl_tke_max(ig,klcl(ig))+(pbl_tke_max(ig,klcl(ig)+1) &
- &                   -pbl_tke_max(ig,klcl(ig)))*interp(ig)
+                     -pbl_tke_max(ig,klcl(ig)))*interp(ig)
      if (therm_tke_max0(ig)>=20.) therm_tke_max0(ig)=20.
      if (env_tke_max0(ig)>=20.) env_tke_max0(ig)=20.
@@ -316 +316 @@
   alp_bl_det(ig)=0.5*coef_m*rhobarz0(ig)*(w0(ig)**3)*fraca0(ig)*(1.-2.*fraca0(ig))/((1.-fraca0(ig))**2)
   alp_bl_fluct_m(ig)=1.5*rhobarz0(ig)*fraca0(ig)*(w_conv(ig)+coef_m*w0(ig))* &
- &                   (w0(ig)**2)
+                     (w0(ig)**2)
   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)
+                      +3.*rhobarz0(ig)*w_conv(ig)*pbl_tke_max0(ig)
     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)
+                     (w0(ig)**2)
     else
     alp_bl_conv(ig)=0.