Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
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
217
218
219
220
221
222
223
224
225
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
283
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
326
327
328
329
330
331
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
403
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
476
477
478
479
480
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
523
524
525
526
527
528
529
530
531
532
533
534
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
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
!OPTION! -pvctl noloopfusion
SUBROUTINE APL_AROME(YDGEOMETRY,YDSURF, YDCFU, YDXFU, YDMODEL, KBL, KGPCOMP, KIDIA , KFDIA , KLON ,&
& KTDIA , KLEV , KSTEP ,&
& KMAXDRAFT, KSGST, KNFRRC, PDT, LDXFUMSE, PINDX, PINDY ,&
& PGEMU,PGELAM,POROG,PGM,PMU0,PMU0LU,PMU0M,PMU0N,PCLON, PSLON,PVO3ABC,PLSM,&
& PAESEA , PAELAN , PAESOO , PAEDES , PAESUL, PAEVOL,&
& PGP2DSDT, PGP2DSPP, &
!---------------------------------------------------------------------
! - INPUT A M
& PAPHIM,PAPHIFM,PAPRSM, PAPRSFM, PRDELPM, PDELPM, PTM, PQVM ,&
& PCPM , PRM ,PALPHM , PLNPRM,&
& PQCM , PQIM ,PQRM , PQSM, PQGM, PQHM,&
& PLIMAM , &
& PTKEM , PEFB1 ,PEFB2 , PEFB3,&
& PSIGM,PSVM,&
& PUM , PVM, PWM, PEDR,&
& PFORCEU,PFORCEV,PFORCET,PFORCEQ,&
!---------------------------------------------------------------------
! - INOUT A S
& PGPAR, PEMTD, PEMTU, PTRSO,&
& PGDEOSI, PGUEOSI, PGMU0, PGMU0_MIN, PGMU0_MAX,&
& PGDEOTI, PGDEOTI2, PGUEOTI, PGUEOTI2, PGEOLT, PGEOXT,&
& PGRPROX, PGMIXP, PGFLUXC, PGRSURF,&
& PTURB3D,&
! - OUT A S
& PQLRAD, PQIRAD, PRH, PCLFS, PSIGS,&
& PTENDT, PTENDR, PTENDU, PTENDV,PTENDW,&
& PTENDLIMA, &
& PTENDTKE, PTENDEFB1, PTENDEFB2, PTENDEFB3,&
& PTENDEXT,PFRTH, PFRSO,PFRTHDS, PFRSODS, PFRSOPS, PFRSDNI,&
& PFRSOPT, PFRTHC, PFRSOC, &
!---------------------------------------------------------------------
! - IN FOR RADIATION IF NO SURFACE SCHEME
& PALBIN,PEMIS,&
! - INOUT for easy diag
& PEZDIAG,&
! - INOUT for CFU XFU
& PCLCH,PCLCL,PCLCM,PCLCT,PFPLSL,PFPLSN,PFPLSG,PFPLSH,PSTRTU,PSTRTV,PFCS,PFCLL,&
& PFCLN,PUCLS,PVCLS,PNUCLS,PNVCLS,PTCLS,PQCLS,PHUCLS,PUGST,PVGST,PFEVL,PFEVN, PPBLH,PSPSG,PSPSGR,&
& PSDUR,PDIAGH,PFLASH,PSFORC,PTPWCLS,PDPRECIPS,PDPRECIPS2,PVISICLD,PVISIHYDRO,PMXCLWC,&
! daand: radflex
& YDPROCSET ,YDDDH)
!**** *APL_AROME * - CALL OF PHYSICAL PARAMETERISATIONS FOR ALARO/AROME
! Sujet.
! ------
! - APPEL DES SOUS-PROGRAMMES DE PARAMETRISATION
!** Interface.
! ----------
! *CALL* *APL_AROME*
!-----------------------------------------------------------------------
! - ARGUMENTS D'ENTREE.
! - INPUT ARGUMENTS.
! -------------------
! - NOM DES PARAMETRES DE DIMENSIONNEMENT DE LA PHYSIQUE.
! - DIMENSIONS.
! KBL : NUMERO DE BLOC NPROMA
! KBL : NPROMA-PACKETS NUMBER
! KGPCOMP : NOMBRE TOTAL DE POINTS DE GRILLE SUR LE DOMAINE
! KGPCOMP : TOTAL GRID POINTS NUMBER IN THE DOMAIN
! KIDIA, KFDIA : BORNES BOUCLES HORIZONTALES (IST,IEND DANS CPG).
! KIDIA, KFDIA : START/END OF HORIZONTAL LOOP (IST,IEND IN *CPG*).
! KLON : DIMENSION HORIZONTALE (NPROMA DANS CPG).
! KLON : HORIZONTAL DIMENSION (NPROMA IN *CPG*).
! KTDIA : DEBUT BOUCLE VERTICALE DANS LA PHYSIQUE.
! KTDIA : START OF THE VERTICAL LOOP IN THE PHYSICS (IF SOME LEVELS ARE
! SKIPPED AT THE TOP OF THE MODEL).
! KLEV : FIN BOUCLE VERTICE ET DIMENSION VERTICALE (NFLEVG DANS CPG).
! KLEV : END OF VERTICAL LOOP AND VERTICAL DIMENSION(NFLEVG IN *CPG*).
! KSTEP : TIME STEP NUMBER (starting with zero)
! KMAXDRAFT : MAX NUMBER OF DRAFTS (FOR DIMENSIONNING)
! KSGST : NUMBER OF SUBGRID SURFACE TEMPERATURES AND FLUXES (NTSSG IN *CPG*)
! KNFRRC : FREQUENCY FOR CLEAR SKY RADIATION CALCULATION
! PDT : TIME STEP (in s)
! LDXFUMSE : T if CDCONF=X in order not to increment surfex timer in that case
!-----------------------------------------------------------------------
! PGEMU : SINE OF GEOGRAPHICAL LATITUDE
! PGELAM : LONGITUDE
! POROG : g * OROGRAPHY
! PGM : MAP FACTOR (used in ALARO convection only)
! PMU0 : COSINUS LOCAL INSTANTANE DE L'ANGLE ZENITHAL SOLAIRE.
! PMU0LU : COSINUS LOCAL INSTANTANE DE L'ANGLE ZENITHAL LUNAIRE.
! PMU0 : LOCAL COSINE OF INSTANTANEOUS SOLAR ZENITH ANGLE.
! PMU0M : COSINUS LOCAL MOYEN DE L'ANGLE ZENITHAL.
! PMU0M : LOCAL COSINE OF AVERAGED SOLAR ZENITH ANGLE.
! PMU0N : COSINUS LOCAL AU PAS DE TEMPS SUIVANT DE L'ANGLE ZENITHAL SOLAIRE.
! PMU0N : NEXT TIME STEP COSINUS LOCAL INSTANTANE DE L'ANGLE ZENITHAL SOLAIRE.
! PCLON : cosine of geographical longitude.
! PSLON : sine of geographical longitude.
! PVO3ABC : OZONE COEFFICIENTS FOR ARPEGE PROFILES
! PLSM : -ATMOSPHERIC MODEL- LAND-SEA MASK (! MAY BE DIFFERENT FROM
! THE SURFACE ONE)
! PAESEA : MARINE AEROSOLS (IF NVCLIA >= 4)
! PAELAN : CONTINENTAL AEROSOLS (IF NVCLIA >= 4)
! PAESOO : SOOT AEROSOLS (IF NVCLIA >= 4)
! PAEDES : DESERT AEROSOLS (IF NVCLIA >= 4)
! PAESUL : SULFATE AEROSOLS (IF LAEROSUL=.T.)
! PAEVOL : VOLCANO AEROSOLS (IF LAEROVOL=.T.)
! PGP2DSDT : STOCHASTIC PHYSICS PATTERNS
! FIELDS WITH SUBSCRIPT M FOR TIME T-DT IN 3TL OR T IN 2TL
! PAPHIM : GEOPOTENTIAL ON HALF-LEVELS
! PAPHIFM : GEOPOTENTIAL ON FULL-LEVELS
! PAPRSM : PRESSURE ON HALF LEVELS
! PAPRSFM : PRESSURE ON FULL LEVELS.
! PRDELPM : INVERSE OF PDELP
! PDELPM : LAYER THICKNESS IN PRESSURE UNITS
! PTM : TEMPERATURE.
! PQVM : SPECIFIC HUMIDITY OF WATER VAPOR
! PCPM : SPECIFIC HEAT AT CONSTANT PRESSURE FOR AIR
! PRM : GAS CONSTANT FOR AIR
! PALPHM : "alpha" on layers
! PLNPRM : "delta" on layers
! PQCM : SPECIFIC HUMIDITY OF CLOUD WATER
! PQIM : SPECIFIC HUMIDITY OF ICE
! PQRM : SPECIFIC HUMIDITY OF RAIN
! PQSM : SPECIFIC HUMIDITY OF SNOW
! PQGM : SPECIFIC HUMIDITY OF GRAUPEL
! PQHM : SPECIFIC HUMIDITY OF HAIL
! PTKEM : TURBULENT KINETIC ENERGY
! PSVM : PASSIVE SCALARS
! PSIGM : SIGMA FOR SUBGRIDCOND
! PUM : ZONAL WIND
! PVM : MERIDIAN WIND
! PWM : VERTICAL VELOCITY (m/s)
!-----------------------------------------------------------------------
! - INOUT
! PGPAR : BUFFER FOR 2D FIELDS - CONTAINS PRECIP, ALBEDO, EMISS, TS
! : SURFACE FLUXES
! PEMTD : DOWNWARD LONGWAVE EMISSIVITY
! PEMTU : UPWARD LONGWAVE EMISSIVITY
! PTRSO : SHORTWAVE TRANSMISSIVITY
! ACRANEB2 intermittency storage
! PGDEOSI : DESCENDING INCREMENTAL OPTICAL DEPTHS, SOLAR
! PGUEOSI : ASCENDING INCREMENTAL OPTICAL DEPTHS, SOLAR
! PGMU0 : COSINE OF SOLAR ZENITH ANGLE, APPROXIMATE ACTUAL VALUE
! PGMU0_MIN : COSINE OF SOLAR ZENITH ANGLE, MIN VALUE
! PGMU0_MAX : COSINE OF SOLAR ZENITH ANGLE, MAX VALUE
! PGDEOTI : DESCENDING INCREMENTAL OPTICAL DEPTHS, dB/dT(T0) WEIGHTS
! PGDEOTI2 : DESCENDING INCREMENTAL OPTICAL DEPTHS, B WEIGHTS WITH
! LINEAR T_e CORRECTION
! PGUEOTI : ASCENDING INCREMENTAL OPTICAL DEPTHS, dB/dT(T0) WEIGHTS
! PGUEOTI2 : ASCENDING INCREMENTAL OPTICAL DEPTHS, B WEIGHTS WITH
! LINEAR T_e CORRECTION
! PGEOLT : LOCAL OPTICAL DEPTHS, dB/dT(T0) WEIGHTS
! PGEOXT : MAXIMUM OPTICAL DEPTHS FOR EBL-EAL, dB/dT(T0) WEIGHTS
! PGRPROX : CORRECTION TERM FOR ADJACENT EXCHANGES
! PGMIXP : NON-STATISTICAL WEIGHTS FOR BRACKETING
! PGFLUXC : OUT OF BRACKET PART OF CLEARSKY EBL, RESP. EBL-EAL FLUX
! PGRSURF : CORRECTIVE RATIO FOR SURFACE CTS CONTRIBUTION
! PTURB3D : MATRICE DE GRADIENTS HORIZONTAUX
!-----------------------------------------------------------------------
! - OUTPUT (SUBSCRIPT S FOR T+DT)
! PCLFS : CLOUD FRACTION
! PQLRAD : SPECIFIC HUMIDITY OF CLOUD WATER FOR RTTOV
! PQIRAD : SPECIFIC HUMIDITY OF ICE FOR RTTOV
! PSIGS : SIGMA FOR SUBGRIDCOND
! PTENDT : TEMPERATURE TENDENCY
! PTENDR : HYDROMETEORE TENDENCIES
! PTENDU : ZONAL WIND TENDENCY
! PTENDV : MERIDIAN WIND TENDENCY
! PTENDW : VERTICAL VELOCITY TENDENCY
! PTENDTKE : TKE TENDENCY
! PTENDEXT : PASSIVE SCALARS TENDENCY
! PFRTH : LONGWAVE RADIATIVE FLUX
! PFRSO : SHORTWAVE RADIATIVE FLUX
! PFRTHDS : LONGWAVE DOWNWARD SURFACE RADIATIVE FLUX
! PFRSOPS : SHORTWAVE DOWNWARD SURFACE RADIATIVE FLUX DIRECT
! PFRSDNI : SHORTWAVE DIRECT NORMAL IRRADIANCE
! PFRSODS : SHORTWAVE DOWNWARD SURFACE RADIATIVE FLUX GLOBAL
! PFRSOPT : SHORTWAVE DOWNWARD TOP RADIATIVE FLUX DIRECT
! - 2D (0:1)
! PFRTHC : LONGWAVE CLEAR SKY NET RADIATIVE FLUX
! PFRSOC : SHORTWAVE CLEAR SKY NET RADIATIVE FLUX
! variables used in input for radiation in case no surface scheme is used
! PALBIN : MODEL SURFACE SHORTWAVE ALBEDO.
! PEMIS : MODEL SURFACE LONGWAVE EMISSIVITY.
! Part of GFL strcture dedicated to easy diagnostics (to be used as a print...)
! PEZDIAG : MULPITPLE ARRAY TO BE FILLED BY THE USER BY 3D FIELDS
! (NGFL_EZDIAG ONES)
! output for CFU XFU
! PCLCH : HIGH CLOUD COVER (DIAGNOSTIC).
! PCLCL : LOW CLOUD COVER (DIAGNOSTIC).
! PCLCM : MEDIUM CLOUD COVER (DIAGNOSTIC).
! PCLCT : TOTAL CLOUD COVER (DIAGNOSTIC).
! PFPLSL : RESOLVED PRECIPITATION AS RAIN.
! PFPLSN : RESOLVED PRECIPITATION AS SNOW
! PFPLSG : RESOLVED PRECIPITATION AS GRAUPEL
! PFPLSH : RESOLVED PRECIPITATION AS HAIL
! PSTRTU : TURBULENT FLUX OF MOMENTUM "U".
! PSTRTV : TURBULENT FLUX OF MOMENTUM "V".
! PFCS : SENSIBLE HEAT FLUX AT SURFACE LEVEL.
! PFCLL : LATENT HEAT FLUX AT SURFACE LEVEL OVER WATER.
! PFCLN : LATENT HEAT FLUX AT SURFACE LEVEL OVER SNOW.
! PUCLS : SORTIE DIAGNOSTIQUE DU VENT EN X A HUV METEO.
! PUCLS : U-COMPONENT OF WIND AT 10 METERS (DIAGNOSTIC).
! PVCLS : SORTIE DIAGNOSTIQUE DU VENT EN Y A HUV METEO.
! PVCLS : V-COMPONENT OF WIND AT 10 METERS (DIAGNOSTIC).
! PNUCLS : SORTIE DIAGNOSTIQUE DU VENT NEUTRE EN X A HUV METEO.
! PNUCLS : U-COMPONENT OF NEUTRAL WIND AT 10 METERS (DIAGNOSTIC).
! PNVCLS : SORTIE DIAGNOSTIQUE DU VENT NEUTRE EN Y A HUV METEO.
! PNVCLS : V-COMPONENT OF NEUTRAL WIND AT 10 METERS (DIAGNOSTIC).
! PTCLS : SORTIE DIAGNOSTIQUE DE LA TEMPERATURE A HTQ METEO.
! PTCLS : TEMPERATURE AT 2 METERS (DIAGNOSTIC).
! PQCLS : SORTIE DIAGNOSTIQUE DE L'HUMIDITE SPECIFIQUE A HTQ METEO.
! PQCLS : SPECIFIC HUMIDITY AT 2 METERS (DIAGNOSTIC).
! PHUCLS : SORTIE DIAGNOSTIQUE DE L'HUMIDITE RELATIVE A HTQ METEO.
! PHUCLS : RELATIVE HUMIDITY AT 2 METERS (DIAGNOSTIC).
! PUGST : SORTIE DIAGNOSTIQUE DU VENT RAFALE EN X A HUV METEO.
! PUGST : U-COMPONENT OF WIND GUST AT 10 METERS (DIAGNOSTIC).
! PVGST : SORTIE DIAGNOSTIQUE DU VENT RAFALE EN Y A HUV METEO.
! PVGST : V-COMPONENT OF WIND GUST AT 10 METERS (DIAGNOSTIC).
! PFEVL : FLUX DE VAPEUR D'EAU SUR EAU LIQUIDE (OU SOL HUMIDE).
! PFEVL : WATER VAPOUR FLUX OVER LIQUID WATER (OR WET SOIL)
! PFEVN : FLUX DE VAPEUR D'EAU SUR NEIGE (OU GLACE) ET SOL GELE.
! PFEVN : WATER VAPOUR FLUX OVER SNOW (OR ICE) AND FROZEN SOIL.
! PPBLH : PSEUDO-HISTORICAL ARRAY FOR PBL HEIGHT
! PSPSG : SNOW COVER
! PSPSGR : SNOW DENSITY
! PSDUR : SUNSHINE DURATION [s]
! PDIAGH : HAIL DIAGNOSTIC
! PFLASH : LIGHTNING DIAGNOSTICS
! PVISICLD : VISIBILITY DUE TO CLOUD WATER AND CLOUD ICE
! PVISIHYDRO : VISIBILITY DUE TO RAIN AND SNOW AND GRAUPEL
! PMXCLWC : CLOUD WATER LIQUID CONTENT AT HVISI METERS
! PDPRECIPS : PRECIPITATION TYPE
! PDPRECIPS2 : PRECIPITATION TYPE FOR 2NDE PERIOD
!-----------------------------------------------------------------------
! Externes.
! ---------
! Method
! ------
! - convert aladin variables into mesonh variables (level inversion
! and q to r, t to theta)
! - call mesoNH physics and ECMWF radiation scheme
! - convert mesoNH tendencies to aladin tendencies
! Auteur.
! -------
! S.Malardel et Y. Seity
! 10-03-03
! big cleaning (18/06/04) S. Malardel and Y. Seity
! externalisation of surface scheme call + small cleaning (20-07-04) Y.Seity
! Modifications
! -------------
! G. Hello 04-02-06: Add the call of KFB-convection scheme
! for future use in ALARO
! T.Kovacic 04-05-05: Added ZCVTENDPR_ and ZCVTENDPRS_
! M.Hamrud 01-Oct-2003 CY28 Cleaning
! F.Bouyssel 04-05-05: New arguments in ACRADIN
! Y. Seity 30-Sept-2005 Add MNH Chemistry scheme
! R. Zaaboul 15-feb-2006 add surface scheme call
! T.Kovacic 2006-03-23: calls to subroutines for budgets
! and new arguments PFRTH and PFRSO
! Y. Seity 2007-05-07: add CFU and XFU calculations
! and call aro_ground_diag
! S.Ivatek-S 2007-04-17: Over dimensioning of PGPAR by NGPAR+1 just
! (KLON,NGPAR) is used boundary checking bf
! T.Kovacic 2007-03-16: Fourth dim. in APFT
! JJMorcrette, ECMWF, 20080325: dummy arguments for RADACT to allow for
! using a new sulphate climatology in the ECMWF model
! Y. Seity 2008-06-15: correct calculations of PFRTHDS, PFRSODS and PFCLL
! Y. Seity 2008-09-29: phasing Chemistry corrections
! O.Riviere 2008-10-01: introduction of new data flow for DDH in Arome
! Y. Seity 2009-05-03: new version of EDKF and implementation of EDMF
! Y. Seity 2009-10-03: add missed deallocations
! S. Riette 2009-03-25: Arguments modification for AROCLDIA to add HTKERAF
! K. Yessad (Jul 2009): remove CDLOCK + some cleanings
! A. Alias 2009-09-01: Sulfate and Volcano aerosols added (call radaer)
! S. Riette 2010-01-19: ZUM__, ZVM__ and ZDEPTH_HEIGHT_ are given
! ARO_GROUND_DIAG in 3D.
! Y. Seity 2010-03-09: add PFEVN and PFEVL
! Y. Bouteloup 2010-03-26 : Add PQLRAD et PQIRAD
! Y. Seity : Test TKE > 0.
! Y. Seity : Optimized version of EDKF + diag HCLS
! Y. Seity : 2010-09 Save Ts at the end of apl_arome for ICMSH+0000
! L. Bengtsson (2010): Introduce cloud diagnostics based on geop.
! height (LWMOCLOUD), AND cloud-overlap assumptions
! from C. Wittman 2009 (LACPANMX + WMXOV)
! S. Riette: 2010-12 aro_ground_diag interface modified
! Y. Seity: 2010-12 add hail diagnostic
! R. El Khatib 30-Jun-2010 NEC directive noloopfusion to preserve critical regions
! P.Marguinaud 2010-06-29 : KSURFEXCTL flag (disable SURFEX)
! 2010-12 B. Decharme : modify the radiative coupling with surfex (SW per band in ACRADIN and RADHEAT)
! 2011-02 A. Voldoire : add ZAERINDS to CALL RADAER and ACRADIN
! for sulfate indirect effect computation
! 2011-06: M. Jerczynski - some cleaning to meet norms
! S. Riette: 2011-10 : Modifications for DUAL-MF scheme (according to Wim de Rooy's apl_arome version)
! Ice in EDKF
! Y. Seity : 2012-03 : add LMNHLEV option to revert/or not arrays for MesoNH parameterisations
! F. Bouttier: 2012-07 add SPPT stochastic physics
! JJMorcrette, ECMWF, 20120815 additional dummy due to changes in RADACT
! P. Marguinaud : 2012-09 : Add control threshold for orography
! Y. Seity : 2013-01 Cleaning LMNHLEV and remove JPVEXT points
! Y. Seity : 2013-02 Cleaning (add compute_neb)
! L. Bengtsson: 2013-02: add LOLSMC and LOTOWNC options to compute (or not) cloud sedimentation
! using different cloud droplet number conc. depending on land/sea/town.
! 2013-11, D. Degrauwe: Introduction of radflex interface, export
! upper-air precipitation fluxes PFPR.
! 2013-11, J. Masek: Inclusion of ACRANEB2 radiation scheme.
! S. Riette: 2013-11: subgrid precipitation
! K. Yessad (July 2014): Move some variables.
! 2014-09, C. Wastl: Adaptations for orographic shadowing
! 2014-11, Y. Seity: add TKE budgets for DDH
! 2016-03, E. Bazile: Phasing MUSC for surf_ideal_flux
! 2016-04, J. Masek: LRNUEXP cloud overlap option (COMPUTE_NEB replaced
! by ACNPART), passing of sushine duration, fix of
! E. Gleeson for ACRANEB2 with SURFEX.
! 2016-09, J. Masek: Proper calculation of sunshine duration in ACRANEB2.
! 2016-10, P. Marguinaud : Port to single precision
! S. Riette 2016-11: Changes in ICE3/ICE4
! 2018-09, E. Gleeson: Corrected misplaced arguments in ACRANEB2 call.
! 2019-09-24 J.M. Piriou arguments for convective gusts.
! R. El Khatib 30-Oct-2018 substantial rewrite for optimization and coding standards respect.
! 2018-10, I. Etchevers : add Visibilities
! 2019-01, I. Etchevers, Y. Seity : add Precipitation Type
! 2019-09, J. Masek: Corrected dimensioning of dummy argument PGMU0.
! Modified call to ACRANEB2 (clearsky fluxes).
! 2019-10, I. Etchevers : Visibilities in ACVISIH, AROCLDIA=>ACCLDIA
! 2019-10, Y.Bouteloup and M. Bouzghaiam : Radiation modifications. Remove acradin.F90 direct
! call to recmwf.F90 and add interface to ecrad (in recmwf !)
! 2020-10, J. Masek: Modified call to ACCLDIA.
! 2020-12, F. Meier add call to latent heat nudging if LNUDGLH is TRUE
! 2020-12, U. Andrae : Introduce SPP for HARMONIE-AROME
! R. El Khatib 24-Aug-2021 NPROMICRO specific cache-blocking factor for microphysics
! End modifications
!-------------------------------------------------------------------------------
USE GEOMETRY_MOD , ONLY : GEOMETRY
USE SURFACE_FIELDS_MIX , ONLY : TSURF
USE YOMXFU , ONLY : TXFU
USE YOMCFU , ONLY : TCFU
USE TYPE_MODEL , ONLY : MODEL
USE PARKIND1 , ONLY : JPIM ,JPRB ,JPRD
USE YOMHOOK , ONLY : LHOOK ,DR_HOOK
! AROME SPECIFIC
! OTHERS
USE YOESW , ONLY : RSUN2
USE YOMCST , ONLY : RG ,RCPD ,RD ,RATM ,RTT ,&
& RCW ,RCPV ,RLVTT ,RCS ,RLSTT ,RGAMW ,&
& RBETW ,RALPW ,RGAMS ,RBETS ,RALPS ,RGAMD ,&
& RBETD ,RALPD ,RETV ,RV ,RKAPPA ,RHOUR
USE YOMLUN , ONLY : NULOUT
USE YOMCT0 , ONLY : LTWOTL, LSFORCS
USE YOMVERT , ONLY : VP00
USE YOMRIP0 , ONLY : NINDAT
USE YOMNUDGLH , ONLY : LNUDGLH, NSTARTNUDGLH, NSTOPNUDGLH, NINTNUDGLH, NTAUNUDGLH, &
& RAMPLIFY,RMAXNUDGLH,RMINNUDGLH,LNUDGLHCOMPT,NTIMESPLITNUDGLH
USE YOMNSV , ONLY : NSV_CO2
USE DDH_MIX , ONLY : ADD_FIELD_3D, NEW_ADD_FIELD_3D, TYP_DDH ! for new diag data flow
USE YOMSPSDT , ONLY : YSPPT_CONFIG, YSPPT
USE SPP_MOD , ONLY : YSPP_CONFIG, YSPP
USE YOMLSFORC , ONLY : LMUSCLFA, NMUSCLFA, REMIS_FORC, RALB_FORC
! daand: radflex
USE INTFLEX_MOD, ONLY : LINTFLEX, LRADFLEX,&
& TYPE_INTPROC, TYPE_INTPROCSET,&
& NEWINTFIELD, NEWINTPROC
USE YOMMP0 , ONLY : MYPROC
! --------------------------------------------------------------------------
IMPLICIT NONE
TYPE(GEOMETRY) ,INTENT(IN) :: YDGEOMETRY
TYPE(TSURF) ,INTENT(INOUT) :: YDSURF
TYPE(TXFU) ,INTENT(INOUT) :: YDXFU
TYPE(TCFU) ,INTENT(INOUT) :: YDCFU
TYPE(MODEL) ,INTENT(INOUT) :: YDMODEL
INTEGER(KIND=JPIM),INTENT(IN) :: KLON
INTEGER(KIND=JPIM),INTENT(IN) :: KLEV
INTEGER(KIND=JPIM),INTENT(IN) :: KBL
INTEGER(KIND=JPIM),INTENT(IN) :: KGPCOMP
INTEGER(KIND=JPIM),INTENT(IN) :: KIDIA
INTEGER(KIND=JPIM),INTENT(IN) :: KFDIA
INTEGER(KIND=JPIM),INTENT(IN) :: KTDIA
INTEGER(KIND=JPIM),INTENT(IN) :: KSTEP
INTEGER(KIND=JPIM),INTENT(IN) :: KMAXDRAFT
INTEGER(KIND=JPIM),INTENT(IN) :: KSGST
INTEGER(KIND=JPIM),INTENT(IN) :: KNFRRC
REAL(KIND=JPRB) ,INTENT(IN) :: PDT
LOGICAL ,INTENT(IN) :: LDXFUMSE
REAL(KIND=JPRB) ,INTENT(IN) :: PINDX(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PINDY(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PGEMU(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PGELAM(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PGM(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: POROG(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PMU0(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PMU0LU(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PMU0M(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PMU0N(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PCLON(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PSLON(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PVO3ABC(KLON,3)
REAL(KIND=JPRB) ,INTENT(IN) :: PLSM(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PAESEA(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PAELAN(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PAESOO(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PAEDES(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PAESUL(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PAEVOL(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PGP2DSDT(KLON,YSPPT%YGPSDT(1)%NG2D,YSPPT%N2D)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGP2DSPP(KLON,YSPP%N2D)
REAL(KIND=JPRB) ,INTENT(IN) :: PAPHIM(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PAPHIFM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PAPRSM(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PAPRSFM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PRDELPM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PDELPM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PTM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQVM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PCPM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PRM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PALPHM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PLNPRM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQCM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQIM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQRM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQSM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQGM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PQHM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PLIMAM(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NLIMA)
REAL(KIND=JPRB) ,INTENT(IN), TARGET :: PTKEM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PEFB1(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PEFB2(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PEFB3(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PSVM(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
REAL(KIND=JPRB) ,INTENT(IN) :: PSIGM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PUM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PVM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PWM(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PFORCEU(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PFORCEV(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PFORCET(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(IN) :: PFORCEQ(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGPAR(KLON,YDMODEL%YRML_PHY_MF%YRPARAR%NGPAR+1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PEMTD(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PEMTU(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PTRSO(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGDEOSI(KLON,0:KLEV,2)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGUEOSI(KLON,0:KLEV,2)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGMU0(KLON,0:YDMODEL%YRML_PHY_MF%YRPHY%NSORAYFR-1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGMU0_MIN(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGMU0_MAX(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGDEOTI(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGDEOTI2(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGUEOTI(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGUEOTI2(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGEOLT(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGEOXT(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGRPROX(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGMIXP(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGFLUXC(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PGRSURF(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PTURB3D(KLON,YDMODEL%YRML_PHY_MF%YRARPHY%NGRADIENTS,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PCLFS(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PQLRAD(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PEDR(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PQIRAD(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PSIGS(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDT(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDR(KLON,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDU(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDV(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDW(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDLIMA(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NLIMA)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDTKE(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PRH(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PTENDEFB1(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PTENDEFB2(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PTENDEFB3(KLON,KLEV)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTENDEXT(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
REAL(KIND=JPRB) ,INTENT(IN) :: PALBIN(KLON)
! daand: radflex; made target
REAL(KIND=JPRB) ,INTENT(INOUT), TARGET :: PFRTH(KLON,0:KLEV,KSGST+1)
REAL(KIND=JPRB) ,INTENT(OUT), TARGET :: PFRSO(KLON,0:KLEV,KSGST+1)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRTHDS(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRSOPS(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRSDNI(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRSODS(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRSOPT(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRSOC(KLON,0:1)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFRTHC(KLON,0:1)
REAL(KIND=JPRB) ,INTENT(IN) :: PEMIS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PEZDIAG(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EZDIAG)
REAL(KIND=JPRB) ,INTENT(OUT) :: PCLCH(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PCLCL(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PCLCM(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PCLCT(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFPLSL(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFPLSN(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFPLSG(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFPLSH(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PSTRTU(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PSTRTV(KLON,0:KLEV)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFCS(KLON,KSGST+1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFCLL(KLON,KSGST+1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFCLN(KLON,KSGST+1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFEVL(KLON,KSGST+1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PFEVN(KLON,KSGST+1)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PUCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PVCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PNUCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PNVCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PTCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PQCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PHUCLS(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PUGST(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PVGST(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PPBLH(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PSPSG(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PSPSGR(KLON)
REAL(KIND=JPRB) ,INTENT(INOUT) :: PSDUR(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PDIAGH(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PTPWCLS(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PFLASH(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PDPRECIPS(KLON,YDMODEL%YRML_PHY_MF%YRPHY%YRDPRECIPS%NDTPREC)
REAL(KIND=JPRB) ,INTENT(OUT) :: PDPRECIPS2(KLON,YDMODEL%YRML_PHY_MF%YRPHY%YRDPRECIPS%NDTPREC2)
REAL(KIND=JPRB) ,INTENT(OUT) :: PVISICLD(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PVISIHYDRO(KLON)
REAL(KIND=JPRB) ,INTENT(OUT) :: PMXCLWC(KLON)
REAL(KIND=JPRB) ,INTENT(IN) :: PSFORC(KLON,YDMODEL%YRML_PHY_MF%YRPHYDS%NSFORC)
! daand: radflex
TYPE(TYPE_INTPROCSET), INTENT(INOUT) :: YDPROCSET
TYPE(TYP_DDH) ,INTENT(INOUT) :: YDDDH
!*
! ------------------------------------------------------------------
! 3D arrays de reference dans mesoNH. En 1D, thetavref=thetavM, mais la question
! concernant la facon d initialiser cette variable dans le 3D reste ouverte (idem pour RHODREF)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! READ ME, PLEASE !
! CODING CONVENTIONS FOR ARPEGE vs MNH PHYSICS
! The horizontal representation in MNH physics is such that 2 dimensions are needed, while a single dimension
! is used in ARPEGE and AROME. The remapping from 1 to 2 dimensions will be made inplicitly through the
! subroutines interfaces (this is a fortran property). Therefore there is non need to add an explicit dimension
! sized to 1.
! Local 3D arrays with extra levels for Meso-NH turbulence scheme :
! - first dimension is KFDIA not KDLON in order to limit array copies
! - suffixed with two underscore to be easily identified
! These arrays are passed in argument as ZXXX__(:,1:KLEV) except for aro_turb_mnh where they are passed as ZXXX__.
! Local 3D arrays with regular number of levels for Meso-NH interfaces :
! - first dimension is KFDIA not KDLON in order to limit array copies due to non-contiguous data.
! - suffixed with one underscore to be easily identified.
! Local 4D arrays with regular number of levels for Meso-NH interfaces :
! - first dimension is KFDIA not KDLON in order to limit array copies
! - suffixed with one underscore to be easily identified
! Local 2D arrays for Meso-NH interfaces :
! - first dimension is KFDIA not KDLON as a convention as well (no arrays copies to fear)
! - suffixed with one underscore to be easily identified
! Other arrays, which can be dummy arguments, or local but used as argument to IFS/ARPEGE physics
! should remained dimensionned KDLON and should not be suffixed with undersores.
! DO NOT USE ARRAY SYNTAX FOR COMPUTATIONAL LOOPS !!
! - They make the code less performant because memory cache is poorly used
! - They can make the code even less readable if the indexes are removed
! AVOID ARRAYS COPIES, OR MAKE THEM FAST !
! - if you do need to initialize or copy an array, do it as follows with explicit array syntax in first dimension
! because the compiler will be able to use an optimized function to initialize/copy a segment of memory,
! and may be able to address simultaneously several cach lines :
! 1D array :
! Z(KIDIA:KFDIA)=value
! 2D arrays :
! DO JLEV=1,KLEV
! ZX(KIDIA:KFDIA,JLEV)=xval
! ZY(KIDIA:KFDIA,JLEV)=yval
! ENDDO
! - if you need the bakup of an array, use a swapp mechanism, as what is done here for instance for
! ZSVM_ and ZSVMIN_ :
! for the developer ZSVMIN_ is always the bakup and ZSVM_ is always the current array.
! ZSVMSWAP_ or ZSVMSAVE_ are used for the swapp mechanism, they should never be used by the developer.
! - do not initialize if not necessary. To avoid useless initialization use the mechanizm 'INIT0' coded below :
! IF (INIT0 == 0) THEN
! ZVALUE=HUGE(1._JPRB)
! ELSE
! ZVALUE=default_value
! ENDIF
! IF (INIT0 >= 0) THEN
! Z(:)=ZVALUE
! ENDIF
! INIT0= 0 : initialize to HUGE (testing/debugging)
! INIT0= 1 : initialize to realistic value (discouraged !)
! INIT0=-1 : no initialization (optimized code) - this is the default.
! - pointers can advantageously be used to avoid copies or to avoid the initialization to zero of a sum.
! Example :
! DO JI=IFIRST,ILAST
! IF (JI == IFIRST) THEN
! ! Fill the sum at the first iteration
! ZARG => ZSUM(:)
! ELSE
! ! increment
! ZARG => ZINC(:)
! ENDIF
! CALL COMPUTE(ZARG)
! IF (JI > IFIRST) THEN
! ! Add increment
! ZSUM(KIDIA:KFDIA)=ZSUM(KIDIA:KFDIA)+ZINC(KIDIA:KFDIA)
! ENDIF
! ENDDO
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
REAL(KIND=JPRB) :: ZRHODJM__(KFDIA,0:KLEV+1), ZRHODREFM__(KFDIA,0:KLEV+1), ZPABSM__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB) :: ZUM__(KFDIA,0:KLEV+1), ZVM__(KFDIA,0:KLEV+1), ZTHM__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB) :: ZUS__(KFDIA,0:KLEV+1), ZVS__(KFDIA,0:KLEV+1), ZWS__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB) :: ZTKES_OUT__(KFDIA,0:KLEV+1), ZMF_UP__(KFDIA,0:KLEV+1), ZTHVREFM__(KFDIA,0:KLEV+1) ! thetav de l etat
REAL(KIND=JPRB) :: ZTENDU_TURB__(KFDIA,0:KLEV+1), ZTENDV_TURB__(KFDIA,0:KLEV+1), ZTENDTHL_TURB__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB) :: ZTENDRT_TURB__(KFDIA,0:KLEV+1), ZTKEM__(KFDIA,0:KLEV+1), ZSRCS__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB) :: ZHLC_HRC__(KFDIA,0:KLEV+1), ZHLC_HCF__(KFDIA,0:KLEV+1), &
& ZHLI_HRI__(KFDIA,0:KLEV+1), ZHLI_HCF__(KFDIA,0:KLEV+1)
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
922
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
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
REAL(KIND=JPRB) :: ZSIGS__(KFDIA,0:KLEV+1), ZEDR__(KFDIA,0:KLEV+1)
! THE DDH budgets
REAL(KIND=JPRB) :: ZDP__(KFDIA,0:KLEV+1), ZTP__(KFDIA,0:KLEV+1), ZTPMF__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB) :: ZTDIFF__(KFDIA,0:KLEV+1), ZTDISS__(KFDIA,0:KLEV+1)
! length scales for momentum and heat for mnh level definitions in case LHARATU=TRUE
REAL(KIND=JPRB) :: ZLENGTHM__(KFDIA,0:KLEV+1), ZLENGTHH__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB), POINTER :: ZTHS__(:,:)
! horizontal gradients and diagnostics
REAL(KIND=JPRB) :: ZTURB3D__(KFDIA,0:KLEV+1,YDMODEL%YRML_PHY_MF%YRARPHY%NGRADIENTS)
! WARNING ! Don't use ZTHSWAP__ or ZTHSAVE__ below because they may be swapped !
! Use only the pointer ZTHS__, and possibly ZTHSIN_ if you need the backup of input data.
REAL(KIND=JPRB), TARGET :: ZTHSWAP__(KFDIA,0:KLEV+1), ZTHSAVE__(KFDIA,0:KLEV+1)
REAL(KIND=JPRB), TARGET :: ZFLXZTHVMF_SUM__(KFDIA,0:KLEV+1), ZWM__(KFDIA,0:KLEV+1)
! Updraft characteristics for Meso-NH world (input of ARO_SHALLOW_MF)
REAL(KIND=JPRB) :: ZTHETAL_UP_(KFDIA,KLEV), ZTHETAV_UP_(KFDIA,KLEV), ZZFRAC_UP_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZRT_UP_(KFDIA,KLEV), ZRC_UP_(KFDIA,KLEV), ZRI_UP_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZZU_UP_(KFDIA,KLEV), ZZV_UP_(KFDIA,KLEV), ZZW_UP_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZZRV_UP_(KFDIA,KLEV), ZTKES_(KFDIA,KLEV), ZZZ_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZDZZ_(KFDIA,KLEV), ZZZ_F_(KFDIA,KLEV), ZDZZ_F_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZCIT_(KFDIA,KLEV), ZMFM_(KFDIA,KLEV), ZEXNREFM_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZSIGM_(KFDIA,KLEV), ZNEBMNH_(KFDIA,KLEV), ZEVAP_(KFDIA,KLEV)
! additions for MF scheme (Pergaud et al)
REAL(KIND=JPRB) :: ZSIGMF_(KFDIA,KLEV), ZRC_MF_(KFDIA,KLEV), ZRI_MF_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZCF_MF_(KFDIA,KLEV), ZAERD_(KFDIA,KLEV), ZCVTENDT_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZCVTENDRV_(KFDIA,KLEV), ZCVTENDRC_(KFDIA,KLEV), ZCVTENDRI_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZMFS_(KFDIA,KLEV), ZTHLS_(KFDIA,KLEV), ZRTS_(KFDIA,KLEV)
REAL(KIND=JPRB) :: ZMFUS_(KFDIA,KLEV), ZMFVS_(KFDIA,KLEV), ZDEPTH_HEIGHT_(KFDIA,KLEV)
REAL(KIND=JPRB), TARGET :: ZFLXZTHVMF_(KFDIA,KLEV)
REAL(KIND=JPRB), POINTER :: ZARG_FLXZTHVMF_(:,:)
! WARNING ! Don't use ZRSWAP_ or ZRSAVE_ below because they may be swapped !
! Use only the pointer ZRS_, and possibly ZRSIN_ if you need the backup of input data.
REAL(KIND=JPRB), POINTER :: ZRSIN_(:,:,:), ZRS_(:,:,:)
REAL(KIND=JPRB), TARGET :: ZRSWAP_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
REAL(KIND=JPRB), TARGET :: ZRSAVE_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
REAL(KIND=JPRB), POINTER :: ZPTRWNU_(:,:), ZTHSIN_(:,:)
REAL(KIND=JPRB), TARGET :: ZWNU_(KFDIA,KLEV)
! WARNING ! Don't use ZSVSWAP_ or ZSVSAVE_ below because they may be swapped !
! Use only the pointer ZSVS_, and possibly ZSVSIN_ if you need the backup of input data.
REAL(KIND=JPRB), POINTER :: ZSVSIN_(:,:,:), ZSVS_(:,:,:)
REAL(KIND=JPRB), TARGET :: ZSVSWAP_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
REAL(KIND=JPRB), TARGET :: ZSVSAVE_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
REAL(KIND=JPRB) :: ZSVXXX_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
! WARNING ! Don't use ZSVMSWAP_ or ZSVMSAVE_ below because they may be swapped !
! Use only the pointer ZSVM_, and possibly ZSVMIN_ if you need the backup of input data.
REAL(KIND=JPRB), POINTER :: ZSVMIN_(:,:,:), ZSVM_(:,:,:)
REAL(KIND=JPRB), TARGET :: ZSVMSWAP_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
REAL(KIND=JPRB), TARGET :: ZSVMSAVE_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
REAL(KIND=JPRB) :: ZSVMB_(KFDIA,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT)
! WARNING ! Don't use ZLIMASWAP_ or ZLIMASAVE_ below because they may be swapped !
! Use only the pointer ZLIMAS_, and possibly ZLIMASIN_ if you need the backup of input data.
REAL(KIND=JPRB), POINTER :: ZLIMAS_(:,:,:), ZLIMASIN_(:,:,:)
REAL(KIND=JPRB), TARGET :: ZLIMASWAP_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NLIMA)
REAL(KIND=JPRB), TARGET :: ZLIMASAVE_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NLIMA)
REAL(KIND=JPRB) :: ZLIMAM_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NLIMA)
!INTEGER(KIND=JPIM) :: KSV_TURB !CPtoclean?
!CPtoclean REAL(KIND=JPRB) :: ZTURBM(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT+YDMODEL%YRML_GCONF%YGFL%NLIMA)
!CPtoclean REAL(KIND=JPRB) :: ZTURBS(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT+YDMODEL%YRML_GCONF%YGFL%NLIMA)
!not (yet ?) used. REK
!REAL(KIND=JPRB) :: ZSFTURB(KLON,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT+YDMODEL%YRML_GCONF%YGFL%NLIMA) ! surf. flux of SV (=0)
!REAL(KIND=JPRB) :: ZTENDSV_TURB2(KLON,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT+YDMODEL%YRML_GCONF%YGFL%NLIMA) ! SV (=0)
REAL(KIND=JPRB) :: ZSFSVLIMA_(KFDIA,YDMODEL%YRML_GCONF%YGFL%NLIMA) ! surf. flux of LIMA vars
REAL(KIND=JPRB) :: ZTENDSV_TURBLIMA_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NLIMA) ! LIMA
! For radiation scheme
REAL(KIND=JPRB) :: ZRM_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
REAL(KIND=JPRB) :: ZPFPR_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
REAL(KIND=JPRB) :: ZPEZDIAG_(KFDIA,KLEV,YDMODEL%YRML_GCONF%YGFL%NGFL_EZDIAG)
REAL(KIND=JPRB) :: ZSFSV_(KFDIA,YDMODEL%YRML_GCONF%YGFL%NGFL_EXT) ! surf. flux of scalars
! Single scattering albedo of dust (points,lev,wvl) :
REAL(KIND=JPRB) :: ZPIZA_DST_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NSWB_MNH)
! Assymetry factor for dust (points,lev,wvl) :
REAL(KIND=JPRB) :: ZCGA_DST_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NSWB_MNH)
! tau/tau_{550} dust (points,lev,wvl) :
REAL(KIND=JPRB) :: ZTAUREL_DST_(KFDIA,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NSWB_MNH)
! surface flux of theta and surface flux of vapor ; surface flux of CO2
REAL(KIND=JPRB) :: ZSFTH_(KFDIA), ZSFRV_(KFDIA), ZSFCO2_(KFDIA)
REAL(KIND=JPRB) :: ZACPRG_(KFDIA), ZINPRG_NOTINCR_(KFDIA), ZINPRG_(KFDIA)
REAL(KIND=JPRB) :: ZACPRR_(KFDIA), ZINPRR_NOTINCR_(KFDIA), ZINPRR_(KFDIA)
REAL(KIND=JPRB) :: ZACPRS_(KFDIA), ZINPRS_NOTINCR_(KFDIA), ZINPRS_(KFDIA)
REAL(KIND=JPRB) :: ZCFBTH_(KFDIA), ZINPRH_NOTINCR_(KFDIA), ZINPRH_(KFDIA)
REAL(KIND=JPRB) :: ZZS_(KFDIA), ZSSO_STDEV_(KFDIA), ZALB_UV_(KIDIA)
REAL(KIND=JPRB) :: ZLAT_(KIDIA), ZLON_(KIDIA), ZZENITH_(KIDIA)
REAL(KIND=JPRB) :: ZGZ0_(KFDIA), ZGZ0H_(KFDIA), ZTOWNS_(KFDIA)
REAL(KIND=JPRB) :: ZCFAQ_(KFDIA), ZCFBQ_(KFDIA), ZCFATH_(KFDIA)
REAL(KIND=JPRB) :: ZCFAU_(KFDIA), ZCFBU_(KFDIA), ZCFBV_(KFDIA)
REAL(KIND=JPRB) :: ZBUDTH_ (KFDIA),ZBUDSO_(KFDIA), ZFCLL_(KFDIA)
REAL(KIND=JPRB) :: ZCD_(KFDIA), ZSEA_(KFDIA), ZTOWN_(KFDIA)
REAL(KIND=JPRB) :: ZZTOP_(KFDIA), ZCVTENDPR_(KFDIA), ZCVTENDPRS_(KFDIA)
! surface flux of x and y component of wind. are they really necessary ? REK
REAL(KIND=JPRB) :: ZSFU_(KFDIA), ZSFV_(KFDIA)
! Arpege-style dimensionning :
! --------------------------
!Variables used in case LHARATU=TRUE
! length scales for momentum and heat and TKE
REAL(KIND=JPRB) :: ZLENGTH_M(KLON,KLEV),ZLENGTH_H(KLON,KLEV)
REAL(KIND=JPRB) :: ZTKEEDMF(KLON,KLEV)
REAL(KIND=JPRB) :: ZTKEEDMFS(KLON,KLEV)
REAL(KIND=JPRB) :: ZEMIS (KLON)
REAL(KIND=JPRB) :: ZQICE(KLON,KLEV), ZQLIQ(KLON,KLEV)
REAL(KIND=JPRB) :: ZAER(KLON,KLEV,6)
REAL(KIND=JPRB) :: ZAERINDS(KLON,KLEV)
REAL(KIND=JPRB) :: ZQSAT(KLON,KLEV)
REAL(KIND=JPRB) :: ZFRSOFS(KLON)
REAL(KIND=JPRB) :: ZLH(KLON,KLEV), ZLSCPE(KLON,KLEV), ZGEOSLC(KLON,KLEV)
REAL(KIND=JPRB) :: ZQDM(KLON,KLEV), ZQV(KLON,KLEV)
REAL(KIND=JPRB) :: ZQCO2(KLON,KLEV)
REAL(KIND=JPRB), TARGET :: ZTKEM4SLDDH(KLON,KLEV)
REAL(KIND=JPRB), POINTER :: ZTKEM(:,:)
REAL(KIND=JPRB) :: ZQW(KLON,KLEV), ZTW(KLON,KLEV)
REAL(KIND=JPRB) :: ZTENT(KLON,KLEV)
REAL(KIND=JPRB) :: ZTENDT(KLON,KLEV) ! array to save heating profile for LHN
REAL(KIND=JPRB) :: ZMAXTEND,ZMINTEND
REAL(KIND=JPRB) :: ZDZZ(KLON,KLEV)
REAL(KIND=JPRB) :: ZTPW(KLON,KLEV)
! POUR GROUND
REAL(KIND=JPRB) :: ZZS_FSWDIR(KLON,YDMODEL%YRML_PHY_RAD%YRERAD%NSW), ZZS_FSWDIF(KLON,YDMODEL%YRML_PHY_RAD%YRERAD%NSW)
REAL(KIND=JPRB) :: ZTRSODIR (KLON,YDMODEL%YRML_PHY_RAD%YRERAD%NSW), ZTRSODIF (KLON,YDMODEL%YRML_PHY_RAD%YRERAD%NSW)
REAL(KIND=JPRB) :: ZTRSOD(KLON)
REAL(KIND=JPRB) :: ZSUDU(KLON), ZSDUR(KLON), ZDSRP(KLON)
REAL(KIND=JPRB) :: ZCEMTR(KLON,2), ZCTRSO(KLON,2)
REAL(KIND=JPRB) :: ZALBD(KLON,YDMODEL%YRML_PHY_RAD%YRERAD%NSW), ZALBP(KLON,YDMODEL%YRML_PHY_RAD%YRERAD%NSW)
REAL(KIND=JPRB) :: ZALBD1(KLON), ZALBP1(KLON)
REAL(KIND=JPRB) :: ZAPHIM(KLON,0:KLEV), ZAPHIFM(KLON,KLEV)
REAL(KIND=JPRB) :: ZTM(KLON,KLEV), ZQVM(KLON,KLEV), ZQIM(KLON,KLEV)
REAL(KIND=JPRB) :: ZQCM(KLON,KLEV),ZQHM(KLON,KLEV), ZQHGM(KLON,KLEV)
REAL(KIND=JPRB) :: ZQRM(KLON,KLEV), ZQSM(KLON,KLEV), ZQGM(KLON,KLEV)
REAL(KIND=JPRB) :: ZUPGENL(KLON,KLEV)
REAL(KIND=JPRB) :: ZUPGENN(KLON,KLEV)
REAL(KIND=JPRB) :: ZCLFR(KLON)
REAL(KIND=JPRB) :: ZCPM(KLON,KLEV), ZRHM(KLON,KLEV)
! Variables concerning updraft rain/snow for EDMF
REAL(KIND=JPRB) :: ZTENDTUP(KLON,KLEV), ZTENDQVUP(KLON,KLEV)
! specific to new data flow for diagnostics
REAL(KIND=JPRB) :: ZTENDTBAK(KLON,KLEV), ZTENDRBAK(KLON,KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
REAL(KIND=JPRB) :: ZTMPAF(KLON,KLEV)
! daand: radflex
REAL(KIND=JPRB) :: ZFPR(KLON,0:KLEV,YDMODEL%YRML_PHY_MF%YRPARAR%NRR)
! Target should not be necessary. REK
REAL(KIND=JPRB), TARGET :: ZCON1(KLON,KLEV)
REAL(KIND=JPRB), TARGET :: ZCON2(KLON,KLEV)
REAL(KIND=JPRB), TARGET :: ZCON3(KLON,KLEV)
! Ajout pour MF Dual Scheme (KNMI et al)
! Updraft characteristics in Arpege/IFS world
REAL(KIND=JPRB) :: ZMF_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZTHETAL_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZQT_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZTHTV_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZQC_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZQI_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZU_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZV_UP(KLON,0:KLEV,KMAXDRAFT)
REAL(KIND=JPRB) :: ZTSURF(KLON), ZTN(KLON), ZQS(KLON)
REAL(KIND=JPRB) :: ZFRSOLU(KLON), ZFRSODS(KLON)
REAL(KIND=JPRB) :: ZFSDNN(KLON), ZFSDNV(KLON)
REAL(KIND=JPRB) :: ZSURFPREP(KLON), ZSURFSNOW(KLON)
REAL(KIND=JPRB) :: ZQO3(KLON,0:KLEV)
REAL(KIND=JPRB) :: ZZS_FTH_(KLON), ZZS_FRV_(KLON), ZZS_FU_(KLON), ZZS_FV_(KLON)
! Surface forcing arrays for MUSC
REAL(KIND=JPRB) :: ZRHODREFM(KLON), ZTHETAS(KLON)
! ACRANEB2 local variables
REAL(KIND=JPRB) :: ZNEB0 (KLON,KLEV) ! protected cloud fractions
REAL(KIND=JPRB) :: ZCLCT_RAD(KLON) ! total cloud cover for radiation
REAL(KIND=JPRB) :: ZDECRD (KLON) ! decorrelation depth
! Stochastic physics pattern & dummy tendencies for calling sppten
! Bof. REK
REAL(KIND=JPRB) :: ZDUMMY(KLON,KLEV)
REAL(KIND=JPRB) :: ZDUMMY1(KLON)
REAL(KIND=JPRB) :: ZROZ(KLON,KLEV)
! Can we remove ? REK
REAL(KIND=JPRB) :: ZEPSM(0,0,0) ! Dissipation of TKE (eps) at time t-dt
REAL(KIND=JPRB) :: ZEPSS(0,0,0) ! Dissipation of TKE at time t+dt
! Integers
INTEGER(KIND=JPIM) :: JLEV, JLON, JRR, JGFL, JGR, ISPLITR
INTEGER(KIND=JPIM) :: IJN ! max. number of day/night slices within NRPOMA
INTEGER(KIND=JPIM) :: IKL !ordering of vert levels 1:MNH -1:AROME
INTEGER(KIND=JPIM) :: IOFF_MFSHAL, IEZDIAG_CHEM
INTEGER(KIND=JPIM) :: IKA,IKB,IKU,IKT,IKTE,IKTB ! vertical points as in mpa
INTEGER(KIND=JPIM) :: JSG, JK, JR, JSW
INTEGER(KIND=JPIM) :: IDRAFT,JDRAFT,INDRAFT
INTEGER(KIND=JPIM) :: ISURFEX
INTEGER(KIND=JPIM) :: IDAY,IYEAR,IMONTH
INTEGER(KIND=JPIM) :: INIT0 ! Kind of safety/debugging initialization :
! 0 = initialize to HUGE (debugging)
! 1 = initialize to realistic value (discouraged)
! -1 = no initialization (optimized code) - this is the default.
INTEGER(KIND=JPIM) :: ICLPH(KLON) !PBL top level
INTEGER(KIND=JPIM) :: JLHSTEP,ISTEP
! Real
REAL(KIND=JPRB) :: ZRHO
REAL(KIND=JPRB) :: ZAEO, ZAEN, ZSALBCOR
REAL(KIND=JPRB) :: ZDT, ZDT2, ZINVDT, ZINVG, ZRSCP, ZINVATM, Z_WMAX, Z_WMIN
! pas de temps pour la surface externalise
REAL(KIND=JPRB) :: ZDTMSE,ZRHGMT,ZSTATI
REAL(KIND=JPRB) :: ZDELTA
REAL(KIND=JPRB) :: ZEPSNEB
! default values for initialization :
REAL(KIND=JPRB) :: ZVALUE, ZVALUE_ONE, ZVALUE_T, ZVALUE_P, ZVALUE_L, ZVALUE_EPSILON
REAL(KIND=JPRB) :: ZVETAH(0:KLEV)
! Boolean
LOGICAL :: LLMSE, LLMSE_PARAM, LLMSE_DIAG
LOGICAL :: LLAROME
LOGICAL :: LLRAD
LOGICAL :: LLSWAP_THS, LLSWAP_RS, LLSWAP_SVS, LLSWAP_SVM, LLSWAP_LIMAS ! logical to swap or not pointers in and out
LOGICAL :: LLHN(KLON,KLEV)
LOGICAL :: LNUDGLHNREAD
! Characters
CHARACTER(LEN=11) :: CLNAME
CHARACTER(LEN=2),DIMENSION(7):: CLVARNAME=(/"QV","QL","QR","QI","QS","QG","QH"/)
! daand: radflex
REAL(KIND=JPRB), POINTER :: ZFRSO(:,:), ZFRTH(:,:)
TYPE(TYPE_INTPROC), POINTER :: YLRADPROC
REAL(KIND=JPRB) :: ZCAPE(KLON), ZDCAPE(KLON)
!
! Phaser team note from CY43T1:
! there was a USE MODD_CTURB for accessing XTKEMIN here, but that created a forbidden
! dependence of APL_AROME (in "ifsarp") to the Méso-NH/Arome interfaces (in "mpa").
! There should be no USE MODD_* in APL_*.
! We decided to change the variable here to a local one, with the classical initial value for TKE.
!
REAL(KIND=JPRB), PARAMETER :: PPTKEMIN = 1.E-6
! Perturbed radiation-cloud interaction coef
REAL(KIND=JPRB), DIMENSION (KLON) :: ZRADGR,ZRADSN
REAL(KIND=JPRB) :: ZMU,ZVAL
INTEGER(KIND=JPIM) :: JKO,JKE
REAL(KIND=JPRB) :: ZHOOK_HANDLE
! --------------------------------------------------------------------------
#include "abor1.intfb.h"
#include "recmwf.intfb.h"
#include "acraneb2.intfb.h"
#include "actqsat.intfb.h"
#include "acnpart.intfb.h"
#include "bri2acconv.intfb.h"
#include "gpgeo.intfb.h"
#include "gprcp.intfb.h"
#include "radheat.intfb.h"
#include "suozon.intfb.h"
#include "radaer.intfb.h"
#include "radozc.intfb.h"
#include "accldia.intfb.h"
#include "vdfhghthl.intfb.h"
#include "sppten.intfb.h"
#include "surf_ideal_flux.intfb.h"
#include "ecr1d.intfb.h"
#include "apl_arome2intflex.intfb.h"
#include "aro_rain_ice.h"
#include "nudglhprecip.intfb.h"
#include "nudglh.intfb.h"
#include "nudglhclimprof.intfb.h"
#include "nudglhprep.intfb.h"
#include "aro_turb_mnh.h"
#include "aro_adjust.h"
#include "aro_mnhc.h"
#include "aro_mnhdust.h"
#include "aro_startbu.h"
#include "aro_convbu.h"
#include "aro_ground_param.h"
#include "aro_ground_diag.h"
#include "aro_shallow_mf.h"
#include "aro_rainaero.h"
#include "aro_lima.h"
#include "diagflash.intfb.h"
#include "dprecips.intfb.h"
#include "ppwetpoint.intfb.h"
#include "acvisih.intfb.h"
#include "aro_ground_diag_2isba.h"
! ------------------------------------------------------------------
#include "fcttrm.func.h"
! --------------------------------------------------------------------------
IF (LHOOK) CALL DR_HOOK('APL_AROME',0,ZHOOK_HANDLE)
ASSOCIATE(YDDIM=>YDGEOMETRY%YRDIM,YDDIMV=>YDGEOMETRY%YRDIMV,YDGEM=>YDGEOMETRY%YRGEM, YDMP=>YDGEOMETRY%YRMP, &
& YDVETA=>YDGEOMETRY%YRVETA, YDVFE=>YDGEOMETRY%YRVFE, YDSTA=>YDGEOMETRY%YRSTA, YDLAP=>YDGEOMETRY%YRLAP, &
& YDCSGLEG=>YDGEOMETRY%YRCSGLEG, YDVSPLIP=>YDGEOMETRY%YRVSPLIP, YDVSLETA=>YDGEOMETRY%YRVSLETA, &
& YDHSLMER=>YDGEOMETRY%YRHSLMER, YDCSGEOM=>YDGEOMETRY%YRCSGEOM, YDCSGEOM_NB=>YDGEOMETRY%YRCSGEOM_NB, &
& YDGSGEOM=>YDGEOMETRY%YRGSGEOM, YDGSGEOM_NB=>YDGEOMETRY%YRGSGEOM_NB, YDSPGEOM=>YDGEOMETRY%YSPGEOM, &
& YDPHY=>YDMODEL%YRML_PHY_MF%YRPHY, &
& YLDDH=>YDMODEL%YRML_DIAG%YRLDDH,YDTOPH=>YDMODEL%YRML_PHY_MF%YRTOPH, &
& YDERDI=>YDMODEL%YRML_PHY_RAD%YRERDI,YDRIP=>YDMODEL%YRML_GCONF%YRRIP,YDARPHY=>YDMODEL%YRML_PHY_MF%YRARPHY, &
& YDERAD=>YDMODEL%YRML_PHY_RAD%YRERAD,YDPHY3=>YDMODEL%YRML_PHY_MF%YRPHY3,YDPHYDS=>YDMODEL%YRML_PHY_MF%YRPHYDS, &
& YDPHY0=>YDMODEL%YRML_PHY_MF%YRPHY0, YDVISI=>YDMODEL%YRML_PHY_MF%YRPHY%YRDVISI,&
& YDPHY2=>YDMODEL%YRML_PHY_MF%YRPHY2,&
& YGFL=>YDMODEL%YRML_GCONF%YGFL,YDMSE=>YDMODEL%YRML_PHY_MF%YRMSE, &
& YDPARAR=>YDMODEL%YRML_PHY_MF%YRPARAR,YDPRECIPS=>YDMODEL%YRML_PHY_MF%YRPHY%YRDPRECIPS, &