Changeset 429


Ignore:
Timestamp:
Nov 14, 2002, 1:08:23 PM (22 years ago)
Author:
lmdz
Message:

pb comgeomphy. En fait, diagphy n'est pas utilise par la version chimiste
mais on est obligé de l'inclure pour que cvsweb se porte bien
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/phylmd/diagphy.F

    r374 r429  
    1       SUBROUTINE diagphy(tit,iprt
     1      SUBROUTINE diagphy(airephy,tit,iprt
    22     $    , tops, topl, sols, soll, sens
    33     $    , evap, rain_fall, snow_fall, ts
     
    1212C
    1313C Arguments:
     14C airephy-------input-R-  grid area
    1415C tit---------input-A15- Comment to be added in PRINT (CHARACTER*15)
    1516C iprt--------input-I-  PRINT level ( <=0 : no PRINT)
     
    4546#include "dimensions.h"
    4647#include "dimphy.h"
    47 #include "comgeophy.h"
    4848#include "YOMCST.h"
    4949#include "YOETHF.h"
    5050C
    5151C     Input variables
     52      real airephy(klon)
    5253      CHARACTER*15 tit
    5354      INTEGER iprt
     
    127128     $    tit, pas, fs_bound, d_etp_tot, fq_bound, d_qt_tot
    128129C
     130      IF (iprt.ge.1) write(6,6668)
     131     $    tit, pas, d_etp_tot+d_ec_tot-fs_bound, d_qt_tot-fq_bound
     132C
    129133      IF (iprt.ge.2) write(6,6667)
    130134     $    tit, pas, stops,stopl,ssols,ssoll,ssens,slat,evap_tot
     
    133137      return
    134138
    135  6666 format('Flux Budget ',a15,1i6,2f8.2,2(1pE13.5))
    136  6667 format('Boundary Flux ',a15,1i6,6f8.2,2(1pE13.5))
     139 6666 format('Phys. Flux Budget ',a15,1i6,2f8.2,2(1pE13.5))
     140 6667 format('Phys. Boundary Flux ',a15,1i6,6f8.2,2(1pE13.5))
     141 6668 format('Phys. Total Budget ',a15,1i6,f8.2,2(1pE13.5))
    137142
    138143      end
    139144
    140145C======================================================================
    141       SUBROUTINE diagetpq(tit,iprt,idiag,idiag2,dtime
     146      SUBROUTINE diagetpq(airephy,tit,iprt,idiag,idiag2,dtime
    142147     e  ,t,q,ql,qs,u,v,paprs,pplay
    143148     s  , d_h_vcol, d_qt, d_qw, d_ql, d_qs, d_ec)
     
    160165c======================================================================
    161166C Arguments:
     167C airephy-------input-R-  grid area
    162168C tit-----imput-A15- Comment added in PRINT (CHARACTER*15)
    163169C iprt----input-I-  PRINT level ( <=1 : no PRINT)
     
    192198C
    193199C     other (COMMON...)
    194 C     RG : gravity
    195 C     airephy(klon) : mesh surface
     200C     RCPD, RCPV, ....
    196201C
    197202C J.L. Dufresne, July 2002
     
    202207#include "dimensions.h"
    203208#include "dimphy.h"
    204 #include "comgeophy.h"
    205209#include "YOMCST.h"
    206210#include "YOETHF.h"
    207211C
    208212c     Input variables
     213      real airephy(klon)
    209214      CHARACTER*15 tit
    210215      INTEGER iprt,idiag, idiag2
     
    292297C         Cinetic Energy
    293298          zec_col(i) =  zec_col(i)
    294      $        +0.5*SQRT(u(i,k)**2+v(i,k)**2)*zairm(i,k)
     299     $        +0.5*(u(i,k)**2+v(i,k)**2)*zairm(i,k)
    295300C         Air enthalpy
    296301          zh_dair_col(i) = zh_dair_col(i)
     
    376381      IF (iprt.ge.2) THEN
    377382        WRITE(6,9000) tit,pas(idiag),d_qt,d_qw,d_ql,d_qs
    378  9000   format('Watter Mass Budget (kg/m2/s)',A15
     383 9000   format('Phys. Watter Mass Budget (kg/m2/s)',A15
    379384     $      ,1i6,10(1pE14.6))
    380385        WRITE(6,9001) tit,pas(idiag), d_h_vcol
    381  9001   format('Enthalpy Budget (W/m2) ',A15,1i6,10(F8.2))
     386 9001   format('Phys. Enthalpy Budget (W/m2) ',A15,1i6,10(F8.2))
    382387        WRITE(6,9002) tit,pas(idiag), d_ec
    383  9002   format('Cinetic Energy Budget (W/m2) ',A15,1i6,10(F8.2))
     388 9002   format('Phys. Cinetic Energy Budget (W/m2) ',A15,1i6,10(F8.2))
    384389      END IF
    385390C
Note: See TracChangeset for help on using the changeset viewer.