Changeset 40
- Timestamp:
- Feb 9, 2000, 11:26:34 AM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/phylmd/condsurf.F
r2 r40 43 43 REAL lmt_rug(klon) ! longeur de rugosite du sol 44 44 REAL lmt_bils(klon) 45 c46 c Couplage OASIS:47 #include "oasis.h"48 REAL cpl_sst(iim,jjm+1), cpl_sic(iim,jjm+1)49 REAL cpl_alb_sst(iim,jjm+1), cpl_alb_sic(iim,jjm+1)50 45 c 51 46 c Variables locales: … … 299 294 PRINT*, 'SST, ALB, RUG, etc. sont lus pour jour: ', jour 300 295 c 301 IF (ok_oasis) THEN302 C303 CALL fromcpl(jour,(jjm+1)*iim,304 . cpl_sst,cpl_sic,cpl_alb_sst,cpl_alb_sic)305 DO i = 1, iim-1 ! un seul point pour le pole nord306 cpl_sst(i,1) = cpl_sst(iim,1)307 cpl_sic(i,1) = cpl_sic(iim,1)308 cpl_alb_sst(i,1) = cpl_alb_sst(iim,1)309 cpl_alb_sic(i,1) = cpl_alb_sic(iim,1)310 ENDDO311 DO i = 2, iim ! un seul point pour le pole sud312 cpl_sst(i,jjm+1) = cpl_sst(1,jjm+1)313 cpl_sic(i,jjm+1) = cpl_sic(1,jjm+1)314 cpl_alb_sst(i,jjm+1) = cpl_alb_sst(1,jjm+1)315 cpl_alb_sic(i,jjm+1) = cpl_alb_sic(1,jjm+1)316 ENDDO317 c318 ig = 1319 IF (pctsrf(ig,is_oce).GT.epsfra .OR.320 . pctsrf(ig,is_sic).GT.epsfra) THEN321 pctsrf(ig,is_oce) = pctsrf(ig,is_oce)322 . - (cpl_sic(1,1)-pctsrf(ig,is_sic))323 pctsrf(ig,is_sic) = cpl_sic(1,1)324 lmt_sst(ig) = cpl_sst(1,1)325 ENDIF326 DO j = 2, jjm327 DO i = 1, iim328 ig = ig + 1329 IF (pctsrf(ig,is_oce).GT.epsfra .OR.330 . pctsrf(ig,is_sic).GT.epsfra) THEN331 pctsrf(ig,is_oce) = pctsrf(ig,is_oce)332 . - (cpl_sic(i,j)-pctsrf(ig,is_sic))333 pctsrf(ig,is_sic) = cpl_sic(i,j)334 lmt_sst(ig) = cpl_sst(i,j)335 ENDIF336 ENDDO337 ENDDO338 ig = ig + 1339 IF (pctsrf(ig,is_oce).GT.epsfra .OR.340 . pctsrf(ig,is_sic).GT.epsfra) THEN341 pctsrf(ig,is_oce) = pctsrf(ig,is_oce)342 . - (cpl_sic(1,jjm+1)-pctsrf(ig,is_sic))343 pctsrf(ig,is_sic) = cpl_sic(1,jjm+1)344 lmt_sst(ig) = cpl_sst(1,jjm+1)345 ENDIF346 c347 ENDIF ! ok_oasis348 c349 296 RETURN 350 297 END
Note: See TracChangeset
for help on using the changeset viewer.