| 1 | ! |
|---|
| 2 | ! $Id: calcaerosolstrato_rrtm.f90 5338 2024-11-20 09:05:48Z evignon $ |
|---|
| 3 | ! |
|---|
| 4 | SUBROUTINE calcaerosolstrato_rrtm(pplay,t_seri,paprs,debut) |
|---|
| 5 | |
|---|
| 6 | USE clesphys_mod_h |
|---|
| 7 | USE iniprint_mod_h |
|---|
| 8 | USE phys_state_var_mod, ONLY: tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, tau_aero_lw_rrtm |
|---|
| 9 | USE phys_local_var_mod, ONLY: mdw, tausum_aero, tausum_strat, tau_strat_550, tau_strat_1020, stratomask |
|---|
| 10 | USE aero_mod |
|---|
| 11 | USE dimphy |
|---|
| 12 | ! USE temps_mod |
|---|
| 13 | USE yomcst_mod_h |
|---|
| 14 | |
|---|
| 15 | !USE paramet_mod_h |
|---|
| 16 | IMPLICIT NONE |
|---|
| 17 | |
|---|
| 18 | |
|---|
| 19 | |
|---|
| 20 | |
|---|
| 21 | ! Variable input |
|---|
| 22 | REAL,DIMENSION(klon,klev),INTENT(IN) :: t_seri ! Temperature |
|---|
| 23 | REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le mileu de chaque couche (en Pa) |
|---|
| 24 | LOGICAL,INTENT(IN) :: debut ! le flag de l'initialisation de la physique |
|---|
| 25 | REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa) |
|---|
| 26 | |
|---|
| 27 | ! Stratospheric aerosols optical properties |
|---|
| 28 | REAL, DIMENSION(klon,klev,nbands_sw_rrtm) :: tau_strat, piz_strat, cg_strat |
|---|
| 29 | REAL, DIMENSION(klon,klev,nwave_sw+nwave_lw) :: tau_strat_wave |
|---|
| 30 | REAL, DIMENSION(klon,klev,nbands_lw_rrtm) :: tau_lw_abs_rrtm |
|---|
| 31 | |
|---|
| 32 | INTEGER k, band, wave, i |
|---|
| 33 | REAL zrho, zdz |
|---|
| 34 | |
|---|
| 35 | !--calculate optical properties of the aerosol size distribution from tr_seri |
|---|
| 36 | tau_strat=0.0 |
|---|
| 37 | piz_strat=0.0 |
|---|
| 38 | cg_strat=0.0 |
|---|
| 39 | tau_strat_wave=0.0 |
|---|
| 40 | tau_lw_abs_rrtm=0.0 |
|---|
| 41 | |
|---|
| 42 | CALL miecalc_aer(tau_strat, piz_strat, cg_strat, tau_strat_wave, tau_lw_abs_rrtm, paprs, debut) |
|---|
| 43 | |
|---|
| 44 | !!--test CK: deactivate radiative effect of aerosol |
|---|
| 45 | ! tau_strat=0.0 |
|---|
| 46 | ! piz_strat=0.0 |
|---|
| 47 | ! cg_strat=0.0 |
|---|
| 48 | ! tau_strat_wave=0.0 |
|---|
| 49 | ! tau_lw_abs_rrtm=0.0 |
|---|
| 50 | |
|---|
| 51 | !--test CK: deactivate SW radiative effect of aerosol (but leave LW) |
|---|
| 52 | ! tau_strat=0.0 |
|---|
| 53 | ! piz_strat=0.0 |
|---|
| 54 | ! cg_strat=0.0 |
|---|
| 55 | |
|---|
| 56 | ! DO wave=1, nwave_sw |
|---|
| 57 | ! tau_strat_wave(:,:,wave)=0.0 |
|---|
| 58 | ! ENDDO |
|---|
| 59 | |
|---|
| 60 | !--test CK: deactivate LW radiative effect of aerosol (but leave SW) |
|---|
| 61 | ! tau_lw_abs_rrtm=0.0 |
|---|
| 62 | |
|---|
| 63 | ! DO wave=nwave_sw+1, nwave_sw+nwave_lw |
|---|
| 64 | ! tau_strat_wave(:,:,wave)=0.0 |
|---|
| 65 | ! ENDDO |
|---|
| 66 | |
|---|
| 67 | !--total vertical aod at the 5 SW + 1 LW wavelengths |
|---|
| 68 | DO wave=1, nwave_sw+nwave_lw |
|---|
| 69 | DO k=1, klev |
|---|
| 70 | tausum_aero(:,wave,id_STRAT_phy)=tausum_aero(:,wave,id_STRAT_phy)+tau_strat_wave(:,k,wave) |
|---|
| 71 | ENDDO |
|---|
| 72 | ENDDO |
|---|
| 73 | |
|---|
| 74 | !--weighted average for cg, piz and tau, adding strat aerosols on top of tropospheric ones |
|---|
| 75 | DO band=1, nbands_sw_rrtm |
|---|
| 76 | !--no stratospheric aerosol in index 1 |
|---|
| 77 | cg_aero_sw_rrtm(:,:,1,band) = cg_aero_sw_rrtm(:,:,2,band) |
|---|
| 78 | piz_aero_sw_rrtm(:,:,1,band) = piz_aero_sw_rrtm(:,:,2,band) |
|---|
| 79 | tau_aero_sw_rrtm(:,:,1,band) = tau_aero_sw_rrtm(:,:,2,band) |
|---|
| 80 | |
|---|
| 81 | !--tropospheric and stratospheric aerosol in index 2 |
|---|
| 82 | cg_aero_sw_rrtm(:,:,2,band) = ( cg_aero_sw_rrtm(:,:,2,band)*piz_aero_sw_rrtm(:,:,2,band)*tau_aero_sw_rrtm(:,:,2,band) + & |
|---|
| 83 | cg_strat(:,:,band)*piz_strat(:,:,band)*tau_strat(:,:,band) ) / & |
|---|
| 84 | MAX( piz_aero_sw_rrtm(:,:,2,band)*tau_aero_sw_rrtm(:,:,2,band) + & |
|---|
| 85 | piz_strat(:,:,band)*tau_strat(:,:,band), 1.e-15 ) |
|---|
| 86 | piz_aero_sw_rrtm(:,:,2,band)= ( piz_aero_sw_rrtm(:,:,2,band)*tau_aero_sw_rrtm(:,:,2,band) + & |
|---|
| 87 | piz_strat(:,:,band)*tau_strat(:,:,band) ) / & |
|---|
| 88 | MAX( tau_aero_sw_rrtm(:,:,2,band) + tau_strat(:,:,band), 1.e-15 ) |
|---|
| 89 | tau_aero_sw_rrtm(:,:,2,band)= tau_aero_sw_rrtm(:,:,2,band) + tau_strat(:,:,band) |
|---|
| 90 | ENDDO |
|---|
| 91 | |
|---|
| 92 | DO band=1, nbands_lw_rrtm |
|---|
| 93 | !--no stratospheric aerosols in index 1 |
|---|
| 94 | tau_aero_lw_rrtm(:,:,1,band) = tau_aero_lw_rrtm(:,:,2,band) |
|---|
| 95 | !--tropospheric and stratospheric aerosol in index 2 |
|---|
| 96 | tau_aero_lw_rrtm(:,:,2,band) = tau_aero_lw_rrtm(:,:,2,band) + tau_lw_abs_rrtm(:,:,band) |
|---|
| 97 | ENDDO |
|---|
| 98 | |
|---|
| 99 | WHERE (tau_aero_sw_rrtm .LT. 1.e-14) piz_aero_sw_rrtm=1.0 |
|---|
| 100 | WHERE (tau_aero_sw_rrtm .LT. 1.e-14) tau_aero_sw_rrtm=1.e-15 |
|---|
| 101 | WHERE (tau_aero_lw_rrtm .LT. 1.e-14) tau_aero_lw_rrtm=1.e-15 |
|---|
| 102 | |
|---|
| 103 | tausum_strat(:,:)=0.0 |
|---|
| 104 | DO i=1,klon |
|---|
| 105 | DO k=1,klev |
|---|
| 106 | IF (stratomask(i,k).GT.0.5) THEN |
|---|
| 107 | tausum_strat(i,1)=tausum_strat(i,1)+tau_strat_wave(i,k,2) !--550 nm |
|---|
| 108 | tausum_strat(i,2)=tausum_strat(i,2)+tau_strat_wave(i,k,5) !--1020 nm |
|---|
| 109 | tausum_strat(i,3)=tausum_strat(i,3)+tau_strat_wave(i,k,6) !--10 um |
|---|
| 110 | ENDIF |
|---|
| 111 | ENDDO |
|---|
| 112 | ENDDO |
|---|
| 113 | |
|---|
| 114 | DO i=1,klon |
|---|
| 115 | DO k=1,klev |
|---|
| 116 | zrho=pplay(i,k)/t_seri(i,k)/RD !air density in kg/m3 |
|---|
| 117 | zdz=(paprs(i,k)-paprs(i,k+1))/zrho/RG !thickness of layer in m |
|---|
| 118 | tau_strat_550(i,k)=tau_strat_wave(i,k,2)/zdz |
|---|
| 119 | tau_strat_1020(i,k)=tau_strat_wave(i,k,5)/zdz |
|---|
| 120 | ENDDO |
|---|
| 121 | ENDDO |
|---|
| 122 | |
|---|
| 123 | END SUBROUTINE calcaerosolstrato_rrtm |
|---|
