Changeset 5253 for LMDZ6/trunk/libf/phylmdiso
- Timestamp:
- Oct 22, 2024, 2:29:31 PM (6 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r5252 r5253 100 100 ! 101 101 ! 102 #ifdef CPP_Dust103 102 USE phytracr_spl_mod, ONLY: phytracr_spl, phytracr_spl_out_init 104 103 USE phys_output_write_spl_mod 105 #else106 104 USE phytrac_mod, ONLY : phytrac_init, phytrac 107 105 USE phys_output_write_mod 108 #endif109 106 110 107 … … 433 430 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra 434 431 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 436 433 437 434 … … 1643 1640 iflag_phytrac = 1 ! by default we do want to call phytrac 1644 1641 CALL getin_p('iflag_phytrac',iflag_phytrac) 1645 #ifdef CPP_Dust 1646 IF (iflag_phytrac.EQ.0) THEN 1642 IF (CPPKEY_DUST) THEN 1643 IF (iflag_phytrac.EQ.0) THEN 1647 1644 WRITE(lunout,*) 'In order to run with SPLA, iflag_phytrac will be forced to 1' 1648 1645 iflag_phytrac = 1 1649 1646 ENDIF 1650 #endif 1647 END IF 1651 1648 nvm_lmdz = 13 1652 1649 CALL getin_p('NVM',nvm_lmdz) … … 2047 2044 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2048 2045 2049 #ifdef CPP_Dust 2046 IF (CPPKEY_DUST) THEN 2050 2047 ! Quand on utilise SPLA, on force iflag_phytrac=1 2051 2048 CALL phytracr_spl_out_init() … … 2056 2053 ptconvth, d_t, qx, d_qx, d_tr_dyn, zmasse, & 2057 2054 flag_aerosol, flag_aerosol_strat, ok_cdnc) 2058 #else 2055 ELSE 2059 2056 ! phys_output_write écrit des variables traceurs seulement si iflag_phytrac == 1 2060 2057 ! donc seulement dans ce cas on doit appeler phytrac_init() … … 2068 2065 ptconvth, d_u, d_t, qx, d_qx, zmasse, & 2069 2066 flag_aerosol, flag_aerosol_strat, ok_cdnc, t, u1, v1) 2070 #endif 2067 END IF 2071 2068 2072 2069 … … 5718 5715 !--new aerosol properties SW and LW 5719 5716 ! 5720 #ifdef CPP_Dust 5717 IF (CPPKEY_DUST) THEN 5721 5718 !--SPL aerosol model 5722 5719 CALL splaerosol_optic_rrtm( ok_alw, pplay, paprs, t_seri, rhcl, & … … 5724 5721 tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, & 5725 5722 tausum_aero, tau3d_aero) 5726 #else 5723 ELSE 5727 5724 !--climatologies or INCA aerosols 5728 5725 CALL readaerosol_optic_rrtm( debut, aerosol_couple, ok_alw, ok_volcan, & … … 5732 5729 tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, & 5733 5730 tausum_aero, drytausum_aero, tau3d_aero) 5734 #endif 5731 END IF 5735 5732 5736 5733 IF (flag_aerosol .EQ. 7) THEN … … 6902 6899 ENDIF 6903 6900 6904 #ifdef CPP_Dust 6901 IF (CPPKEY_DUST) THEN 6905 6902 ! Avec SPLA, iflag_phytrac est forcé =1 6906 6903 CALL phytracr_spl ( debut,lafin , jD_cur,jH_cur,iflag_con, & ! I … … 6919 6916 d_tr_dyn,tr_seri) 6920 6917 6921 #else 6918 ELSE 6922 6919 IF (iflag_phytrac == 1 ) THEN 6923 6920 CALL phytrac ( & … … 6962 6959 ENDIF ! (iflag_phytrac=1) 6963 6960 6964 #endif 6961 END IF 6965 6962 !ENDIF ! (iflag_phytrac=1) 6966 6963 … … 7361 7358 !On effectue les sorties: 7362 7359 7363 #ifdef CPP_Dust 7360 IF (CPPKEY_DUST) THEN 7364 7361 CALL phys_output_write_spl(itap, pdtphys, paprs, pphis, & 7365 7362 pplay, lmax_th, aerosol_couple, & … … 7368 7365 ptconvth, d_t, qx, d_qx, d_tr_dyn, zmasse, & 7369 7366 flag_aerosol, flag_aerosol_strat, ok_cdnc) 7370 #else 7367 ELSE 7371 7368 CALL phys_output_write(itap, pdtphys, paprs, pphis, & 7372 7369 pplay, lmax_th, aerosol_couple, & … … 7375 7372 ptconvth, d_u, d_t, qx, d_qx, zmasse, & 7376 7373 flag_aerosol, flag_aerosol_strat, ok_cdnc,t, u1, v1) 7377 #endif 7374 END IF 7378 7375 7379 7376 #ifndef CPP_XIOS
Note: See TracChangeset
for help on using the changeset viewer.