--- trunk/Sources/phylmd/Radlwsw/lwu.f	2015/05/26 17:46:03	139
+++ trunk/Sources/phylmd/Radlwsw/lwu.f	2016/03/11 18:47:26	178
@@ -68,7 +68,7 @@
 
     INTEGER jl, jk, jkj, jkjr, jkjp, ig1
     INTEGER jki, jkip1, ja, jj
-    INTEGER jkl, jkp1, jkk, jkjpn
+    INTEGER jkl, jkk, jkjpn
     INTEGER jae1, jae2, jae3, jae, jjpn
     INTEGER ir, jc, jcp1
     DOUBLE PRECISION zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup
@@ -91,28 +91,28 @@
     DOUBLE PRECISION AT(8, 3), BT(8, 3)
     DOUBLE PRECISION:: OCT(4) = (/- 0.326D-3, - 0.102D-5, 0.137D-2, - 0.535D-5/)
 
-    DATA RAER / .038520, .037196, .040532, .054934, .038520, &
-         .12613, .18313, .10357, .064106, .126130, &
-         .012579, .013649, .018652, .025181, .012579, &
-         .011890, .016142, .021105, .028908, .011890, &
-         .013792, .026810, .052203, .066338, .013792 /
-
-    DATA (AT(1, IR), IR = 1, 3) / 0.298199E-02, - .394023E-03, 0.319566E-04 /
-    DATA (BT(1, IR), IR = 1, 3) / - 0.106432E-04, 0.660324E-06, 0.174356E-06 /
-    DATA (AT(2, IR), IR = 1, 3) / 0.143676E-01, 0.366501E-02, -.160822E-02 /
-    DATA (BT(2, IR), IR = 1, 3) / -0.553979E-04, - .101701E-04, 0.920868E-05 /
-    DATA (AT(3, IR), IR = 1, 3) / 0.197861E-01, 0.315541E-02, - .174547E-02 /
-    DATA (BT(3, IR), IR = 1, 3) / - 0.877012E-04, 0.513302E-04, 0.523138E-06 /
-    DATA (AT(4, IR), IR = 1, 3) / 0.289560E-01, - .208807E-02, - .121943E-02 /
-    DATA (BT(4, IR), IR = 1, 3) / - 0.165960E-03, 0.157704E-03, - .146427E-04 /
-    DATA (AT(5, IR), IR = 1, 3) / 0.103800E-01, 0.436296E-02, - .161431E-02 /
-    DATA (BT(5, IR), IR = 1, 3) / - .276744E-04, - .327381E-04, 0.127646E-04 /
-    DATA (AT(6, IR), IR = 1, 3) / 0.868859E-02, - .972752E-03, 0.000000E-00 /
-    DATA (BT(6, IR), IR = 1, 3) / - .278412E-04, - .713940E-06, 0.117469E-05 /
-    DATA (AT(7, IR), IR = 1, 3) / 0.250073E-03, 0.455875E-03, 0.109242E-03 /
-    DATA (BT(7, IR), IR = 1, 3) / 0.199846E-05, - .216313E-05, 0.175991E-06 /
-    DATA (AT(8, IR), IR = 1, 3) / 0.307423E-01, 0.110879E-02, - .322172E-03 /
-    DATA (BT(8, IR), IR = 1, 3) / - 0.108482E-03, 0.258096E-05, - .814575E-06 /
+    DATA RAER / .038520d0, .037196d0, .040532d0, .054934d0, .038520d0, &
+         .12613d0, .18313d0, .10357d0, .064106d0, .126130d0, &
+         .012579d0, .013649d0, .018652d0, .025181d0, .012579d0, &
+         .011890d0, .016142d0, .021105d0, .028908d0, .011890d0, &
+         .013792d0, .026810d0, .052203d0, .066338d0, .013792d0 /
+
+    DATA (AT(1, IR), IR = 1, 3) / 0.298199D-02, - .394023D-03, 0.319566D-04 /
+    DATA (BT(1, IR), IR = 1, 3) / - 0.106432D-04, 0.660324D-06, 0.174356D-06 /
+    DATA (AT(2, IR), IR = 1, 3) / 0.143676D-01, 0.366501D-02, -.160822D-02 /
+    DATA (BT(2, IR), IR = 1, 3) / -0.553979D-04, - .101701D-04, 0.920868D-05 /
+    DATA (AT(3, IR), IR = 1, 3) / 0.197861D-01, 0.315541D-02, - .174547D-02 /
+    DATA (BT(3, IR), IR = 1, 3) / - 0.877012D-04, 0.513302D-04, 0.523138D-06 /
+    DATA (AT(4, IR), IR = 1, 3) / 0.289560D-01, - .208807D-02, - .121943D-02 /
+    DATA (BT(4, IR), IR = 1, 3) / - 0.165960D-03, 0.157704D-03, - .146427D-04 /
+    DATA (AT(5, IR), IR = 1, 3) / 0.103800D-01, 0.436296D-02, - .161431D-02 /
+    DATA (BT(5, IR), IR = 1, 3) / - .276744D-04, - .327381D-04, 0.127646D-04 /
+    DATA (AT(6, IR), IR = 1, 3) / 0.868859D-02, - .972752D-03, 0.000000D-00 /
+    DATA (BT(6, IR), IR = 1, 3) / - .278412D-04, - .713940D-06, 0.117469D-05 /
+    DATA (AT(7, IR), IR = 1, 3) / 0.250073D-03, 0.455875D-03, 0.109242D-03 /
+    DATA (BT(7, IR), IR = 1, 3) / 0.199846D-05, - .216313D-05, 0.175991D-06 /
+    DATA (AT(8, IR), IR = 1, 3) / 0.307423D-01, 0.110879D-02, - .322172D-03 /
+    DATA (BT(8, IR), IR = 1, 3) / - 0.108482D-03, 0.258096D-05, - .814575D-06 /
 
     !-----------------------------------------------------------------------
 
@@ -159,7 +159,6 @@
     end DO
 
     DO JK = 1, KFLEV
-       JKP1 = JK + 1
        JKL = KFLEV + 1 - JK
        DO JL = 1, KDLON
           ZXWV(JL) = MAX(PWV(JL, JK), ZEPSCQ)