Changeset 2064 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_closure.F90

Timestamp:

Jan 10, 2019, 3:44:36 PM (6 years ago)

Author:

aboissinot

Message:

Fix a bug in thermcell_closure arguments which implies wrong fluxes normalization.

File:

• trunk/LMDZ.GENERIC/libf/phystd/thermcell_closure.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/thermcell_closure.F90

@@ -3 +3 @@
 !
       SUBROUTINE thermcell_closure(ngrid,nlay,ptimestep,rho,zlev,             &
-      &                    lalim,alim_star,f_star,zmax,wmax,f,lev_out)
-      
+      &                            lalim,alim_star,f_star,                    &
+      &                            zmax,wmax,f,lev_out)
       
 !==============================================================================
@@ -77 +77 @@
       DO ig=1,ngrid
          IF (lalim(ig)>llmax) THEN
-            llmax=lalim(ig)
+            llmax = lalim(ig)
          ENDIF
       ENDDO
@@ -88 +88 @@
          DO ig=1,ngrid
             IF (k<lalim(ig)) THEN
-               alim_star2(ig)=alim_star2(ig)+alim_star(ig,k)**2  &
-               &                    /(rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k)))
-               alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,k)
+               alim_star2(ig) = alim_star2(ig) + alim_star(ig,k)**2           &
+               &              / (rho(ig,k) * (zlev(ig,k+1) - zlev(ig,k)))
+               alim_star_tot(ig) = alim_star_tot(ig) + alim_star(ig,k)
             ENDIF
          ENDDO
@@ -97 +97 @@
       DO ig=1,ngrid
          IF (alim_star2(ig)>1.e-10) THEN
-            f(ig) = wmax(ig) * alim_star_tot(ig)/                             &
-            &     (max(500.,zmax(ig))*r_aspect_thermals*alim_star2(ig))
+            f(ig) = wmax(ig) * alim_star_tot(ig)                              &
+            &     / (max(500.,zmax(ig)) * r_aspect_thermals * alim_star2(ig))
          ENDIF
       ENDDO