Changeset 110 for trunk/libf/dyn3d/top_bound.F

Timestamp:

Apr 15, 2011, 3:44:44 PM (14 years ago)

Author:

slebonnois

Message:

SL: corrections de bugs suite a compilations venus et titan de la version 109.

File:

• trunk/libf/dyn3d/top_bound.F (modified) (3 diffs)

trunk/libf/dyn3d/top_bound.F

@@ -1 @@
-      SUBROUTINE top_bound( vcov,ucov,teta,masse, du,dv,dh )
+      SUBROUTINE top_bound( vcov,ucov,teta,phi,masse, du,dv,dh )
       IMPLICIT NONE
 c
@@ -35 +35 @@
 
       REAL ucov(iip1,jjp1,llm),vcov(iip1,jjm,llm),teta(iip1,jjp1,llm)
+      REAL phi(iip1,jjp1,llm)                  ! geopotentiel
       REAL masse(iip1,jjp1,llm)
       REAL dv(iip1,jjm,llm),du(iip1,jjp1,llm),dh(iip1,jjp1,llm)
@@ -77 +78 @@
          print*,'TOP_BOUND mode',mode_top_bound
          print*,'Coeffs pour la couche eponge a l equateur'
-         print*,'p (Pa)  z(km)  tau (s)   dt*rdamp'
+         print*,'p (Pa)  z(km)  tau (s)'
          do l=1,llm
            if (rdamp(l).ne.0.) then
             zkm        = phi(iip1/2,jjp1/2,l)/(1000*g)
-          print*,presnivs(l),zkm,
-     .          1./rdamp(l),
-     .          dt*rdamp(l)
+          print*,presnivs(l),zkm,1./rdamp(l)
            endif
          enddo