- Timestamp:
- 2017-05-30T10:23:55+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90
r7521 r8096 410 410 SELECT CASE ( cd_type ) 411 411 412 CASE( 'T', 'W' )412 CASE( 'T', 'W','U','V' ) 413 413 DO jk = 1, jpk 414 414 zsurf (:,:,jk) = p_e12(:,:) * p_e3(:,:,jk) * p_mask(:,:,jk) … … 1171 1171 ! check 1172 1172 !========================================================================== 1173 !CALL FLUSH(narea+1000-1) 1173 1174 !WRITE(narea+1000-1,*)"nfipproc(ji,jj),narea :",nfipproc(iproci,iprocj),narea 1174 1175 !WRITE(narea+1000-1,*)"proc i,j ",iproci,iprocj … … 1190 1191 !WRITE(narea+1000-1,*)"glamt gphit ",glamt(1,1),gphit(jpi,jpj),glamt(jpi,jpj),gphit(jpi,jpj) 1191 1192 !WRITE(narea+1000-1,*)"min max tmask ",MINVAL(tmask),MAXVAL(tmask) 1192 1193 !CALL FLUSH(narea+1000-1) 1193 1194 !========================================================================== 1194 1195 ! coarsened domain: dimensions along I … … 1233 1234 mis2_crs(ji)=ijis 1234 1235 mie2_crs(ji)=ijie 1236 !WRITE(narea+1000-1,*)"glo crs",ji,ijis,ijie,ijis-nimpp+1,ijie-nimpp+1 1235 1237 ENDDO 1236 1238 … … 1254 1256 ijis=ji 1255 1257 ii_start = mis2_crs(ijis)-nimpp+1 1258 WRITE(narea+1000-1,*)"start ",ijis,mis2_crs(ijis),ii_start ; CALL FLUSH(narea+1000-1) 1256 1259 1257 1260 !---------------------------------------------------------------------------------------------- … … 1292 1295 1293 1296 ji=jpiglo_crs 1294 DO WHILE( mie2_crs(ji) - nimpp + 1 .GT. jpi )1297 DO WHILE( mie2_crs(ji) - nimpp + 1 .GT. nlci ) 1295 1298 ji=ji-1 1296 1299 IF( ji==1 )EXIT 1297 1300 END DO 1298 1301 ijie=ji 1302 !WRITE(narea+1000-1,*)"end ",ijie ; CALL FLUSH(narea+1000-1) 1299 1303 1300 1304 !------------------------------------------------------------------------------- … … 1302 1306 !------------------------------------------------------------------------------- 1303 1307 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 1306 1313 1307 1314 !------------------------------------------------------------------------------- … … 1310 1317 DO ji = 1, jpi_crs 1311 1318 mig_crs(ji) = ji + nimpp_crs - 1 1319 WRITE(narea+1000-1,*)"crs loctoglo",ji,mig_crs(ji) ; CALL FLUSH(narea+1000-1) 1312 1320 ENDDO 1313 1321 DO ji = 1, jpiglo_crs … … 1320 1328 !--------------------------------------------------------------------------------------- 1321 1329 DO ji = 1, nlei_crs 1330 IF( mig_crs(ji) .GT. jpiglo_crs )WRITE(narea+1000-1,*)"BUG1 " ; CALL FLUSH(narea+1000-1) 1322 1331 mis_crs(ji) = mis2_crs(mig_crs(ji)) - nimpp + 1 1323 1332 mie_crs(ji) = mie2_crs(mig_crs(ji)) - nimpp + 1 1324 IF( iproci == jpni .AND. ji == nlei_crs )THEN1325 mie_crs(ji) = nlei1326 mie2_crs(mig_crs(ji)) = mig(nlei)1327 ENDIF1333 !IF( iproci == jpni .AND. ji == nlei_crs )THEN 1334 ! mie_crs(ji) = nlei 1335 ! mie2_crs(mig_crs(ji)) = mig(nlei) 1336 !ENDIF 1328 1337 nfactx(ji) = mie_crs(ji)-mis_crs(ji)+1 1329 1338 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 1330 1369 1331 1370 !========================================================================== … … 1454 1493 ENDIF 1455 1494 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) 1456 1499 !------------------------------------------------------------------------------- 1457 1500 ! J-7 local to global and global to local indices for CRS grid … … 1521 1564 ! check 1522 1565 !========================================================================== 1566 !WRITE(narea+1000-1,*)"mpp_crs ",nimpp_crs,njmpp_crs 1523 1567 !WRITE(narea+1000-1,*)"loc crs jpi nldi,nlei,nlci ",jpi_crs, nldi_crs ,nlei_crs ,nlci_crs 1524 1568 !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.