Ignore:
Timestamp:
Jun 7, 2007, 12:01:52 PM (17 years ago)
Author:
Laurent Fairhead
Message:

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le

LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/iniphysiq.F

    r766 r776  
    99     $           plat,plon,parea,pcu,pcv,
    1010     $           prad,pg,pr,pcpp)
    11       USE dimphy,klon2=>klon,klon=>klon2
     11      USE dimphy
     12      USE mod_grid_phy_lmdz
     13      USE mod_phys_lmdz_para
    1214      USE comgeomphy
    1315      IMPLICIT NONE
     
    5254 
    5355      INTEGER ngrid,nlayer
    54       REAL plat(ngrid),plon(ngrid),parea(klon),pcu(klon),pcv(klon)
     56      REAL plat(ngrid),plon(ngrid),parea(klon_glo)
     57      REAL pcu(klon_glo),pcv(klon_glo)
    5558      INTEGER pdayref
    5659      INTEGER :: ibegin,iend
     
    6669      ENDIF
    6770
    68       IF (ngrid.NE.klon) THEN
     71      IF (ngrid.NE.klon_glo) THEN
    6972         PRINT*,'STOP in inifis'
    7073         PRINT*,'Probleme de dimensions :'
    7174         PRINT*,'ngrid     = ',ngrid
    72          PRINT*,'klon   = ',klon
     75         PRINT*,'klon   = ',klon_glo
    7376         STOP
    7477      ENDIF
    7578c$OMP PARALLEL PRIVATE(ibegin,iend)
    7679c$OMP+         SHARED(parea,pcu,pcv,plon,plat)
    77       ibegin=klon_begin+klon_omp_begin(omp_rank)-1
    78       iend=ibegin+klon_omp-1
    7980
    80       airephy(1:klon_omp)=parea(ibegin:iend)
    81       cuphy(1:klon_omp)=pcu(ibegin:iend)
    82       cvphy(1:klon_omp)=pcv(ibegin:iend)
    83       rlond(1:klon_omp) = plon(ibegin:iend)
    84       rlatd(1:klon_omp) = plat(ibegin:iend)
     81      airephy(1:klon_omp)=parea(klon_omp_begin:klon_omp_end)
     82      cuphy(1:klon_omp)=pcu(klon_omp_begin:klon_omp_end)
     83      cvphy(1:klon_omp)=pcv(klon_omp_begin:klon_omp_end)
     84      rlond(1:klon_omp) = plon(klon_omp_begin:klon_omp_end)
     85      rlatd(1:klon_omp) = plat(klon_omp_begin:klon_omp_end)
    8586
    8687      call suphec
Note: See TracChangeset for help on using the changeset viewer.