Ignore:
Timestamp:
Oct 23, 2007, 3:31:10 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Correction un peu trop rapide
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phytherm/thermcell_flux.F90

    r855 r856  
    459459      return
    460460      end
     461
     462      subroutine printflux(ngrid,klev,lunout,igout,f,lmax,lalim, &
     463    &    ptimestep,masse,entr,detr,fm,descr)
     464
     465     implicit none
     466
     467      integer ngrid,klev,lunout,igout,l,lm
     468
     469      integer lmax(klev),lalim(klev)
     470      real ptimestep,masse(ngrid,klev),entr(ngrid,klev),detr(ngrid,klev)
     471      real fm(ngrid,klev+1),f(ngrid)
     472
     473      character*3 descr
     474
     475      lm=lmax(igout)+5
     476      if(lm.gt.klev) lm=klev
     477
     478      print*,'Impression jusque lm=',lm
     479
     480         write(lunout,*) 'Dans thermcell_flux '//descr
     481         write(lunout,*) 'flux base ',f(igout)
     482         write(lunout,*) 'lmax ',lmax(igout)
     483         write(lunout,*) 'lalim ',lalim(igout)
     484         write(lunout,*) 'ig= ',igout
     485         write(lunout,'(a3,4a14)') 'l','M','E','D','F'
     486         write(lunout,'(i4,4e14.4)') (l,masse(igout,l)/ptimestep, &
     487     &     entr(igout,l),detr(igout,l) &
     488     &    ,fm(igout,l+1),l=1,lm)
     489
     490
     491      do l=lmax(igout)+1,klev
     492          if (abs(entr(igout,l))+abs(detr(igout,l))+abs(fm(igout,l)).gt.0.) then
     493          print*,'cas 1 : igout,l,lmax(igout)',igout,l,lmax(igout)
     494          print*,'entr(igout,l)',entr(igout,l)
     495          print*,'detr(igout,l)',detr(igout,l)
     496          print*,'fm(igout,l)',fm(igout,l)
     497          stop
     498          endif
     499      enddo
     500
     501      return
     502      end
     503
Note: See TracChangeset for help on using the changeset viewer.