--- trunk/libf/phylmd/Radlwsw/sw1s.f	2010/03/03 13:23:49	24
+++ trunk/libf/phylmd/Radlwsw/sw1s.f	2013/07/08 18:12:18	71
@@ -42,59 +42,59 @@
 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 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 PRMU(KDLON)
-      REAL*8 PSEC(KDLON)
-      REAL*8 PTAU(KDLON,2,KFLEV)
-      REAL*8 PUD(KDLON,5,KFLEV+1)
+      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 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 PRMU(KDLON)
+      DOUBLE PRECISION PSEC(KDLON)
+      DOUBLE PRECISION PTAU(KDLON,2,KFLEV)
+      DOUBLE PRECISION PUD(KDLON,5,KFLEV+1)
 C
-      REAL*8 PFD(KDLON,KFLEV+1)
-      REAL*8 PFU(KDLON,KFLEV+1)
+      DOUBLE PRECISION PFD(KDLON,KFLEV+1)
+      DOUBLE PRECISION PFU(KDLON,KFLEV+1)
 C
 C* LOCAL VARIABLES:
 C
       INTEGER IIND(4)
 C      
-      REAL*8 ZCGAZ(KDLON,KFLEV) 
-      REAL*8 ZDIFF(KDLON)
-      REAL*8 ZDIRF(KDLON)        
-      REAL*8 ZPIZAZ(KDLON,KFLEV)
-      REAL*8 ZRAYL(KDLON)
-      REAL*8 ZRAY1(KDLON,KFLEV+1)
-      REAL*8 ZRAY2(KDLON,KFLEV+1)
-      REAL*8 ZREFZ(KDLON,2,KFLEV+1)
-      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 ZRMUE(KDLON,KFLEV+1)
-      REAL*8 ZRMU0(KDLON,KFLEV+1)
-      REAL*8 ZR(KDLON,4)
-      REAL*8 ZTAUAZ(KDLON,KFLEV)
-      REAL*8 ZTRA1(KDLON,KFLEV+1)
-      REAL*8 ZTRA2(KDLON,KFLEV+1)
-      REAL*8 ZW(KDLON,4)
+      DOUBLE PRECISION ZCGAZ(KDLON,KFLEV) 
+      DOUBLE PRECISION ZDIFF(KDLON)
+      DOUBLE PRECISION ZDIRF(KDLON)        
+      DOUBLE PRECISION ZPIZAZ(KDLON,KFLEV)
+      DOUBLE PRECISION ZRAYL(KDLON)
+      DOUBLE PRECISION ZRAY1(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZRAY2(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZREFZ(KDLON,2,KFLEV+1)
+      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 ZRMUE(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZRMU0(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZR(KDLON,4)
+      DOUBLE PRECISION ZTAUAZ(KDLON,KFLEV)
+      DOUBLE PRECISION ZTRA1(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZTRA2(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZW(KDLON,4)
 C
       INTEGER jl, jk, k, jaj, ikm1, ikl
 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 /