--- trunk/libf/phylmd/Radlwsw/lwvb.f	2010/03/03 13:23:49	24
+++ trunk/libf/phylmd/Radlwsw/lwvb.f	2013/07/08 18:12:18	71
@@ -46,45 +46,45 @@
 C
       INTEGER KUAER,KTRAER, KLIM
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
-      REAL*8 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
-      REAL*8 PADJU(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
-      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) ! SPECTRAL SURFACE PLANCK FUNCTION
-      REAL*8 PBSUI(KDLON) ! SURFACE PLANCK FUNCTION
-      REAL*8 PBTOP(KDLON,Ninter) ! SPECTRAL T.O.A. PLANCK FUNCTION
-      REAL*8 PDISD(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
-      REAL*8 PDISU(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
-      REAL*8 PEMIS(KDLON) ! SURFACE EMISSIVITY
-      REAL*8 PPMB(KDLON,KFLEV+1) ! PRESSURE MB
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
-      REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
-      REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
-      REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
+      DOUBLE PRECISION PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
+      DOUBLE PRECISION PADJD(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
+      DOUBLE PRECISION PADJU(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
+      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) ! SPECTRAL SURFACE PLANCK FUNCTION
+      DOUBLE PRECISION PBSUI(KDLON) ! SURFACE PLANCK FUNCTION
+      DOUBLE PRECISION PBTOP(KDLON,Ninter) ! SPECTRAL T.O.A. PLANCK FUNCTION
+      DOUBLE PRECISION PDISD(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
+      DOUBLE PRECISION PDISU(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
+      DOUBLE PRECISION PEMIS(KDLON) ! SURFACE EMISSIVITY
+      DOUBLE PRECISION PPMB(KDLON,KFLEV+1) ! PRESSURE MB
+      DOUBLE PRECISION PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      DOUBLE PRECISION PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
+      DOUBLE PRECISION PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
+      DOUBLE PRECISION PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
+      DOUBLE PRECISION PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
 C
-      REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
-      REAL*8 PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
+      DOUBLE PRECISION PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
+      DOUBLE PRECISION PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZBGND(KDLON)
-      REAL*8 ZFD(KDLON)
-      REAL*8  ZFN10(KDLON)
-      REAL*8 ZFU(KDLON)
-      REAL*8  ZTT(KDLON,NTRA)
-      REAL*8 ZTT1(KDLON,NTRA)
-      REAL*8 ZTT2(KDLON,NTRA)
-      REAL*8  ZUU(KDLON,NUA) 
-      REAL*8 ZCNSOL(KDLON)
-      REAL*8 ZCNTOP(KDLON)
+      DOUBLE PRECISION ZBGND(KDLON)
+      DOUBLE PRECISION ZFD(KDLON)
+      DOUBLE PRECISION  ZFN10(KDLON)
+      DOUBLE PRECISION ZFU(KDLON)
+      DOUBLE PRECISION  ZTT(KDLON,NTRA)
+      DOUBLE PRECISION ZTT1(KDLON,NTRA)
+      DOUBLE PRECISION ZTT2(KDLON,NTRA)
+      DOUBLE PRECISION  ZUU(KDLON,NUA) 
+      DOUBLE PRECISION ZCNSOL(KDLON)
+      DOUBLE PRECISION ZCNTOP(KDLON)
 C
       INTEGER jk, jl, ja
       INTEGER jstra, jstru
       INTEGER ind1, ind2, ind3, ind4, in, jlim
-      REAL*8 zctstr
+      DOUBLE PRECISION zctstr
 C-----------------------------------------------------------------------
 C
 C*         1.    INITIALIZATION