--- trunk/libf/phylmd/Radlwsw/swde.f	2010/03/03 13:23:49	24
+++ trunk/libf/phylmd/Radlwsw/swde.f	2013/07/08 18:12:18	71
@@ -32,27 +32,28 @@
 C     ------------------------------------------------------------------
 C* ARGUMENTS:
 C
-      REAL*8 PGG(KDLON)   ! ASSYMETRY FACTOR
-      REAL*8 PREF(KDLON)  ! REFLECTIVITY OF THE UNDERLYING LAYER
-      REAL*8 PRMUZ(KDLON) ! COSINE OF SOLAR ZENITH ANGLE
-      REAL*8 PTO1(KDLON)  ! OPTICAL THICKNESS
-      REAL*8 PW(KDLON)    ! SINGLE SCATTERING ALBEDO
-      REAL*8 PRE1(KDLON)  ! LAYER REFLECTIVITY (NO UNDERLYING-LAYER REFLECTION)
-      REAL*8 PRE2(KDLON)  ! LAYER REFLECTIVITY
-      REAL*8 PTR1(KDLON)  ! LAYER TRANSMISSIVITY (NO UNDERLYING-LAYER REFLECTION)
-      REAL*8 PTR2(KDLON)  ! LAYER TRANSMISSIVITY
+      DOUBLE PRECISION PGG(KDLON)   ! ASSYMETRY FACTOR
+      DOUBLE PRECISION PREF(KDLON)  ! REFLECTIVITY OF THE UNDERLYING LAYER
+      DOUBLE PRECISION PRMUZ(KDLON) ! COSINE OF SOLAR ZENITH ANGLE
+      DOUBLE PRECISION PTO1(KDLON)  ! OPTICAL THICKNESS
+      DOUBLE PRECISION PW(KDLON)    ! SINGLE SCATTERING ALBEDO
+      DOUBLE PRECISION PRE1(KDLON)  ! LAYER REFLECTIVITY (NO UNDERLYING-LAYER REFLECTION)
+      DOUBLE PRECISION PRE2(KDLON)  ! LAYER REFLECTIVITY
+      DOUBLE PRECISION PTR1(KDLON)  ! LAYER TRANSMISSIVITY (NO UNDERLYING-LAYER REFLECTION)
+      DOUBLE PRECISION PTR2(KDLON)  ! LAYER TRANSMISSIVITY
 C
 C* LOCAL VARIABLES:
 C
       INTEGER jl
-      REAL*8 ZFF, ZGP, ZTOP, ZWCP, ZDT, ZX1, ZWM
-      REAL*8 ZRM2, ZRK, ZX2, ZRP, ZALPHA, ZBETA, ZARG
-      REAL*8 ZEXMU0, ZARG2, ZEXKP, ZEXKM, ZXP2P, ZXM2P, ZAP2B, ZAM2B
-      REAL*8 ZA11, ZA12, ZA13, ZA21, ZA22, ZA23
-      REAL*8 ZDENA, ZC1A, ZC2A, ZRI0A, ZRI1A
-      REAL*8 ZRI0B, ZRI1B
-      REAL*8 ZB21, ZB22, ZB23, ZDENB, ZC1B, ZC2B
-      REAL*8 ZRI0C, ZRI1C, ZRI0D, ZRI1D
+      DOUBLE PRECISION ZFF, ZGP, ZTOP, ZWCP, ZDT, ZX1, ZWM
+      DOUBLE PRECISION ZRM2, ZRK, ZX2, ZRP, ZALPHA, ZBETA, ZARG
+      DOUBLE PRECISION ZEXMU0, ZARG2, ZEXKP, ZEXKM, ZXP2P, ZXM2P, ZAP2B
+      DOUBLE PRECISION ZAM2B
+      DOUBLE PRECISION ZA11, ZA12, ZA13, ZA21, ZA22, ZA23
+      DOUBLE PRECISION ZDENA, ZC1A, ZC2A, ZRI0A, ZRI1A
+      DOUBLE PRECISION ZRI0B, ZRI1B
+      DOUBLE PRECISION ZB21, ZB22, ZB23, ZDENB, ZC1B, ZC2B
+      DOUBLE PRECISION ZRI0C, ZRI1C, ZRI0D, ZRI1D
 C     ------------------------------------------------------------------
 C
 C*         1.      DELTA-EDDINGTON CALCULATIONS