Changeset 8096


Ignore:
Timestamp:
2017-05-30T10:23:55+02:00 (3 years ago)
Author:
cbricaud
Message:

correction for mpp indices computation in CRS space

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90

    r7521 r8096  
    410410            SELECT CASE ( cd_type ) 
    411411             
    412                CASE( 'T', 'W' ) 
     412               CASE( 'T', 'W','U','V' ) 
    413413                  DO jk = 1, jpk 
    414414                     zsurf   (:,:,jk) =  p_e12(:,:) * p_e3(:,:,jk) *  p_mask(:,:,jk)  
     
    11711171         ! check 
    11721172         !========================================================================== 
     1173         !CALL FLUSH(narea+1000-1) 
    11731174         !WRITE(narea+1000-1,*)"nfipproc(ji,jj),narea :",nfipproc(iproci,iprocj),narea 
    11741175         !WRITE(narea+1000-1,*)"proc i,j ",iproci,iprocj 
     
    11901191         !WRITE(narea+1000-1,*)"glamt gphit ",glamt(1,1),gphit(jpi,jpj),glamt(jpi,jpj),gphit(jpi,jpj) 
    11911192         !WRITE(narea+1000-1,*)"min max tmask ",MINVAL(tmask),MAXVAL(tmask) 
    1192  
     1193         !CALL FLUSH(narea+1000-1) 
    11931194         !========================================================================== 
    11941195         ! coarsened domain: dimensions along I 
     
    12331234               mis2_crs(ji)=ijis 
    12341235               mie2_crs(ji)=ijie 
     1236               !WRITE(narea+1000-1,*)"glo crs",ji,ijis,ijie,ijis-nimpp+1,ijie-nimpp+1 
    12351237            ENDDO 
    12361238 
     
    12541256         ijis=ji 
    12551257         ii_start = mis2_crs(ijis)-nimpp+1 
     1258         WRITE(narea+1000-1,*)"start ",ijis,mis2_crs(ijis),ii_start ; CALL FLUSH(narea+1000-1) 
    12561259 
    12571260         !---------------------------------------------------------------------------------------------- 
     
    12921295 
    12931296         ji=jpiglo_crs 
    1294          DO WHILE( mie2_crs(ji) - nimpp + 1 .GT. jpi ) 
     1297         DO WHILE( mie2_crs(ji) - nimpp + 1 .GT. nlci ) 
    12951298            ji=ji-1 
    12961299            IF( ji==1 )EXIT 
    12971300         END DO 
    12981301         ijie=ji 
     1302         !WRITE(narea+1000-1,*)"end ",ijie ; CALL FLUSH(narea+1000-1) 
    12991303 
    13001304         !------------------------------------------------------------------------------- 
     
    13021306         !------------------------------------------------------------------------------- 
    13031307         nlei_crs=ijie-nimpp_crs+1 
    1304          nlci_crs=nlei_crs+jpreci 
    1305          IF( iproci == jpni )nlei_crs=nlci_crs 
     1308         !nlci_crs=nlei_crs ! cbr ???? +jpreci 
     1309         !IF( iproci == jpni ) THEN ; nlci_crs=nlei_crs ! cbr ???? +jpreci 
     1310         !ELSE                      ; nlci_crs=nlei_crs+1 
     1311         !ENDIF 
     1312         !cbr???? IF( iproci == jpni )nlei_crs=nlci_crs 
    13061313 
    13071314         !------------------------------------------------------------------------------- 
     
    13101317         DO ji = 1, jpi_crs 
    13111318            mig_crs(ji) = ji + nimpp_crs - 1 
     1319            WRITE(narea+1000-1,*)"crs loctoglo",ji,mig_crs(ji) ; CALL FLUSH(narea+1000-1) 
    13121320         ENDDO 
    13131321         DO ji = 1, jpiglo_crs 
     
    13201328         !--------------------------------------------------------------------------------------- 
    13211329         DO ji = 1, nlei_crs 
     1330            IF( mig_crs(ji) .GT. jpiglo_crs )WRITE(narea+1000-1,*)"BUG1 " ; CALL FLUSH(narea+1000-1) 
    13221331            mis_crs(ji) = mis2_crs(mig_crs(ji)) - nimpp + 1 
    13231332            mie_crs(ji) = mie2_crs(mig_crs(ji)) - nimpp + 1 
    1324             IF( iproci == jpni  .AND. ji == nlei_crs )THEN 
    1325                mie_crs(ji) = nlei 
    1326                mie2_crs(mig_crs(ji)) = mig(nlei) 
    1327             ENDIF 
     1333            !IF( iproci == jpni  .AND. ji == nlei_crs )THEN 
     1334            !   mie_crs(ji) = nlei 
     1335            !   mie2_crs(mig_crs(ji)) = mig(nlei) 
     1336            !ENDIF 
    13281337            nfactx(ji)  = mie_crs(ji)-mis_crs(ji)+1 
    13291338         ENDDO 
     1339 
     1340         !--------- 
     1341         !cbr  
     1342         IF( iproci == 1 ) THEN 
     1343            nldi_crs=1 
     1344            mis_crs(1) = 1 
     1345            mie_crs(1) = 1 
     1346            mis2_crs(1) = 1 
     1347            mie2_crs(1) = 1 
     1348         ENDIF 
     1349 
     1350         IF( iproci == jpni ) THEN 
     1351            nlei_crs=jpiglo_crs-nimpp_crs+1 
     1352            nlci_crs=nlei_crs 
     1353            mis_crs(nlei_crs) = 1 
     1354            mie_crs(nlei_crs) = 1 
     1355            mis2_crs(nlei_crs) = 1 
     1356            mie2_crs(nlei_crs) = 1 
     1357            nfactx(nlei_crs)=0 
     1358         ELSE 
     1359            nlci_crs=nlei_crs+1 
     1360         ENDIF 
     1361 
     1362         !WRITE(narea+1000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs            ,nlei_crs             ,nlci_crs 
     1363         !CALL FLUSH(narea+1000-1) 
     1364         !WRITE(narea+1000-1,*)"glo crs jpi nldi,nlei      ",jpi_crs, nldi_crs+nimpp_crs-1,nlei_crs+nimpp_crs-1 
     1365         !CALL FLUSH(narea+1000-1) 
     1366         !DO ji = 1, jpi_crs 
     1367         !   WRITE(narea+1000-1,*)"crs i",ji,ji+nimpp_crs-1,mis_crs(ji),mie_crs(ji),mis_crs(ji)+nimpp-1,mie_crs(ji)+nimpp-1,nfactx(ji) 
     1368         !ENDDO 
    13301369 
    13311370         !========================================================================== 
     
    14541493         ENDIF 
    14551494 
     1495         !WRITE(narea+1000-1,*)"loc crs jpj nldj,nlej,nlcj ",jpj_crs, nldj_crs            ,nlej_crs             ,nlcj_crs 
     1496         !CALL FLUSH(narea+1000-1) 
     1497         !WRITE(narea+1000-1,*)"glo crs jpj nldj,nlej      ",jpj_crs, nldj_crs+njmpp_crs-1,nlej_crs+njmpp_crs-1 
     1498         !CALL FLUSH(narea+1000-1) 
    14561499         !------------------------------------------------------------------------------- 
    14571500         ! J-7 local to global and global to local indices for CRS grid 
     
    15211564         ! check 
    15221565         !========================================================================== 
     1566         !WRITE(narea+1000-1,*)"mpp_crs ",nimpp_crs,njmpp_crs 
    15231567         !WRITE(narea+1000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs            ,nlei_crs             ,nlci_crs 
    15241568         !WRITE(narea+1000-1,*)"glo crs jpi nldi,nlei      ",jpi_crs, nldi_crs+nimpp_crs-1,nlei_crs+nimpp_crs-1 
Note: See TracChangeset for help on using the changeset viewer.