Ignore:
Timestamp:
Oct 22, 2024, 2:29:31 PM (6 weeks ago)
Author:
abarral
Message:

Wrap uses of cpp key DUST

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r5252 r5253  
    100100!
    101101!
    102 #ifdef CPP_Dust
    103102    USE phytracr_spl_mod, ONLY: phytracr_spl, phytracr_spl_out_init
    104103    USE phys_output_write_spl_mod
    105 #else
    106104    USE phytrac_mod, ONLY : phytrac_init, phytrac
    107105    USE phys_output_write_mod
    108 #endif
    109106
    110107
     
    433430       reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
    434431       USE output_physiqex_mod, ONLY: output_physiqex
    435        USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA, CPPKEY_REPROBUS, CPPKEY_STRATAER
     432       USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA, CPPKEY_REPROBUS, CPPKEY_STRATAER, CPPKEY_DUST
    436433
    437434
     
    16431640       iflag_phytrac = 1 ! by default we do want to call phytrac
    16441641       CALL getin_p('iflag_phytrac',iflag_phytrac)
    1645 #ifdef CPP_Dust
    1646        IF (iflag_phytrac.EQ.0) THEN 
     1642IF (CPPKEY_DUST) THEN
     1643       IF (iflag_phytrac.EQ.0) THEN
    16471644         WRITE(lunout,*) 'In order to run with SPLA, iflag_phytrac will be forced to 1'
    16481645         iflag_phytrac = 1
    16491646       ENDIF
    1650 #endif
     1647END IF
    16511648       nvm_lmdz = 13
    16521649       CALL getin_p('NVM',nvm_lmdz)
     
    20472044!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    20482045
    2049 #ifdef CPP_Dust
     2046IF (CPPKEY_DUST) THEN
    20502047       ! Quand on utilise SPLA, on force iflag_phytrac=1
    20512048       CALL phytracr_spl_out_init()
     
    20562053                                ptconvth, d_t, qx, d_qx, d_tr_dyn, zmasse,      &
    20572054                                flag_aerosol, flag_aerosol_strat, ok_cdnc)
    2058 #else
     2055ELSE
    20592056       ! phys_output_write écrit des variables traceurs seulement si iflag_phytrac == 1
    20602057       ! donc seulement dans ce cas on doit appeler phytrac_init()
     
    20682065                              ptconvth, d_u, d_t, qx, d_qx, zmasse,           &
    20692066                              flag_aerosol, flag_aerosol_strat, ok_cdnc, t, u1, v1)
    2070 #endif
     2067END IF
    20712068
    20722069
     
    57185715                   !--new aerosol properties SW and LW
    57195716                   !
    5720 #ifdef CPP_Dust
     5717IF (CPPKEY_DUST) THEN
    57215718                   !--SPL aerosol model
    57225719                   CALL splaerosol_optic_rrtm( ok_alw, pplay, paprs, t_seri, rhcl, &
     
    57245721                        tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,  &
    57255722                        tausum_aero, tau3d_aero)
    5726 #else
     5723ELSE
    57275724                   !--climatologies or INCA aerosols
    57285725                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, ok_alw, ok_volcan, &
     
    57325729                        tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,  &
    57335730                        tausum_aero, drytausum_aero, tau3d_aero)
    5734 #endif
     5731END IF
    57355732
    57365733                   IF (flag_aerosol .EQ. 7) THEN
     
    69026899    ENDIF
    69036900
    6904 #ifdef CPP_Dust
     6901IF (CPPKEY_DUST) THEN
    69056902    !  Avec SPLA, iflag_phytrac est forcé =1
    69066903    CALL       phytracr_spl ( debut,lafin , jD_cur,jH_cur,iflag_con,       &  ! I
     
    69196916                      d_tr_dyn,tr_seri)
    69206917
    6921 #else
     6918ELSE
    69226919    IF (iflag_phytrac == 1 ) THEN
    69236920      CALL phytrac ( &
     
    69626959    ENDIF    ! (iflag_phytrac=1)
    69636960
    6964 #endif
     6961END IF
    69656962    !ENDIF    ! (iflag_phytrac=1)
    69666963
     
    73617358    !On effectue les sorties:
    73627359
    7363 #ifdef CPP_Dust
     7360IF (CPPKEY_DUST) THEN
    73647361  CALL phys_output_write_spl(itap, pdtphys, paprs, pphis,  &
    73657362       pplay, lmax_th, aerosol_couple,                 &
     
    73687365       ptconvth, d_t, qx, d_qx, d_tr_dyn, zmasse,      &
    73697366       flag_aerosol, flag_aerosol_strat, ok_cdnc)
    7370 #else
     7367ELSE
    73717368    CALL phys_output_write(itap, pdtphys, paprs, pphis,  &
    73727369         pplay, lmax_th, aerosol_couple,                 &
     
    73757372         ptconvth, d_u, d_t, qx, d_qx, zmasse,           &
    73767373         flag_aerosol, flag_aerosol_strat, ok_cdnc,t, u1, v1)
    7377 #endif
     7374END IF
    73787375
    73797376#ifndef CPP_XIOS
Note: See TracChangeset for help on using the changeset viewer.