source: LMDZ5/branches/testing/libf/phylmd/radlwsw_m.F90 @ 1750

Last change on this file since 1750 was 1707, checked in by Laurent Fairhead, 12 years ago

Version testing basée sur la r1706

Testing release based on r1706
File size: 18.4 KB
RevLine 
[1687]1module radlwsw_m
2
3  IMPLICIT NONE
4
5contains
6
7SUBROUTINE radlwsw( &
8   dist, rmu0, fract, &
9   paprs, pplay,tsol,alb1, alb2, &
10   t,q,wo,&
11   cldfra, cldemi, cldtaupd,&
12   ok_ade, ok_aie, flag_aerosol,&
13   tau_aero, piz_aero, cg_aero,&
14   cldtaupi, new_aod, &
15   qsat, flwc, fiwc, &
16   heat,heat0,cool,cool0,radsol,albpla,&
17   topsw,toplw,solsw,sollw,&
18   sollwdown,&
19   topsw0,toplw0,solsw0,sollw0,&
20   lwdn0, lwdn, lwup0, lwup,&
21   swdn0, swdn, swup0, swup,&
22   topswad_aero, solswad_aero,&
23   topswai_aero, solswai_aero, &
24   topswad0_aero, solswad0_aero,&
25   topsw_aero, topsw0_aero,&
26   solsw_aero, solsw0_aero, &
27   topswcf_aero, solswcf_aero)
28
29
30
31  USE DIMPHY
32  USE assert_m, ONLY : assert
33  USE infotrac, ONLY : type_trac
34#ifdef REPROBUS
35  USE CHEM_REP, ONLY : solaireTIME, ok_SUNTIME, ndimozon
36#endif
37
38  !======================================================================
39  ! Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
40  ! Objet: interface entre le modele et les rayonnements
41  ! Arguments:
42  ! dist-----input-R- distance astronomique terre-soleil
43  ! rmu0-----input-R- cosinus de l'angle zenithal
44  ! fract----input-R- duree d'ensoleillement normalisee
45  ! co2_ppm--input-R- concentration du gaz carbonique (en ppm)
46  ! paprs----input-R- pression a inter-couche (Pa)
47  ! pplay----input-R- pression au milieu de couche (Pa)
48  ! tsol-----input-R- temperature du sol (en K)
49  ! alb1-----input-R- albedo du sol(entre 0 et 1) dans l'interval visible
50  ! alb2-----input-R- albedo du sol(entre 0 et 1) dans l'interval proche infra-rouge   
51  ! t--------input-R- temperature (K)
52  ! q--------input-R- vapeur d'eau (en kg/kg)
53  ! cldfra---input-R- fraction nuageuse (entre 0 et 1)
54  ! cldtaupd---input-R- epaisseur optique des nuages dans le visible (present-day value)
55  ! cldemi---input-R- emissivite des nuages dans l'IR (entre 0 et 1)
56  ! ok_ade---input-L- apply the Aerosol Direct Effect or not?
57  ! ok_aie---input-L- apply the Aerosol Indirect Effect or not?
58  ! flag_aerosol-input-I- aerosol flag from 0 to 6
59  ! tau_ae, piz_ae, cg_ae-input-R- aerosol optical properties (calculated in aeropt.F)
60  ! cldtaupi-input-R- epaisseur optique des nuages dans le visible
61  !                   calculated for pre-industrial (pi) aerosol concentrations, i.e. with smaller
62  !                   droplet concentration, thus larger droplets, thus generally cdltaupi cldtaupd
63  !                   it is needed for the diagnostics of the aerosol indirect radiative forcing     
64  !
65  ! heat-----output-R- echauffement atmospherique (visible) (K/jour)
66  ! cool-----output-R- refroidissement dans l'IR (K/jour)
67  ! radsol---output-R- bilan radiatif net au sol (W/m**2) (+ vers le bas)
68  ! albpla---output-R- albedo planetaire (entre 0 et 1)
69  ! topsw----output-R- flux solaire net au sommet de l'atm.
70  ! toplw----output-R- ray. IR montant au sommet de l'atmosphere
71  ! solsw----output-R- flux solaire net a la surface
72  ! sollw----output-R- ray. IR montant a la surface
73  ! solswad---output-R- ray. solaire net absorbe a la surface (aerosol dir)
74  ! topswad---output-R- ray. solaire absorbe au sommet de l'atm. (aerosol dir)
75  ! solswai---output-R- ray. solaire net absorbe a la surface (aerosol ind)
76  ! topswai---output-R- ray. solaire absorbe au sommet de l'atm. (aerosol ind)
77  !
78  ! ATTENTION: swai and swad have to be interpreted in the following manner:
79  ! ---------
80  ! ok_ade=F & ok_aie=F -both are zero
81  ! ok_ade=T & ok_aie=F -aerosol direct forcing is F_{AD} = topsw-topswad
82  !                        indirect is zero
83  ! ok_ade=F & ok_aie=T -aerosol indirect forcing is F_{AI} = topsw-topswai
84  !                        direct is zero
85  ! ok_ade=T & ok_aie=T -aerosol indirect forcing is F_{AI} = topsw-topswai
86  !                        aerosol direct forcing is F_{AD} = topswai-topswad
87  !
88 
89  !======================================================================
90 
91  ! ====================================================================
92  ! Adapte au modele de chimie INCA par Celine Deandreis & Anne Cozic -- 2009
93  ! 1 = ZERO   
94  ! 2 = AER total   
95  ! 3 = NAT   
96  ! 4 = BC   
97  ! 5 = SO4   
98  ! 6 = POM   
99  ! 7 = DUST   
100  ! 8 = SS   
101  ! 9 = NO3   
102  !
103  ! ====================================================================
104  include "YOETHF.h"
105  include "YOMCST.h"
106  include "clesphys.h"
107  include "iniprint.h"
108
109! Input arguments
110  REAL,    INTENT(in)  :: dist
111  REAL,    INTENT(in)  :: rmu0(KLON), fract(KLON)
112  REAL,    INTENT(in)  :: paprs(KLON,KLEV+1), pplay(KLON,KLEV)
113  REAL,    INTENT(in)  :: alb1(KLON), alb2(KLON), tsol(KLON)
114  REAL,    INTENT(in)  :: t(KLON,KLEV), q(KLON,KLEV)
115
116  REAL, INTENT(in):: wo(:, :, :) ! dimension(KLON,KLEV, 1 or 2)
117  ! column-density of ozone in a layer, in kilo-Dobsons
118  ! "wo(:, :, 1)" is for the average day-night field,
119  ! "wo(:, :, 2)" is for daylight time.
120
121  LOGICAL, INTENT(in)  :: ok_ade, ok_aie                                 ! switches whether to use aerosol direct (indirect) effects or not
122  INTEGER, INTENT(in)  :: flag_aerosol                                   ! takes value 0 (no aerosol) or 1 to 6 (aerosols)
123  REAL,    INTENT(in)  :: cldfra(KLON,KLEV), cldemi(KLON,KLEV), cldtaupd(KLON,KLEV)
124  REAL,    INTENT(in)  :: tau_aero(KLON,KLEV,9,2)                        ! aerosol optical properties (see aeropt.F)
125  REAL,    INTENT(in)  :: piz_aero(KLON,KLEV,9,2)                        ! aerosol optical properties (see aeropt.F)
126  REAL,    INTENT(in)  :: cg_aero(KLON,KLEV,9,2)                         ! aerosol optical properties (see aeropt.F)
127  REAL,    INTENT(in)  :: cldtaupi(KLON,KLEV)                            ! cloud optical thickness for pre-industrial aerosol concentrations
128  LOGICAL, INTENT(in)  :: new_aod                                        ! flag pour retrouver les resultats exacts de l'AR4 dans le cas ou l'on ne travaille qu'avec les sulfates
129  REAL,    INTENT(in)  :: qsat(klon,klev) ! Variable pour iflag_rrtm=1
130  REAL,    INTENT(in)  :: flwc(klon,klev) ! Variable pour iflag_rrtm=1
131  REAL,    INTENT(in)  :: fiwc(klon,klev) ! Variable pour iflag_rrtm=1
132
133! Output arguments
134  REAL,    INTENT(out) :: heat(KLON,KLEV), cool(KLON,KLEV)
135  REAL,    INTENT(out) :: heat0(KLON,KLEV), cool0(KLON,KLEV)
136  REAL,    INTENT(out) :: radsol(KLON), topsw(KLON), toplw(KLON)
137  REAL,    INTENT(out) :: solsw(KLON), sollw(KLON), albpla(KLON)
138  REAL,    INTENT(out) :: topsw0(KLON), toplw0(KLON), solsw0(KLON), sollw0(KLON)
139  REAL,    INTENT(out) :: sollwdown(KLON)
140  REAL,    INTENT(out) :: swdn(KLON,kflev+1),swdn0(KLON,kflev+1)
141  REAL,    INTENT(out) :: swup(KLON,kflev+1),swup0(KLON,kflev+1)
142  REAL,    INTENT(out) :: lwdn(KLON,kflev+1),lwdn0(KLON,kflev+1)
143  REAL,    INTENT(out) :: lwup(KLON,kflev+1),lwup0(KLON,kflev+1)
144  REAL,    INTENT(out) :: topswad_aero(KLON), solswad_aero(KLON)         ! output: aerosol direct forcing at TOA and surface
145  REAL,    INTENT(out) :: topswai_aero(KLON), solswai_aero(KLON)         ! output: aerosol indirect forcing atTOA and surface
146  REAL, DIMENSION(klon), INTENT(out)    :: topswad0_aero
147  REAL, DIMENSION(klon), INTENT(out)    :: solswad0_aero
148  REAL, DIMENSION(kdlon,9), INTENT(out) :: topsw_aero
149  REAL, DIMENSION(kdlon,9), INTENT(out) :: topsw0_aero
150  REAL, DIMENSION(kdlon,9), INTENT(out) :: solsw_aero
151  REAL, DIMENSION(kdlon,9), INTENT(out) :: solsw0_aero
152  REAL, DIMENSION(kdlon,3), INTENT(out) :: topswcf_aero
153  REAL, DIMENSION(kdlon,3), INTENT(out) :: solswcf_aero
154
155! Local variables
156  REAL(KIND=8) ZFSUP(KDLON,KFLEV+1)
157  REAL(KIND=8) ZFSDN(KDLON,KFLEV+1)
158  REAL(KIND=8) ZFSUP0(KDLON,KFLEV+1)
159  REAL(KIND=8) ZFSDN0(KDLON,KFLEV+1)
160  REAL(KIND=8) ZFLUP(KDLON,KFLEV+1)
161  REAL(KIND=8) ZFLDN(KDLON,KFLEV+1)
162  REAL(KIND=8) ZFLUP0(KDLON,KFLEV+1)
163  REAL(KIND=8) ZFLDN0(KDLON,KFLEV+1)
164  REAL(KIND=8) zx_alpha1, zx_alpha2
165  INTEGER k, kk, i, j, iof, nb_gr
166  REAL(KIND=8) PSCT
167  REAL(KIND=8) PALBD(kdlon,2), PALBP(kdlon,2)
168  REAL(KIND=8) PEMIS(kdlon), PDT0(kdlon), PVIEW(kdlon)
169  REAL(KIND=8) PPSOL(kdlon), PDP(kdlon,KLEV)
170  REAL(KIND=8) PTL(kdlon,kflev+1), PPMB(kdlon,kflev+1)
171  REAL(KIND=8) PTAVE(kdlon,kflev)
172  REAL(KIND=8) PWV(kdlon,kflev), PQS(kdlon,kflev)
173
174  real(kind=8) POZON(kdlon, kflev, size(wo, 3)) ! mass fraction of ozone
175  ! "POZON(:, :, 1)" is for the average day-night field,
176  ! "POZON(:, :, 2)" is for daylight time.
177
178  REAL(KIND=8) PAER(kdlon,kflev,5)
179  REAL(KIND=8) PCLDLD(kdlon,kflev)
180  REAL(KIND=8) PCLDLU(kdlon,kflev)
181  REAL(KIND=8) PCLDSW(kdlon,kflev)
182  REAL(KIND=8) PTAU(kdlon,2,kflev)
183  REAL(KIND=8) POMEGA(kdlon,2,kflev)
184  REAL(KIND=8) PCG(kdlon,2,kflev)
185  REAL(KIND=8) zfract(kdlon), zrmu0(kdlon), zdist
186  REAL(KIND=8) zheat(kdlon,kflev), zcool(kdlon,kflev)
187  REAL(KIND=8) zheat0(kdlon,kflev), zcool0(kdlon,kflev)
188  REAL(KIND=8) ztopsw(kdlon), ztoplw(kdlon)
189  REAL(KIND=8) zsolsw(kdlon), zsollw(kdlon), zalbpla(kdlon)
190  REAL(KIND=8) zsollwdown(kdlon)
191  REAL(KIND=8) ztopsw0(kdlon), ztoplw0(kdlon)
192  REAL(KIND=8) zsolsw0(kdlon), zsollw0(kdlon)
193  REAL(KIND=8) zznormcp
194  REAL(KIND=8) tauaero(kdlon,kflev,9,2)                     ! aer opt properties
195  REAL(KIND=8) pizaero(kdlon,kflev,9,2)
196  REAL(KIND=8) cgaero(kdlon,kflev,9,2)
197  REAL(KIND=8) PTAUA(kdlon,2,kflev)                         ! present-day value of cloud opt thickness (PTAU is pre-industrial value), local use
198  REAL(KIND=8) POMEGAA(kdlon,2,kflev)                       ! dito for single scatt albedo
199  REAL(KIND=8) ztopswadaero(kdlon), zsolswadaero(kdlon)     ! Aerosol direct forcing at TOAand surface
200  REAL(KIND=8) ztopswad0aero(kdlon), zsolswad0aero(kdlon)   ! Aerosol direct forcing at TOAand surface
201  REAL(KIND=8) ztopswaiaero(kdlon), zsolswaiaero(kdlon)     ! dito, indirect
202  REAL(KIND=8) ztopsw_aero(kdlon,9), ztopsw0_aero(kdlon,9)
203  REAL(KIND=8) zsolsw_aero(kdlon,9), zsolsw0_aero(kdlon,9)
204  REAL(KIND=8) ztopswcf_aero(kdlon,3), zsolswcf_aero(kdlon,3)     
205  real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
206
207  call assert(size(wo, 1) == klon, size(wo, 2) == klev, "radlwsw wo")
208  ! initialisation
209  tauaero(:,:,:,:)=0.
210  pizaero(:,:,:,:)=0.
211  cgaero(:,:,:,:)=0.
212 
213  !
214  !-------------------------------------------
215  nb_gr = KLON / kdlon
216  IF (nb_gr*kdlon .NE. KLON) THEN
217      PRINT*, "kdlon mauvais:", KLON, kdlon, nb_gr
218      CALL abort
219  ENDIF
220  IF (kflev .NE. KLEV) THEN
221      PRINT*, "kflev differe de KLEV, kflev, KLEV"
222      CALL abort
223  ENDIF
224  !-------------------------------------------
225  DO k = 1, KLEV
226    DO i = 1, KLON
227      heat(i,k)=0.
228      cool(i,k)=0.
229      heat0(i,k)=0.
230      cool0(i,k)=0.
231    ENDDO
232  ENDDO
233  !
234  zdist = dist
235  !
236  PSCT = solaire/zdist/zdist
237
238  IF (type_trac == 'repr') THEN
239#ifdef REPROBUS
240     if(ok_SUNTIME) PSCT = solaireTIME/zdist/zdist
241     print*,'Constante solaire: ',PSCT*zdist*zdist
242#endif
243  END IF
244
245  DO j = 1, nb_gr
246    iof = kdlon*(j-1)
247    DO i = 1, kdlon
248      zfract(i) = fract(iof+i)
249      zrmu0(i) = rmu0(iof+i)
250      PALBD(i,1) = alb1(iof+i)
251      PALBD(i,2) = alb2(iof+i)
252      PALBP(i,1) = alb1(iof+i)
253      PALBP(i,2) = alb2(iof+i)
254      PEMIS(i) = 1.0
255      PVIEW(i) = 1.66
256      PPSOL(i) = paprs(iof+i,1)
257      zx_alpha1 = (paprs(iof+i,1)-pplay(iof+i,2))/(pplay(iof+i,1)-pplay(iof+i,2))
258      zx_alpha2 = 1.0 - zx_alpha1
259      PTL(i,1) = t(iof+i,1) * zx_alpha1 + t(iof+i,2) * zx_alpha2
260      PTL(i,KLEV+1) = t(iof+i,KLEV)
261      PDT0(i) = tsol(iof+i) - PTL(i,1)
262    ENDDO
263    DO k = 2, kflev
264      DO i = 1, kdlon
265        PTL(i,k) = (t(iof+i,k)+t(iof+i,k-1))*0.5
266      ENDDO
267    ENDDO
268    DO k = 1, kflev
269      DO i = 1, kdlon
270        PDP(i,k) = paprs(iof+i,k)-paprs(iof+i,k+1)
271        PTAVE(i,k) = t(iof+i,k)
272        PWV(i,k) = MAX (q(iof+i,k), 1.0e-12)
273        PQS(i,k) = PWV(i,k)
274        POZON(i,k, :) = wo(iof+i, k, :) * RG * dobson_u * 1e3 &
275             / (paprs(iof+i, k) - paprs(iof+i, k+1))
276        PCLDLD(i,k) = cldfra(iof+i,k)*cldemi(iof+i,k)
277        PCLDLU(i,k) = cldfra(iof+i,k)*cldemi(iof+i,k)
278        PCLDSW(i,k) = cldfra(iof+i,k)
279        PTAU(i,1,k) = MAX(cldtaupi(iof+i,k), 1.0e-05)! 1e-12 serait instable
280        PTAU(i,2,k) = MAX(cldtaupi(iof+i,k), 1.0e-05)! pour 32-bit machines
281        POMEGA(i,1,k) = 0.9999 - 5.0e-04 * EXP(-0.5 * PTAU(i,1,k))
282        POMEGA(i,2,k) = 0.9988 - 2.5e-03 * EXP(-0.05 * PTAU(i,2,k))
283        PCG(i,1,k) = 0.865
284        PCG(i,2,k) = 0.910
285        !-
286        ! Introduced for aerosol indirect forcings.
287        ! The following values use the cloud optical thickness calculated from
288        ! present-day aerosol concentrations whereas the quantities without the
289        ! "A" at the end are for pre-industial (natural-only) aerosol concentrations
290        !
291        PTAUA(i,1,k) = MAX(cldtaupd(iof+i,k), 1.0e-05)! 1e-12 serait instable
292        PTAUA(i,2,k) = MAX(cldtaupd(iof+i,k), 1.0e-05)! pour 32-bit machines
293        POMEGAA(i,1,k) = 0.9999 - 5.0e-04 * EXP(-0.5 * PTAUA(i,1,k))
294        POMEGAA(i,2,k) = 0.9988 - 2.5e-03 * EXP(-0.05 * PTAUA(i,2,k))
295      ENDDO
296    ENDDO
297
298    IF (type_trac == 'repr') THEN
299#ifdef REPROBUS
300       ndimozon = size(wo, 3)
301       CALL RAD_INTERACTIF(POZON,iof)
302#endif
303    END IF
304
305    !
306    DO k = 1, kflev+1
307      DO i = 1, kdlon
308        PPMB(i,k) = paprs(iof+i,k)/100.0
309      ENDDO
310    ENDDO
311    !
312    DO kk = 1, 5
313      DO k = 1, kflev
314        DO i = 1, kdlon
315          PAER(i,k,kk) = 1.0E-15
316        ENDDO
317      ENDDO
318    ENDDO
319    DO k = 1, kflev
320      DO i = 1, kdlon
321        tauaero(i,k,:,1)=tau_aero(iof+i,k,:,1)
322        pizaero(i,k,:,1)=piz_aero(iof+i,k,:,1)
323        cgaero(i,k,:,1) =cg_aero(iof+i,k,:,1)
324        tauaero(i,k,:,2)=tau_aero(iof+i,k,:,2)
325        pizaero(i,k,:,2)=piz_aero(iof+i,k,:,2)
326        cgaero(i,k,:,2) =cg_aero(iof+i,k,:,2)
327      ENDDO
328    ENDDO
329
330!
331!===== iflag_rrtm ================================================
332!     
333    IF (iflag_rrtm == 0) THEN
334       ! Old radiation scheme, used for AR4 runs
335       ! average day-night ozone for longwave
336       CALL LW_LMDAR4(&
337            PPMB, PDP,&
338            PPSOL,PDT0,PEMIS,&
339            PTL, PTAVE, PWV, POZON(:, :, 1), PAER,&
340            PCLDLD,PCLDLU,&
341            PVIEW,&
342            zcool, zcool0,&
343            ztoplw,zsollw,ztoplw0,zsollw0,&
344            zsollwdown,&
345            ZFLUP, ZFLDN, ZFLUP0,ZFLDN0)
346
347       ! daylight ozone, if we have it, for short wave
348       IF (.NOT. new_aod) THEN
349          ! use old version
350          CALL SW_LMDAR4(PSCT, zrmu0, zfract,&
351               PPMB, PDP, &
352               PPSOL, PALBD, PALBP,&
353               PTAVE, PWV, PQS, POZON(:, :, size(wo, 3)), PAER,&
354               PCLDSW, PTAU, POMEGA, PCG,&
355               zheat, zheat0,&
356               zalbpla,ztopsw,zsolsw,ztopsw0,zsolsw0,&
357               ZFSUP,ZFSDN,ZFSUP0,ZFSDN0,&
358               tauaero(:,:,5,:), pizaero(:,:,5,:), cgaero(:,:,5,:),&
359               PTAUA, POMEGAA,&
360               ztopswadaero,zsolswadaero,&
361               ztopswaiaero,zsolswaiaero,&
362               ok_ade, ok_aie, flag_aerosol)
363         
364       ELSE ! new_aod=T         
365          CALL SW_AEROAR4(PSCT, zrmu0, zfract,&
366               PPMB, PDP,&
367               PPSOL, PALBD, PALBP,&
368               PTAVE, PWV, PQS, POZON(:, :, size(wo, 3)), PAER,&
369               PCLDSW, PTAU, POMEGA, PCG,&
370               zheat, zheat0,&
371               zalbpla,ztopsw,zsolsw,ztopsw0,zsolsw0,&
372               ZFSUP,ZFSDN,ZFSUP0,ZFSDN0,&
373               tauaero, pizaero, cgaero, &
374               PTAUA, POMEGAA,&
375               ztopswadaero,zsolswadaero,&
376               ztopswad0aero,zsolswad0aero,&
377               ztopswaiaero,zsolswaiaero, &
378               ztopsw_aero,ztopsw0_aero,&
379               zsolsw_aero,zsolsw0_aero,&
380               ztopswcf_aero,zsolswcf_aero, &
381               ok_ade, ok_aie, flag_aerosol)
382       ENDIF
383
384    ELSE 
385!===== iflag_rrtm=1, on passe dans SW via RECMWFL ===============
386       WRITE(lunout,*) "Option iflag_rrtm=T ne fonctionne pas encore !!!"
387       CALL abort_gcm('radlwsw','iflag_rrtm=T not valid',1)
388
389    ENDIF ! iflag_rrtm
390!======================================================================
391
392    DO i = 1, kdlon
393      radsol(iof+i) = zsolsw(i) + zsollw(i)
394      topsw(iof+i) = ztopsw(i)
395      toplw(iof+i) = ztoplw(i)
396      solsw(iof+i) = zsolsw(i)
397      sollw(iof+i) = zsollw(i)
398      sollwdown(iof+i) = zsollwdown(i)
399      DO k = 1, kflev+1
400        lwdn0 ( iof+i,k)   = ZFLDN0 ( i,k)
401        lwdn  ( iof+i,k)   = ZFLDN  ( i,k)
402        lwup0 ( iof+i,k)   = ZFLUP0 ( i,k)
403        lwup  ( iof+i,k)   = ZFLUP  ( i,k)
404      ENDDO
405      topsw0(iof+i) = ztopsw0(i)
406      toplw0(iof+i) = ztoplw0(i)
407      solsw0(iof+i) = zsolsw0(i)
408      sollw0(iof+i) = zsollw0(i)
409      albpla(iof+i) = zalbpla(i)
410
411      DO k = 1, kflev+1
412        swdn0 ( iof+i,k)   = ZFSDN0 ( i,k)
413        swdn  ( iof+i,k)   = ZFSDN  ( i,k)
414        swup0 ( iof+i,k)   = ZFSUP0 ( i,k)
415        swup  ( iof+i,k)   = ZFSUP  ( i,k)
416      ENDDO
417    ENDDO
418    !-transform the aerosol forcings, if they have
419    ! to be calculated
420    IF (ok_ade) THEN
421        DO i = 1, kdlon
422          topswad_aero(iof+i) = ztopswadaero(i)
423          topswad0_aero(iof+i) = ztopswad0aero(i)
424          solswad_aero(iof+i) = zsolswadaero(i)
425          solswad0_aero(iof+i) = zsolswad0aero(i)
426! MS the following lines seem to be wrong, why is iof on right hand side???
427!          topsw_aero(iof+i,:) = ztopsw_aero(iof+i,:)
428!          topsw0_aero(iof+i,:) = ztopsw0_aero(iof+i,:)
429!          solsw_aero(iof+i,:) = zsolsw_aero(iof+i,:)
430!          solsw0_aero(iof+i,:) = zsolsw0_aero(iof+i,:)
431          topsw_aero(iof+i,:) = ztopsw_aero(i,:)
432          topsw0_aero(iof+i,:) = ztopsw0_aero(i,:)
433          solsw_aero(iof+i,:) = zsolsw_aero(i,:)
434          solsw0_aero(iof+i,:) = zsolsw0_aero(i,:)
435          topswcf_aero(iof+i,:) = ztopswcf_aero(i,:)
436          solswcf_aero(iof+i,:) = zsolswcf_aero(i,:)         
437        ENDDO
438    ELSE
439        DO i = 1, kdlon
440          topswad_aero(iof+i) = 0.0
441          solswad_aero(iof+i) = 0.0
442          topswad0_aero(iof+i) = 0.0
443          solswad0_aero(iof+i) = 0.0
444          topsw_aero(iof+i,:) = 0.
445          topsw0_aero(iof+i,:) =0.
446          solsw_aero(iof+i,:) = 0.
447          solsw0_aero(iof+i,:) = 0.
448        ENDDO
449    ENDIF
450    IF (ok_aie) THEN
451        DO i = 1, kdlon
452          topswai_aero(iof+i) = ztopswaiaero(i)
453          solswai_aero(iof+i) = zsolswaiaero(i)
454        ENDDO
455    ELSE
456        DO i = 1, kdlon
457          topswai_aero(iof+i) = 0.0
458          solswai_aero(iof+i) = 0.0
459        ENDDO
460    ENDIF
461    DO k = 1, kflev
462      DO i = 1, kdlon
463        !        scale factor to take into account the difference between
464        !        dry air and watter vapour scpecifi! heat capacity
465        zznormcp=1.0+RVTMP2*PWV(i,k)
466        heat(iof+i,k) = zheat(i,k)/zznormcp
467        cool(iof+i,k) = zcool(i,k)/zznormcp
468        heat0(iof+i,k) = zheat0(i,k)/zznormcp
469        cool0(iof+i,k) = zcool0(i,k)/zznormcp
470      ENDDO
471    ENDDO
472
473 ENDDO ! j = 1, nb_gr
474
475END SUBROUTINE radlwsw
476
477end module radlwsw_m
Note: See TracBrowser for help on using the repository browser.