Changeset 972 for LMDZ4/trunk/libf/phylmd/thermcell_closure.F90

Timestamp:

Jun 19, 2008, 12:24:22 PM (16 years ago)

Author:

lmdzadmin

Message:

Version thermique FH/CRio
Ajout tests cas physiques non pris en comptes et ajout/enleve prints
Nouvelle routine thermcell_flux2.F90
IM

File:

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

LMDZ4/trunk/libf/phylmd/thermcell_closure.F90

@@ -1 +1 @@
       SUBROUTINE thermcell_closure(ngrid,nlay,r_aspect,ptimestep,rho,  &
-     &   zlev,lalim,alim_star,zmax_sec,wmax_sec,zmax,wmax,f,f0,lev_out)
+     &   zlev,lalim,alim_star,zmax_sec,wmax_sec,zmax,wmax,f,lev_out)
 
 !-------------------------------------------------------------------------
@@ -24 +24 @@
 
       REAL f(ngrid)
-      REAL f0(ngrid)
 
       do ig=1,ngrid
@@ -52 +51 @@
              f(ig)=wmax_sec(ig)/(max(500.,zmax_sec(ig))*r_aspect  &
      &             *alim_star2(ig))
-             f(ig)=f(ig)+(f0(ig)-f(ig))*exp((-ptimestep/  &
-     &                     zmax_sec(ig))*wmax_sec(ig))
+!            f(ig)=f(ig)+(f0(ig)-f(ig))*exp((-ptimestep/  &
+!    &                     zmax_sec(ig))*wmax_sec(ig))
              else
              f(ig)=wmax(ig)/zdenom
-            f(ig)=f(ig)+(f0(ig)-f(ig))*exp((-ptimestep/  &
-     &                     zmax(ig))*wmax(ig))
+!            f(ig)=f(ig)+(f0(ig)-f(ig))*exp((-ptimestep/  &
+!     &                     zmax(ig))*wmax(ig))
              endif
          endif
-         f0(ig)=f(ig)
+!         f0(ig)=f(ig)
       enddo
       if (prt_level.ge.1) print*,'apres fermeture'