Changeset 409

Timestamp:

Oct 15, 2002, 5:00:59 PM (22 years ago)

Author:

lmdzadmin

Message:

correction bug JLD
IM/LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/diagphy.F (modified) (5 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/diagphy.F

@@ -128 +128 @@
      $    tit, pas, fs_bound, d_etp_tot, fq_bound, d_qt_tot
 C
+      IF (iprt.ge.1) write(6,6668) 
+     $    tit, pas, d_etp_tot+d_ec_tot-fs_bound, d_qt_tot-fq_bound
+C
       IF (iprt.ge.2) write(6,6667) 
      $    tit, pas, stops,stopl,ssols,ssoll,ssens,slat,evap_tot
@@ -134 +137 @@
       return
 
- 6666 format('Flux Budget ',a15,1i6,2f8.2,2(1pE13.5))
- 6667 format('Boundary Flux ',a15,1i6,6f8.2,2(1pE13.5))
+ 6666 format('Phys. Flux Budget ',a15,1i6,2f8.2,2(1pE13.5))
+ 6667 format('Phys. Boundary Flux ',a15,1i6,6f8.2,2(1pE13.5))
+ 6668 format('Phys. Total Budget ',a15,1i6,f8.2,2(1pE13.5))
 
       end
@@ -194 +198 @@
 C
 C     other (COMMON...)
-C     RG : gravity
+C     RCPD, RCPV, ....
 C
 C J.L. Dufresne, July 2002
@@ -293 +297 @@
 C         Cinetic Energy
           zec_col(i) =  zec_col(i)
-     $        +0.5*SQRT(u(i,k)**2+v(i,k)**2)*zairm(i,k)
+     $        +0.5*(u(i,k)**2+v(i,k)**2)*zairm(i,k)
 C         Air enthalpy
           zh_dair_col(i) = zh_dair_col(i) 
@@ -377 +381 @@
       IF (iprt.ge.2) THEN
         WRITE(6,9000) tit,pas(idiag),d_qt,d_qw,d_ql,d_qs
- 9000   format('Watter Mass Budget (kg/m2/s)',A15
+ 9000   format('Phys. Watter Mass Budget (kg/m2/s)',A15
      $      ,1i6,10(1pE14.6))
         WRITE(6,9001) tit,pas(idiag), d_h_vcol
- 9001   format('Enthalpy Budget (W/m2) ',A15,1i6,10(F8.2))
+ 9001   format('Phys. Enthalpy Budget (W/m2) ',A15,1i6,10(F8.2))
         WRITE(6,9002) tit,pas(idiag), d_ec
- 9002   format('Cinetic Energy Budget (W/m2) ',A15,1i6,10(F8.2))
+ 9002   format('Phys. Cinetic Energy Budget (W/m2) ',A15,1i6,10(F8.2))
       END IF 
 C