Changeset 2695 for LMDZ5/trunk/libf/phylmd/StratAer/cond_evap_tstep_mod.F90
- Timestamp:
- Nov 1, 2016, 11:19:45 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/StratAer/cond_evap_tstep_mod.F90
r2690 r2695 18 18 USE aerophys 19 19 USE infotrac 20 USE YOMCST 20 USE YOMCST, ONLY : RPI 21 21 22 22 IMPLICIT NONE … … 139 139 USE aerophys 140 140 USE infotrac 141 USE YOMCST, ONLY : RPI 141 142 142 143 IMPLICIT NONE 143 144 include "YOMCST.h"145 144 146 145 ! input variables … … 175 174 ff(IK,:)=0.0 176 175 DO k=1, nbtr_bin 177 IF (k.LE.(nbtr_bin-1).AND.Vbin_wet(k).LE.Vnew.AND.Vnew.LT.Vbin_wet(k+1)) THEN 178 ff(IK,k)= Vbin_wet(k)/Vnew*(Vbin_wet(k+1)-Vnew)/(Vbin_wet(k+1)-Vbin_wet(k)) 176 IF (k.LE.(nbtr_bin-1)) THEN 177 IF (Vbin_wet(k).LE.Vnew.AND.Vnew.LT.Vbin_wet(k+1)) THEN 178 ff(IK,k)= Vbin_wet(k)/Vnew*(Vbin_wet(k+1)-Vnew)/(Vbin_wet(k+1)-Vbin_wet(k)) 179 ENDIF 179 180 ENDIF 180 181 IF (k.EQ.1.AND.Vnew.LE.Vbin_wet(k)) THEN 181 182 ff(IK,k)= 1. 182 183 ENDIF 183 IF (k.GT.1.AND.Vbin_wet(k-1).LT.Vnew.AND.Vnew.LT.Vbin_wet(k)) THEN 184 ff(IK,k)= 1.-ff(IK,k-1) 184 IF (k.GT.1) THEN 185 IF (Vbin_wet(k-1).LT.Vnew.AND.Vnew.LT.Vbin_wet(k)) THEN 186 ff(IK,k)= 1.-ff(IK,k-1) 187 ENDIF 185 188 ENDIF 186 189 IF (k.EQ.nbtr_bin.AND.Vnew.GE.Vbin_wet(k)) THEN
Note: See TracChangeset
for help on using the changeset viewer.