apl_arome.F90

 & YDSTOPH=>YDMODEL%YRML_PHY_STOCH%YRSTOPH)

ASSOCIATE(MINPRR=>YDPARAR%MINPRR, MINPRS=>YDPARAR%MINPRS, MVQS=>YDPARAR%MVQS, &
 & MINPRG=>YDPARAR%MINPRG, LOTOWNC=>YDPARAR%LOTOWNC, LFPREC3D=>YDPARAR%LFPREC3D, &
 & NGPAR=>YDPARAR%NGPAR, CSUBG_PR_PDF=>YDPARAR%CSUBG_PR_PDF, &
 & NRRI=>YDPARAR%NRRI, NRRL=>YDPARAR%NRRL, CSUBG_AUCV_RC=>YDPARAR%CSUBG_AUCV_RC, &
 & CSUBG_AUCV_RI=>YDPARAR%CSUBG_AUCV_RI, CCONDENS=>YDPARAR%CCONDENS, &
 & CSUBG_MF_PDF=>YDPARAR%CSUBG_MF_PDF, &
 & LTOTPREC=>YDPARAR%LTOTPREC, CSUBG_RC_RR_ACCR=>YDPARAR%CSUBG_RC_RR_ACCR, CSUBG_RR_EVAP=>YDPARAR%CSUBG_RR_EVAP, &
 & NPRINTFR=>YDPARAR%NPRINTFR, CMF_CLOUD=>YDPARAR%CMF_CLOUD, &
 & MALBDIR=>YDPARAR%MALBDIR, NSWB_MNH=>YDPARAR%NSWB_MNH, &
 & XSW_BANDS=>YDPARAR%XSW_BANDS, MACPRG=>YDPARAR%MACPRG, MSWDIR=>YDPARAR%MSWDIR, &
 & LMIXUV=>YDPARAR%LMIXUV, MSWDIF=>YDPARAR%MSWDIF, LOLSMC=>YDPARAR%LOLSMC, &
 & NDIAGWMAX=>YDPARAR%NDIAGWMAX, MACPRS=>YDPARAR%MACPRS, MACPRR=>YDPARAR%MACPRR, &
 & LSQUALL=>YDPARAR%LSQUALL, VSIGQSAT=>YDPARAR%VSIGQSAT, &
 & MALBSCA=>YDPARAR%MALBSCA,&
 & RADSN=>YDPARAR%RADSN, LOSEDIC=>YDPARAR%LOSEDIC, LDIAGWMAX=>YDPARAR%LDIAGWMAX, &
 & CSEDIM=>YDPARAR%CSEDIM, CLAMBDA3=>YDPARAR%CLAMBDA3, &
 & NPTP=>YDPARAR%NPTP, NSPLITR=>YDPARAR%NSPLITR, NSPLITG=>YDPARAR%NSPLITG, NSV=>YDPARAR%NSV, &
 & CFRAC_ICE_SHALLOW_MF=>YDPARAR%CFRAC_ICE_SHALLOW_MF, CFRAC_ICE_ADJUST=>YDPARAR%CFRAC_ICE_ADJUST, &
 & MVTS=>YDPARAR%MVTS, NREFROI2=>YDPARAR%NREFROI2, NREFROI1=>YDPARAR%NREFROI1, &
 & MVEMIS=>YDPARAR%MVEMIS, LOWARM=>YDPARAR%LOWARM, LOCND2=>YDPARAR%LOCND2, &
 & LGRSN=>YDPARAR%LGRSN, LOSIGMAS=>YDPARAR%LOSIGMAS, NRR=>YDPARAR%NRR, &
 & LOSUBG_COND=>YDPARAR%LOSUBG_COND, RADGR=>YDPARAR%RADGR, &
 & CMF_UPDRAFT=>YDPARAR%CMF_UPDRAFT, LHARATU=>YDPARAR%LHARATU, &
 & XMINLM=>YDPHY0%XMINLM, XMAXLM=>YDPHY0%XMAXLM, AERCS1=>YDPHY0%AERCS1, &
 & AERCS3=>YDPHY0%AERCS3, AERCS5=>YDPHY0%AERCS5, &
 & RDECRD1=>YDPHY0%RDECRD1, RDECRD2=>YDPHY0%RDECRD2, &
 & RDECRD3=>YDPHY0%RDECRD3, RDECRD4=>YDPHY0%RDECRD4, &
 & LMPA=>YDARPHY%LMPA, LUSECHEM=>YDARPHY%LUSECHEM, LKFBCONV=>YDARPHY%LKFBCONV, &
 & LMFSHAL=>YDARPHY%LMFSHAL, LMICRO=>YDARPHY%LMICRO, &
 & CCOUPLING=>YDARPHY%CCOUPLING, LTURB=>YDARPHY%LTURB, LGRADHPHY=>YDARPHY%LGRADHPHY, &
 & NSURFEX_ITER=>YDARPHY%NSURFEX_ITER, LRDUST=>YDARPHY%LRDUST, &
 & NGRADIENTS=>YDARPHY%NGRADIENTS, &
 & LRDEPOS=>YDARPHY%LRDEPOS, LSURFEX_CRITICAL=>YDARPHY%LSURFEX_CRITICAL, &
 & LRCO2=>YDARPHY%LRCO2, LMSE=>YDARPHY%LMSE, &
 & LSURFEX_KFROM=>YDARPHY%LSURFEX_KFROM, &
 & NSURFEXCTL=>YDMSE%NSURFEXCTL, XZSEPS=>YDMSE%XZSEPS, &
 & NDLUNG=>YDDIM%NDLUNG, NDGUNG=>YDDIM%NDGUNG, NPROMA=>YDDIM%NPROMA, &
 & NDLUXG=>YDDIM%NDLUXG, NDGUXG=>YDDIM%NDGUXG, &
 & NSFORC=>YDPHYDS%NSFORC, &
 & NGFL_EXT=>YGFL%NGFL_EXT, YLRAD=>YGFL%YLRAD, YIRAD=>YGFL%YIRAD, &
 & NGFL_EZDIAG=>YGFL%NGFL_EZDIAG, &
 & NLIMA=>YGFL%NLIMA, CMICRO=>YDPARAR%CMICRO,NPROMICRO=>YDPARAR%NPROMICRO, &
 & YSD_VAD=>YDSURF%YSD_VAD, &
 & QCO2=>YDPHY3%QCO2, &
 & NTEND_DIAG_POS=>YDPHY%NTEND_DIAG_POS, NTEND_DIAG_FREQ_RESET=>YDPHY%NTEND_DIAG_FREQ_RESET, &
 & NRAY=>YDPHY%NRAY, LRAYFM=>YDPHY%LRAYFM, &
 & LO3ABC=>YDPHY%LO3ABC, LRAY=>YDPHY%LRAY, &
 & LAEROVOL=>YDPHY%LAEROVOL, LRSTAER=>YDPHY%LRSTAER, LRNUEXP=>YDPHY%LRNUEXP, &
 & AMAGSTOPH_CASBS=> YDSTOPH%AMAGSTOPH_CASBS, LFORCENL=>YDSTOPH%LFORCENL, &
 & NFORCESTART=>YDSTOPH%NFORCESTART, NFORCEEND=>YDSTOPH%NFORCEEND, &
 & NTRADI=>YDTOPH%NTRADI, NTQSAT=>YDTOPH%NTQSAT, NTNEBU=>YDTOPH%NTNEBU,&
 & NAER=>YDERAD%NAER, &
 & LHLRADUPD=>YDPHY%LHLRADUPD, &
 & TSPHY=>YDPHY2%TSPHY,&
 & NMODE=>YDERAD%NMODE, &
 & NOZOCL=>YDERAD%NOZOCL, &
 & NRADFR=>YDERAD%NRADFR, &
 & NSW=>YDERAD%NSW, &
 & RCARDI=>YDERDI%RCARDI, &
 & LFLEXDIA=>YLDDH%LFLEXDIA, LDDH_OMP=>YLDDH%LDDH_OMP, LRSLDDH=>YLDDH%LRSLDDH, &
 & RDECLI=>YDRIP%RDECLI, &
 & RCODEC=>YDRIP%RCODEC, &
 & RHGMT=>YDRIP%RHGMT, &
 & RSIDEC=>YDRIP%RSIDEC, &
 & RSOVR=>YDRIP%RSOVR, &
 & RSTATI=>YDRIP%RSTATI, &
 & TSTEP=>YDRIP%TSTEP, &
 & STPREH=>YDSTA%STPREH, &
 & LXXDIAGH=>YDXFU%LXXDIAGH,&
 & LFLASH =>YDCFU%LFLASH,&
 & LDPRECIPS=>YDPHY%LDPRECIPS,LDPRECIPS2=>YDPHY%LDPRECIPS2,&
 & NDTPREC=>YDPRECIPS%NDTPREC,NDTPREC2=>YDPRECIPS%NDTPREC2,&
 & NDTPRECCUR=>YDPRECIPS%NDTPRECCUR,NDTPRECCUR2=>YDPRECIPS%NDTPRECCUR2,&
 & NGPTOT=>YDGEM%NGPTOT,&
 & NGPBLKS=>YDDIM%NGPBLKS)
!     --------------------------------------------------------------------------

!    ------------------------------------------------------------------
!     1 - Initialisations
!    - --------------------------------------------------------------------

INIT0=-1

IF (INIT0 == 0) THEN
  ZVALUE=HUGE(1._JPRB)
  ZVALUE_ONE=HUGE(1._JPRB)
  ZVALUE_T=HUGE(1._JPRB)
  ZVALUE_P=HUGE(1._JPRB)
  ZVALUE_L=HUGE(1._JPRB)
  ZVALUE_EPSILON=HUGE(1._JPRB)
ELSE
  ZVALUE=0._JPRB
  ZVALUE_ONE=1._JPRB
  ZVALUE_T=293._JPRB
  ZVALUE_P=101325._JPRB
  ZVALUE_L=0.01_JPRB
  ZVALUE_EPSILON=1E-12_JPRB
ENDIF

LLSWAP_THS=.TRUE. ! it can be as well true as false, actually : this is just to start up the swapp process
LLSWAP_RS=.TRUE. ! it can be as well true as false, actually : this is just to start up the swapp process
LLSWAP_SVS=.TRUE. ! it can be as well true as false, actually : this is just to start up the swapp process
LLSWAP_SVM=.TRUE. ! it can be as well true as false, actually : this is just to start up the swapp process
LLSWAP_LIMAS=.TRUE. ! it can be as well true as false, actually : this is just to start up the swapp process

!        1.0 numerical safety

IF (JPRD == JPRB) THEN
  ZEPSNEB=1.E-12
ELSE
  ZEPSNEB=1.E-06
ENDIF

NSV=0

!         1.3 time step initialisation
!             the mesoNH physics (turb and microphysics) is written 
!             for leap frog scheme
!             !!! be carefull for 2TL or 3TL 

IF (LTWOTL) THEN
  ZDT=PDT/2._JPRB
ELSE
  IF (KSTEP/=0) THEN
    ZDT=PDT/2._JPRB
  ELSE
    ZDT=PDT
  ENDIF
ENDIF

ZINVDT=1/PDT

ZINVG=1._JPRB/RG 

! initialisation de ZDTMSE
IF (LDXFUMSE) THEN
  ZDTMSE=0.01_JPRB
  ZSTATI=REAL(RSTATI,JPRB)-ZDTMSE/2._JPRB
ELSE
  ZDTMSE=PDT
  ZSTATI=REAL(RSTATI,JPRB)
ENDIF

IF(LTWOTL) THEN
  ZRHGMT=REAL(RHGMT,JPRB)-ZDTMSE/2._JPRB
ELSE
  ZRHGMT=REAL(RHGMT,JPRB)
ENDIF


LLMSE=LMSE.AND.(NSURFEXCTL >= 2)
LLMSE_PARAM=LLMSE
LLMSE_DIAG=LLMSE.AND.(NSURFEXCTL >= 3)


!  Vertical points
IKA=KLEV
IKB=KLEV
IKU=1
IKT=KLEV
IKTE=KLEV
IKTB=1
IKL=-1

!  SETUP

IF (INIT0 >= 0) THEN

  ZUM__(:,:)=ZVALUE
  ZUS__(:,:)=ZVALUE
  ZVM__(:,:)=ZVALUE
  ZVS__(:,:)=ZVALUE
  ZWM__(:,:)=ZVALUE
  ZTHSWAP__(:,:)=ZVALUE
  ZTHSAVE__(:,:)=ZVALUE
  ZSRCS__(:,:)=ZVALUE
  ZSIGS__(:,:)=ZVALUE
  ZTHM__(:,:)=ZVALUE_T
  ZRHODREFM__(:,:)=ZVALUE_ONE
  ZPABSM__(:,:)=ZVALUE_P
  ZTURB3D__(:,:,:)=ZVALUE
  ZLENGTHM__(:,:)=ZVALUE_L
  ZLENGTHH__(:,:)=ZVALUE_L

  ZZZ_(:,:)=ZVALUE
  ZMFM_(:,:)=ZVALUE
  ZSIGM_(:,:)=ZVALUE
  ZNEBMNH_(:,:)=ZVALUE
  ZEVAP_(:,:)=ZVALUE

  ZRSWAP_(:,:,:)=ZVALUE
  ZRSAVE_(:,:,:)=ZVALUE

  ZRM_(:,:,:)=ZVALUE

  ZLIMASWAP_(:,:,:)=ZVALUE
  ZLIMASAVE_(:,:,:)=ZVALUE

  ZLIMAM_(:,:,:)=ZVALUE

  ZFLXZTHVMF_(:,:)=ZVALUE
  ZSIGMF_(:,:)=ZVALUE
  ZRC_MF_(:,:)=ZVALUE
  ZRI_MF_(:,:)=ZVALUE
  ZCF_MF_(:,:)=ZVALUE

  ZSVSWAP_(:,:,:)=ZVALUE
  ZSVSAVE_(:,:,:)=ZVALUE

  ZSVMSWAP_(:,:,:)=ZVALUE
  ZSVMSAVE_(:,:,:)=ZVALUE
  ZSVMB_(:,:)=ZVALUE

  ZPIZA_DST_(:,:,:)  = ZVALUE
  ZCGA_DST_(:,:,:)   = ZVALUE
  ZTAUREL_DST_(:,:,:) = ZVALUE_EPSILON

  ZAERD_(:,:)=ZVALUE

  ZMFS_(:,:)=ZVALUE

  ZFPR(:,:,:)=ZVALUE

  IF(LKFBCONV) THEN
    ZCVTENDRV_(:,:)=ZVALUE
    ZCVTENDRC_(:,:)=ZVALUE
    ZCVTENDRI_(:,:)=ZVALUE
    ZCVTENDT_(:,:)=ZVALUE
  ENDIF

  ZACPRG_(:)=ZVALUE
  ZINPRR_NOTINCR_(:)=ZVALUE
  ZINPRS_NOTINCR_(:)=ZVALUE
  ZINPRG_NOTINCR_(:)=ZVALUE
  ZINPRH_NOTINCR_(:)=ZVALUE

  ZTHLS_(:,:)=ZVALUE
  ZMFUS_(:,:)=ZVALUE
  ZMFVS_(:,:)=ZVALUE

  ZTKEEDMF(:,:)=ZVALUE

  ZSFTH_(:)=ZVALUE
  ZSFRV_(:)=ZVALUE
  ZSFU_(:)=ZVALUE
  ZSFV_(:)=ZVALUE

  ZSFCO2_(:)=ZVALUE
  ZEMIS(:)=ZVALUE_ONE

  ZQICE(:,:)=ZVALUE
  ZQLIQ(:,:)=ZVALUE
  ZQO3(:,:)=ZVALUE

  ZAER(:,:,:)=ZVALUE
  ZAERINDS(:,:)=ZVALUE

  ZQSAT(:,:)=ZVALUE
  ZLH(:,:)=ZVALUE
  ZLSCPE(:,:)=ZVALUE
  ZGEOSLC(:,:)=ZVALUE

  ZFRSOFS(:)=ZVALUE

  ZQW(:,:)=ZVALUE
  PRH(:,:)=ZVALUE
  ZTW(:,:)=ZVALUE

  ZTRSODIF(:,:)=ZVALUE
  ZTRSODIR(:,:)=ZVALUE
  ZZS_FSWDIR(:,:)=ZVALUE
  ZZS_FSWDIF(:,:)=ZVALUE

  ZSDUR(:)=ZVALUE
  ZDSRP(:)=ZVALUE
  ZALBD(:,:)=ZVALUE
  ZALBP(:,:)=ZVALUE
  ZALBD1(:)=ZVALUE
  ZALBP1(:)=ZVALUE

  ZQV(:,:)=ZVALUE
  ZSFSV_(:,:)=ZVALUE

  PFRSOC(:,:)=ZVALUE
  PFRTHC(:,:)=ZVALUE
  PFRSOPS(:)=ZVALUE
  PDIAGH(:)=ZVALUE

  ZTENDSV_TURBLIMA_(:,:,:)=ZVALUE

ENDIF

!  INITIALIZE (CUMULATED) TENDENCIES

DO JLEV=1,KLEV
  PTENDT(KIDIA:KFDIA,JLEV)=0.0_JPRB
  PTENDU(KIDIA:KFDIA,JLEV)=0.0_JPRB
  PTENDV(KIDIA:KFDIA,JLEV)=0.0_JPRB
  PTENDW(KIDIA:KFDIA,JLEV)=0.0_JPRB
  PTENDTKE(KIDIA:KFDIA,JLEV)=0.0_JPRB
ENDDO
DO JRR=1,NRR
  DO JLEV=1,KLEV
    PTENDR(KIDIA:KFDIA,JLEV,JRR)=0.0_JPRB
  ENDDO
ENDDO
DO JGFL=1,NGFL_EXT
  DO JLEV=1,KLEV
    PTENDEXT(KIDIA:KFDIA,JLEV,JGFL)=0.0_JPRB
  ENDDO
ENDDO
DO JGFL=1,NLIMA
  DO JLEV=1,KLEV
    PTENDLIMA(KIDIA:KFDIA,JLEV,JGFL)=0.0_JPRB
  ENDDO
ENDDO

!  INITIALIZE CUMULATED STUFF

! Small array, OK. REK
ZINPRH_(KIDIA:KFDIA)=0._JPRB
ZINPRR_(KIDIA:KFDIA)=0._JPRB
ZACPRR_(KIDIA:KFDIA)=0._JPRB
ZINPRS_(KIDIA:KFDIA)=0._JPRB
ZACPRS_(KIDIA:KFDIA)=0._JPRB
ZINPRG_(KIDIA:KFDIA)=0._JPRB


DO JLEV = 1,KLEV
  DO JLON = KIDIA,KFDIA
    ZZZ_F_(JLON,JLEV)=PAPHIFM(JLON,JLEV)*ZINVG
    ZTENDT(JLON,JLEV)=0._JPRB
  ENDDO
ENDDO

! adhoc solution to avoid negative tke values
! when SL advective ddh is activated 
IF (LRSLDDH) THEN
  DO JLEV=1, KLEV
    DO JLON = KIDIA, KFDIA
      ZTKEM4SLDDH(JLON,JLEV)=MAX(PTKEM(JLON,JLEV),PPTKEMIN)
    ENDDO
  ENDDO
  ZTKEM => ZTKEM4SLDDH(:,:)
ELSE
  ZTKEM => PTKEM(:,:)
  !test TKE > 0.
  IF (MINVAL(ZTKEM(KIDIA:KFDIA,1:KLEV)) <= 0._JPRB) THEN
    CALL ABOR1('TKE < 0 under APL_AROME check YTKE_NL%NREQIN')
  ENDIF
ENDIF


!test invalid combinations
IF (LHARATU .AND. CMF_UPDRAFT == 'EDKF') THEN
  CALL ABOR1('Combination LHARATU and EDKF not valid!')
ENDIF

!initialisation of first useful field for EZDIAG use in Chemistry/Dust
IOFF_MFSHAL=1
IF(LFPREC3D) IOFF_MFSHAL=2

!    ------------------------------------------------------------------
!     2 - CHANGEMENTS DE VARIABLES ET INVERSION DES NIVEAUX 
!     --------------------------------------------------------------------

IF (LMICRO.OR.LTURB.OR.LLMSE.OR.LKFBCONV) THEN

  ZRSCP=RD/RCPD
  ZINVATM=1/RATM
  !initialisation de ZZZ_
  DO JLEV = 1,KLEV
   !initialisation de qdm (utile localement pour calculer rho  
   !et convertir q en r 
    IF (NRR==7) THEN
      DO JLON = KIDIA,KFDIA
        ZQDM(JLON,JLEV)=1._JPRB-PQVM(JLON,JLEV)-PQCM(JLON,JLEV)-PQRM(JLON,JLEV)&
         & -PQIM(JLON,JLEV)-PQSM(JLON,JLEV)-PQGM(JLON,JLEV)-PQHM(JLON,JLEV)   
      ENDDO
    ELSE
      DO JLON = KIDIA,KFDIA
        ZQDM(JLON,JLEV)=1._JPRB-PQVM(JLON,JLEV)-PQCM(JLON,JLEV)-PQRM(JLON,JLEV)&
         & -PQIM(JLON,JLEV)-PQSM(JLON,JLEV)-PQGM(JLON,JLEV)
      ENDDO
    ENDIF
    DO JLON = KIDIA,KFDIA
   !initialisation de ZRHODREFM__ (=qd*zrho) 
      ZRHO=PAPRSFM(JLON,JLEV)/(PRM(JLON,JLEV)*PTM(JLON,JLEV))
      ZRHODREFM__(JLON,JLEV)=ZRHO*ZQDM(JLON,JLEV)
      ZRHODJM__(JLON,JLEV)=PDELPM(JLON,JLEV)*ZINVG
  !initialisation de ZEXNREFM_
      ZEXNREFM_(JLON,JLEV)=(PAPRSFM(JLON,JLEV)*ZINVATM)**(ZRSCP)
 ! vent horizontal et TKE
      ZPABSM__(JLON,JLEV)=PAPRSFM(JLON,JLEV)
      ZUM__(JLON,JLEV)= PUM(JLON,JLEV)
      ZVM__(JLON,JLEV)= PVM(JLON,JLEV)
      ZWM__(JLON,JLEV)= PWM(JLON,JLEV)
      ZTKEM__(JLON,JLEV)= ZTKEM(JLON,JLEV)
      ZZZ_(JLON,JLEV)=PAPHIM(JLON,JLEV)*ZINVG
    ENDDO
  ENDDO
  !initialise sigma for subgrid condensation coming
  !from previous time step turbulence scheme
  IF (LOSIGMAS) THEN
    DO JLEV = 1, KLEV 
      ZSIGM_(KIDIA:KFDIA,JLEV)= PSIGM(KIDIA:KFDIA,JLEV)
    ENDDO
  ENDIF
  !initialise convective mas flux for subgrid condensation coming 
  !from previous time step convection scheme
  IF ( LKFBCONV.AND.LOSUBG_COND.AND..NOT.LOSIGMAS) THEN
    DO JLEV = 1, KLEV 
      ZMFM_(KIDIA:KFDIA,JLEV)=PSIGM(KIDIA:KFDIA,JLEV) 
    ENDDO
  ENDIF
!!! initialisation des variables d etat MNH �t

  !initialisation de ZRM_ pour les hydrometeores (ri=qi/qd)
  DO JLEV = 1, KLEV
    DO JLON = KIDIA,KFDIA
      ZTHM__(JLON,JLEV)=PTM(JLON,JLEV)/ZEXNREFM_(JLON,JLEV)
      ZRM_(JLON,JLEV,1)=PQVM(JLON,JLEV)/ZQDM(JLON,JLEV)  
      ZRM_(JLON,JLEV,2)=PQCM(JLON,JLEV)/ZQDM(JLON,JLEV)  
      ZRM_(JLON,JLEV,3)=PQRM(JLON,JLEV)/ZQDM(JLON,JLEV)  
      ZRM_(JLON,JLEV,4)=PQIM(JLON,JLEV)/ZQDM(JLON,JLEV)  
      ZRM_(JLON,JLEV,5)=PQSM(JLON,JLEV)/ZQDM(JLON,JLEV)  
      ZRM_(JLON,JLEV,6)=PQGM(JLON,JLEV)/ZQDM(JLON,JLEV)  
    ENDDO
  ENDDO
  
  IF (NRR==7) THEN
    DO JLEV = 1, KLEV
      DO JLON = KIDIA,KFDIA
        ZRM_(JLON,JLEV,7)=PQHM(JLON,JLEV)/ZQDM(JLON,JLEV)  
      ENDDO
    ENDDO
  ENDIF
   
  IF (NRR==6) THEN
    !initialisation de ZTHVREFM__
    DO JLEV = 1, KLEV 
      DO JLON = KIDIA,KFDIA
        ZTHVREFM__(JLON,JLEV)=ZTHM__(JLON,JLEV)*&
         & (1._JPRB+ZRM_(JLON,JLEV,1)*(RV/RD))/&
         & (1._JPRB+ZRM_(JLON,JLEV,1)+ZRM_(JLON,JLEV,2) +&
         & ZRM_(JLON,JLEV,3)+ZRM_(JLON,JLEV,4)+&
         & ZRM_(JLON,JLEV,5)+ZRM_(JLON,JLEV,6))
      ENDDO
    ENDDO
  ELSEIF (NRR==7) THEN
    DO JLEV = 1, KLEV 
      DO JLON = KIDIA,KFDIA
        ZTHVREFM__(JLON,JLEV)=ZTHM__(JLON,JLEV)*&
         & (1._JPRB+ZRM_(JLON,JLEV,1)*(RV/RD))/&
         & (1._JPRB+ZRM_(JLON,JLEV,1)+ZRM_(JLON,JLEV,2) +&
         & ZRM_(JLON,JLEV,3)+ZRM_(JLON,JLEV,4)+&
         & ZRM_(JLON,JLEV,5)+ZRM_(JLON,JLEV,6)+&
         &  ZRM_(JLON,JLEV,7) )
      ENDDO
    ENDDO
  ENDIF

!!! initialisation des variables d etat MNH a t+dt
!!! division pas le pas de temps
!!!(la multiplication par rhodj est faite plus tard, si necessaire, 
!!! suivant les parametrisations)   

  ! initialise pointers :
  CALL SWAP_THS
  ! vent horizontal
  DO JLEV = 1, KLEV 
    DO JLON=KIDIA,KFDIA
      ZUS__(JLON,JLEV)= PUM(JLON,JLEV)*ZINVDT
      ZVS__(JLON,JLEV)= PVM(JLON,JLEV)*ZINVDT
      ZWS__(JLON,JLEV)= PWM(JLON,JLEV)*ZINVDT
      ZTKES_(JLON,JLEV)= ZTKEM(JLON,JLEV)*ZINVDT
      ZTHS__(JLON,JLEV)=ZTHM__(JLON,JLEV)*ZINVDT
    ENDDO
  ENDDO

  !initialisation de ZRS_ pour les hydrometeores
  ! initialise pointers :
  CALL SWAP_RS
  DO JRR=1,NRR 
    DO JLEV = 1, KLEV
      DO JLON=KIDIA,KFDIA
        ZRS_(JLON,JLEV,JRR)=ZRM_(JLON,JLEV,JRR)*ZINVDT
      ENDDO
    ENDDO
  ENDDO

!!! Initialisations temporaires d'arguments non-utilises
  !initialisation de ZCIT_
  ZCIT_(KIDIA:KFDIA,1:IKT)=0.0_JPRB
  
  !initialisation des tableaux de precipitations inst. and cumulated 
  !and surface fluxes for turbulence
  IF (LLMSE.OR.LSFORCS) THEN
    ZACPRR_(KIDIA:KFDIA)=PGPAR(KIDIA:KFDIA,MACPRR)
    ZACPRS_(KIDIA:KFDIA)=PGPAR(KIDIA:KFDIA,MACPRS)
    ZACPRG_(KIDIA:KFDIA)=PGPAR(KIDIA:KFDIA,MACPRG)
    ZINPRR_NOTINCR_(KIDIA:KFDIA)=PGPAR(KIDIA:KFDIA,MINPRR)
    ZINPRS_NOTINCR_(KIDIA:KFDIA)=PGPAR(KIDIA:KFDIA,MINPRS)
    ZINPRG_NOTINCR_(KIDIA:KFDIA)=PGPAR(KIDIA:KFDIA,MINPRG)
  ENDIF

  !initialisation des scalaires passifs
  ! initialise pointers :
  CALL SWAP_SVM
  CALL SWAP_SVS
  DO JGFL=1,NGFL_EXT
    DO JLEV = 1, KLEV
      DO JLON=KIDIA,KFDIA
        ZSVM_(JLON,JLEV,JGFL)=PSVM(JLON,JLEV,JGFL)
        ZSVS_(JLON,JLEV,JGFL)=PSVM(JLON,JLEV,JGFL)*ZINVDT
      ENDDO
    ENDDO
  ENDDO

  !initialisation des concentrations LIMA
  ! initialise pointers :
  CALL SWAP_LIMAS
  DO JGFL=1,NLIMA
    DO JLEV = 1, KLEV
      DO JLON=KIDIA,KFDIA
        ZLIMAM_(JLON,JLEV,JGFL)=PLIMAM(JLON,JLEV,JGFL)
        ZLIMAS_(JLON,JLEV,JGFL)=PLIMAM(JLON,JLEV,JGFL)*ZINVDT
      ENDDO
    ENDDO
  ENDDO

ENDIF

! daand: radflex
ZFRSO => PFRSO(:,:,1)
ZFRTH => PFRTH(:,:,1)

!    ------------------------------------------------------------------
!     3 - PRINTS FOR DIAGNOSTICS IF NEEDED
!    ------------------------------------------------------------------
IF (LDIAGWMAX) THEN
  IF (MOD(KSTEP+1,NDIAGWMAX)==0) THEN
  ! calcul de wmax
    DO JLEV = 1 , KLEV
      Z_WMAX=0._JPRB
      Z_WMIN=0._JPRB
      DO JLON=KIDIA,KFDIA
        IF (PWM(JLON,JLEV)>Z_WMAX) THEN
          Z_WMAX=PWM(JLON,JLEV)
        ENDIF
        IF (PWM(JLON,JLEV)<Z_WMIN) THEN
          Z_WMIN=PWM(JLON,JLEV)
        ENDIF
      ENDDO
    ENDDO
  ENDIF
ENDIF

IF (LFLEXDIA) THEN
  !save tendencies
  ZTENDTBAK(KIDIA:KFDIA,1:KLEV)=PTENDT(KIDIA:KFDIA,1:KLEV)
  DO JR=1,NRR
    ZTENDRBAK(KIDIA:KFDIA,1:KLEV,JR)=PTENDR(KIDIA:KFDIA,1:KLEV,JR)
  ENDDO
ENDIF


!    ------------------------------------------------------------------
!     4 - ADJUSTMENT (CALLED IF THE MICROPHYSICS IS SWITCH ON) 
!    ------------------------------------------------------------------

IF (LMICRO) THEN

  ! Swap pointers because input values of THS and RS should be saved
  CALL SWAP_THS
  CALL SWAP_RS

  IF (LMFSHAL .AND. (CMF_CLOUD=='DIRE'.OR.CMF_CLOUD=='BIGA')) THEN
    IOFF_MFSHAL=IOFF_MFSHAL+3
    IF (KSTEP==0) THEN
      DO JLEV = 1, KLEV 
        ZRC_MF_(KIDIA:KFDIA,JLEV)=0._JPRB
        ZRI_MF_(KIDIA:KFDIA,JLEV)=0._JPRB
        ZCF_MF_(KIDIA:KFDIA,JLEV)=0._JPRB
      ENDDO
    ELSE
      DO JLEV = 1, KLEV 
        ZRC_MF_(KIDIA:KFDIA,JLEV)=PEZDIAG(KIDIA:KFDIA,JLEV,1)
        ZRI_MF_(KIDIA:KFDIA,JLEV)=PEZDIAG(KIDIA:KFDIA,JLEV,3)
        ZCF_MF_(KIDIA:KFDIA,JLEV)=PEZDIAG(KIDIA:KFDIA,JLEV,2)
      ENDDO
    ENDIF
    PEZDIAG(KIDIA:KFDIA,1:KLEV,1:3)=0._JPRB
  ENDIF

  IF (MOD(KSTEP+1,NPRINTFR)==0) THEN
    WRITE(NULOUT,*)'avant aro_adjust sous apl_arome'
    WRITE(NULOUT,*)'JLEV   ZZZ_F_   RHODJM   EXNREFM       PABSM       THM      SIGM         MFM    '
    DO JLEV=1,KLEV+1 
      WRITE(NULOUT,'(I2,X,7F10.3)') JLEV,ZZZ_F_(NPTP,JLEV),ZRHODJM__(NPTP,JLEV),&
       & ZEXNREFM_(NPTP,JLEV), ZPABSM__(NPTP,JLEV), ZTHM__(NPTP,JLEV), ZSIGM_(NPTP,JLEV), ZMFM_(NPTP,JLEV)
    ENDDO 
    WRITE (NULOUT,*)'JLEV  rhoQv  rhoQc   rhoQr   rhoQi   rhoQs   rhoQg'
    DO JLEV=1,KLEV
      WRITE(NULOUT,'(I2,X,6E11.4)') JLEV,ZRM_(NPTP,JLEV,1),&
       & ZRM_(NPTP,JLEV,2), ZRM_(NPTP,JLEV,3),ZRM_(NPTP,JLEV,4),ZRM_(NPTP,JLEV,5), ZRM_(NPTP,JLEV,6)  
    ENDDO
    WRITE (NULOUT,*)'JLEV  ZRC_MF_  ZRI_MF_  ZCF_MF_ ZTHSIN_  ZSRCS__ ZNEBMNH_'
    DO JLEV=1,KLEV
      WRITE(NULOUT,'(I2,X,6E11.4)') JLEV,ZRC_MF_(NPTP,JLEV),&
       & ZRI_MF_(NPTP,JLEV),ZCF_MF_(NPTP,JLEV), ZTHSIN_(NPTP,JLEV),ZSRCS__(NPTP,JLEV), ZNEBMNH_(NPTP,JLEV)
    ENDDO
  ENDIF

  ! for now a copy is needed (see below, inside). I don't like than :-( REK
  ZTHS__(KIDIA:KFDIA,1:KLEV)=ZTHSIN_(KIDIA:KFDIA,1:KLEV)
  ZRS_(KIDIA:KFDIA,1:KLEV,1:NRR)=ZRSIN_(KIDIA:KFDIA,1:KLEV,1:NRR)

  IF (CMICRO == 'LIMA') THEN

    CALL SWAP_LIMAS
    ! for now a copy is needed (see below, inside). I don't like than :-( REK
    ZLIMAS_(KIDIA:KFDIA,1:KLEV,1:NLIMA)=ZLIMASIN_(KIDIA:KFDIA,1:KLEV,1:NLIMA)

    CALL ARO_ADJUST_LIMA (KLEV,IKU,IKL,KFDIA,KLEV,NRR,NLIMA,KSTEP+1,&
     & LOSUBG_COND, LOSIGMAS, LOCND2, &
     & ZDT,VSIGQSAT,ZZZ_F_,&
     & ZRHODJM__(:,1:KLEV),&
     & ZRHODREFM__(:,1:KLEV),&
     & ZEXNREFM_,&
     & ZPABSM__(:,1:KLEV),&
     & ZTHM__(:,1:KLEV),&
     & ZRM_,&
     & ZLIMAM_,&
     & ZSIGM_,&
     & ZMFM_,ZRC_MF_,&
     & ZRI_MF_,ZCF_MF_,&
     & ZTHS__(:,1:KLEV),ZRS_,&
     & ZLIMAS_,&
     & ZSRCS__(:,1:KLEV),ZNEBMNH_,&
     & YDDDH, YDMODEL%YRML_DIAG%YRLDDH, YDMODEL%YRML_DIAG%YRMDDH)
  ELSE

    CALL ARO_ADJUST (KLEV,IKU,IKL,KFDIA,KLEV,NRR,&
     & NGFL_EZDIAG, &
     & CFRAC_ICE_ADJUST, CCONDENS, CLAMBDA3, LOSUBG_COND, &
     & LOSIGMAS, CMICRO, LOCND2, CSUBG_MF_PDF, &
     & ZDT,VSIGQSAT,ZZZ_F_,&
     & ZRHODJM__(:,1:KLEV),&
     & ZEXNREFM_,&
     & ZRHODREFM__(:,1:KLEV),&
     & ZPABSM__(:,1:KLEV),&
     & ZTHM__(:,1:KLEV),&
     & ZRM_,ZSIGM_,&
     & ZMFM_,ZRC_MF_,&
     & ZRI_MF_,ZCF_MF_,&
     & ZTHS__(:,1:KLEV),ZRS_,&
     & ZSRCS__(:,1:KLEV),ZNEBMNH_,&
     & ZHLC_HRC__(:,1:KLEV), ZHLC_HCF__(:,1:KLEV),&
     & ZHLI_HRI__(:,1:KLEV), ZHLI_HCF__(:,1:KLEV),&
     & PGP2DSPP,PEZDIAG,&
     & YDDDH, YDMODEL%YRML_DIAG%YRLDDH, YDMODEL%YRML_DIAG%YRMDDH)

  ENDIF
  
  IF (MOD(KSTEP+1,NPRINTFR)==0) THEN
    WRITE(NULOUT,*)'apres aro_adjust sous apl_arome'
    WRITE(NULOUT,*)'JLEV   ZZZ_F_   RHODJM   EXNREFM       PABSM       THM      SIGM         MFM    '   
    DO JLEV=1,KLEV+1 
      WRITE(NULOUT,'(I2,X,7F10.3)') JLEV,ZZZ_F_(NPTP,JLEV),ZRHODJM__(NPTP,JLEV),&
       & ZEXNREFM_(NPTP,JLEV), ZPABSM__(NPTP,JLEV), ZTHM__(NPTP,JLEV), ZSIGM_(NPTP,JLEV), ZMFM_(NPTP,JLEV)
    ENDDO 
    WRITE (NULOUT,*)'JLEV  rhoQv  rhoQc   rhoQr   rhoQi   rhoQs   rhoQg'
    DO JLEV=1,KLEV
      WRITE(NULOUT,'(I2,X,6E11.4)') JLEV,ZRS_(NPTP,JLEV,1),&
       & ZRS_(NPTP,JLEV,2), ZRS_(NPTP,JLEV,3),ZRS_(NPTP,JLEV,4),ZRS_(NPTP,JLEV,5), ZRS_(NPTP,JLEV,6)  
    ENDDO
    WRITE (NULOUT,*)'JLEV  ZRC_MF_  ZRI_MF_  ZCF_MF_ ZTHS__ ZSRCS__ ZNEBMNH_'
    DO JLEV=1,KLEV
      WRITE(NULOUT,'(I2,X,6E11.4)')JLEV,ZRC_MF_(NPTP,JLEV),&
       & ZRI_MF_(NPTP,JLEV),ZCF_MF_(NPTP,JLEV), ZTHS__(NPTP,JLEV),ZSRCS__(NPTP,JLEV), ZNEBMNH_(NPTP,JLEV)
    ENDDO
  ENDIF

  DO JLEV=1,KLEV
    PCLFS(KIDIA:KFDIA,JLEV)=ZNEBMNH_(KIDIA:KFDIA,JLEV) 
  ENDDO

  !adjusted zthm and zrm
  DO JLEV = 1, KLEV
    DO JLON = KIDIA,KFDIA
      ZTHM__(JLON,JLEV)=ZTHS__(JLON,JLEV)*PDT
    ENDDO
  ENDDO

  DO JRR=1,NRR
    DO JLEV = 1, KLEV
      DO JLON = KIDIA,KFDIA
        ZRM_(JLON,JLEV,JRR)=ZRS_(JLON,JLEV,JRR)*PDT
      ENDDO
    ENDDO
  ENDDO


  !initialisation de qdm utile pour 
  !convertir tendance de r en tendance de q 
  IF (NRR==6) THEN
    DO JLEV=1,KLEV
      DO JLON= KIDIA, KFDIA
        ZQDM(JLON,JLEV)=1._JPRB/(1._JPRB+ZRM_(JLON,JLEV,1)+&
        &ZRM_(JLON,JLEV,2)+ZRM_(JLON,JLEV,3)+ZRM_(JLON,JLEV,4)+ZRM_(JLON,JLEV,5)+&
        &ZRM_(JLON,JLEV,6) )
      ENDDO
    ENDDO
  ELSEIF (NRR==7) THEN
    DO JLEV=1,KLEV
      DO JLON= KIDIA, KFDIA
        ZQDM(JLON,JLEV)=1._JPRB/(1._JPRB+ZRM_(JLON,JLEV,1)+&
        &ZRM_(JLON,JLEV,2)+ZRM_(JLON,JLEV,3)+ZRM_(JLON,JLEV,4)+ZRM_(JLON,JLEV,5)+&
        &ZRM_(JLON,JLEV,6)+ZRM_(JLON,JLEV,7) )
      ENDDO
    ENDDO
  ENDIF 
  !reinitialisation des qi
  DO JLEV = 1, KLEV
    DO JLON = KIDIA,KFDIA
      ZQVM(JLON,JLEV)=ZRM_(JLON,JLEV,1)*ZQDM(JLON,JLEV)
      ZQCM(JLON,JLEV)=ZRM_(JLON,JLEV,2)*ZQDM(JLON,JLEV)
      ZQRM(JLON,JLEV)=ZRM_(JLON,JLEV,3)*ZQDM(JLON,JLEV)
      ZQIM(JLON,JLEV)=ZRM_(JLON,JLEV,4)*ZQDM(JLON,JLEV)
      ZQSM(JLON,JLEV)=ZRM_(JLON,JLEV,5)*ZQDM(JLON,JLEV)
      ZQGM(JLON,JLEV)=ZRM_(JLON,JLEV,6)*ZQDM(JLON,JLEV)
    ENDDO
  ENDDO

  IF (NRR==7) THEN
    DO JLEV = 1,KLEV
      DO JLON = KIDIA,KFDIA
        ZQHM(JLON,JLEV)=ZRM_(JLON,JLEV,7)*ZQDM(JLON,JLEV)
      ENDDO
    ENDDO
  ELSE
    ZQHM(KIDIA:KFDIA,1:KLEV)=0._JPRB
  ENDIF

  ! Tendances des variables LIMA
  DO JGFL=1,NLIMA
    DO JLEV = 1, KLEV
      ! Réinitialisation des variables LIMA
      DO JLON=KIDIA,KFDIA
        ZLIMAM_(JLON,JLEV,JGFL)=ZLIMAS_(JLON,JLEV,JGFL)*PDT
        PTENDLIMA(JLON,JLEV,JGFL)=PTENDLIMA(JLON,JLEV,JGFL)+(ZLIMAS_(JLON,JLEV,JGFL)-ZLIMASIN_(JLON,JLEV,JGFL))  
      ENDDO
    ENDDO
  ENDDO

  !modif de R et CP
  ZQHGM(KIDIA:KFDIA,:)=ZQHM(KIDIA:KFDIA,:)+ZQGM(KIDIA:KFDIA,:)
  CALL GPRCP(KLON,KIDIA,KFDIA,KLEV,PQ=ZQVM,PQI=ZQIM,PQL=ZQCM,PQR=ZQRM,PQS=ZQSM,PQG=ZQHGM,PCP=ZCPM,PR=ZRHM)  

  DO JLEV = 1,KLEV
    DO JLON = KIDIA,KFDIA
      ZTM(JLON,JLEV)=ZTHM__(JLON,JLEV)*ZEXNREFM_(JLON,JLEV)
  !reinitialisation de ZRHODREFM__ (=qd*zrho)
      ZRHO=PAPRSFM(JLON,JLEV)/(ZRHM(JLON,JLEV)*ZTM(JLON,JLEV))
      ZRHODREFM__(JLON,JLEV)=ZRHO*ZQDM(JLON,JLEV)
    ENDDO
  ENDDO

  !geopotentiel calculation
 
  ZAPHIM(KIDIA:KFDIA,KLEV)=PAPHIM(KIDIA:KFDIA,KLEV)
  CALL GPGEO(KLON,KIDIA,KFDIA,KLEV,ZAPHIM,ZAPHIFM,ZTM,ZRHM,PLNPRM,PALPHM,YDGEOMETRY%YRVERT_GEOM)
   
  !calcul de l'altitude
  DO JLEV = 1, KLEV
    DO JLON = KIDIA,KFDIA
      ZZZ_(JLON,JLEV)=ZAPHIM(JLON,JLEV)*ZINVG
      !initialisation de ZZZ_F_
      ZZZ_F_(JLON,JLEV)=ZAPHIFM(JLON,JLEV)*ZINVG
      ! tendency of T
      PTENDT(JLON,JLEV)=PTENDT(JLON,JLEV)+(ZTHS__(JLON,JLEV)-ZTHSIN_(JLON,JLEV))*ZEXNREFM_(JLON,JLEV)
      ZTENDT(JLON,JLEV)=ZTHS__(JLON,JLEV)-ZTHSIN_(JLON,JLEV)
    ENDDO
  ENDDO
  !inversion niveaux tendances des ri et conversion en qi en multipliant par qd 
  DO JR=1,NRR
    DO JLEV=1,KLEV
      DO JLON=KIDIA,KFDIA
        PTENDR(JLON,JLEV,JR)=PTENDR(JLON,JLEV,JR)+(ZRS_(JLON,JLEV,JR)-ZRSIN_(JLON,JLEV,JR))*ZQDM(JLON,JLEV)
      ENDDO
    ENDDO
  ENDDO
  !initialisation de ZDZZ_
  DO JLON = KIDIA,KFDIA
    ZDZZ_(JLON,1)=ZAPHIM(JLON,0)*ZINVG-ZZZ_(JLON,1)
  ENDDO
  DO JLEV = 2, KLEV
    DO JLON = KIDIA,KFDIA
      ZDZZ_(JLON,JLEV)=ZZZ_(JLON,JLEV+IKL)-ZZZ_(JLON,JLEV)
    ENDDO
  ENDDO

ELSE

  ZTM (KIDIA:KFDIA,1:KLEV)=PTM(KIDIA:KFDIA,1:KLEV)
  ZRHM(KIDIA:KFDIA,1:KLEV)=PRM(KIDIA:KFDIA,1:KLEV)
  ZQVM(KIDIA:KFDIA,1:KLEV)=PQVM(KIDIA:KFDIA,1:KLEV)
  ZQIM(KIDIA:KFDIA,1:KLEV)=PQIM(KIDIA:KFDIA,1:KLEV)
  ZQCM(KIDIA:KFDIA,1:KLEV)=PQCM(KIDIA:KFDIA,1:KLEV)
  ZQRM(KIDIA:KFDIA,1:KLEV)=PQRM(KIDIA:KFDIA,1:KLEV)
  ZQSM(KIDIA:KFDIA,1:KLEV)=PQSM(KIDIA:KFDIA,1:KLEV)
  ZQGM(KIDIA:KFDIA,1:KLEV)=PQGM(KIDIA:KFDIA,1:KLEV)
  IF (NRR==7) THEN
    ZQHM(KIDIA:KFDIA,1:KLEV)=PQHM(KIDIA:KFDIA,1:KLEV)
  ELSE
    ZQHM(KIDIA:KFDIA,1:KLEV)=0._JPRB
  ENDIF
  ZCPM(KIDIA:KFDIA,KTDIA:KLEV)=PCPM(KIDIA:KFDIA,KTDIA:KLEV)
  ZAPHIM(KIDIA:KFDIA,0:KLEV)=PAPHIM(KIDIA:KFDIA,0:KLEV)
  ZAPHIFM(KIDIA:KFDIA,1:KLEV)=PAPHIFM(KIDIA:KFDIA,1:KLEV)
  ZZZ_(KIDIA:KFDIA,1:KLEV)=PAPHIM(KIDIA:KFDIA,1:KLEV)*ZINVG
  !initialisation of PCLFS outside LMICRO to be zero in case LMICRO=F
  PCLFS(KIDIA:KFDIA,1:KLEV)=0._JPRB

ENDIF ! ADJUSTMENT LMICRO

IF (LFLEXDIA) THEN
  DO JLEV = 1, KLEV
    DO JLON=KIDIA,KFDIA 
      ZTMPAF(JLON,JLEV)=(PTENDT(JLON,JLEV)-ZTENDTBAK(JLON,JLEV))*PDELPM(JLON,JLEV)*ZINVG*ZCPM(JLON,JLEV)
    ENDDO
  ENDDO
  IF (LDDH_OMP) THEN
    CALL NEW_ADD_FIELD_3D(YDMODEL%YRML_DIAG%YRMDDH,ZTMPAF,'TCTADJU',YDDDH)
  ELSE
    CALL ADD_FIELD_3D(YLDDH,ZTMPAF,'TCTADJU','T','ARP',.TRUE.,.TRUE.)
  ENDIF
  DO JR=1,NRR
    CLNAME='T'//CLVARNAME(JR)//'ADJU'
    DO JLEV = 1, KLEV
      DO JLON=KIDIA,KFDIA 
        ZTMPAF(JLON,JLEV)=(PTENDR(JLON,JLEV,JR)-ZTENDRBAK(JLON,JLEV,JR))*PDELPM(JLON,JLEV)*ZINVG
      ENDDO
    ENDDO
    IF (LDDH_OMP) THEN
      CALL NEW_ADD_FIELD_3D(YDMODEL%YRML_DIAG%YRMDDH,ZTMPAF,CLNAME,YDDDH)
    ELSE
      CALL ADD_FIELD_3D(YLDDH,ZTMPAF,CLNAME,'T','ARP',.TRUE.,.TRUE.)
    ENDIF
    DO JLEV = 1, KLEV
      DO JLON=KIDIA,KFDIA 
        ZTMPAF(JLON,JLEV)=PCLFS(JLON,JLEV)*PDELPM(JLON,JLEV)
      ENDDO
    ENDDO
    IF (LDDH_OMP) THEN
      CALL NEW_ADD_FIELD_3D(YDMODEL%YRML_DIAG%YRMDDH,ZTMPAF,'VNT',YDDDH)
    ELSE
      CALL ADD_FIELD_3D(YLDDH,ZTMPAF,'VNT','V','ARP',.TRUE.,.TRUE.)
    ENDIF
  ENDDO
! specific to new data flow for diagnostics
  IF (LDDH_OMP) THEN
    DO JLEV = 1, KLEV
      ZCON1(KIDIA:KFDIA,JLEV) = 1.0_JPRB
      ZCON2(KIDIA:KFDIA,JLEV) = ZQDM(KIDIA:KFDIA,JLEV)
    ENDDO
    DO JLEV = 1, KLEV
      DO JLON=KIDIA,KFDIA 
        ZCON3(JLON,JLEV) = PCPM(JLON,JLEV)*ZEXNREFM_(JLON,JLEV)
      ENDDO
    ENDDO
    ! missing interface !!! REK
    CALL ARO_SUINTBUDGET_OMP(KLON,KLEV,KSTEP,ZCON1,ZCON2,ZCON3,YDDDH)
  ELSE
    ! missing interface !!! REK
    CALL ARO_SUINTBUDGET(KLON,KLEV,KSTEP,KFDIA,ZQDM,ZEXNREFM_,PCPM)
  ENDIF
ENDIF


DO JLEV = 1, KLEV-1
  DO JLON = KIDIA,KFDIA
    ZDZZ_F_(JLON,JLEV)=ZZZ_F_(JLON,JLEV)-ZZZ_F_(JLON,JLEV-IKL)
  ENDDO
ENDDO
DO JLON = KIDIA,KFDIA
  ZDZZ_F_(JLON,KLEV)=ZZZ_F_(JLON,KLEV)-POROG(JLON)*ZINVG
ENDDO


!     --------------------------------------------------------------------
!     5 - COMPUTE DUST PROPERTIES FOR RADIATION IF LRDUST=T
!     --------------------------------------------------------------------
IF (LRDUST) THEN
  PEZDIAG(KIDIA:KFDIA,1:KLEV,IOFF_MFSHAL:NGFL_EZDIAG)=0.0_JPRB
  ! input dust scalar concentration in ppp from
  IEZDIAG_CHEM=NGFL_EZDIAG-IOFF_MFSHAL+1
  ! Swapp because IN and OUT will be needed simultaneously
  CALL SWAP_SVM
  ! input dust scalar concentration in ppp from
  CALL ARO_MNHDUST (IKL,KFDIA,KLEV,NGFL_EXT, PDT,ZSVMIN_,ZZZ_,ZDZZ_,&
   & ZPABSM__(:,1:KLEV),ZTHM__(:,1:KLEV),ZRHODREFM__(:,1:KLEV),&
   & NSWB_MNH,KSTEP+1,ZSVM_,ZPIZA_DST_,ZCGA_DST_,ZTAUREL_DST_,ZAERD_,IEZDIAG_CHEM,&
   & ZPEZDIAG_(:,:,IOFF_MFSHAL:NGFL_EZDIAG)           )
  PEZDIAG(KIDIA:KFDIA,1:KLEV,IOFF_MFSHAL:NGFL_EZDIAG)=ZPEZDIAG_(KIDIA:KFDIA,1:KLEV,IOFF_MFSHAL:NGFL_EZDIAG)
! return to tendency
  DO JGFL=1, NGFL_EXT
    DO JLEV = 1,KLEV
      DO JLON = KIDIA,KFDIA
        PTENDEXT(JLON,JLEV,JGFL)=PTENDEXT(JLON,JLEV,JGFL)+(ZSVM_(JLON,JLEV,JGFL)-ZSVMIN_(JLON,JLEV,JGFL))*ZINVDT
      ENDDO
    ENDDO
  ENDDO
ENDIF ! LRDUST

IF (LSFORCS) THEN   ! <== Surface forcing for MUSC

  ZTSURF(KIDIA:KFDIA) = PTM(KIDIA:KFDIA,KLEV)
  ZTN(KIDIA:KFDIA)    = PTM(KIDIA:KFDIA,KLEV)
  ZQS(KIDIA:KFDIA)    = PQVM(KIDIA:KFDIA,KLEV)
  DO JLON=KIDIA,KFDIA
    ZRHODREFM(JLON) = PAPRSFM(JLON,KLEV)/(PTM(JLON,KLEV)*PRM(JLON,KLEV))
    ZTHETAS(JLON)   = ZTSURF(JLON)*(RATM/PAPRSM(JLON,KLEV))**RKAPPA
  ENDDO

  LLAROME=.TRUE.
  CALL SURF_IDEAL_FLUX(YDRIP,YDPHY0,YDPHYDS, LLAROME, KIDIA , KFDIA  , KLON, PAPHIFM(:,KLEV), &
   & ZRHODREFM, PSFORC,ZTN,ZTSURF,PLSM,PQVM(:,KLEV), PUM(:,KLEV), PVM(:,KLEV), ZTHETAS, &
   & ZSFTH_, ZSFRV_, ZSFU_, ZSFV_)

!* Compute PBL-diagnostics
   
   ZCAPE(:)=0._JPRB
   ZDCAPE(:)=0._JPRB
  CALL ACCLDIA(YDXFU,YDPHY,YDMODEL%YRML_PHY_MF%YRPHY2,YDTOPH,KIDIA,KFDIA,KLON,KLEV,PUCLS,&
   & PVCLS,PUM(:,1:KLEV),PVM(:,1:KLEV),ZCAPE,ZDCAPE,ZTKEM(:,1:KLEV),PAPHIFM(:,1:KLEV),POROG,&
   & PUGST,PVGST,PPBLH,ICLPH)

  PPBLH(KIDIA:KFDIA)=MIN(XMAXLM,MAX(XMINLM,PPBLH(KIDIA:KFDIA))) 

ENDIF    ! <== End of surface forcing for MUSC

!     --------------------------------------------------------------------
!     6 - RADIATION LRAYFM (IFS) or LRAY (ACRANEB2) 
!     --------------------------------------------------------------------
IF (LRAYFM.OR.LRAY) THEN
  ! prepare some input for both radiation schemes at every time step

  ! test de coherence sur le nombre de bandes spectrales entre ce qui sort de
  ! la surface et ce qu'attend le rayonnement
  IF( NSWB_MNH /= NSW) THEN
    CALL ABOR1 (' NSWB_MNH must be equal to NSW !')
  ENDIF

  ! compute saturated specific humidity
  CALL ACTQSAT ( YDPHY,KIDIA,KFDIA,KLON,NTQSAT,KLEV, PAPRSFM, ZCPM, ZQVM, ZTM,&
   & ZGEOSLC, ZLH, ZLSCPE, ZQSAT, ZQW, PRH, ZTW)  

  IF (YSPP_CONFIG%LSPP.AND.YSPP_CONFIG%LPERT_RADGR) THEN
   IF (YSPP_CONFIG%LLNN_MEAN1.OR.YSPP_CONFIG%LLNN_MEAN1_RADGR) THEN
    ZMU = -0.5_JPRB * (YSPP_CONFIG%CMPERT_RADGR * YSPP_CONFIG%SDEV)**2
   ELSE
    ZMU = 0._JPRB
   ENDIF
   DO JLON=KIDIA,KFDIA
    ZVAL = RADGR*EXP(ZMU+YSPP_CONFIG%CMPERT_RADGR*PGP2DSPP(JLON,YSPP%MP_RADGR))
    ZRADGR(JLON) = MAX(YSPP_CONFIG%CLIP_RADGR(1),MIN(ZVAL,YSPP_CONFIG%CLIP_RADGR(2)))
   ENDDO
   IF ( YSPP_CONFIG%IEZDIAG_POS > 0 ) THEN
    JKO=2*YSPP%MP_RADGR-1
    JKE=2*YSPP%MP_RADGR
    DO JLON=KIDIA,KFDIA
     PEZDIAG(JLON,JKO,YSPP_CONFIG%IEZDIAG_POS) = PGP2DSPP(JLON,YSPP%MP_RADGR)
     PEZDIAG(JLON,JKE,YSPP_CONFIG%IEZDIAG_POS) = ZRADGR(JLON)
    ENDDO
   ENDIF
  ELSE
   DO JLON=KIDIA,KFDIA
    ZRADGR(JLON) = RADGR
   ENDDO
  ENDIF

  IF (YSPP_CONFIG%LSPP.AND.YSPP_CONFIG%LPERT_RADSN) THEN
   IF (YSPP_CONFIG%LLNN_MEAN1.OR.YSPP_CONFIG%LLNN_MEAN1_RADSN) THEN
    ZMU = -0.5_JPRB * (YSPP_CONFIG%CMPERT_RADSN * YSPP_CONFIG%SDEV)**2
   ELSE
    ZMU = 0._JPRB
   ENDIF
   DO JLON=KIDIA,KFDIA
    ZVAL = RADSN*EXP(ZMU+YSPP_CONFIG%CMPERT_RADSN*PGP2DSPP(JLON,YSPP%MP_RADSN))
    ZRADSN(JLON) = MAX(YSPP_CONFIG%CLIP_RADSN(1),MIN(ZVAL,YSPP_CONFIG%CLIP_RADSN(2)))
   ENDDO
   IF ( YSPP_CONFIG%IEZDIAG_POS > 0 ) THEN
    JKO=2*YSPP%MP_RADSN-1
    JKE=2*YSPP%MP_RADSN
    DO JLON=KIDIA,KFDIA
     PEZDIAG(JLON,JKO,YSPP_CONFIG%IEZDIAG_POS) = PGP2DSPP(JLON,YSPP%MP_RADSN)