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 |
|
|
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 |