Changeset 1299 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/fisrtilp_tr.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_tr.F
r766 r1299 1 1 ! 2 ! $ Header$2 ! $Id$ 3 3 ! 4 4 c … … 140 140 PRINT*, 'fisrtilp, evap_prec:', evap_prec 141 141 PRINT*, 'fisrtilp, cpartiel:', cpartiel 142 IF (ABS(dtime/ FLOAT(ninter)-360.0).GT.0.001) THEN142 IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN 143 143 PRINT*, 'fisrtilp: Ce n est pas prevu, voir Z.X.Li', dtime 144 144 PRINT*, 'Je prefere un sous-intervalle de 6 minutes' … … 335 335 zfice(i) = zfice(i)**nexpo 336 336 zneb(i) = MAX(rneb(i,k), seuil_neb) 337 radliq(i,k) = zoliq(i)/ FLOAT(ninter+1)337 radliq(i,k) = zoliq(i)/REAL(ninter+1) 338 338 ENDIF 339 339 ENDDO … … 342 342 DO i = 1, klon 343 343 IF (rneb(i,k).GT.0.0) THEN 344 zchau(i) = ct*dtime/ FLOAT(ninter) * zoliq(i)344 zchau(i) = ct*dtime/REAL(ninter) * zoliq(i) 345 345 . * (1.0-EXP(-(zoliq(i)/zneb(i)/cl)**2)) *(1.-zfice(i)) 346 346 zrhol(i) = zrho(i) * zoliq(i) / zneb(i) 347 zfroi(i) = dtime/ FLOAT(ninter)/zdz(i)*zoliq(i)347 zfroi(i) = dtime/REAL(ninter)/zdz(i)*zoliq(i) 348 348 . *fallv(zrhol(i)) * zfice(i) 349 349 ztot(i) = zchau(i) + zfroi(i) … … 351 351 ztot(i) = MIN(MAX(ztot(i),0.0),zoliq(i)) 352 352 zoliq(i) = MAX(zoliq(i)-ztot(i), 0.0) 353 radliq(i,k) = radliq(i,k) + zoliq(i)/ FLOAT(ninter+1)353 radliq(i,k) = radliq(i,k) + zoliq(i)/REAL(ninter+1) 354 354 ENDIF 355 355 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.