--- trunk/libf/phylmd/Radlwsw/lwv.f	2013/06/24 15:39:52	70
+++ trunk/libf/phylmd/Radlwsw/lwv.f	2013/07/08 18:12:18	71
@@ -42,33 +42,33 @@
 C* ARGUMENTS:
       INTEGER KUAER,KTRAER, KLIM
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
-      REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
-      REAL*8 PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
-      REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
-      REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
-      REAL*8 PBTOP(KDLON,Ninter) ! T.O.A. SPECTRAL PLANCK FUNCTION
-      REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
-      REAL*8 PEMIS(KDLON) ! SURFACE EMISSIVITY
-      REAL*8 PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
-      REAL*8 PTAVE(KDLON,KFLEV) ! TEMPERATURE
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGASUR(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGBSUR(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGATOP(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGBTOP(KDLON,8,2) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
+      DOUBLE PRECISION PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
+      DOUBLE PRECISION PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
+      DOUBLE PRECISION PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
+      DOUBLE PRECISION PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
+      DOUBLE PRECISION PBTOP(KDLON,Ninter) ! T.O.A. SPECTRAL PLANCK FUNCTION
+      DOUBLE PRECISION PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
+      DOUBLE PRECISION PEMIS(KDLON) ! SURFACE EMISSIVITY
+      DOUBLE PRECISION PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
+      DOUBLE PRECISION PTAVE(KDLON,KFLEV) ! TEMPERATURE
+      DOUBLE PRECISION PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGASUR(KDLON,8,2) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGBSUR(KDLON,8,2) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGATOP(KDLON,8,2) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGBTOP(KDLON,8,2) ! PADE APPROXIMANTS
 C
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
-      REAL*8 PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
-      REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
+      DOUBLE PRECISION PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
+      DOUBLE PRECISION PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
+      DOUBLE PRECISION PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
 C-----------------------------------------------------------------------
 C LOCAL VARIABLES:
-      REAL*8 ZADJD(KDLON,KFLEV+1)
-      REAL*8 ZADJU(KDLON,KFLEV+1)
-      REAL*8 ZDBDT(KDLON,Ninter,KFLEV)
-      REAL*8 ZDISD(KDLON,KFLEV+1)
-      REAL*8 ZDISU(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZADJD(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZADJU(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZDBDT(KDLON,Ninter,KFLEV)
+      DOUBLE PRECISION ZDISD(KDLON,KFLEV+1)
+      DOUBLE PRECISION ZDISU(KDLON,KFLEV+1)
 C
       INTEGER jk, jl
 C-----------------------------------------------------------------------