source: LMDZ6/branches/contrails/libf/phylmd/rrtm/srtm_srtm_224gp.F90 @ 5424

Last change on this file since 5424 was 5294, checked in by Laurent Fairhead, 7 weeks ago

Keeping clesphys.h was not the right solution
LF
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:keywords set to Author Date Id Revi
File size: 13.1 KB
Line 
1!
2! $Id: srtm_srtm_224gp.F90 5294 2024-10-29 18:35:00Z fhourdin $
3!
4SUBROUTINE SRTM_SRTM_224GP &
5 & ( KIDIA , KFDIA  , KLON  , KLEV  , KSW , KOVLP ,&
6 &   PAER  , PALBD  , PALBP , PAPH  , PAP ,&
7 &   PTS   , PTH    , PT    ,&
8 &   PQ    , PCCO2  , POZN  , PRMU0 ,&
9 &   PFRCL , PTAUC  , PASYC , POMGC ,&
10 &   PALBT , PFSUX  , PFSUC &
11 & ) 
12
13!-- interface to RRTM_SW
14!     JJMorcrette 030225
15
16USE PARKIND1  ,ONLY : JPIM     ,JPRB
17USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
18
19USE PARSRTM  , ONLY : JPLAY
20!USE YOERDI   , ONLY : RCH4   , RN2O   
21USE YOERAD   , ONLY : NAER
22USE YOESRTAER, ONLY : RSRTAUA, RSRPIZA, RSRASYA
23USE YOMPHY3  , ONLY : RII0
24USE YOMCST   , ONLY : RI0
25! Temporary fix waiting for cleaner interface (or not)
26USE clesphys_mod_h, ONLY: NSW, rch4, rn2o
27
28
29IMPLICIT NONE
30
31!!#include "clesphys.h"
32
33!-- Input arguments
34
35INTEGER(KIND=JPIM),INTENT(IN)    :: KLON
36INTEGER(KIND=JPIM)               :: KLEV! UNDETERMINED INTENT
37INTEGER(KIND=JPIM)               :: KSW! UNDETERMINED INTENT
38INTEGER(KIND=JPIM),INTENT(IN)    :: KIDIA
39INTEGER(KIND=JPIM),INTENT(IN)    :: KFDIA
40INTEGER(KIND=JPIM),INTENT(IN)    :: KOVLP
41REAL(KIND=JPRB)   ,INTENT(IN)    :: PAER(KLON,6,KLEV)    ! top to bottom
42REAL(KIND=JPRB)   ,INTENT(IN)    :: PALBD(KLON,KSW)
43REAL(KIND=JPRB)   ,INTENT(IN)    :: PALBP(KLON,KSW)
44REAL(KIND=JPRB)   ,INTENT(IN)    :: PAPH(KLON,KLEV+1)
45REAL(KIND=JPRB)   ,INTENT(IN)    :: PAP(KLON,KLEV)
46REAL(KIND=JPRB)   ,INTENT(IN)    :: PTS(KLON)
47REAL(KIND=JPRB)   ,INTENT(IN)    :: PTH(KLON,KLEV+1)
48REAL(KIND=JPRB)   ,INTENT(IN)    :: PT(KLON,KLEV)
49REAL(KIND=JPRB)   ,INTENT(IN)    :: PQ(KLON,KLEV)
50REAL(KIND=JPRB)   ,INTENT(IN)    :: PCCO2
51REAL(KIND=JPRB)   ,INTENT(IN)    :: POZN(KLON,KLEV)
52REAL(KIND=JPRB)   ,INTENT(IN)    :: PRMU0(KLON)
53REAL(KIND=JPRB)   ,INTENT(IN)    :: PFRCL(KLON,KLEV)     ! bottom to top
54REAL(KIND=JPRB)   ,INTENT(IN)    :: PTAUC(KLON,KSW,KLEV) ! bottom to top
55REAL(KIND=JPRB)   ,INTENT(IN)    :: PASYC(KLON,KSW,KLEV) ! bottom to top
56REAL(KIND=JPRB)   ,INTENT(IN)    :: POMGC(KLON,KSW,KLEV) ! bottom to top
57REAL(KIND=JPRB)                  :: PALBT(KLON,KSW) ! Argument NOT used
58REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSUX(KLON,2,KLEV+1)
59REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSUC(KLON,2,KLEV+1)
60!INTEGER_M :: KMOL, KCLDATM, KNFLAG, KCEFLAG, KIQFLAG, KSTR 
61
62!-- Output arguments
63
64!-----------------------------------------------------------------------
65
66!-- dummy integers
67
68INTEGER(KIND=JPIM) :: ICLDATM, INFLAG, ICEFLAG, I_LIQFLAG, I_NMOL, I_NSTR
69
70INTEGER(KIND=JPIM) :: IK, IMOL, J1, J2, JAE, JL, JK, JSW
71
72!-- dummy reals
73
74REAL(KIND=JPRB) :: Z_PZ(0:JPLAY)   , Z_TZ(0:JPLAY)   , Z_PAVEL(JPLAY)  , Z_TAVEL(JPLAY)
75REAL(KIND=JPRB) :: Z_COLDRY(JPLAY) , Z_COLMOL(JPLAY) , Z_WKL(35,JPLAY)
76REAL(KIND=JPRB) :: Z_CO2MULT(JPLAY), Z_COLCH4(JPLAY) , Z_COLCO2(JPLAY) , Z_COLH2O(JPLAY)
77REAL(KIND=JPRB) :: Z_COLN2O(JPLAY) , Z_COLO2(JPLAY)  , Z_COLO3(JPLAY)
78REAL(KIND=JPRB) :: Z_FORFAC(JPLAY) , Z_FORFRAC(JPLAY), Z_SELFFAC(JPLAY), Z_SELFFRAC(JPLAY)
79REAL(KIND=JPRB) :: Z_FAC00(JPLAY)  , Z_FAC01(JPLAY)  , Z_FAC10(JPLAY)  , Z_FAC11(JPLAY)
80REAL(KIND=JPRB) :: Z_TBOUND        , Z_ONEMINUS    , ZRMU0 , ZADJI0
81REAL(KIND=JPRB) :: ZALBD(KSW)    , ZALBP(KSW)    , ZFRCL(JPLAY)
82REAL(KIND=JPRB) :: ZTAUC(JPLAY,KSW), ZASYC(JPLAY,KSW), ZOMGC(JPLAY,KSW)
83REAL(KIND=JPRB) :: ZTAUA(JPLAY,KSW), ZASYA(JPLAY,KSW), ZOMGA(JPLAY,KSW)
84
85REAL(KIND=JPRB) :: ZBBCD(JPLAY+1), ZBBCU(JPLAY+1), ZBBFD(JPLAY+1), ZBBFU(JPLAY+1)
86REAL(KIND=JPRB) :: ZUVCD(JPLAY+1), ZUVCU(JPLAY+1), ZUVFD(JPLAY+1), ZUVFU(JPLAY+1)
87REAL(KIND=JPRB) :: ZVSCD(JPLAY+1), ZVSCU(JPLAY+1), ZVSFD(JPLAY+1), ZVSFU(JPLAY+1)
88REAL(KIND=JPRB) :: ZNICD(JPLAY+1), ZNICU(JPLAY+1), ZNIFD(JPLAY+1), ZNIFU(JPLAY+1)
89
90INTEGER(KIND=JPIM) :: I_LAYTROP, I_LAYSWTCH, I_LAYLOW
91INTEGER(KIND=JPIM) :: INDFOR(JPLAY), INDSELF(JPLAY)
92INTEGER(KIND=JPIM) :: JP(JPLAY), JT(JPLAY), JT1(JPLAY)
93
94REAL(KIND=JPRB) :: Z_AMD                  ! Effective molecular weight of dry air (g/mol)
95REAL(KIND=JPRB) :: Z_AMW                  ! Molecular weight of water vapor (g/mol)
96REAL(KIND=JPRB) :: Z_AMCO2                ! Molecular weight of carbon dioxide (g/mol)
97REAL(KIND=JPRB) :: Z_AMO                  ! Molecular weight of ozone (g/mol)
98REAL(KIND=JPRB) :: Z_AMCH4                ! Molecular weight of methane (g/mol)
99REAL(KIND=JPRB) :: Z_AMN2O                ! Molecular weight of nitrous oxide (g/mol)
100REAL(KIND=JPRB) :: Z_AMC11                ! Molecular weight of CFC11 (g/mol) - CFCL3
101REAL(KIND=JPRB) :: Z_AMC12                ! Molecular weight of CFC12 (g/mol) - CF2CL2
102REAL(KIND=JPRB) :: Z_AVGDRO               ! Avogadro's number (molecules/mole)
103REAL(KIND=JPRB) :: Z_GRAVIT               ! Gravitational acceleration (cm/sec2)
104REAL(KIND=JPRB) :: Z_AMM
105
106! Atomic weights for conversion from mass to volume mixing ratios; these
107!  are the same values used in ECRT to assure accurate conversion to vmr
108data Z_AMD   /  28.970_JPRB    /
109data Z_AMW   /  18.0154_JPRB   /
110data Z_AMCO2 /  44.011_JPRB    /
111data Z_AMO   /  47.9982_JPRB   /
112data Z_AMCH4 /  16.043_JPRB    /
113data Z_AMN2O /  44.013_JPRB    /
114data Z_AMC11 / 137.3686_JPRB   /
115data Z_AMC12 / 120.9140_JPRB   /
116data Z_AVGDRO/ 6.02214E23_JPRB /
117data Z_GRAVIT/ 9.80665E02_JPRB /
118
119REAL(KIND=JPRB) :: ZCLEAR, ZCLOUD, ZEPSEC, ZTOTCC
120
121INTEGER(KIND=JPIM) :: IOVLP
122REAL(KIND=JPRB) :: ZHOOK_HANDLE
123
124
125#include "srtm_setcoef.intfb.h"
126#include "srtm_spcvrt.intfb.h"
127
128
129!-----------------------------------------------------------------------
130!-- calculate information needed ny the radiative transfer routine
131
132IF (LHOOK) CALL DR_HOOK('SRTM_SRTM_224GP',0,ZHOOK_HANDLE)
133ZEPSEC  = 1.E-06_JPRB
134Z_ONEMINUS=1.0_JPRB -  ZEPSEC
135ZADJI0 = RII0 / RI0
136!-- overlap: 1=max-ran, 2=maximum, 3=random
137IOVLP=3
138
139!print *,'Entering srtm_srtm_224gp'
140
141ICLDATM = 1
142INFLAG    = 2
143ICEFLAG    = 3
144I_LIQFLAG = 1
145I_NMOL    = 6
146I_NSTR    = 2
147
148DO JL = KIDIA, KFDIA
149  ZRMU0=PRMU0(JL)
150  IF (ZRMU0 > 0.0_JPRB) THEN
151
152!- coefficients related to the cloud optical properties (original RRTM_SW)
153
154!  print *,'just before SRTM_CLDPROP'
155
156!  DO JK=1,KLEV
157!    CLDFRAC(JK) = PFRCL (JL,JK)
158!    CLDDAT1(JK) = PSCLA1(JL,JK)
159!    CLDDAT2(JK) = PSCLA2(JL,JK)
160!    CLDDAT3(JK) = PSCLA3(JL,JK)
161!    CLDDAT4(JK) = PSCLA4(JL,JK)
162!    DO JMOM=0,16
163!      CLDDATMOM(JMOM,JK)=PSCLMOM(JL,JMOM,JK)
164!    ENDDO
165!    print 9101,JK,CLDFRAC(JK),CLDDAT1(JK),CLDDAT2(JK),CLDDAT3(JK)&
166!    &,CLDDAT4(JK),(CLDDATMOM(JMOM,JK),JMOM=0,NSTR)
167    9101 format(1x,'srtm_srtm_224gp Cld :',I3,f7.4,7E12.5)
168!  ENDDO
169
170!  CALL SRTM_CLDPROP &
171!    &( KLEV, ICLDATM, INFLAG, ICEFLAG, LIQFLAG, NSTR &
172!    &, CLDFRAC, CLDDAT1, CLDDAT2, CLDDAT3, CLDDAT4, CLDDATMOM &
173!    &, TAUCLDORIG, TAUCLOUD, SSACLOUD, XMOM &
174!    &)
175
176!- coefficients for the temperature and pressure dependence of the
177! molecular absorption coefficients
178
179    DO J1=1,35
180      DO J2=1,KLEV
181        Z_WKL(J1,J2)=0.0_JPRB
182      ENDDO
183    ENDDO
184
185    Z_TBOUND=PTS(JL)
186    Z_PZ(0) = paph(JL,klev+1)/100._JPRB
187    Z_TZ(0) = pth (JL,klev+1)
188
189    ZCLEAR=1.0_JPRB
190    ZCLOUD=0.0_JPRB
191    ZTOTCC=0.0_JPRB
192    DO JK = 1, KLEV
193      Z_PAVEL(JK) = pap(JL,KLEV-JK+1) /100._JPRB
194      Z_TAVEL(JK) = pt (JL,KLEV-JK+1)
195      Z_PZ(JK)    = paph(JL,KLEV-JK+1)/100._JPRB
196      Z_TZ(JK)    = pth (JL,KLEV-JK+1)
197      Z_WKL(1,JK) = pq(JL,KLEV-JK+1)  *Z_AMD/Z_AMW
198      Z_WKL(2,JK) = pcco2             *Z_AMD/Z_AMCO2
199      Z_WKL(3,JK) = pozn(JL,KLEV-JK+1)*Z_AMD/Z_AMO
200      Z_WKL(4,JK) = rn2o              *Z_AMD/Z_AMN2O
201      Z_WKL(6,JK) = rch4              *Z_AMD/Z_AMCH4
202      Z_AMM = (1-Z_WKL(1,JK))*Z_AMD + Z_WKL(1,JK)*Z_AMW
203      Z_COLDRY(JK) = (Z_PZ(JK-1)-Z_PZ(JK))*1.E3_JPRB*Z_AVGDRO/(Z_GRAVIT*Z_AMM*(1+Z_WKL(1,JK)))
204!    print 9200,JK,PAVEL(JK),TAVEL(JK),(WKL(JA,JK),JA=1,4),WKL(6,JK),COLDRY(JK)
205      9200 format(1x,'SRTM ',I3,2F7.1,6E13.5)
206
207      IF (KOVLP == 1) THEN
208        ZCLEAR=ZCLEAR*(1.0_JPRB-MAX(PFRCL(JL,JK),ZCLOUD)) &
209         & /(1.0_JPRB-MIN(ZCLOUD,1.0_JPRB-ZEPSEC)) 
210        ZCLOUD=PFRCL(JL,JK)
211        ZTOTCC=1.0_JPRB-ZCLEAR
212      ELSEIF (KOVLP == 2) THEN
213        ZCLOUD=MAX(ZCLOUD,PFRCL(JL,JK))
214        ZCLEAR=1.0_JPRB-ZCLOUD
215        ZTOTCC=ZCLOUD
216      ELSEIF (KOVLP == 3) THEN
217        ZCLEAR=ZCLEAR*(1.0_JPRB-PFRCL(JL,JK))
218        ZCLOUD=1.0_JPRB-ZCLEAR
219        ZTOTCC=ZCLOUD
220      ENDIF
221
222    ENDDO
223
224!  print *,'ZTOTCC ZCLEAR : ',ZTOTCC,' ',ZCLEAR
225
226    DO IMOL=1,I_NMOL
227      DO JK=1,KLEV
228        Z_WKL(IMOL,JK)=Z_COLDRY(JK)* Z_WKL(IMOL,JK)
229      ENDDO
230    ENDDO
231
232!    IF (ZTOTCC == 0.0_JPRB) THEN
233!      DO JK=1,KLEV
234!        ZFRCL(JK)=0.0_JPRB   
235!      ENDDO
236!    ELSE
237!      DO JK=1,KLEV
238!        ZFRCL(JK)=PFRCL(JL,JK)/ZTOTCC
239!      ENDDO
240!    ENDIF
241
242!  print *,'just before SRTM_SETCOEF'
243
244    ZFRCL(1:KLEV)=PFRCL(JL,1:KLEV)
245    ZCLEAR=0._JPRB
246    ZCLOUD=1._JPRB
247
248    CALL SRTM_SETCOEF &
249     & ( KLEV   , I_NMOL,&
250     & Z_PAVEL  , Z_TAVEL   , Z_PZ     , Z_TZ     , Z_TBOUND,&
251     & Z_COLDRY , Z_WKL,&
252     & I_LAYTROP, I_LAYSWTCH, I_LAYLOW,&
253     & Z_CO2MULT, Z_COLCH4  , Z_COLCO2 , Z_COLH2O , Z_COLMOL  , Z_COLN2O  , Z_COLO2 , Z_COLO3,&
254     & Z_FORFAC , Z_FORFRAC , INDFOR , Z_SELFFAC, Z_SELFFRAC, INDSELF,&
255     & Z_FAC00  , Z_FAC01   , Z_FAC10  , Z_FAC11,&
256     & JP     , JT      , JT1     &
257     & ) 
258 
259!  print *,'just after SRTM_SETCOEF'
260
261!- call the radiation transfer routine
262 
263    DO JSW=1,KSW
264      ZALBD(JSW)=PALBD(JL,JSW)
265      ZALBP(JSW)=PALBP(JL,JSW)
266      DO JK=1,KLEV
267        ZTAUC(JK,JSW) = PTAUC(JL,JSW,JK)
268        ZASYC(JK,JSW) = PASYC(JL,JSW,JK)
269        ZOMGC(JK,JSW) = POMGC(JL,JSW,JK)
270!      print 9002,JSW,JK,ZFRCL(JK),ZTAUC(JK,JSW),ZASYC(JK,JSW),ZOMGC(JK,JSW)
271        9002  format(1x,'srtm_224gp ClOPropECmodel ',2I3,f8.4,3E12.5)
272      ENDDO
273    ENDDO
274
275!- mixing of aerosols
276 
277!  print *,'Aerosol optical properties computations'
278!  DO JSW=1,KSW
279!    print 9012,JSW,(JAE,RSRTAUA(JSW,JAE),RSRPIZA(JSW,JAE),RSRASYA(JSW,JAE),JAE=1,6)
280    9012 format(I3,(/,I3,3E13.5))
281!  ENDDO
282
283!  DO JK=1,KLEV
284!    print 9013,JK,(PAER(JL,JAE,JK),JAE=1,6)
285    9013 format(1x,I3,6E12.5)
286!  ENDDO
287
288    IF (NAER == 0) THEN
289      DO JSW=1,KSW
290        DO JK=1,KLEV
291          ZTAUA(JK,JSW)= 0.0_JPRB
292          ZASYA(JK,JSW)= 0.0_JPRB
293          ZOMGA(JK,JSW)= 1.0_JPRB
294        ENDDO
295      ENDDO
296    ELSE
297      DO JSW=1,KSW
298        DO JK=1,KLEV
299          IK=KLEV+1-JK
300          ZTAUA(JK,JSW)=0.0_JPRB
301          ZASYA(JK,JSW)=0.0_JPRB
302          ZOMGA(JK,JSW)=0.0_JPRB
303          DO JAE=1,6
304            ZTAUA(JK,JSW)=ZTAUA(JK,JSW)+RSRTAUA(JSW,JAE)*PAER(JL,JAE,IK)
305            ZOMGA(JK,JSW)=ZOMGA(JK,JSW)+RSRTAUA(JSW,JAE)*PAER(JL,JAE,IK) &
306             & *RSRPIZA(JSW,JAE) 
307            ZASYA(JK,JSW)=ZASYA(JK,JSW)+RSRTAUA(JSW,JAE)*PAER(JL,JAE,IK) &
308             & *RSRPIZA(JSW,JAE)*RSRASYA(JSW,JAE) 
309          ENDDO
310          IF (ZOMGA(JK,JSW) /= 0.0_JPRB) THEN
311            ZASYA(JK,JSW)=ZASYA(JK,JSW)/ZOMGA(JK,JSW)
312          ENDIF
313          IF (ZTAUA(JK,JSW) /= 0.0_JPRB) THEN
314            ZOMGA(JK,JSW)=ZOMGA(JK,JSW)/ZTAUA(JK,JSW)
315          ENDIF
316!      print 9003,JSW,JK,ZTAUA(JK,JSW),ZOMGA(JK,JSW),ZASYA(JK,JSW)
3179003  format(1x,'Aerosols ',2I3,3F10.4)
318        ENDDO
319      ENDDO
320    ENDIF
321
322    DO JK=1,KLEV+1
323      ZBBCU(JK)=0.0_JPRB
324      ZBBCD(JK)=0.0_JPRB
325      ZBBFU(JK)=0.0_JPRB
326      ZBBFD(JK)=0.0_JPRB
327      ZUVCU(JK)=0.0_JPRB
328      ZUVCD(JK)=0.0_JPRB
329      ZUVFU(JK)=0.0_JPRB
330      ZUVFD(JK)=0.0_JPRB
331      ZVSCU(JK)=0.0_JPRB
332      ZVSCD(JK)=0.0_JPRB
333      ZVSFU(JK)=0.0_JPRB
334      ZVSFD(JK)=0.0_JPRB
335      ZNICU(JK)=0.0_JPRB
336      ZNICD(JK)=0.0_JPRB
337      ZNIFU(JK)=0.0_JPRB
338      ZNIFD(JK)=0.0_JPRB
339    ENDDO
340
341!  print *,'just before calling STRM_SPCVRT for JL=',JL,' and ZRMU0=',ZRMU0
342
343    CALL SRTM_SPCVRT &
344     & ( KLEV   , I_NMOL    , KSW    , Z_ONEMINUS,&
345     & Z_PAVEL  , Z_TAVEL   , Z_PZ     , Z_TZ     , Z_TBOUND  , ZALBD   , ZALBP,&
346     & ZFRCL  , ZTAUC   , ZASYC  , ZOMGC  , ZTAUA   , ZASYA   , ZOMGA , ZRMU0,&
347     & Z_COLDRY , Z_WKL,&
348     & I_LAYTROP, I_LAYSWTCH, I_LAYLOW,&
349     & Z_CO2MULT, Z_COLCH4  , Z_COLCO2 , Z_COLH2O , Z_COLMOL  , Z_COLN2O  , Z_COLO2 , Z_COLO3,&
350     & Z_FORFAC , Z_FORFRAC , INDFOR , Z_SELFFAC, Z_SELFFRAC, INDSELF,&
351     & Z_FAC00  , Z_FAC01   , Z_FAC10  , Z_FAC11,&
352     & JP     , JT      , JT1,&
353     & ZBBFD  , ZBBFU   , ZUVFD  , ZUVFU  , ZVSFD   , ZVSFU   , ZNIFD , ZNIFU,&
354     & ZBBCD  , ZBBCU   , ZUVCD  , ZUVCU  , ZVSCD   , ZVSCU   , ZNICD , ZNICU &
355     & ) 
356
357!  print *,'SRTM_SRTM_224GP before potential scaling'
358!    IF (IOVLP == 3) THEN
359!      DO JK=1,KLEV+1
360!!      print 9004,JK,ZBBCU(JK),ZBBCD(JK),ZBBFU(JK),ZBBFD(JK)
361        9004 format(1x,'Clear-sky and total fluxes U & D ',I3,4F10.3)
362!        PFSUC(JL,1,JK)=ZBBCU(JK)
363!        PFSUC(JL,2,JK)=ZBBCD(JK)
364!        PFSUX(JL,1,JK)=ZBBFU(JK)
365!        PFSUX(JL,2,JK)=ZBBFD(JK)
366!      ENDDO
367!    ELSE
368!    print *,'SRTM_SRTM_224GP after potential scaling'
369      DO JK=1,KLEV+1
370        PFSUC(JL,1,JK)=ZADJI0 * ZBBCU(JK)
371        PFSUC(JL,2,JK)=ZADJI0 * ZBBCD(JK)
372        PFSUX(JL,1,JK)=ZADJI0 * ( (1.0_JPRB-ZCLEAR)*ZBBFU(JK)+ZCLEAR*ZBBCU(JK) )
373        PFSUX(JL,2,JK)=ZADJI0 * ( (1.0_JPRB-ZCLEAR)*ZBBFD(JK)+ZCLEAR*ZBBCD(JK) )
374      ENDDO
375!    ENDIF
376
377!  DO JK=1,KLEV+1
378!    print 9005,JK,PFSUC(JL,1,JK),PFSUC(JL,2,JK),PFSUX(JL,1,JK),PFSUX(JL,2,JK)
379    9005 format(1x,'Clear-sky and total fluxes U & D ',I3,4F10.3)
380!  ENDDO
381 
382  ELSE
383    DO JK=1,KLEV+1
384      PFSUC(JL,1,JK)=0.0_JPRB
385      PFSUC(JL,2,JK)=0.0_JPRB
386      PFSUX(JL,1,JK)=0.0_JPRB
387      PFSUX(JL,2,JK)=0.0_JPRB
388    ENDDO
389  ENDIF
390ENDDO
391
392!PRINT *,'OUT OF SRTM_224GP'
393
394!-----------------------------------------------------------------------
395IF (LHOOK) CALL DR_HOOK('SRTM_SRTM_224GP',1,ZHOOK_HANDLE)
396END SUBROUTINE SRTM_SRTM_224GP
397
Note: See TracBrowser for help on using the repository browser.