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., |