Changeset 3607 for LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/radlwsw_m.F90
- Timestamp:
- Nov 29, 2019, 9:08:29 AM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/radlwsw_m.F90
r3604 r3607 16 16 t,q,wo,& 17 17 cldfra, cldemi, cldtaupd,& 18 ok_ade, ok_aie, ok_volcan, flag_ volc_surfstrat, flag_aerosol,&18 ok_ade, ok_aie, ok_volcan, flag_aerosol,& 19 19 flag_aerosol_strat, flag_aer_feedback, & 20 20 tau_aero, piz_aero, cg_aero,& … … 102 102 ! ok_aie---input-L- apply the Aerosol Indirect Effect or not? 103 103 ! ok_volcan-input-L- activate volcanic diags (SW heat & LW cool rate, SW & LW flux) 104 ! flag_volc_surfstrat-input-I- activate volcanic surf cooling or strato heating (or nothing)105 104 ! flag_aerosol-input-I- aerosol flag from 0 to 6 106 105 ! flag_aerosol_strat-input-I- use stratospheric aerosols flag (0, 1, 2) … … 201 200 LOGICAL, INTENT(in) :: ok_volcan ! produce volcanic diags (SW/LW heat flux and rate) 202 201 LOGICAL :: lldebug 203 INTEGER, INTENT(in) :: flag_volc_surfstrat ! allow to impose volcanic cooling rate at surf or heating in strato204 202 INTEGER, INTENT(in) :: flag_aerosol ! takes value 0 (no aerosol) or 1 to 6 (aerosols) 205 203 INTEGER, INTENT(in) :: flag_aerosol_strat ! use stratospheric aerosols … … 321 319 REAL(KIND=8) zswadaero(kdlon,kflev+1) ! SW Aerosol direct forcing 322 320 REAL(KIND=8) zlwadaero(kdlon,kflev+1) ! LW Aerosol direct forcing 323 !-- VolMIP324 REAL(KIND=8) volmip_solsw(kdlon) ! SW clear sky in the case of VOLMIP325 321 !-LW by CK 326 322 REAL(KIND=8) ztoplwadaero(kdlon), zsollwadaero(kdlon) ! LW Aerosol direct forcing at TOAand surface … … 879 875 ZTOPLWAIAERO,ZSOLLWAIAERO, & 880 876 ZLWADAERO, & !--NL 881 volmip_solsw, flag_volc_surfstrat, & !--VOLMIP882 877 ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat, & 883 878 flag_aer_feedback) ! flags aerosols … … 1030 1025 ! print*,'OK2' 1031 1026 1032 !!--add VOLMIP (surf cool or strat heat activate)1033 IF (flag_volc_surfstrat > 0) THEN1034 DO i = 1, kdlon1035 zsolsw(i) = volmip_solsw(i)*fract(i)1036 ENDDO1037 ENDIF1038 1039 1027 ! extrait de SW_AR4 1040 1028 ! DO k = 1, KFLEV
Note: See TracChangeset
for help on using the changeset viewer.