/[lmdze]/trunk/libf/phylmd/Radlwsw/lwb.f
ViewVC logotype

Diff of /trunk/libf/phylmd/Radlwsw/lwb.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 70 by guez, Wed Mar 3 13:23:49 2010 UTC revision 71 by guez, Mon Jul 8 18:12:18 2013 UTC
# Line 57  C--------------------------------------- Line 57  C---------------------------------------
57  C  C
58  C ARGUMENTS:  C ARGUMENTS:
59  C  C
60        REAL*8 PDT0(KDLON)        DOUBLE PRECISION PDT0(KDLON)
61        REAL*8 PTAVE(KDLON,KFLEV)        DOUBLE PRECISION PTAVE(KDLON,KFLEV)
62        REAL*8 PTL(KDLON,KFLEV+1)        DOUBLE PRECISION PTL(KDLON,KFLEV+1)
63  C  C
64        REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION        DOUBLE PRECISION PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION
65        REAL*8 PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION        DOUBLE PRECISION PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION
66        REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION        DOUBLE PRECISION PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
67        REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION        DOUBLE PRECISION PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
68        REAL*8 PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION        DOUBLE PRECISION PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION
69        REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT        DOUBLE PRECISION PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
70        REAL*8 PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS        DOUBLE PRECISION PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
71        REAL*8 PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS        DOUBLE PRECISION PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
72        REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS        DOUBLE PRECISION PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
73        REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS        DOUBLE PRECISION PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
74        REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS        DOUBLE PRECISION PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
75        REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS        DOUBLE PRECISION PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
76  C  C
77  C-------------------------------------------------------------------------  C-------------------------------------------------------------------------
78  C*  LOCAL VARIABLES:  C*  LOCAL VARIABLES:
79        INTEGER INDB(KDLON),INDS(KDLON)        INTEGER INDB(KDLON),INDS(KDLON)
80        REAL*8 ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)        DOUBLE PRECISION ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
81        REAL*8 ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)        DOUBLE PRECISION ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)
82  c  c
83        INTEGER jk, jl, ic, jnu, jf, jg        INTEGER jk, jl, ic, jnu, jf, jg
84        INTEGER jk1, jk2        INTEGER jk1, jk2
85        INTEGER k, j, ixtox, indto, ixtx, indt        INTEGER k, j, ixtox, indto, ixtx, indt
86        INTEGER indsu, indtp        INTEGER indsu, indtp
87        REAL*8 zdsto1, zdstox, zdst1, zdstx        DOUBLE PRECISION zdsto1, zdstox, zdst1, zdstx
88  c  c
89  C* Quelques parametres:  C* Quelques parametres:
90        REAL*8 TSTAND        DOUBLE PRECISION TSTAND
91        PARAMETER (TSTAND=250.0)        PARAMETER (TSTAND=250.0)
92        REAL*8 TSTP        DOUBLE PRECISION TSTP
93        PARAMETER (TSTP=12.5)        PARAMETER (TSTP=12.5)
94        INTEGER MXIXT        INTEGER MXIXT
95        PARAMETER (MXIXT=10)        PARAMETER (MXIXT=10)
96  C  C
97  C* Used Data Block:  C* Used Data Block:
98        REAL*8 TINTP(11)        DOUBLE PRECISION TINTP(11)
99        SAVE TINTP        SAVE TINTP
100        REAL*8 GA(11,16,3), GB(11,16,3)        DOUBLE PRECISION GA(11,16,3), GB(11,16,3)
101        SAVE GA, GB        SAVE GA, GB
102        REAL*8 XP(6,6)        DOUBLE PRECISION XP(6,6)
103        SAVE XP        SAVE XP
104  c  c
105        DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,        DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
Legend:
Removed from v.70  
changed lines
  Added in v.71
  ViewVC Help
Powered by ViewVC 1.1.21