--- trunk/libf/phylmd/Radlwsw/lwvn.f 2010/03/03 13:23:49 24 +++ trunk/libf/phylmd/Radlwsw/lwvn.f 2013/07/08 18:12:18 71 @@ -38,32 +38,32 @@ C INTEGER KUAER,KTRAER C - REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS - REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-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 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS - REAL*8 PADJU(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS - REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX - REAL*8 PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT + DOUBLE PRECISION PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS + DOUBLE PRECISION PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-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 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS + DOUBLE PRECISION PADJU(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS + DOUBLE PRECISION PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX + DOUBLE PRECISION PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT C C* LOCAL ARRAYS: C - REAL*8 ZGLAYD(KDLON) - REAL*8 ZGLAYU(KDLON) - REAL*8 ZTT(KDLON,NTRA) - REAL*8 ZTT1(KDLON,NTRA) - REAL*8 ZTT2(KDLON,NTRA) - REAL*8 ZUU(KDLON,NUA) + DOUBLE PRECISION ZGLAYD(KDLON) + DOUBLE PRECISION ZGLAYU(KDLON) + DOUBLE PRECISION ZTT(KDLON,NTRA) + DOUBLE PRECISION ZTT1(KDLON,NTRA) + DOUBLE PRECISION ZTT2(KDLON,NTRA) + DOUBLE PRECISION ZUU(KDLON,NUA) C INTEGER jk, jl, ja, im12, ind, inu, ixu, jg INTEGER ixd, ibs, idd, imu, jk1, jk2, jnu - REAL*8 zwtr + DOUBLE PRECISION zwtr c C* Data Block: c - REAL*8 WG1(2) + DOUBLE PRECISION WG1(2) SAVE WG1 DATA (WG1(jk),jk=1,2) /1.0, 1.0/ C-----------------------------------------------------------------------