--- trunk/libf/phylmd/Radlwsw/lwc.f 2010/03/03 13:23:49 24 +++ trunk/phylmd/Radlwsw/lwc.f 2013/11/15 18:45:49 76 @@ -54,30 +54,31 @@ C----------------------------------------------------------------------- C* ARGUMENTS: INTEGER klim - REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES - REAL*8 PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION - REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION - REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) !CLEAR-SKY ENERGY EXCHANGE - REAL*8 PCTS(KDLON,KFLEV) ! CLEAR-SKY LAYER COOLING-TO-SPACE + DOUBLE PRECISION PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES + DOUBLE PRECISION PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION + DOUBLE PRECISION PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION + DOUBLE PRECISION PCNTRB(KDLON,KFLEV+1,KFLEV+1) !CLEAR-SKY ENERGY EXCHANGE + DOUBLE PRECISION PCTS(KDLON,KFLEV) ! CLEAR-SKY LAYER COOLING-TO-SPACE c - REAL*8 PCLDLD(KDLON,KFLEV) - REAL*8 PCLDLU(KDLON,KFLEV) - REAL*8 PEMIS(KDLON) + DOUBLE PRECISION PCLDLD(KDLON,KFLEV) + DOUBLE PRECISION PCLDLU(KDLON,KFLEV) + DOUBLE PRECISION PEMIS(KDLON) C - REAL*8 PFLUX(KDLON,2,KFLEV+1) + DOUBLE PRECISION PFLUX(KDLON,2,KFLEV+1) C----------------------------------------------------------------------- C* LOCAL VARIABLES: INTEGER IMX(KDLON), IMXP(KDLON) C - REAL*8 ZCLEAR(KDLON),ZCLOUD(KDLON),ZDNF(KDLON,KFLEV+1,KFLEV+1) + DOUBLE PRECISION ZCLEAR(KDLON),ZCLOUD(KDLON) + DOUBLE PRECISION ZDNF(KDLON,KFLEV+1,KFLEV+1) S , ZFD(KDLON), ZFN10(KDLON), ZFU(KDLON) S , ZUPF(KDLON,KFLEV+1,KFLEV+1) - REAL*8 ZCLM(KDLON,KFLEV+1,KFLEV+1) + DOUBLE PRECISION ZCLM(KDLON,KFLEV+1,KFLEV+1) C INTEGER jk, jl, imaxc, imx1, imx2, jkj, jkp1, jkm1 INTEGER jk1, jk2, jkc, jkcp1, jcloud INTEGER imxm1, imxp1 - REAL*8 zcfrac + DOUBLE PRECISION zcfrac C ------------------------------------------------------------------ C C* 1. INITIALIZATION