/[lmdze]/trunk/libf/phylmd/Thermcell/thermcell.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/Thermcell/thermcell.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 10 by guez, Fri Apr 18 14:45:53 2008 UTC
# Line 45  c   ---------- Line 45  c   ----------
45        REAL pu(ngrid,nlay),pduadj(ngrid,nlay)        REAL pu(ngrid,nlay),pduadj(ngrid,nlay)
46        REAL pv(ngrid,nlay),pdvadj(ngrid,nlay)        REAL pv(ngrid,nlay),pdvadj(ngrid,nlay)
47        REAL po(ngrid,nlay),pdoadj(ngrid,nlay)        REAL po(ngrid,nlay),pdoadj(ngrid,nlay)
48        REAL pplay(ngrid,nlay)        REAL, intent(in):: pplay(ngrid,nlay)
49        real, intent(in):: pplev(ngrid,nlay+1)        real, intent(in):: pplev(ngrid,nlay+1)
50        real pphi(ngrid,nlay)        real pphi(ngrid,nlay)
51    
# Line 722  c    s         ,'   FM=',fm(ig,l) Line 722  c    s         ,'   FM=',fm(ig,l)
722              endif              endif
723              if(.not.masse(ig,l).ge.1.e-10              if(.not.masse(ig,l).ge.1.e-10
724       s         .or..not.masse(ig,l).le.1.e4) then       s         .or..not.masse(ig,l).le.1.e4) then
 c     print*,'WARN!!! masse exagere ig=',ig,'   l=',l  
 c    s         ,'   M=',masse(ig,l)  
 c     print*,'rho(ig,l),pplay(ig,l),zpspsk(ig,l),RD,zh(ig,l)',  
 c    s                 rho(ig,l),pplay(ig,l),zpspsk(ig,l),RD,zh(ig,l)  
 c     print*,'zlev(ig,l+1),zlev(ig,l)'  
 c    s                ,zlev(ig,l+1),zlev(ig,l)  
 c     print*,'pphi(ig,l-1),pphi(ig,l),pphi(ig,l+1)'  
 c    s                ,pphi(ig,l-1),pphi(ig,l),pphi(ig,l+1)  
725              endif              endif
726              if(.not.entr(ig,l).ge.0..or..not.entr(ig,l).le.10.) then              if(.not.entr(ig,l).ge.0..or..not.entr(ig,l).le.10.) then
727  c     print*,'WARN!!! entr exagere ig=',ig,'   l=',l  c     print*,'WARN!!! entr exagere ig=',ig,'   l=',l
Legend:
Removed from v.3  
changed lines
  Added in v.10
  ViewVC Help
Powered by ViewVC 1.1.21