Changeset 40 for LMDZ.3.3


Ignore:
Timestamp:
Feb 9, 2000, 11:26:34 AM (25 years ago)
Author:
lmdz
Message:

L'appel au coupleur est enleve (pour pouvoir changer la frequence de couplage). L.Li
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/phylmd/condsurf.F

    r2 r40  
    4343      REAL lmt_rug(klon) ! longeur de rugosite du sol
    4444      REAL lmt_bils(klon)
    45 c
    46 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)
    5045c
    5146c Variables locales:
     
    299294      PRINT*, 'SST, ALB, RUG, etc. sont lus pour jour: ', jour
    300295c
    301       IF (ok_oasis) THEN
    302 C
    303          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 nord
    306             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          ENDDO
    311          DO i = 2, iim ! un seul point pour le pole sud
    312             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          ENDDO
    317 c
    318          ig = 1
    319          IF (pctsrf(ig,is_oce).GT.epsfra .OR.
    320      .       pctsrf(ig,is_sic).GT.epsfra) THEN
    321             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          ENDIF
    326          DO j = 2, jjm
    327          DO i = 1, iim
    328          ig = ig + 1
    329          IF (pctsrf(ig,is_oce).GT.epsfra .OR.
    330      .       pctsrf(ig,is_sic).GT.epsfra) THEN
    331            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          ENDIF
    336          ENDDO
    337          ENDDO
    338          ig = ig + 1
    339          IF (pctsrf(ig,is_oce).GT.epsfra .OR.
    340      .       pctsrf(ig,is_sic).GT.epsfra) THEN
    341             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          ENDIF
    346 c
    347       ENDIF ! ok_oasis
    348 c
    349296      RETURN
    350297      END
Note: See TracChangeset for help on using the changeset viewer.