4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE diagphy(airephy, tit, iprt, tops, topl, sols, soll, sens, evap, & |
SUBROUTINE diagphy(airephy, tit, ip_ebil, tops, topl, sols, soll, sens, & |
8 |
rain_fall, snow_fall, ts, d_etp_tot, d_qt_tot, d_ec_tot) |
evap, rain_fall, snow_fall, ts, d_etp_tot, d_qt_tot, d_ec_tot) |
9 |
|
|
10 |
! From LMDZ4/libf/phylmd/diagphy.F, version 1.1.1.1 2004/05/19 12:53:08 |
! From LMDZ4/libf/phylmd/diagphy.F, version 1.1.1.1, 2004/05/19 12:53:08 |
11 |
|
|
12 |
! Purpose: compute the thermal flux and the water mass flux at |
! Purpose: compute the thermal flux and the water mass flux at |
13 |
! the atmospheric boundaries. Print them and print the atmospheric |
! the atmospheric boundaries. Print them and print the atmospheric |
18 |
USE dimphy, ONLY: klon |
USE dimphy, ONLY: klon |
19 |
USE suphec_m, ONLY: rcpd, rcpv, rcs, rcw, rlstt, rlvtt |
USE suphec_m, ONLY: rcpd, rcpv, rcs, rcw, rlstt, rlvtt |
20 |
|
|
21 |
! Arguments: |
real, intent(in):: airephy(:) ! (klon) grid area |
|
|
|
|
! Input variables |
|
|
real, intent(in):: airephy(klon) ! grid area |
|
22 |
CHARACTER(len=15), intent(in):: tit ! comment to be added in PRINT |
CHARACTER(len=15), intent(in):: tit ! comment to be added in PRINT |
23 |
INTEGER, intent(in):: iprt ! PRINT level (<=0 : no PRINT) |
INTEGER, intent(in):: ip_ebil ! PRINT level (<=0 : no PRINT) |
24 |
real, intent(in):: tops(klon) ! SW rad. at TOA (W/m2), positive up |
real, intent(in):: tops(klon) ! SW rad. at TOA (W/m2), positive up |
25 |
real, intent(in):: topl(klon) ! LW rad. at TOA (W/m2), positive down |
real, intent(in):: topl(klon) ! LW rad. at TOA (W/m2), positive down |
26 |
|
|
68 |
|
|
69 |
!------------------------------------------------------------------ |
!------------------------------------------------------------------ |
70 |
|
|
71 |
IF (iprt >= 1) then |
IF (ip_ebil >= 1) then |
72 |
pas=pas+1 |
pas=pas+1 |
73 |
stops=0. |
stops=0. |
74 |
stopl=0. |
stopl=0. |
121 |
print 6666, tit, pas, fs_bound, d_etp_tot, fq_bound, d_qt_tot |
print 6666, tit, pas, fs_bound, d_etp_tot, fq_bound, d_qt_tot |
122 |
print 6668, tit, pas, d_etp_tot+d_ec_tot-fs_bound, d_qt_tot - fq_bound |
print 6668, tit, pas, d_etp_tot+d_ec_tot-fs_bound, d_qt_tot - fq_bound |
123 |
|
|
124 |
IF (iprt >= 2) print 6667, tit, pas, stops, stopl, ssols, ssoll, ssens, & |
IF (ip_ebil >= 2) print 6667, tit, pas, stops, stopl, ssols, ssoll, & |
125 |
slat, evap_tot, rain_fall_tot + snow_fall_tot |
ssens, slat, evap_tot, rain_fall_tot + snow_fall_tot |
126 |
end IF |
end IF |
127 |
|
|
128 |
6666 format('Physics flux budget ', a15, 1i6, 2f8.2, 2(1pE13.5)) |
6666 format('Physics flux budget ', a15, 1i6, 2f8.2, 2(1pE13.5)) |