/[lmdze]/trunk/Sources/phylmd/diagphy.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/diagphy.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 98 by guez, Tue May 13 17:23:16 2014 UTC
# Line 5  module diagphy_m Line 5  module diagphy_m
5  contains  contains
6    
7    SUBROUTINE diagphy(airephy, tit, iprt, tops, topl, sols, soll, sens, evap, &    SUBROUTINE diagphy(airephy, tit, iprt, tops, topl, sols, soll, sens, evap, &
8         rain_fall, snow_fall, ts, d_etp_tot, d_qt_tot, d_ec_tot, fs_bound, &         rain_fall, snow_fall, ts, d_etp_tot, d_qt_tot, d_ec_tot)
        fq_bound)  
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    
# Line 22  contains Line 21  contains
21      ! Arguments:      ! Arguments:
22    
23      ! Input variables      ! Input variables
24      real airephy(klon)      real, intent(in):: airephy(klon) ! grid area
25      ! airephy-------input-R- grid area      CHARACTER(len=15), intent(in):: tit ! comment to be added in PRINT
26      CHARACTER(len=15) tit      INTEGER, intent(in):: iprt ! PRINT level (<=0 : no PRINT)
27      ! tit---------input-A15- Comment to be added in PRINT (CHARACTER*15)      real, intent(in):: tops(klon) ! SW rad. at TOA (W/m2), positive up
28      INTEGER iprt      real, intent(in):: topl(klon) ! LW rad. at TOA (W/m2), positive down
29      ! iprt--------input-I- PRINT level (<=0 : no PRINT)  
30      real tops(klon), sols(klon)      real, intent(in):: sols(klon)
31      ! tops(klon)--input-R- SW rad. at TOA (W/m2), positive up.      ! net SW flux above surface (W/m2), positive up (i.e. -1 * flux
32      ! sols(klon)--input-R- Net SW flux above surface (W/m2), positive up      ! absorbed by the surface)
     ! (i.e. -1 * flux absorbed by the surface)  
33    
34      real, intent(in):: soll(klon)      real, intent(in):: soll(klon)
35      ! net longwave flux above surface (W/m2), positive up (i. e. flux emited      ! net longwave flux above surface (W/m2), positive up (i. e. flux
36      ! - flux absorbed by the surface)      ! emited - flux absorbed by the surface)
37    
38      real, intent(in):: topl(klon) !LW rad. at TOA (W/m2), positive down      real, intent(in):: sens(klon)
39      real sens(klon)      ! sensible Flux at surface (W/m2), positive down
40      ! sens(klon)--input-R- Sensible Flux at surface (W/m2), positive down  
41      real evap(klon)      real, intent(in):: evap(klon)
42      ! evap(klon)--input-R- Evaporation + sublimation water vapour mass flux      ! evaporation + sublimation water vapour mass flux (kg/m2/s),
43      ! (kg/m2/s), positive up      ! positive up
44    
45      real, intent(in):: rain_fall(klon)      real, intent(in):: rain_fall(klon)
46      ! liquid water mass flux (kg/m2/s), positive down      ! liquid water mass flux (kg/m2/s), positive down
47    
48      real snow_fall(klon)      real, intent(in):: snow_fall(klon)
49      ! snow_fall(klon)      ! solid water mass flux (kg/m2/s), positive down
50      ! --input-R- Solid water mass flux (kg/m2/s), positive down  
51      REAL ts(klon)      REAL, intent(in):: ts(klon) ! surface temperature (K)
52      ! ts(klon)----input-R- Surface temperature (K)  
53      REAL d_etp_tot, d_qt_tot, d_ec_tot      REAL, intent(in):: d_etp_tot
54      ! d_etp_tot---input-R- Heat flux equivalent to atmospheric enthalpy      ! heat flux equivalent to atmospheric enthalpy change (W/m2)
     ! change (W/m2)  
     ! d_qt_tot----input-R- Mass flux equivalent to atmospheric water mass  
     ! change (kg/m2/s)  
     ! d_ec_tot----input-R- Flux equivalent to atmospheric cinetic energy  
     ! change (W/m2)  
   
     ! Output variables  
     REAL fs_bound  
     ! fs_bound---output-R- Thermal flux at the atmosphere boundaries (W/m2)  
     real fq_bound  
     ! fq_bound---output-R- Water mass flux at the atmosphere  
     ! boundaries (kg/m2/s)  
55    
56      ! Local variables:      REAL, intent(in):: d_qt_tot
57        ! Mass flux equivalent to atmospheric water mass change (kg/m2/s)
58    
59        REAL, intent(in):: d_ec_tot
60        ! flux equivalent to atmospheric cinetic energy change (W/m2)
61    
62        ! Local:
63        REAL fs_bound ! thermal flux at the atmosphere boundaries (W/m2)
64        real fq_bound ! water mass flux at the atmosphere boundaries (kg/m2/s)
65      real stops, stopl, ssols, ssoll      real stops, stopl, ssols, ssoll
66      real ssens, sfront, slat      real ssens, sfront, slat
67      real airetot, zcpvap, zcwat, zcice      real airetot, zcpvap, zcwat, zcice
68      REAL rain_fall_tot, snow_fall_tot, evap_tot      REAL rain_fall_tot, snow_fall_tot, evap_tot
   
69      integer i      integer i
70      integer:: pas = 0      integer:: pas = 0
71    
72      !------------------------------------------------------------------      !------------------------------------------------------------------
73    
74      pas=pas+1      IF (iprt >= 1) then
75      stops=0.         pas=pas+1
76      stopl=0.         stops=0.
77      ssols=0.         stopl=0.
78      ssoll=0.         ssols=0.
79      ssens=0.         ssoll=0.
80      sfront = 0.         ssens=0.
81      evap_tot = 0.         sfront = 0.
82      rain_fall_tot = 0.         evap_tot = 0.
83      snow_fall_tot = 0.         rain_fall_tot = 0.
84      airetot=0.         snow_fall_tot = 0.
85           airetot=0.
86      ! Pour les chaleur specifiques de la vapeur d'eau, de l'eau et de  
87      ! la glace, on travaille par difference a la chaleur specifique de         ! Pour les chaleurs spécifiques de la vapeur d'eau, de l'eau et de
88      ! l' air sec. En effet, comme on travaille a niveau de pression         ! la glace, on travaille par différence à la chaleur spécifique de
89      ! donne, toute variation de la masse d'un constituant est         ! l'air sec. En effet, comme on travaille à niveau de pression
90      ! totalement compense par une variation de masse d'air.         ! donné, toute variation de la masse d'un constituant est
91           ! totalement compensée par une variation de masse d'air.
92      zcpvap=RCPV-RCPD  
93      zcwat=RCW-RCPD         zcpvap=RCPV-RCPD
94      zcice=RCS-RCPD         zcwat=RCW-RCPD
95           zcice=RCS-RCPD
96      do i=1, klon  
97         stops=stops+tops(i)*airephy(i)         do i=1, klon
98         stopl=stopl+topl(i)*airephy(i)            stops=stops+tops(i)*airephy(i)
99         ssols=ssols+sols(i)*airephy(i)            stopl=stopl+topl(i)*airephy(i)
100         ssoll=ssoll+soll(i)*airephy(i)            ssols=ssols+sols(i)*airephy(i)
101         ssens=ssens+sens(i)*airephy(i)            ssoll=ssoll+soll(i)*airephy(i)
102         sfront = sfront &            ssens=ssens+sens(i)*airephy(i)
103              + (evap(i)*zcpvap-rain_fall(i)*zcwat-snow_fall(i)*zcice) * ts(i) &            sfront = sfront + (evap(i) * zcpvap - rain_fall(i) * zcwat &
104              * airephy(i)                 - snow_fall(i) * zcice) * ts(i) * airephy(i)
105         evap_tot = evap_tot + evap(i)*airephy(i)            evap_tot = evap_tot + evap(i)*airephy(i)
106         rain_fall_tot = rain_fall_tot + rain_fall(i)*airephy(i)            rain_fall_tot = rain_fall_tot + rain_fall(i)*airephy(i)
107         snow_fall_tot = snow_fall_tot + snow_fall(i)*airephy(i)            snow_fall_tot = snow_fall_tot + snow_fall(i)*airephy(i)
108         airetot=airetot+airephy(i)            airetot=airetot+airephy(i)
109      enddo         enddo
110      stops=stops/airetot         stops=stops/airetot
111      stopl=stopl/airetot         stopl=stopl/airetot
112      ssols=ssols/airetot         ssols=ssols/airetot
113      ssoll=ssoll/airetot         ssoll=ssoll/airetot
114      ssens=ssens/airetot         ssens=ssens/airetot
115      sfront = sfront/airetot         sfront = sfront/airetot
116      evap_tot = evap_tot /airetot         evap_tot = evap_tot /airetot
117      rain_fall_tot = rain_fall_tot/airetot         rain_fall_tot = rain_fall_tot/airetot
118      snow_fall_tot = snow_fall_tot/airetot         snow_fall_tot = snow_fall_tot/airetot
119    
120      slat = RLVTT * rain_fall_tot + RLSTT * snow_fall_tot         slat = RLVTT * rain_fall_tot + RLSTT * snow_fall_tot
121      ! Heat flux at atm. boundaries         fs_bound = stops-stopl - (ssols+ssoll)+ssens+sfront + slat
122      fs_bound = stops-stopl - (ssols+ssoll)+ssens+sfront + slat         fq_bound = evap_tot - rain_fall_tot -snow_fall_tot
123      ! Water flux at atm. boundaries  
124      fq_bound = evap_tot - rain_fall_tot -snow_fall_tot         print 6666, tit, pas, fs_bound, d_etp_tot, fq_bound, d_qt_tot
125           print 6668, tit, pas, d_etp_tot+d_ec_tot-fs_bound, d_qt_tot - fq_bound
126      IF (iprt >= 1) print 6666, tit, pas, fs_bound, d_etp_tot, fq_bound, d_qt_tot  
127           IF (iprt >= 2) print 6667, tit, pas, stops, stopl, ssols, ssoll, ssens, &
128      IF (iprt >= 1) print 6668, tit, pas, d_etp_tot+d_ec_tot-fs_bound, &              slat, evap_tot, rain_fall_tot + snow_fall_tot
129           d_qt_tot-fq_bound      end IF
130    
131      IF (iprt >= 2) print 6667, tit, pas, stops, stopl, ssols, ssoll, ssens, &  6666 format('Physics flux budget ', a15, 1i6, 2f8.2, 2(1pE13.5))
132           slat, evap_tot, rain_fall_tot+snow_fall_tot  6667 format('Physics boundary flux ', a15, 1i6, 6f8.2, 2(1pE13.5))
133    6668 format('Physics total budget ', a15, 1i6, f8.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))  
134    
135    end SUBROUTINE diagphy    end SUBROUTINE diagphy
136    
Legend:
Removed from v.82  
changed lines
  Added in v.98
  ViewVC Help
Powered by ViewVC 1.1.21