--- trunk/libf/phylmd/Radlwsw/sw2s.f	2010/03/03 13:23:49	24
+++ trunk/phylmd/Radlwsw/sw2s.f	2013/11/15 18:45:49	76
@@ -48,87 +48,87 @@
 C
       INTEGER KNU
 c-OB
-      real*8 flag_aer
-      real*8 tauae(kdlon,kflev,2)
-      real*8 pizae(kdlon,kflev,2)
-      real*8 cgae(kdlon,kflev,2)
-      REAL*8 PAER(KDLON,KFLEV,5)
-      REAL*8 PAKI(KDLON,2)
-      REAL*8 PALBD(KDLON,2)
-      REAL*8 PALBP(KDLON,2)
-      REAL*8 PCG(KDLON,2,KFLEV)
-      REAL*8 PCLD(KDLON,KFLEV)
-      REAL*8 PCLDSW(KDLON,KFLEV)
-      REAL*8 PCLEAR(KDLON)
-      REAL*8 PDSIG(KDLON,KFLEV)
-      REAL*8 POMEGA(KDLON,2,KFLEV)
-      REAL*8 POZ(KDLON,KFLEV)
-      REAL*8 PQS(KDLON,KFLEV)
-      REAL*8 PRMU(KDLON)
-      REAL*8 PSEC(KDLON)
-      REAL*8 PTAU(KDLON,2,KFLEV)
-      REAL*8 PUD(KDLON,5,KFLEV+1)
-      REAL*8 PWV(KDLON,KFLEV)
+      double precision flag_aer
+      double precision tauae(kdlon,kflev,2)
+      double precision pizae(kdlon,kflev,2)
+      double precision cgae(kdlon,kflev,2)
+      DOUBLE PRECISION PAER(KDLON,KFLEV,5)
+      DOUBLE PRECISION PAKI(KDLON,2)
+      DOUBLE PRECISION PALBD(KDLON,2)
+      DOUBLE PRECISION PALBP(KDLON,2)
+      DOUBLE PRECISION PCG(KDLON,2,KFLEV)
+      DOUBLE PRECISION PCLD(KDLON,KFLEV)
+      DOUBLE PRECISION PCLDSW(KDLON,KFLEV)
+      DOUBLE PRECISION PCLEAR(KDLON)
+      DOUBLE PRECISION PDSIG(KDLON,KFLEV)
+      DOUBLE PRECISION POMEGA(KDLON,2,KFLEV)
+      DOUBLE PRECISION POZ(KDLON,KFLEV)
+      DOUBLE PRECISION PQS(KDLON,KFLEV)
+      DOUBLE PRECISION PRMU(KDLON)
+      DOUBLE PRECISION PSEC(KDLON)
+      DOUBLE PRECISION PTAU(KDLON,2,KFLEV)
+      DOUBLE PRECISION PUD(KDLON,5,KFLEV+1)
+      DOUBLE PRECISION PWV(KDLON,KFLEV)
 C
-      REAL*8 PFDOWN(KDLON,KFLEV+1)
-      REAL*8 PFUP(KDLON,KFLEV+1)
+      DOUBLE PRECISION PFDOWN(KDLON,KFLEV+1)
+      DOUBLE PRECISION PFUP(KDLON,KFLEV+1)
 C
 C* LOCAL VARIABLES:
 C
       INTEGER IIND2(2), IIND3(3)
-      REAL*8 ZCGAZ(KDLON,KFLEV)
-      REAL*8 ZFD(KDLON,KFLEV+1)
-      REAL*8 ZFU(KDLON,KFLEV+1) 
-      REAL*8 ZG(KDLON)
-      REAL*8 ZGG(KDLON)
-      REAL*8 ZPIZAZ(KDLON,KFLEV)
-      REAL*8 ZRAYL(KDLON)
-      REAL*8 ZRAY1(KDLON,KFLEV+1)
-      REAL*8 ZRAY2(KDLON,KFLEV+1)
-      REAL*8 ZREF(KDLON)
-      REAL*8 ZREFZ(KDLON,2,KFLEV+1)
-      REAL*8 ZRE1(KDLON)
-      REAL*8 ZRE2(KDLON)
-      REAL*8 ZRJ(KDLON,6,KFLEV+1)
-      REAL*8 ZRJ0(KDLON,6,KFLEV+1)
-      REAL*8 ZRK(KDLON,6,KFLEV+1)
-      REAL*8 ZRK0(KDLON,6,KFLEV+1)
-      REAL*8 ZRL(KDLON,8)
-      REAL*8 ZRMUE(KDLON,KFLEV+1)
-      REAL*8 ZRMU0(KDLON,KFLEV+1)
-      REAL*8 ZRMUZ(KDLON)
-      REAL*8 ZRNEB(KDLON)
-      REAL*8 ZRUEF(KDLON,8)
-      REAL*8 ZR1(KDLON) 
-      REAL*8 ZR2(KDLON,2)
-      REAL*8 ZR3(KDLON,3)
-      REAL*8 ZR4(KDLON)
-      REAL*8 ZR21(KDLON)
-      REAL*8 ZR22(KDLON)
-      REAL*8 ZS(KDLON)
-      REAL*8 ZTAUAZ(KDLON,KFLEV)
-      REAL*8 ZTO1(KDLON)
-      REAL*8 ZTR(KDLON,2,KFLEV+1)
-      REAL*8 ZTRA1(KDLON,KFLEV+1)
-      REAL*8 ZTRA2(KDLON,KFLEV+1)
-      REAL*8 ZTR1(KDLON)
-      REAL*8 ZTR2(KDLON)
-      REAL*8 ZW(KDLON)   
-      REAL*8 ZW1(KDLON)
-      REAL*8 ZW2(KDLON,2)
-      REAL*8 ZW3(KDLON,3)
-      REAL*8 ZW4(KDLON)
-      REAL*8 ZW5(KDLON)
+      DOUBLE PRECISION ZCGAZ(KDLON,KFLEV)
+      DOUBLE PRECISION ZFD(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZFU(KDLON,KFLEV+1) 
+      DOUBLE PRECISION ZG(KDLON)
+      DOUBLE PRECISION ZGG(KDLON)
+      DOUBLE PRECISION ZPIZAZ(KDLON,KFLEV)
+      DOUBLE PRECISION ZRAYL(KDLON)
+      DOUBLE PRECISION ZRAY1(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZRAY2(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZREF(KDLON)
+      DOUBLE PRECISION ZREFZ(KDLON,2,KFLEV+1)
+      DOUBLE PRECISION ZRE1(KDLON)
+      DOUBLE PRECISION ZRE2(KDLON)
+      DOUBLE PRECISION ZRJ(KDLON,6,KFLEV+1)
+      DOUBLE PRECISION ZRJ0(KDLON,6,KFLEV+1)
+      DOUBLE PRECISION ZRK(KDLON,6,KFLEV+1)
+      DOUBLE PRECISION ZRK0(KDLON,6,KFLEV+1)
+      DOUBLE PRECISION ZRL(KDLON,8)
+      DOUBLE PRECISION ZRMUE(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZRMU0(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZRMUZ(KDLON)
+      DOUBLE PRECISION ZRNEB(KDLON)
+      DOUBLE PRECISION ZRUEF(KDLON,8)
+      DOUBLE PRECISION ZR1(KDLON) 
+      DOUBLE PRECISION ZR2(KDLON,2)
+      DOUBLE PRECISION ZR3(KDLON,3)
+      DOUBLE PRECISION ZR4(KDLON)
+      DOUBLE PRECISION ZR21(KDLON)
+      DOUBLE PRECISION ZR22(KDLON)
+      DOUBLE PRECISION ZS(KDLON)
+      DOUBLE PRECISION ZTAUAZ(KDLON,KFLEV)
+      DOUBLE PRECISION ZTO1(KDLON)
+      DOUBLE PRECISION ZTR(KDLON,2,KFLEV+1)
+      DOUBLE PRECISION ZTRA1(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZTRA2(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZTR1(KDLON)
+      DOUBLE PRECISION ZTR2(KDLON)
+      DOUBLE PRECISION ZW(KDLON)   
+      DOUBLE PRECISION ZW1(KDLON)
+      DOUBLE PRECISION ZW2(KDLON,2)
+      DOUBLE PRECISION ZW3(KDLON,3)
+      DOUBLE PRECISION ZW4(KDLON)
+      DOUBLE PRECISION ZW5(KDLON)
 C
       INTEGER jl, jk, k, jaj, ikm1, ikl, jn, jabs, jkm1
       INTEGER jref, jkl, jklp1, jajp, jkki, jkkp4, jn2j, iabs
-      REAL*8 ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
+      DOUBLE PRECISION ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
 C
 C* Prescribed Data:
 C
-      REAL*8 RSUN(2)
+      DOUBLE PRECISION RSUN(2)
       SAVE RSUN
-      REAL*8 RRAY(2,6)
+      DOUBLE PRECISION RRAY(2,6)
       SAVE RRAY
       DATA RSUN(1) / 0.441676 /
       DATA RSUN(2) / 0.558324 /