Changeset 5081 for LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer
- Timestamp:
- Jul 19, 2024, 4:15:44 PM (7 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/cond_evap_tstep_mod.F90
r4950 r5081 227 227 ! H2SO4 mass fraction in aerosol 228 228 WH2=R2SO4*1.0E-2 229 IF(WH2 .EQ.0.0) RETURN229 IF(WH2==0.0) RETURN 230 230 ! ACTIVITY COEFFICIENT(SEE GIAUQUE,1951) 231 231 ! AYERS ET AL (1980) … … 324 324 ff(IK,:)=0.0 325 325 DO k=1, nbtr_bin 326 IF (k .LE.(nbtr_bin-1)) THEN327 IF (Vbin_wet(k) .LE.Vnew.AND.Vnew.LT.Vbin_wet(k+1)) THEN326 IF (k<=(nbtr_bin-1)) THEN 327 IF (Vbin_wet(k)<=Vnew.AND.Vnew<Vbin_wet(k+1)) THEN 328 328 ff(IK,k)= Vbin_wet(k)/Vnew*(Vbin_wet(k+1)-Vnew)/(Vbin_wet(k+1)-Vbin_wet(k)) 329 329 ENDIF 330 330 ENDIF 331 IF (k .EQ.1.AND.Vnew.LE.Vbin_wet(k)) THEN331 IF (k==1.AND.Vnew<=Vbin_wet(k)) THEN 332 332 ff(IK,k)= 1. 333 333 ENDIF 334 IF (k .GT.1) THEN335 IF (Vbin_wet(k-1) .LT.Vnew.AND.Vnew.LT.Vbin_wet(k)) THEN334 IF (k>1) THEN 335 IF (Vbin_wet(k-1)<Vnew.AND.Vnew<Vbin_wet(k)) THEN 336 336 ff(IK,k)= 1.-ff(IK,k-1) 337 337 ENDIF 338 338 ENDIF 339 IF (k .EQ.nbtr_bin.AND.Vnew.GE.Vbin_wet(k)) THEN339 IF (k==nbtr_bin.AND.Vnew>=Vbin_wet(k)) THEN 340 340 ff(IK,k)= 1. 341 341 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.