Changeset 1299 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/fisrtilp.F
- Timestamp:
- Jan 20, 2010, 3:27:21 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/fisrtilp.F
r1279 r1299 1 1 ! 2 ! $ Header$2 ! $Id$ 3 3 ! 4 4 c … … 135 135 PRINT*, 'fisrtilp, evap_prec:', evap_prec 136 136 PRINT*, 'fisrtilp, cpartiel:', cpartiel 137 IF (ABS(dtime/ FLOAT(ninter)-360.0).GT.0.001) THEN137 IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN 138 138 PRINT*, 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime 139 139 PRINT*, 'Je prefere un sous-intervalle de 6 minutes' … … 436 436 zfice(i) = zfice(i)**nexpo 437 437 zneb(i) = MAX(rneb(i,k), seuil_neb) 438 radliq(i,k) = zoliq(i)/ FLOAT(ninter+1)438 radliq(i,k) = zoliq(i)/REAL(ninter+1) 439 439 ENDIF 440 440 ENDDO … … 453 453 zcl =cld_lc_con 454 454 zct =1./cld_tau_con 455 zfroi = dtime/ FLOAT(ninter)/zdz(i)*zoliq(i)455 zfroi = dtime/REAL(ninter)/zdz(i)*zoliq(i) 456 456 . *fallvc(zrhol(i)) * zfice(i) 457 457 else 458 458 zcl =cld_lc_lsc 459 459 zct =1./cld_tau_lsc 460 zfroi = dtime/ FLOAT(ninter)/zdz(i)*zoliq(i)460 zfroi = dtime/REAL(ninter)/zdz(i)*zoliq(i) 461 461 . *fallvs(zrhol(i)) * zfice(i) 462 462 endif 463 zchau = zct *dtime/ FLOAT(ninter) * zoliq(i)463 zchau = zct *dtime/REAL(ninter) * zoliq(i) 464 464 . *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl )**2)) *(1.-zfice(i)) 465 465 ztot = zchau + zfroi … … 468 468 ztot = MIN(ztot,zoliq(i)) 469 469 zoliq(i) = MAX(zoliq(i)-ztot , 0.0) 470 radliq(i,k) = radliq(i,k) + zoliq(i)/ FLOAT(ninter+1)470 radliq(i,k) = radliq(i,k) + zoliq(i)/REAL(ninter+1) 471 471 ENDIF 472 472 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.