Changeset 855 for LMDZ4/trunk/libf/phytherm

Timestamp:

Oct 23, 2007, 3:03:13 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Bidouillage interne :-)
LF

Location:

LMDZ4/trunk/libf/phytherm

Files:

• thermcell_flux.F90 (modified) (14 diffs)
• thermcell_flux1.F90 (deleted)
• thermcell_main.F90 (modified) (1 diff)

LMDZ4/trunk/libf/phytherm/thermcell_flux.F90

@@ +1 @@
+!
+! $Header$
+!
+
+
       SUBROUTINE thermcell_flux(ngrid,klev,ptimestep,masse, &
      &       lalim,lmax,alim_star,  &
@@ -152 +157 @@
     &    ptimestep,masse,entr,detr,fm,'2  ')
 
+
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! FH Version en cours de test;
+! par rapport a thermcell_flux, on fait une grande boucle sur "l"
+! et on modifie le flux avec tous les contrôles appliques d'affilee
+! pour la meme couche
+! Momentanement, on duplique le calcule du flux pour pouvoir comparer
+! les flux avant et apres modif
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+      do l=1,klev
+
+         do ig=1,ngrid
+            if (l.lt.lmax(ig)) then
+               fm(ig,l+1)=fm(ig,l)+entr(ig,l)-detr(ig,l)
+            elseif(l.eq.lmax(ig)) then
+               fm(ig,l+1)=0.
+               detr(ig,l)=fm(ig,l)+entr(ig,l)
+            else
+               fm(ig,l+1)=0.
+            endif
+         enddo
+
+
 !-------------------------------------------------------------------------
 ! Verification de la positivite des flux de masse
 !-------------------------------------------------------------------------
 
-      do l=1,klev
+!     do l=1,klev
          do ig=1,ngrid
             if (fm(ig,l+1).lt.0.) then
@@ -165 +195 @@
             endif
          enddo
-      enddo
-
-      if (lev_out.ge.10) &
-    &    call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
-    &    ptimestep,masse,entr,detr,fm,'3  ')
+!     enddo
+
+      if (lev_out.ge.10) &
+     &   write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
+     &     entr(igout,l),detr(igout,l),fm(igout,l+1)
 
 !-------------------------------------------------------------------------
@@ -175 +205 @@
 !-------------------------------------------------------------------------
 
-      do l=1,klev
+
+      if (1.eq.0) then
+!     do l=1,klev
          do ig=1,ngrid
           if (l.ge.lalim(ig).and.l.le.lmax(ig) &
@@ -189 +221 @@
           endif
         enddo
-      enddo
-
-      if (lev_out.ge.10) &
-    &    call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
-    &    ptimestep,masse,entr,detr,fm,'4  ')
+!     enddo
+
+      if (lev_out.ge.10) &
+     &   write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
+     &     entr(igout,l),detr(igout,l),fm(igout,l+1)
+      else
+         print*,'Test sur les fractions croissantes inhibe dans thermcell_flux2'
+      endif
 
 
@@ -200 +235 @@
 !-------------------------------------------------------------------------
 
-      do l=1,klev
+!     do l=1,klev
          do ig=1,ngrid
             if ((fm(ig,l+1).gt.fm(ig,l)).and.(l.gt.lalim(ig))) then
@@ -209 +244 @@
             endif
          enddo
-      enddo
-
-      if (lev_out.ge.10) &
-    &    call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
-    &    ptimestep,masse,entr,detr,fm,'5  ')
+!     enddo
+
+      if (lev_out.ge.10) &
+     &   write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
+     &     entr(igout,l),detr(igout,l),fm(igout,l+1)
 
 !-------------------------------------------------------------------------
@@ -221 +256 @@
       if(1.eq.1) then
 
-      do l=1,klev
+!     do l=1,klev
          do ig=1,ngrid
             if (entr(ig,l)<0.) then
@@ -256 +291 @@
             endif
          enddo
-      enddo
-      endif
-
-
-      if (lev_out.ge.10) &
-    &    call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
-    &    ptimestep,masse,entr,detr,fm,'6  ')
+!     enddo
+      endif
+
+
+      if (lev_out.ge.10) &
+     &   write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
+     &     entr(igout,l),detr(igout,l),fm(igout,l+1)
 
 !-------------------------------------------------------------------------
@@ -268 +303 @@
 !-------------------------------------------------------------------------
 
-      do l=1,klev
+!     do l=1,klev
          do ig=1,ngrid
             if (fm(ig,l+1).lt.0.) then
@@ -283 +318 @@
             endif
         enddo
-     enddo
-
-      if (lev_out.ge.10) &
-    &    call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
-    &    ptimestep,masse,entr,detr,fm,'7  ')
+!    enddo
+
+      if (lev_out.ge.10) &
+     &   write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
+     &     entr(igout,l),detr(igout,l),fm(igout,l+1)
 
 !-----------------------------------------------------------------------
@@ -308 +343 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-     if(1.eq.0) then
-
-     do l=1,klev
+!    do l=1,klev
         do ig=1,ngrid
            if (zw2(ig,l+1).gt.1.e-10) then
@@ -327 +360 @@
            endif
         enddo
-     enddo
+!    enddo
 !
+
+
+      if (lev_out.ge.10) &
+     &   write(lunout,'(i4,4e14.4)') l,masse(igout,l)/ptimestep, &
+     &     entr(igout,l),detr(igout,l),fm(igout,l+1)
+
+! Fin de la grande boucle sur les niveaux verticaux
+      enddo
+
       if (lev_out.ge.10) &
     &    call printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
     &    ptimestep,masse,entr,detr,fm,'8  ')
 
-     endif
 
 !-----------------------------------------------------------------------
@@ -418 +459 @@
       return
       end
-
-      subroutine printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
-    &    ptimestep,masse,entr,detr,fm,descr)
-
-     implicit none
-
-      integer ngrid,klev,lunout,igout,l,lm
-
-      integer lmax(klev),lalim(klev)
-      real ptimestep,masse(ngrid,klev),entr(ngrid,klev),detr(ngrid,klev)
-      real fm(ngrid,klev+1),f(ngrid)
-
-      character*3 descr
-
-      lm=lmax(igout)+5
-      if(lm.gt.klev) lm=klev
-
-      print*,'Impression jusque lm=',lm
-
-         write(lunout,*) 'Dans thermcell_flux '//descr
-         write(lunout,*) 'flux base ',f(igout)
-         write(lunout,*) 'lmax ',lmax(igout)
-         write(lunout,*) 'lalim ',lalim(igout)
-         write(lunout,*) 'ig= ',igout
-         write(lunout,'(a3,4a14)') 'l','M','E','D','F'
-         write(lunout,'(i4,4e14.4)') (l,masse(igout,l)/ptimestep, &
-     &     entr(igout,l),detr(igout,l) &
-     &    ,fm(igout,l+1),l=1,lm)
-
-
-      do l=lmax(igout)+1,klev
-          if (abs(entr(igout,l))+abs(detr(igout,l))+abs(fm(igout,l)).gt.0.) then
-          print*,'cas 1 : igout,l,lmax(igout)',igout,l,lmax(igout)
-          print*,'entr(igout,l)',entr(igout,l)
-          print*,'detr(igout,l)',detr(igout,l)
-          print*,'fm(igout,l)',fm(igout,l)
-          stop
-          endif
-      enddo
-
-      return
-      end
-

LMDZ4/trunk/libf/phytherm/thermcell_main.F90

@@ -385 +385 @@
 !-------------------------------------------------------------------------------
 
-      CALL thermcell_flux2(ngrid,nlay,ptimestep,masse, &
+      CALL thermcell_flux(ngrid,nlay,ptimestep,masse, &
      &       lalim,lmax,alim_star,  &
      &       entr_star,detr_star,f,rhobarz,zlev,zw2,fm,entr,  &