/[lmdze]/trunk/phylmd/Radlwsw/lwu.f
ViewVC logotype

Diff of /trunk/phylmd/Radlwsw/lwu.f

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

trunk/libf/phylmd/Radlwsw/lwu.f90 revision 71 by guez, Mon Jul 8 18:12:18 2013 UTC trunk/Sources/phylmd/Radlwsw/lwu.f revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC
# Line 4  module LWU_m Line 4  module LWU_m
4    
5  contains  contains
6    
7    SUBROUTINE LWU(PAER, PDP, PPMB, PPSOL, POZ, PTAVE, PVIEW, PWV, PABCU)    SUBROUTINE LWU(PAER, PDP, PPMB, POZ, PTAVE, PVIEW, PWV, PABCU)
8    
9      ! Purpose. Computes absorber amounts including pressure and      ! Purpose. Computes absorber amounts including pressure and
10      ! temperature effects.      ! temperature effects.
# Line 33  contains Line 33  contains
33      DOUBLE PRECISION PAER(KDLON, KFLEV, 5)      DOUBLE PRECISION PAER(KDLON, KFLEV, 5)
34      DOUBLE PRECISION PDP(KDLON, KFLEV)      DOUBLE PRECISION PDP(KDLON, KFLEV)
35      DOUBLE PRECISION PPMB(KDLON, KFLEV + 1)      DOUBLE PRECISION PPMB(KDLON, KFLEV + 1)
     DOUBLE PRECISION PPSOL(KDLON)  
36      DOUBLE PRECISION POZ(KDLON, KFLEV)      DOUBLE PRECISION POZ(KDLON, KFLEV)
37      DOUBLE PRECISION PTAVE(KDLON, KFLEV)      DOUBLE PRECISION PTAVE(KDLON, KFLEV)
38      DOUBLE PRECISION PVIEW(KDLON)      DOUBLE PRECISION PVIEW(KDLON)
# Line 69  contains Line 68  contains
68    
69      INTEGER jl, jk, jkj, jkjr, jkjp, ig1      INTEGER jl, jk, jkj, jkjr, jkjp, ig1
70      INTEGER jki, jkip1, ja, jj      INTEGER jki, jkip1, ja, jj
71      INTEGER jkl, jkp1, jkk, jkjpn      INTEGER jkl, jkk, jkjpn
72      INTEGER jae1, jae2, jae3, jae, jjpn      INTEGER jae1, jae2, jae3, jae, jjpn
73      INTEGER ir, jc, jcp1      INTEGER ir, jc, jcp1
74      DOUBLE PRECISION zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup      DOUBLE PRECISION zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup
# Line 92  contains Line 91  contains
91      DOUBLE PRECISION AT(8, 3), BT(8, 3)      DOUBLE PRECISION AT(8, 3), BT(8, 3)
92      DOUBLE PRECISION:: OCT(4) = (/- 0.326D-3, - 0.102D-5, 0.137D-2, - 0.535D-5/)      DOUBLE PRECISION:: OCT(4) = (/- 0.326D-3, - 0.102D-5, 0.137D-2, - 0.535D-5/)
93    
94      DATA RAER / .038520, .037196, .040532, .054934, .038520, &      DATA RAER / .038520d0, .037196d0, .040532d0, .054934d0, .038520d0, &
95           .12613, .18313, .10357, .064106, .126130, &           .12613d0, .18313d0, .10357d0, .064106d0, .126130d0, &
96           .012579, .013649, .018652, .025181, .012579, &           .012579d0, .013649d0, .018652d0, .025181d0, .012579d0, &
97           .011890, .016142, .021105, .028908, .011890, &           .011890d0, .016142d0, .021105d0, .028908d0, .011890d0, &
98           .013792, .026810, .052203, .066338, .013792 /           .013792d0, .026810d0, .052203d0, .066338d0, .013792d0 /
99    
100      DATA (AT(1, IR), IR = 1, 3) / 0.298199E-02, - .394023E-03, 0.319566E-04 /      DATA (AT(1, IR), IR = 1, 3) / 0.298199D-02, - .394023D-03, 0.319566D-04 /
101      DATA (BT(1, IR), IR = 1, 3) / - 0.106432E-04, 0.660324E-06, 0.174356E-06 /      DATA (BT(1, IR), IR = 1, 3) / - 0.106432D-04, 0.660324D-06, 0.174356D-06 /
102      DATA (AT(2, IR), IR = 1, 3) / 0.143676E-01, 0.366501E-02, -.160822E-02 /      DATA (AT(2, IR), IR = 1, 3) / 0.143676D-01, 0.366501D-02, -.160822D-02 /
103      DATA (BT(2, IR), IR = 1, 3) / -0.553979E-04, - .101701E-04, 0.920868E-05 /      DATA (BT(2, IR), IR = 1, 3) / -0.553979D-04, - .101701D-04, 0.920868D-05 /
104      DATA (AT(3, IR), IR = 1, 3) / 0.197861E-01, 0.315541E-02, - .174547E-02 /      DATA (AT(3, IR), IR = 1, 3) / 0.197861D-01, 0.315541D-02, - .174547D-02 /
105      DATA (BT(3, IR), IR = 1, 3) / - 0.877012E-04, 0.513302E-04, 0.523138E-06 /      DATA (BT(3, IR), IR = 1, 3) / - 0.877012D-04, 0.513302D-04, 0.523138D-06 /
106      DATA (AT(4, IR), IR = 1, 3) / 0.289560E-01, - .208807E-02, - .121943E-02 /      DATA (AT(4, IR), IR = 1, 3) / 0.289560D-01, - .208807D-02, - .121943D-02 /
107      DATA (BT(4, IR), IR = 1, 3) / - 0.165960E-03, 0.157704E-03, - .146427E-04 /      DATA (BT(4, IR), IR = 1, 3) / - 0.165960D-03, 0.157704D-03, - .146427D-04 /
108      DATA (AT(5, IR), IR = 1, 3) / 0.103800E-01, 0.436296E-02, - .161431E-02 /      DATA (AT(5, IR), IR = 1, 3) / 0.103800D-01, 0.436296D-02, - .161431D-02 /
109      DATA (BT(5, IR), IR = 1, 3) / - .276744E-04, - .327381E-04, 0.127646E-04 /      DATA (BT(5, IR), IR = 1, 3) / - .276744D-04, - .327381D-04, 0.127646D-04 /
110      DATA (AT(6, IR), IR = 1, 3) / 0.868859E-02, - .972752E-03, 0.000000E-00 /      DATA (AT(6, IR), IR = 1, 3) / 0.868859D-02, - .972752D-03, 0.000000D-00 /
111      DATA (BT(6, IR), IR = 1, 3) / - .278412E-04, - .713940E-06, 0.117469E-05 /      DATA (BT(6, IR), IR = 1, 3) / - .278412D-04, - .713940D-06, 0.117469D-05 /
112      DATA (AT(7, IR), IR = 1, 3) / 0.250073E-03, 0.455875E-03, 0.109242E-03 /      DATA (AT(7, IR), IR = 1, 3) / 0.250073D-03, 0.455875D-03, 0.109242D-03 /
113      DATA (BT(7, IR), IR = 1, 3) / 0.199846E-05, - .216313E-05, 0.175991E-06 /      DATA (BT(7, IR), IR = 1, 3) / 0.199846D-05, - .216313D-05, 0.175991D-06 /
114      DATA (AT(8, IR), IR = 1, 3) / 0.307423E-01, 0.110879E-02, - .322172E-03 /      DATA (AT(8, IR), IR = 1, 3) / 0.307423D-01, 0.110879D-02, - .322172D-03 /
115      DATA (BT(8, IR), IR = 1, 3) / - 0.108482E-03, 0.258096E-05, - .814575E-06 /      DATA (BT(8, IR), IR = 1, 3) / - 0.108482D-03, 0.258096D-05, - .814575D-06 /
116    
117      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
118    
# Line 160  contains Line 159  contains
159      end DO      end DO
160    
161      DO JK = 1, KFLEV      DO JK = 1, KFLEV
        JKP1 = JK + 1  
162         JKL = KFLEV + 1 - JK         JKL = KFLEV + 1 - JK
163         DO JL = 1, KDLON         DO JL = 1, KDLON
164            ZXWV(JL) = MAX(PWV(JL, JK), ZEPSCQ)            ZXWV(JL) = MAX(PWV(JL, JK), ZEPSCQ)
Legend:
Removed from v.71  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21