- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/lsc_scav_orig.F90
r2630 r5082 177 177 ENDDO 178 178 179 IF (it .gt.1) THEN ! aerosol179 IF (it>1) THEN ! aerosol 180 180 frac_ev=frac_aer 181 181 ELSE ! gas … … 183 183 ENDIF 184 184 185 IF(it .gt.1) then ! aerosol185 IF(it>1) then ! aerosol 186 186 DO k=1, klev 187 187 DO i=1, klon … … 203 203 ! incloud scavenging 204 204 ! if(inscav_fisrt) then 205 if (iflag_lscav .eq.4) then205 if (iflag_lscav == 4) then 206 206 beta=beta_fisrt(i,k)*rneb(i,k) 207 207 else … … 219 219 220 220 ! below-cloud impaction 221 IF(it .eq.1) then221 IF(it==1) then 222 222 d_tr_bcscav(i,k,it)=0. 223 223 ELSE … … 236 236 deltaP(i,k)=max(deltaP(i,k),0.) 237 237 238 if(flxr(i,k+1)+flxs(i,k+1) .gt.1.e-16) then238 if(flxr(i,k+1)+flxs(i,k+1)>1.e-16) then 239 239 beta_ev(i,k)=deltaP(i,k)/(flxr(i,k+1)+flxs(i,k+1)) 240 240 else … … 246 246 !jyg 247 247 248 if(abs(1-(1-frac_ev)*beta_ev(i,k)) .gt.1.e-16) then248 if(abs(1-(1-frac_ev)*beta_ev(i,k))>1.e-16) then 249 249 ! remove tracers from precipitation owing to release by evaporation in his_dh 250 250 ! dxev=frac_ev*beta_ev(i,k)*his_dh(i) *pdtime/(zrho(i,k)*zdz(i,k)) &
Note: See TracChangeset
for help on using the changeset viewer.