Changeset 852 for LMDZ4/trunk/libf/phytherm/thermcell_closure.F90

Timestamp:

Oct 11, 2007, 3:43:42 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour de la physique avec thermiques avec la version de FH d'aout 2007
LF

File:

• LMDZ4/trunk/libf/phytherm/thermcell_closure.F90 (modified) (3 diffs)

LMDZ4/trunk/libf/phytherm/thermcell_closure.F90

@@ -18 +18 @@
       REAL zmax(ngrid),zmax_sec(ngrid)
       REAL wmax(ngrid),wmax_sec(ngrid)
+      real zdenom
 
       REAL alim_star2(ngrid)
@@ -35 +36 @@
      &                    /(rho(ig,k)*(zlev(ig,k+1)-zlev(ig,k)))
              enddo
+             zdenom=max(500.,zmax(ig))*r_aspect*alim_star2(ig)
+             if (zdenom<1.e-14) then
+                print*,'ig=',ig
+                print*,'alim_star2',alim_star2(ig)
+                print*,'zmax',zmax(ig)
+                print*,'r_aspect',r_aspect
+                print*,'zdenom',zdenom
+                print*,'alim_star',alim_star(ig,:)
+                print*,'zmax_sec',zmax_sec(ig)
+                print*,'wmax_sec',wmax_sec(ig)
+                stop
+             endif
              if ((zmax_sec(ig).gt.1.e-10).and.(1.eq.1)) then 
              f(ig)=wmax_sec(ig)/(max(500.,zmax_sec(ig))*r_aspect  &
@@ -41 +54 @@
      &                     zmax_sec(ig))*wmax_sec(ig))
              else
-             f(ig)=wmax(ig)/(max(500.,zmax(ig))*r_aspect*alim_star2(ig))
+             f(ig)=wmax(ig)/zdenom
             f(ig)=f(ig)+(f0(ig)-f(ig))*exp((-ptimestep/  &
      &                     zmax(ig))*wmax(ig))