- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_ts2.F
r3792 r5082 207 207 208 208 DO ig=1,knonv 209 IF (isnoSV(ig) .GT.0) THEN210 IF (isnoSV(ig) .GT.1) THEN209 IF (isnoSV(ig)>0) THEN 210 IF (isnoSV(ig)>1) THEN 211 211 mug(ig)=1./(1.+dzsnSV(ig,isnoSV(ig)-1)/dzsnSV(ig,isnoSV(ig))) !mu 212 212 ELSE … … 217 217 ENDIF 218 218 219 IF (mug(ig) .LE.0.05) THEN219 IF (mug(ig) <= 0.05) THEN 220 220 write(*,*)'Attention mu low', mug(ig) 221 221 ENDIF 222 IF (mug(ig) .GE.0.98) THEN222 IF (mug(ig) >= 0.98) THEN 223 223 write(*,*)'Attention mu high', mug(ig) 224 224 ENDIF … … 331 331 332 332 DO ig = 1,knonv 333 IF (ps__SV(ig) .LT.1.) THEN333 IF (ps__SV(ig)<1.) THEN 334 334 ! write(*,*)'ig',ig,'ps',ps__SV(ig) 335 335 ps__SV(ig)=max(ps__SV(ig),1.e-8) 336 336 ENDIF 337 IF (p1l_SV(ig) .LT.1.) THEN337 IF (p1l_SV(ig)<1.) THEN 338 338 ! write(*,*)'ig',ig,'p1l',p1l_SV(ig) 339 339 p1l_SV(ig)=max(p1l_SV(ig),1.e-8) 340 340 ENDIF 341 IF (TaT_SV(ig) .LT.180.) THEN341 IF (TaT_SV(ig)<180.) THEN 342 342 ! write(*,*)'ig',ig,'TaT',TaT_SV(ig) 343 343 TaT_SV(ig)=max(TaT_SV(ig),180.) 344 344 ENDIF 345 IF (QaT_SV(ig) .LT.1.e-8) THEN345 IF (QaT_SV(ig)<1.e-8) THEN 346 346 ! write(*,*)'ig',ig,'QaT',QaT_SV(ig) 347 347 QaT_SV(ig)=max(QaT_SV(ig),1.e-8) 348 348 ENDIF 349 IF (Tsf_SV(ig) .LT.100.) THEN349 IF (Tsf_SV(ig)<100.) THEN 350 350 ! write(*,*)'ig',ig,'Tsf',Tsf_SV(ig) 351 351 Tsf_SV(ig)=max(Tsf_SV(ig),180.) 352 352 ENDIF 353 IF (Tsf_SV(ig) .GT.500.) THEN353 IF (Tsf_SV(ig)>500.) THEN 354 354 ! write(*,*)'ig',ig,'Tsf',Tsf_SV(ig) 355 355 Tsf_SV(ig)=min(Tsf_SV(ig),400.) … … 359 359 ! Tsrf(ig)=max(Tsrf(ig),TaT_SV(ig)-20.) 360 360 ! ENDIF 361 IF (cdH_SV(ig) .LT.1.e-10) THEN361 IF (cdH_SV(ig)<1.e-10) THEN 362 362 ! IF (ig.le.3) write(*,*)'ig',ig,'cdH',cdH_SV(ig) 363 363 cdH_SV(ig)=.5 … … 380 380 & /LhvH2O / zx_pkh(ig) 381 381 ELSE 382 IF (Tsf_SV(ig) .LT.t_coup) THEN382 IF (Tsf_SV(ig)<t_coup) THEN 383 383 zx_qs = qsats(Tsf_SV(ig)) / ps__SV(ig) 384 384 zx_dq_s_dh = dqsats(Tsf_SV(ig),zx_qs)/LhvH2O & … … 408 408 DO ig = 1,knonv 409 409 zx_sl(ig) = LhvH2O 410 IF (Tsf_SV(ig) .LT.RTT) zx_sl(ig) = LhsH2O410 IF (Tsf_SV(ig) < RTT) zx_sl(ig) = LhsH2O 411 411 zx_k1(ig) = zx_coef(ig) 412 412 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.