Newer
Older
!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
!,! ########################
MODULE MODI_PHYS_PARAM_n
! ########################
!
!
INTERFACE
!
SUBROUTINE PHYS_PARAM_n(KTCOUNT,HFMFILE,OCLOSE_OUT, &
PRAD,PSHADOWS,PKAFR,PGROUND,PMAFL,PDRAG,PTURB,PTRACER,PCHEM, &
PTIME_BU, OMASKkids )
!
INTEGER, INTENT(IN) :: KTCOUNT ! temporal iteration count
CHARACTER (LEN=28),INTENT(IN) :: HFMFILE ! name of the synchronous
! OUTPUT FM-file
LOGICAL, INTENT(IN) :: OCLOSE_OUT! conditional closure of the
! OUTPUT FM-file
! advection schemes
REAL*8,DIMENSION(2), INTENT(INOUT) :: PRAD,PSHADOWS,PKAFR,PGROUND,PTURB,PMAFL,PDRAG,PTRACER ! to store CPU
! time for computing time
REAL*8,DIMENSION(2), INTENT(INOUT) :: PCHEM ! to store CPU time for chemistry
REAL*8,DIMENSION(2), INTENT(INOUT) :: PTIME_BU ! time used in budget&LES budgets
! statistics
LOGICAL, DIMENSION(:,:), INTENT(IN) :: OMASKkids ! kids domains mask
END SUBROUTINE PHYS_PARAM_n
!
END INTERFACE
!
END MODULE MODI_PHYS_PARAM_n
!
! ######################################################################
SUBROUTINE PHYS_PARAM_n(KTCOUNT,HFMFILE,OCLOSE_OUT, &
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
PRAD,PSHADOWS,PKAFR,PGROUND,PMAFL,PDRAG,PTURB,PTRACER,PCHEM, &
PTIME_BU, OMASKkids )
! ######################################################################
!
!!**** *PHYS_PARAM_n * -monitor of the parameterizations used by model _n
!!
!! PURPOSE
!! -------
! The purpose of this routine is to update the sources by adding the
! parameterized terms. This is realized by sequentially calling the
! specialized routines.
!
!!** METHOD
!! ------
!! The first parametrization is the radiation scheme:
!! ----------------
!! * CRAD = 'FIXE'
!! In this case, a temporal interpolation is performed for the downward
!! surface fluxes XFLALWD and XFLASWD.
!! * CRAD = 'ECMWF'
!! Several tests are performed before calling the radiation computations
!! interface with the ECMWF radiation scheme code. A control is made to
!! ensure that:
!! - the full radiation code is called at the first model timestep
!! - there is a priority for calling the full radiation instead of the
!! cloud-only approximation if both must be called at the current
!! timestep
!! - the cloud-only option (approximation) is coherent with the
!! occurence of one cloudy vertical column at least
!! If all the above conditions are fulfilled (GRAD is .TRUE.) then the
!! position of the sun is computed in routine SUNPOS_n and the interfacing
!! routine RADIATIONS is called to update the radiative tendency XDTHRAD
!! and the downward surface fluxes XFLALWD and XFLASWD. Finally, the
!! radiative tendency is integrated as a source term in the THETA prognostic
!! equation.
!!
!! The second parameterization is the soil scheme:
!! -----------
!!
!! externalized surface
!!
!! The third parameterization is the turbulence scheme:
!! -----------------
!! * CTURB='NONE'
!! no turbulent mixing is taken into account
!! * CTURB='TKEL'
!! The turbulent fluxes are computed according to a one and half order
!! closure of the hydrodynamical equations. This scheme is based on a
!! prognostic for the turbulent kinetic energy and a mixing length
!! computation ( the mesh size or a physically based length). Other
!! turbulent moments are diagnosed according to a stationarization of the
!! second order turbulent moments. This turbulent scheme forecasts
!! either a purely vertical turbulent mixing or 3-dimensional mixing
!! according to its internal degrees of freedom.
!!
!!
!! The LAST parameterization is the chemistry scheme:
!! -----------------
!! The chemistry part of MesoNH has two namelists, NAM_SOLVER for the
!! parameters concerning the stiff solver, and NAM_MNHCn concerning the
!! configuration and options of the chemistry module itself.
!! The switch LUSECHEM in NAM_CONF acitvates or deactivates the chemistry.
!! The only variables of MesoNH that are modified by chemistry are the
!! scalar variables. If calculation of chemical surface fluxes is
!! requested, those fluxes are calculated before
!! entering the turbulence scheme, since those fluxes are taken into
!! account by TURB as surface boundary conditions.
!! CAUTION: chemistry has allways to be called AFTER ALL OTHER TERMS
!! that affect the scalar variables (dynamical terms, forcing,
!! parameterizations (like TURB, CONVECTION), since it uses the variables
!! XRSVS as input in case of the time-split option.
!!
!! EXTERNAL
!! --------
!! Subroutine SUNPOS_n : computes the position of the sun
!! Subroutine RADIATIONS : computes the radiative tendency and fluxes
!! Subroutine TSZ0 : computes the surface from temporally
!! interpolated Ts and given z0
!! Subroutine ISBA : computes the surface fluxes from a soil scheme
!! Subroutine TURB : computes the turbulence source terms
!! Subroutine CONVECTION : computes the convection source term
!! Subroutine CH_SURFACE_FLUX_n: computes the surface flux for chemical
!! species
!! Subroutine CH_MONITOR_n : computes the chemistry source terms
!! that are applied to the scalar variables
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!! USE MODD_DYN
!! USE MODD_CONF
!! USE MODD_CONF_n
!! USE MODD_CURVCOR_n
!! USE MODD_DYN_n
!! USE MODD_FIELD_n
!! USE MODD_GR_FIELD_n
!! USE MODD_LSFIELD_n
!! USE MODD_GRID_n
!! USE MODD_LBC_n
!! USE MODD_PARAM_RAD_n
!! USE MODD_RADIATIONS_n
!! USE MODD_REF_n
!! USE MODD_LUNIT_n
!! USE MODD_TIME_n
!! USE MODD_CH_MNHC_n
!!
!! REFERENCE
!! ---------
!! None
!!
!! AUTHOR
!! ------
!! J. Stein * Meteo-France *
!!
!! MODIFICATIONS
!! -------------
!! Original 05/01/95
!! Modifications Feb 14, 1995 (J.Cuxart) add the I/O arguments,
!! the director cosinus and change the names of the surface fluxes
!! Modifications March 21, 1995 (J.M.Carriere) take into account liquid
!! water
!! June 30,1995 (J.Stein) initialize at 0 the surf. fluxes
!! Modifications Sept. 1, 1995 (S.Belair) ISBA scheme
!! Modifications Sept.25, 1995 (J.Stein) switch on the radiation scheme
!! Modifications Sept. 11, 1995 (J.-P. Pinty) radiation scheme
!! Nov. 15, 1995 (J.Stein) cleaning + change the temporal
!! algorithm for the soil scheme-turbulence
!! Jan. 23, 1996 (J.Stein) add a new option for the surface
!! fluxes where Ts and z0 are given
!! March 18, 1996 (J.Stein) add the cloud fraction
!! March 28, 1996 (J.Stein) the soil scheme gives energy
!! fluxes + cleaning
!! June 17, 1996 (Lafore) statistics of computing time
!! August 4, 1996 (K. Suhre) add chemistry
!! Oct. 12, 1996 (J.Stein) use XSRCM in the turbulence
!! scheme
!! Nov. 18, 1996 (J.-P. Pinty) add domain translation
!! change arg. in radiations
!! Fev. 4, 1997 (J.Viviand) change isba's calling for ice
!! Jun. 22, 1997 (J.Stein) change the equation system and use
!! the absolute pressure
!! Jul. 09, 1997 (V.Masson) add directional z0
!! Jan. 24, 1998 (P.Bechtold) add convective transport for tracers
!! Jan. 24, 1998 (J.-P. Pinty) split SW and LW part for radiation
!! Mai. 10, 1999 (P.Bechtold) shallow convection
!! Oct. 20, 1999 (P.Jabouille) domain translation for turbulence
!! Jan. 04, 2000 (V.Masson) removes TSZ0 case
!! Jan. 04, 2000 (V.Masson) modifies albedo computation
! Jul 02, 2000 (F.Solmon/V.Masson) adaptation for patch approach
!! Nov. 15, 2000 (V.Masson) LES routines
!! Nov. 15, 2000 (V.Masson) effect of slopes on surface fluxes
!! Feb. 02, 2001 (P.Tulet) add friction velocities and aerodynamical
!! resistance (patch approach)
!! Jan. 04, 2000 (V.Masson) modify surf_rad_modif computation
!! Mar. 04, 2002 (F.Solmon) new interface for radiation call
!! Nov. 06, 2002 (V.Masson) LES budgets & budget time counters
!! Jan. 2004 (V.Masson) surface externalization
!! Jan. 13, 2004 (J.Escobar) bug correction : compute "GRAD" in parallel
!! Jan. 20, 2005 (P. Tulet) add dust sedimentation
!! Jan. 20, 2005 (P. Tulet) climatologic SSA
!! Jan. 20, 2005 (P. Tulet) add aerosol / dust scavenging
!! Jul. 2005 (N. Asencio) use the two-way result-fields
!! before ground_param call
!! May 2006 Remove EPS
!! Oct. 2007 (J.Pergaud) Add shallow_MF
!! Oct. 2009 (C.Lac) Introduction of different PTSTEP according to the
!! advection schemes
!! Oct. 2009 (V. MAsson) optimization of Pergaud et al massflux scheme
!! Aug. 2010 (V.Masson, C.Lac) Exchange of SBL_DEPTH for
!! reproducibility
!! Oct. 2010 (J.Escobar) init ZTIME_LES_MF ( pb detected with g95 )
!! Feb. 2011 (V.Masson, C.Lac) SBL_DEPTH values on outer pts
!! for RMC01
!! Sept.2011 (J.Escobar) init YINST_SFU ='M'
!!
!! Specific for 2D modeling :
!!
!! 06/2010 (P.Peyrille) add Call to aerozon.f90 if LAERO_FT=T
!! to update
!! aerosols and ozone climatology at each call to
!! phys_param otherwise it is constant to monthly average
Gaelle Tanguy
committed
!! 01/2014 (C.Lac) correction for the nesting of 2D surface
!! fields if the number of the son model does not
!! follow the number of the dad model
ESCOBAR MUNOZ Juan
committed
!! J.Escobar 21/03/2013: for HALOK comment all NHALO=1 test
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
!!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
!
USE MODE_ll
USE MODE_FM
USE MODE_FMWRIT
USE MODD_ARGSLIST_ll, ONLY : LIST_ll
!
USE MODD_CST
USE MODD_DYN
USE MODD_CONF
USE MODD_FRC
USE MODD_PARAMETERS
USE MODD_GRID
USE MODD_NSV
USE MODD_LES
USE MODD_LES_BUDGET
!
USE MODD_CONF_n
USE MODD_CURVCOR_n
USE MODD_DYN_n
USE MODD_FIELD_n
USE MODD_LSFIELD_n
USE MODD_GRID_n
USE MODD_METRICS_n
USE MODD_LBC_n
USE MODD_REF_n
USE MODD_LUNIT_n
USE MODD_OUT_n
USE MODD_PARAM_n
USE MODD_PARAM_RAD_n
USE MODD_PARAM_KAFR_n
USE MODD_RADIATIONS_n
USE MODD_SHADOWS_n
USE MODD_DEEP_CONVECTION_n
USE MODD_TIME_n
USE MODD_TURB_n
USE MODD_CH_MNHC_n, ONLY : LUSECHEM, &! indicates if chemistry is used
LCH_CONV_SCAV, &
LCH_CONV_LINOX
USE MODD_PRECIP_n
USE MODD_PASPOL_n
USE MODD_BUDGET
USE MODD_RAIN_ICE_DESCR, ONLY : XRTMIN
USE MODD_ICE_C1R3_DESCR, ONLY : XRTMIN_C1R3=>XRTMIN
USE MODD_TURB_CLOUD, ONLY : CTURBLEN_CLOUD,NMODEL_CLOUD, &
XCEI,XCEI_MIN,XCEI_MAX,XCOEF_AMPL_SAT
USE MODD_NESTING, ONLY : XWAY,NDAD, NDXRATIO_ALL, NDYRATIO_ALL
USE MODD_PARAM_ICE, ONLY : LSEDIC
USE MODD_PARAM_C2R2, ONLY : LSEDC
USE MODD_SUB_PHYS_PARAM_n
!
USE MODD_PARAM_MFSHALL_n
USE MODI_SHALLOW_MF_PACK
USE MODD_CLOUD_MF_n
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
!
USE MODI_SURF_RAD_MODIF
USE MODI_GROUND_PARAM_n
USE MODI_TURB
USE MODI_SUNPOS_n
USE MODI_RADIATIONS
USE MODI_CONVECTION
USE MODI_TEMPORAL_DIST
USE MODI_CH_MONITOR_n
USE MODI_AER_MONITOR_n
USE MODI_BUDGET
USE MODI_PASPOL
USE MODI_CONDSAMP
USE MODE_FM
USE MODE_MODELN_HANDLER
USE MODI_SEDIM_DUST
USE MODI_SEDIM_SALT
USE MODI_DUST_FILTER
USE MODI_SALT_FILTER
USE MODI_DRAG_VEG
USE MODD_DUST
USE MODD_SALT
USE MODD_PASPOL
USE MODD_CONDSAMP
USE MODD_CH_AEROSOL
USE MODE_DUST_PSD
USE MODE_SALT_PSD
USE MODE_AERO_PSD
USE MODE_MNH_TIMING
USE MODD_TURB_FLUX_AIRCRAFT_BALLOON, ONLY : XTHW_FLUX, XRCW_FLUX, XSVW_FLUX
USE MODD_DRAGTREE
!
USE MODD_TIME, ONLY : TDTEXP ! Ajout PP
USE MODI_AEROZON ! Ajout PP
!
USE MODI_EDDY_FLUX_n ! Ajout PP
USE MODI_EDDYUV_FLUX_n ! Ajout PP
USE MODI_EDDY_FLUX_ONE_WAY_n ! Ajout PP
USE MODI_EDDYUV_FLUX_ONE_WAY_n ! Ajout PP
USE MODD_DEF_EDDY_FLUX_n ! Ajout PP
USE MODD_DEF_EDDYUV_FLUX_n ! Ajout PP
USE MODD_LATZ_EDFLX
USE MODD_PARAM_LIMA, ONLY : MSEDC => LSEDC, XRTMIN_LIMA=>XRTMIN
!
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
IMPLICIT NONE
!
!* 0.1 declarations of arguments
!
INTEGER, INTENT(IN) :: KTCOUNT ! temporal iteration count
CHARACTER (LEN=28),INTENT(IN) :: HFMFILE ! name of the synchronous
! OUTPUT FM-file
LOGICAL, INTENT(IN) :: OCLOSE_OUT! conditional closure of the
! OUTPUT FM-file
! advection schemes
REAL*8,DIMENSION(2), INTENT(INOUT) :: PRAD,PSHADOWS,PKAFR,PGROUND,PTURB,PMAFL,PDRAG,PTRACER ! to store CPU
! time for computing time
! statistics
REAL*8,DIMENSION(2), INTENT(INOUT) :: PCHEM ! to store CPU time for chemistry
REAL*8,DIMENSION(2), INTENT(INOUT) :: PTIME_BU ! time used in budget&LES budgets
LOGICAL, DIMENSION(:,:), INTENT(IN) :: OMASKkids ! kids domains mask
!
!* 0.2 declarations of local variables
!
REAL, DIMENSION(:,:), ALLOCATABLE :: ZSFU ! surface flux of x and
REAL, DIMENSION(:,:), ALLOCATABLE :: ZSFV ! y component of wind
REAL, DIMENSION(:,:), ALLOCATABLE :: ZSFTH ! surface flux of theta
REAL, DIMENSION(:,:), ALLOCATABLE :: ZSFRV ! surface flux of vapor
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZSFSV ! surface flux of scalars
REAL, DIMENSION(:,:), ALLOCATABLE :: ZSFCO2! surface flux of CO2
!
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZDIR_ALB ! direct albedo
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZSCA_ALB ! diffuse albedo
REAL, DIMENSION(:,:), ALLOCATABLE :: ZEMIS ! emissivity
REAL, DIMENSION(:,:), ALLOCATABLE :: ZTSRAD ! surface temperature
REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: ZRGDST,ZSIGDST,ZNDST,ZSVDST
REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: ZRGSLT,ZSIGSLT,ZNSLT,ZSVSLT
REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: ZRGAER,ZSIGAER,ZNAER,ZSVAER
REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: ZSVT
!
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZEXN ! Atmospheric density and Exner
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZSIGMF ! MF contribution to XSIGS
!
REAL, DIMENSION(0:24) :: ZRG_HOUR = (/ 0., 0., 0., 0., 0., 32.04, 114.19, &
228.01, 351.25, 465.49, 557.24, &
616.82, 638.33, 619.43, 566.56, &
474.71, 359.20, 230.87, 115.72, &
32.48, 0., 0., 0., 0., 0. /)
!
REAL, DIMENSION(0:24) :: ZRAT_HOUR = (/ 326.00, 325.93, 325.12, 324.41, &
323.16, 321.95, 322.51, 325.16, &
328.01, 331.46, 335.58, 340.00, &
345.20, 350.32, 354.20, 356.58, &
356.56, 355.33, 352.79, 351.34, &
347.00, 342.00, 337.00, 332.00, &
326.00 /)
!
INTEGER :: IHOUR ! parameters necessary for the temporal
REAL :: ZTIME, ZDT ! interpolation
REAL :: ZTEMP_DIST ! time between 2 instants (in seconds)
!
LOGICAL :: GRAD ! conditionnal call for the full radiation
! computations
REAL :: ZRAD_GLOB_ll ! 'real' global parallel mask of 'GRAD'
INTEGER :: INFO_ll ! error report of parallel routines
LOGICAL :: GCLOUD_ONLY ! conditionnal radiation computations for
! the only cloudy columns
!
REAL*8,DIMENSION(2) :: ZTIME1,ZTIME2,ZTIME3,ZTIME4 ! for computing time analysis
REAL*8,DIMENSION(2) :: ZTIME_LES_MF ! time spent in LES computation in shallow conv.
LOGICAL :: GDCONV ! conditionnal call for the deep convection
! computations
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZRC, ZRI, ZWT ! additional dummies
REAL, DIMENSION(:,:), ALLOCATABLE :: ZDXDY ! grid area
! for rc, ri, w required if main variables not allocated
!
INTEGER :: IIU, IJU, IKU ! dimensional indexes
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
!
INTEGER :: JSV ! Loop index for Scalar Variables
INTEGER :: JSWB ! loop on SW spectral bands
INTEGER :: IIB,IIE,IJB,IJE, IKB, IKE
INTEGER :: IMODEIDX
! index values for the Beginning or the End of the physical
! domain in x and y directions
TYPE(LIST_ll), POINTER :: TZFIELDS_ll ! list of fields to exchange
INTEGER :: IINFO_ll ! return code of parallel routine
!
!* variables for writing in a fm file
!
INTEGER :: IRESP ! IRESP : return-code if a problem appears
!in LFI subroutines at the open of the file
INTEGER :: ILUOUT ! logical unit numbers of output-listing
INTEGER :: IMI ! model index
INTEGER :: JKID ! loop index to look for the KID models
REAL :: ZINIRADIUSI, ZINIRADIUSJ ! ORILAM initial radius
REAL, DIMENSION(NMODE_DST) :: ZINIRADIUS ! DUST initial radius
REAL, DIMENSION(NMODE_SLT) :: ZINIRADIUS_SLT ! Sea Salt initial radius
REAL, DIMENSION(SIZE(XRSVS,1), SIZE(XRSVS,2), SIZE(XRSVS,3), SIZE(XRSVS,4)) :: ZRSVS
REAL, DIMENSION(SIZE(XRSVS,1), SIZE(XRSVS,2), SIZE(XRSVS,3), NSV_AER) :: ZWETDEPAER
LOGICAL :: GCLD ! conditionnal call for dust wet deposition
! * arrays to store the surface fields before radiation and convection scheme
! calls
INTEGER :: IMODSON ! Number of son models of IMI with XWAY=2
INTEGER :: IKIDM ! index loop
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZSAVE_INPRR,ZSAVE_INPRS,ZSAVE_INPRG,ZSAVE_INPRH
REAL, DIMENSION(:,:,:), ALLOCATABLE :: ZSAVE_INPRC,ZSAVE_PRCONV,ZSAVE_PRSCONV
REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: ZSAVE_DIRFLASWD, ZSAVE_SCAFLASWD,ZSAVE_DIRSRFSWD
!
!-----------------------------------------------------------------------------
!
NULLIFY(TZFIELDS_ll)
IMI=GET_CURRENT_MODEL_INDEX()
!
CALL FMLOOK_ll(CLUOUT,CLUOUT,ILUOUT,IRESP)
CALL GET_DIM_EXT_ll ('B',IIU,IJU)
IKU=SIZE(XTHT,3)
IKB = 1 + JPVEXT
IKE = IKU - JPVEXT
CALL GET_INDICE_ll (IIB,IJB,IIE,IJE)
!
ZTIME1 = 0.0
ZTIME2 = 0.0
ZTIME3 = 0.0
ZTIME4 = 0.0
PTIME_BU = 0.
ZTIME_LES_MF = 0.0
ZWETDEPAER(:,:,:,:) = 0.
!
!* allocation of variables used in more than one parameterization
!
ALLOCATE(ZSFU (IIU,IJU)) ! surface schemes + turbulence
ALLOCATE(ZSFV (IIU,IJU))
ALLOCATE(ZSFTH (IIU,IJU))
ALLOCATE(ZSFRV (IIU,IJU))
ALLOCATE(ZSFSV (IIU,IJU,NSV))
ALLOCATE(ZSFCO2(IIU,IJU))
!
!* if XWAY(son)=2 save surface fields before radiation or convective scheme
! calls
!
IMODSON = 0
DO JKID = IMI+1,NMODEL ! min value of the possible kids
IF (IMI == NDAD(JKID) .AND. XWAY(JKID) == 2. .AND. CPROGRAM=='MESONH' &
.AND. (CCONF == 'RESTA' .OR. (CCONF == 'START' .AND. KTCOUNT /= 1))) THEN
IMODSON = IMODSON + 1
END IF
END DO
!
IF (IMODSON /= 0 ) THEN
IF (LUSERC .AND. ( &
(LSEDIC .AND. CCLOUD(1:3) == 'ICE') .OR. &
(LSEDC .AND. (CCLOUD == 'C2R2' .OR. CCLOUD == 'KHKO')) .OR. &
(MSEDC .AND. CCLOUD=='LIMA') &
)) THEN
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
ALLOCATE( ZSAVE_INPRC(SIZE(XINPRC,1),SIZE(XINPRC,2),IMODSON))
ELSE
ALLOCATE( ZSAVE_INPRC(0,0,0))
END IF
IF (LUSERR) THEN
ALLOCATE( ZSAVE_INPRR(SIZE(XINPRR,1),SIZE(XINPRR,2),IMODSON))
ELSE
ALLOCATE( ZSAVE_INPRR(0,0,0))
END IF
IF (LUSERS) THEN
ALLOCATE( ZSAVE_INPRS(SIZE(XINPRS,1),SIZE(XINPRS,2),IMODSON))
ELSE
ALLOCATE( ZSAVE_INPRS(0,0,0))
END IF
IF (LUSERG) THEN
ALLOCATE( ZSAVE_INPRG(SIZE(XINPRG,1),SIZE(XINPRG,2),IMODSON))
ELSE
ALLOCATE( ZSAVE_INPRG(0,0,0))
END IF
IF (LUSERH) THEN
ALLOCATE( ZSAVE_INPRH(SIZE(XINPRH,1),SIZE(XINPRH,2),IMODSON))
ELSE
ALLOCATE( ZSAVE_INPRH(0,0,0))
END IF
IF (CDCONV /= 'NONE') THEN
ALLOCATE( ZSAVE_PRCONV(SIZE(XPRCONV,1),SIZE(XPRCONV,2),IMODSON))
ALLOCATE( ZSAVE_PRSCONV(SIZE(XPRSCONV,1),SIZE(XPRSCONV,2),IMODSON))
ELSE
ALLOCATE( ZSAVE_PRCONV(0,0,0))
ALLOCATE( ZSAVE_PRSCONV(0,0,0))
END IF
IF (CRAD /= 'NONE') THEN
ALLOCATE( ZSAVE_DIRFLASWD(SIZE(XDIRFLASWD,1),SIZE(XDIRFLASWD,2),SIZE(XDIRFLASWD,3),IMODSON))
ALLOCATE( ZSAVE_SCAFLASWD(SIZE(XSCAFLASWD,1),SIZE(XSCAFLASWD,2),SIZE(XSCAFLASWD,3),IMODSON))
ALLOCATE( ZSAVE_DIRSRFSWD(SIZE(XDIRSRFSWD,1),SIZE(XDIRSRFSWD,2),SIZE(XDIRSRFSWD,3),IMODSON))
ELSE
ALLOCATE( ZSAVE_DIRFLASWD(0,0,0,0))
ALLOCATE( ZSAVE_SCAFLASWD(0,0,0,0))
ALLOCATE( ZSAVE_DIRSRFSWD(0,0,0,0))
END IF
ENDIF
!
Gaelle Tanguy
committed
IKIDM=0
DO JKID = IMI+1,NMODEL ! min value of the possible kids
IF (IMI == NDAD(JKID) .AND. XWAY(JKID) == 2. .AND. CPROGRAM=='MESONH' &
.AND. (CCONF == 'RESTA' .OR. (CCONF == 'START' .AND. KTCOUNT /= 1))) THEN
Gaelle Tanguy
committed
! BUG if number of the son does not follow the number of the dad
! IKIDM = JKID-IMI
IKIDM = IKIDM + 1
IF (LUSERC .AND. ( &
(LSEDIC .AND. CCLOUD(1:3) == 'ICE') .OR. &
(LSEDC .AND. (CCLOUD == 'C2R2' .OR. CCLOUD == 'KHKO')) .OR. &
(MSEDC .AND. CCLOUD=='LIMA') &
)) THEN
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
ZSAVE_INPRC(:,:,IKIDM) = XINPRC(:,:)
END IF
IF (LUSERR) THEN
ZSAVE_INPRR(:,:,IKIDM) = XINPRR(:,:)
END IF
IF (LUSERS) THEN
ZSAVE_INPRS(:,:,IKIDM) = XINPRS(:,:)
END IF
IF (LUSERG) THEN
ZSAVE_INPRG(:,:,IKIDM) = XINPRG(:,:)
END IF
IF (LUSERH) THEN
ZSAVE_INPRH(:,:,IKIDM) = XINPRH(:,:)
END IF
IF (CDCONV /= 'NONE') THEN
ZSAVE_PRCONV(:,:,IKIDM) = XPRCONV(:,:)
ZSAVE_PRSCONV(:,:,IKIDM) = XPRSCONV(:,:)
END IF
IF (CRAD /= 'NONE') THEN
ZSAVE_DIRFLASWD(:,:,:,IKIDM) = XDIRFLASWD(:,:,:)
ZSAVE_SCAFLASWD(:,:,:,IKIDM) = XSCAFLASWD(:,:,:)
ZSAVE_DIRSRFSWD(:,:,:,IKIDM) = XDIRSRFSWD(:,:,:)
END IF
ENDIF
END DO
!
!-----------------------------------------------------------------------------
!
!* 1. RADIATION SCHEME
! ----------------
!
!
XTIME_BU_PROCESS = 0.
XTIME_LES_BU_PROCESS = 0.
!
CALL SECOND_MNH2(ZTIME1)
!
!
!* 1.1 Tests to control how the radiation package should be called (at the current timestep)
! -----------------------------------------------------------
!
!
GRAD = .FALSE.
GCLOUD_ONLY = .FALSE.
!
IF (CRAD /='NONE') THEN
!
! test to see if the partial radiations for cloudy must be called
!
IF (CRAD =='ECMW') THEN
CALL TEMPORAL_DIST(TDTCUR%TDATE%YEAR,TDTCUR%TDATE%MONTH, &
TDTCUR%TDATE%DAY, TDTCUR%TIME, &
TDTRAD_CLONLY%TDATE%YEAR,TDTRAD_CLONLY%TDATE%MONTH, &
TDTRAD_CLONLY%TDATE%DAY, TDTRAD_CLONLY%TIME, &
ZTEMP_DIST)
IF( MOD(NINT(ZTEMP_DIST/XTSTEP),NINT(XDTRAD_CLONLY/XTSTEP))==0 ) THEN
TDTRAD_CLONLY = TDTCUR
GRAD = .TRUE.
GCLOUD_ONLY = .TRUE.
END IF
END IF
!
! test to see if the full radiations must be called
!
CALL TEMPORAL_DIST(TDTCUR%TDATE%YEAR,TDTCUR%TDATE%MONTH, &
TDTCUR%TDATE%DAY, TDTCUR%TIME, &
TDTRAD_FULL%TDATE%YEAR,TDTRAD_FULL%TDATE%MONTH, &
TDTRAD_FULL%TDATE%DAY, TDTRAD_FULL%TIME, &
ZTEMP_DIST)
IF( MOD(NINT(ZTEMP_DIST/XTSTEP),NINT(XDTRAD/XTSTEP))==0 ) THEN
TDTRAD_FULL = TDTCUR
GRAD = .TRUE.
GCLOUD_ONLY = .FALSE.
END IF
!
! tests to see if any cloud exists
!
IF (CRAD =='ECMW') THEN
IF (GRAD .AND. NRR.LE.3 ) THEN
IF( MAXVAL(XCLDFR(:,:,:)).LE. 1.E-10 .AND. GCLOUD_ONLY ) THEN
GRAD = .FALSE. ! only the cloudy verticals would be
! refreshed but there is no clouds
END IF
END IF
!
IF (GRAD .AND. NRR.GE.4 ) THEN
IF( CCLOUD(1:3)=='ICE' )THEN
IF( MAXVAL(XRT(:,:,:,2)).LE.XRTMIN(2) .AND. &
MAXVAL(XRT(:,:,:,4)).LE.XRTMIN(4) .AND. GCLOUD_ONLY ) THEN
GRAD = .FALSE. ! only the cloudy verticals would be
! refreshed but there is no cloudwater and ice
END IF
END IF
IF( CCLOUD=='C3R5' )THEN
IF( MAXVAL(XRT(:,:,:,2)).LE.XRTMIN_C1R3(2) .AND. &
MAXVAL(XRT(:,:,:,4)).LE.XRTMIN_C1R3(4) .AND. GCLOUD_ONLY ) THEN
GRAD = .FALSE. ! only the cloudy verticals would be
! refreshed but there is no cloudwater and ice
END IF
END IF
IF( CCLOUD=='LIMA' )THEN
IF( MAXVAL(XRT(:,:,:,2)).LE.XRTMIN_LIMA(2) .AND. &
MAXVAL(XRT(:,:,:,4)).LE.XRTMIN_LIMA(4) .AND. GCLOUD_ONLY ) THEN
GRAD = .FALSE. ! only the cloudy verticals would be
! refreshed but there is no cloudwater and ice
END IF
END IF
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
END IF
END IF
!
END IF
!
! global parallel mask for 'GRAD'
ZRAD_GLOB_ll = 0.0
IF (GRAD) ZRAD_GLOB_ll = 1.0
CALL REDUCESUM_ll(ZRAD_GLOB_ll,INFO_ll)
if (ZRAD_GLOB_ll .NE. 0.0 ) GRAD = .TRUE.
!
!
IF( GRAD ) THEN
ALLOCATE(ZCOSZEN(IIU,IJU))
ALLOCATE(ZSINZEN(IIU,IJU))
ALLOCATE(ZAZIMSOL(IIU,IJU))
!
!
!* 1.2. Astronomical computations
! -------------------------
!
! Ajout PP
IF (.NOT. GCLOUD_ONLY .AND. KTCOUNT /= 1) THEN
IF (LAERO_FT) THEN
CALL AEROZON (XPABST,XTHT,XTSRAD,XLAT,XLON,TDTCUR,TDTEXP, &
NDLON,NFLEV,CAER,NAER,NSTATM, &
XSINDEL,XCOSDEL,XTSIDER,XCORSOL, &
XSTATM,XOZON, XAER)
END IF
END IF
!
CALL SUNPOS_n ( XZENITH, ZCOSZEN, ZSINZEN, ZAZIMSOL )
!
!* 1.3 Call to radiation scheme
! ------------------------
!
SELECT CASE ( CRAD )
!
!* 1.3.1 TOP of Atmposphere radiation
! ----------------------------
CASE('TOPA')
!
XFLALWD (:,:) = 300.
DO JSWB=1,NSWB_MNH
XDIRFLASWD(:,:,JSWB) = XI0 * MAX(COS(XZENITH(:,:)),0.)/FLOAT(NSWB_MNH)
XSCAFLASWD(:,:,JSWB) = 0.
END DO
XDTHRAD(:,:,:) = 0.
!
!* 1.3.1 FIXEd radiative surface fluxes
! ------------------------------
!
CASE('FIXE')
ZTIME = MOD(TDTCUR%TIME +XLON0*240., XDAY)
IHOUR = INT( ZTIME/3600. )
IF (IHOUR < 0) IHOUR=IHOUR + 24
ZDT = ZTIME/3600. - FLOAT(IHOUR)
XDIRFLASWD(:,:,:) =(( ZRG_HOUR(IHOUR+1)-ZRG_HOUR(IHOUR) )*ZDT + ZRG_HOUR(IHOUR)) / FLOAT(NSWB_MNH)
XFLALWD (:,:) = (ZRAT_HOUR(IHOUR+1)-ZRAT_HOUR(IHOUR))*ZDT + ZRAT_HOUR(IHOUR)
DO JSWB=1,NSWB_MNH
WHERE(ZCOSZEN(:,:)<0.) XDIRFLASWD(:,:,JSWB) = 0.
END DO
XSCAFLASWD(:,:,:) = XDIRFLASWD(:,:,:) * 0.2
XDIRFLASWD(:,:,:) = XDIRFLASWD(:,:,:) * 0.8
XDTHRAD(:,:,:) = 0.
!
!
!* 1.3.2 ECMWf radiative surface and atmospheric fluxes
! ----------------------------------------------
!
CASE('ECMW')
IF (LLES_MEAN) GCLOUD_ONLY=.FALSE.
XRADEFF(:,:,:)=0.0
XSWU(:,:,:)=0.0
XSWD(:,:,:)=0.0
XLWU(:,:,:)=0.0
XLWD(:,:,:)=0.0
XDTHRADSW(:,:,:)=0.0
XDTHRADLW(:,:,:)=0.0
CALL RADIATIONS ( OCLOSE_OUT, HFMFILE, CLUOUT, &
LCLEAR_SKY,GCLOUD_ONLY, NCLEARCOL_TM1,CEFRADL, CEFRADI,COPWSW,COPISW,&
COPWLW,COPILW, XFUDG, &
NDLON, NFLEV, NRAD_DIAG, NFLUX, NRAD, NAER,NSWB, NSTATM, NRAD_COLNBR,&
ZCOSZEN, XSEA, XCORSOL, &
XDIR_ALB, XSCA_ALB, XEMIS, XCLDFR, XCCO2, XTSRAD, XSTATM, XTHT, XRT, &
XPABST,XOZON, XAER,XDST_WL, XAER_CLIM, XSVT, &
XDTHRAD, XFLALWD, XDIRFLASWD, XSCAFLASWD, XRHODREF, XZZ , &
XRADEFF, XSWU, XSWD, XLWU, XLWD, XDTHRADSW, XDTHRADLW )
!
WRITE(UNIT=ILUOUT,FMT='(" RADIATIONS called for KTCOUNT=",I6, &
& "with the CLOUD_ONLY option set ",L2)') KTCOUNT,GCLOUD_ONLY
!
WHERE( XDIRFLASWD(:,:,1) + XSCAFLASWD(:,:,1) >0. )
XALBUV(:,:) = ( XDIR_ALB(:,:,1) * XDIRFLASWD(:,:,1) &
+ XSCA_ALB(:,:,1) * XSCAFLASWD(:,:,1) ) &
/ (XDIRFLASWD(:,:,1) + XSCAFLASWD(:,:,1) )
ELSEWHERE
XALBUV(:,:) = XDIR_ALB(:,:,1)
END WHERE
!
END SELECT
!
CALL SECOND_MNH2(ZTIME2)
!
PRAD = PRAD + ZTIME2 - ZTIME1
!
ZTIME1 = ZTIME2
!
CALL SURF_RAD_MODIF (XMAP, XXHAT, XYHAT, &
ZCOSZEN, ZSINZEN, ZAZIMSOL, XZS, XZS_XY, &
XDIRFLASWD, XDIRSRFSWD )
!
!* Azimuthal angle to be sent later to surface processes
! Defined in radian, clockwise, from North
!
XAZIM = ZAZIMSOL
!
CALL SECOND_MNH2(ZTIME2)
!
PSHADOWS = PSHADOWS + ZTIME2 - ZTIME1
!
ZTIME1 = ZTIME2
!
DEALLOCATE(ZCOSZEN)
DEALLOCATE(ZSINZEN)
DEALLOCATE(ZAZIMSOL)
!
END IF
!
!
!* 1.4 control prints
! --------------
!
!* 1.5 Radiative tendency integration
! ------------------------------
!
IF (CRAD /='NONE') THEN
XRTHS(:,:,:) = XRTHS(:,:,:) + XRHODJ(:,:,:)*XDTHRAD(:,:,:)
END IF
!
!* 1.6 budget storage
! --------------
!
IF (CRAD/='NONE' .AND. LBUDGET_TH) CALL BUDGET (XRTHS,4,'RAD_BU_RTH')
!
CALL SECOND_MNH2(ZTIME2)
!
PRAD = PRAD + ZTIME2 - ZTIME1 &
- XTIME_LES_BU_PROCESS - XTIME_BU_PROCESS
!
PTIME_BU = PTIME_BU + XTIME_LES_BU_PROCESS + XTIME_BU_PROCESS
!
!
!-----------------------------------------------------------------------------
!
!* 2. DEEP CONVECTION SCHEME
! ----------------------
!
ZTIME1 = ZTIME2
XTIME_BU_PROCESS = 0.
XTIME_LES_BU_PROCESS = 0.
!
CALL SECOND_MNH2(ZTIME1)
!
IF( CDCONV /= 'NONE' .OR. CSCONV == 'KAFR' ) THEN
!
! test to see if the deep convection scheme should be called
!
GDCONV = .FALSE.
!
CALL TEMPORAL_DIST(TDTCUR%TDATE%YEAR,TDTCUR%TDATE%MONTH, &
TDTCUR%TDATE%DAY, TDTCUR%TIME, &
TDTDCONV%TDATE%YEAR,TDTDCONV%TDATE%MONTH, &
TDTDCONV%TDATE%DAY, TDTDCONV%TIME, &
ZTEMP_DIST)
IF( MOD(NINT(ZTEMP_DIST/XTSTEP),NINT(XDTCONV/XTSTEP))==0 ) THEN
TDTDCONV = TDTCUR
GDCONV = .TRUE.
END IF
!
IF( GDCONV ) THEN
IF (CDCONV == 'KAFR' .OR. CSCONV == 'KAFR' ) THEN
ALLOCATE( ZRC(IIU,IJU,IKU) )
ALLOCATE( ZRI(IIU,IJU,IKU) )
ALLOCATE( ZWT(IIU,IJU,IKU) )
ALLOCATE( ZDXDY(IIU,IJU) )
! Compute grid area
ZDXDY(:,:) = SPREAD(XDXHAT(1:IIU),2,IJU) * SPREAD(XDYHAT(1:IJU),1,IIU)
!
IF( LUSERC .AND. LUSERI ) THEN
ZRC(:,:,:) = XRT(:,:,:,2)
ZRI(:,:,:) = XRT(:,:,:,4)
ELSE IF( LUSERC .AND. (.NOT. LUSERI) ) THEN
ZRC(:,:,:) = XRT(:,:,:,2)
ZRI(:,:,:) = 0.0
ELSE
ZRC(:,:,:) = 0.0
ZRI(:,:,:) = 0.0
END IF
WRITE(UNIT=ILUOUT,FMT='(" CONVECTION called for KTCOUNT=",I6)') &
KTCOUNT
IF ( LFORCING .AND. L1D ) THEN
ZWT(:,:,:) = XWTFRC(:,:,:)
ELSE
ZWT(:,:,:) = XWT(:,:,:)
ENDIF
IF (LDUST) CALL DUST_FILTER(XSVT(:,:,:,NSV_DSTBEG:NSV_DSTEND), XRHODREF(:,:,:))
IF (LSALT) CALL SALT_FILTER(XSVT(:,:,:,NSV_SLTBEG:NSV_SLTEND), XRHODREF(:,:,:))
IF (LCH_CONV_LINOX) THEN
CALL CONVECTION( XDTCONV, CDCONV, CSCONV, LREFRESH_ALL, LDOWN, NICE, &
LSETTADJ, XTADJD, XTADJS, LDIAGCONV, NENSM, &
XPABST, XZZ, ZDXDY, &
XTHT, XRT(:,:,:,1), ZRC, ZRI, XUT, XVT, &
ZWT,XTKET(:,:,IKB), &
NCOUNTCONV, XDTHCONV, XDRVCONV, XDRCCONV, XDRICONV, &
XPRCONV, XPRSCONV, &
XUMFCONV,XDMFCONV,XMFCONV,XPRLFLXCONV,XPRSFLXCONV, &
XCAPE, NCLTOPCONV, NCLBASCONV, &
LCHTRANS, XSVT, XDSVCONV, &
LUSECHEM, LCH_CONV_SCAV, LCH_CONV_LINOX, &
LDUST, LSALT, &
XRHODREF, XIC_RATE, XCG_RATE )
ELSE
CALL CONVECTION( XDTCONV, CDCONV, CSCONV, LREFRESH_ALL, LDOWN, NICE, &
LSETTADJ, XTADJD, XTADJS, LDIAGCONV, NENSM, &
XPABST, XZZ, ZDXDY, &
XTHT, XRT(:,:,:,1), ZRC, ZRI, XUT, XVT, &
ZWT,XTKET(:,:,IKB), &
NCOUNTCONV, XDTHCONV, XDRVCONV, XDRCCONV, XDRICONV, &
XPRCONV, XPRSCONV, &
XUMFCONV,XDMFCONV,XMFCONV,XPRLFLXCONV,XPRSFLXCONV, &
XCAPE, NCLTOPCONV, NCLBASCONV, &
LCHTRANS, XSVT, XDSVCONV, &
LUSECHEM, LCH_CONV_SCAV, LCH_CONV_LINOX, &
LDUST, LSALT, &
XRHODREF )
END IF
!
DEALLOCATE( ZRC )
DEALLOCATE( ZRI )
DEALLOCATE( ZWT )
DEALLOCATE( ZDXDY )
END IF
END IF
!
! Deep convection tendency integration
!
XRTHS(:,:,:) = XRTHS(:,:,:) + XRHODJ(:,:,:) * XDTHCONV(:,:,:)
XRRS(:,:,:,1) = XRRS(:,:,:,1) + XRHODJ(:,:,:) * XDRVCONV(:,:,:)
!
!
! Aerosols size distribution
! Compute Rg and sigma before tracers convection tendency (for orilam, dust and sea
! salt)
!
IF ( LCHTRANS ) THEN ! update tracers for chemical transport
IF (LORILAM) ZRSVS(:,:,:,:) = XRSVS(:,:,:,:) !
IF ((LDUST)) THEN ! dust convective balance
ALLOCATE(ZSIGDST(IIU,IJU,IKU,NMODE_DST))
ALLOCATE(ZRGDST(IIU,IJU,IKU,NMODE_DST))
ALLOCATE(ZNDST(IIU,IJU,IKU,NMODE_DST))
ALLOCATE(ZSVDST(IIU,IJU,IKU,NSV_DST))
!
DO JSV=1,NMODE_DST
IMODEIDX = JPDUSTORDER(JSV)
IF (CRGUNITD=="MASS") THEN
ZINIRADIUS(JSV) = XINIRADIUS(IMODEIDX) * EXP(-3.*(LOG(XINISIG(IMODEIDX)))**2)
ELSE
ZINIRADIUS(JSV) = XINIRADIUS(IMODEIDX)
END IF
ZSIGDST(:,:,:,JSV) = XINISIG(IMODEIDX)
ZRGDST(:,:,:,JSV) = ZINIRADIUS(JSV)
ZNDST(:,:,:,JSV) = XN0MIN(IMODEIDX)
ENDDO
!
DO JSV=NSV_DSTBEG,NSV_DSTEND
ZSVDST(:,:,:,JSV-NSV_DSTBEG+1) = XRSVS(:,:,:,JSV) * XTSTEP / XRHODJ(:,:,:)
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
ENDDO
CALL PPP2DUST(ZSVDST(IIB:IIE,IJB:IJE,IKB:IKE,:), XRHODREF(IIB:IIE,IJB:IJE,IKB:IKE),&
PSIG3D=ZSIGDST(IIB:IIE,IJB:IJE,IKB:IKE,:), PRG3D=ZRGDST(IIB:IIE,IJB:IJE,IKB:IKE,:), &
PN3D=ZNDST(IIB:IIE,IJB:IJE,IKB:IKE,:))
END IF
!
IF ((LSALT)) THEN ! sea salt convective balance
ALLOCATE(ZSIGSLT(IIU,IJU,IKU,NMODE_SLT))
ALLOCATE(ZRGSLT(IIU,IJU,IKU,NMODE_SLT))
ALLOCATE(ZNSLT(IIU,IJU,IKU,NMODE_SLT))
ALLOCATE(ZSVSLT(IIU,IJU,IKU,NSV_SLT))
!
DO JSV=1,NMODE_SLT
IMODEIDX = JPSALTORDER(JSV)
IF (CRGUNITS=="MASS") THEN
ZINIRADIUS_SLT(JSV) = XINIRADIUS_SLT(IMODEIDX) * &
EXP(-3.*(LOG(XINISIG_SLT(IMODEIDX)))**2)
ELSE
ZINIRADIUS_SLT(JSV) = XINIRADIUS_SLT(IMODEIDX)
END IF
ZSIGSLT(:,:,:,JSV) = XINISIG_SLT(IMODEIDX)
ZRGSLT(:,:,:,JSV) = ZINIRADIUS_SLT(JSV)
ZNSLT(:,:,:,JSV) = XN0MIN_SLT(IMODEIDX)
ENDDO
!
DO JSV=NSV_SLTBEG,NSV_SLTEND
ZSVSLT(:,:,:,JSV-NSV_SLTBEG+1) = XRSVS(:,:,:,JSV) * XTSTEP / XRHODJ(:,:,:)
ENDDO
CALL PPP2SALT(ZSVSLT(IIB:IIE,IJB:IJE,IKB:IKE,:), XRHODREF(IIB:IIE,IJB:IJE,IKB:IKE),&
PSIG3D=ZSIGSLT(IIB:IIE,IJB:IJE,IKB:IKE,:), PRG3D=ZRGSLT(IIB:IIE,IJB:IJE,IKB:IKE,:), &
PN3D=ZNSLT(IIB:IIE,IJB:IJE,IKB:IKE,:))
END IF
!
!
! Compute convective tendency for all tracers
!
IF (LCHTRANS) THEN
DO JSV = 1, SIZE(XRSVS,4)
XRSVS(:,:,:,JSV) = XRSVS(:,:,:,JSV) + XRHODJ(:,:,:) * XDSVCONV(:,:,:,JSV)
END DO
IF (LORILAM) THEN
DO JSV = NSV_AERBEG,NSV_AEREND
ZWETDEPAER(:,:,:,JSV-NSV_AERBEG+1) = XDSVCONV(:,:,:,JSV) * XRHODJ(:,:,:)
XRSVS(:,:,:,JSV) = ZRSVS(:,:,:,JSV)
END DO
END IF
END IF
!
IF ((LDUST).AND.(LCHTRANS)) THEN ! dust convective balance
IF (CPROGRAM == "MESONH") THEN
DO JSV=NSV_DSTBEG,NSV_DSTEND
ZSVDST(:,:,:,JSV-NSV_DSTBEG+1) = XRSVS(:,:,:,JSV) * XTSTEP / XRHODJ(:,:,:)
ENDDO
ELSE
DO JSV=NSV_DSTBEG,NSV_DSTEND
ZSVDST(:,:,:,JSV-NSV_DSTBEG+1) = XSVT(:,:,:,JSV)
ENDDO
ENDIF
CALL DUST2PPP(ZSVDST(IIB:IIE,IJB:IJE,IKB:IKE,:), &
XRHODREF(IIB:IIE,IJB:IJE,IKB:IKE), ZSIGDST(IIB:IIE,IJB:IJE,IKB:IKE,:),&
ZRGDST(IIB:IIE,IJB:IJE,IKB:IKE,:))
DO JSV=NSV_DSTBEG,NSV_DSTEND
XRSVS(:,:,:,JSV) = ZSVDST(:,:,:,JSV-NSV_DSTBEG+1) * XRHODJ(:,:,:) / XTSTEP
ENDDO
!
DEALLOCATE(ZSVDST)
DEALLOCATE(ZNDST)
DEALLOCATE(ZRGDST)
DEALLOCATE(ZSIGDST)
END IF
!
IF ((LSALT).AND.(LCHTRANS)) THEN ! sea salt convective balance
IF (CPROGRAM == "MESONH") THEN
DO JSV=NSV_SLTBEG,NSV_SLTEND
ZSVSLT(:,:,:,JSV-NSV_SLTBEG+1) = XRSVS(:,:,:,JSV) * XTSTEP / XRHODJ(:,:,:)
ENDDO
ELSE
DO JSV=NSV_SLTBEG,NSV_SLTEND