--- trunk/libf/phylmd/Radlwsw/lwvd.f	2010/03/03 13:23:49	24
+++ trunk/phylmd/Radlwsw/lwvd.f	2013/11/15 18:45:49	76
@@ -37,27 +37,27 @@
 C
       INTEGER KUAER,KTRAER
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
-      REAL*8 PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-C
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! ENERGY EXCHANGE MATRIX
-      REAL*8 PDISD(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
-      REAL*8 PDISU(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
+      DOUBLE PRECISION PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
+      DOUBLE PRECISION PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT
+      DOUBLE PRECISION PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+C
+      DOUBLE PRECISION PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! ENERGY EXCHANGE MATRIX
+      DOUBLE PRECISION PDISD(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
+      DOUBLE PRECISION PDISU(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZGLAYD(KDLON)
-      REAL*8 ZGLAYU(KDLON)
-      REAL*8 ZTT(KDLON,NTRA)
-      REAL*8 ZTT1(KDLON,NTRA)
-      REAL*8 ZTT2(KDLON,NTRA)
+      DOUBLE PRECISION ZGLAYD(KDLON)
+      DOUBLE PRECISION ZGLAYU(KDLON)
+      DOUBLE PRECISION ZTT(KDLON,NTRA)
+      DOUBLE PRECISION ZTT1(KDLON,NTRA)
+      DOUBLE PRECISION ZTT2(KDLON,NTRA)
 C
       INTEGER jl, jk, ja, ikp1, ikn, ikd1, jkj, ikd2
       INTEGER ikjp1, ikm1, ikj, jlk, iku1, ijkl, iku2
       INTEGER ind1, ind2, ind3, ind4, itt
-      REAL*8 zww, zdzxdg, zdzxmg
+      DOUBLE PRECISION zww, zdzxdg, zdzxmg
 C
 C*         1.    INITIALIZATION
 C                --------------