- Timestamp:
- 2019-05-16T15:23:56+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/BDY/bdyini.F90
r10968 r10986 91 91 ! ----------------------------------------- 92 92 IF ( ln_bdy ) THEN 93 IF(lwp) WRITE(numout,*) 94 IF(lwp) WRITE(numout,*) 'bdy_init : initialization of open boundaries' 95 IF(lwp) WRITE(numout,*) '~~~~~~~~' 93 IF(lwp) THEN 94 WRITE(numout,*) 95 WRITE(numout,*) 'bdy_init : initialization of open boundaries' 96 WRITE(numout,*) '~~~~~~~~' 97 IF(lflush) CALL FLUSH(numout) 98 ENDIF 96 99 ! 97 100 ! Open boundaries definition (arrays and masks) … … 105 108 ! 106 109 ELSE 107 IF(lwp) WRITE(numout,*) 108 IF(lwp) WRITE(numout,*) 'bdy_init : open boundaries not used (ln_bdy = F)' 109 IF(lwp) WRITE(numout,*) '~~~~~~~~' 110 IF(lwp) THEN 111 WRITE(numout,*) 112 WRITE(numout,*) 'bdy_init : open boundaries not used (ln_bdy = F)' 113 WRITE(numout,*) '~~~~~~~~' 114 IF(lflush) CALL FLUSH(numout) 115 ENDIF 110 116 ! 111 117 ENDIF … … 403 409 ENDIF 404 410 ENDIF 411 412 IF(lflush .AND. lwp) CALL FLUSH(numout) 405 413 406 414 ! ------------------------------------------------- … … 568 576 ! 569 577 ibr_max = MAXVAL( nbrdta(:,igrd,ib_bdy) ) 570 IF(lwp) WRITE(numout,*) 571 IF(lwp) WRITE(numout,*) ' Maximum rimwidth in file is ', ibr_max 572 IF(lwp) WRITE(numout,*) ' nn_rimwidth from namelist is ', nn_rimwidth(ib_bdy) 578 IF(lwp) THEN 579 WRITE(numout,*) 580 WRITE(numout,*) ' Maximum rimwidth in file is ', ibr_max 581 WRITE(numout,*) ' nn_rimwidth from namelist is ', nn_rimwidth(ib_bdy) 582 IF(lflush) CALL FLUSH(numout) 583 ENDIF 573 584 IF (ibr_max < nn_rimwidth(ib_bdy)) & 574 585 CALL ctl_stop( 'nn_rimwidth is larger than maximum rimwidth in file',cn_coords_file(ib_bdy) ) … … 1205 1216 IF( i_offset == 1 .and. zefl + zwfl == 2 ) THEN 1206 1217 icount = icount + 1 1207 IF(lwp) WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 1218 IF(lwp) THEN 1219 WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 1220 IF(lflush) CALL FLUSH(numout) 1221 ENDIF 1208 1222 ELSE 1209 1223 idx_bdy(ib_bdy)%flagu(ib,igrd) = -zefl + zwfl … … 1237 1251 ! This error check only works if you are using the bdyXmask arrays 1238 1252 IF( j_offset == 1 .and. znfl + zsfl == 2 ) THEN 1239 IF(lwp) WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 1253 IF(lwp) THEN 1254 WRITE(numout,*) 'Problem with igrd = ',igrd,' at (global) nbi, nbj : ',mig(nbi),mjg(nbj) 1255 IF(lflush) CALL FLUSH(numout) 1256 ENDIF 1240 1257 icount = icount + 1 1241 1258 ELSE … … 1273 1290 !!---------------------------------------------------------------------- 1274 1291 ! 1275 IF (lwp) WRITE(numout,*) ' ' 1276 IF (lwp) WRITE(numout,*) 'bdy_ctl_seg: Check analytical segments' 1277 IF (lwp) WRITE(numout,*) '~~~~~~~~~~~~' 1278 ! 1279 IF(lwp) WRITE(numout,*) 'Number of east segments : ', nbdysege 1280 IF(lwp) WRITE(numout,*) 'Number of west segments : ', nbdysegw 1281 IF(lwp) WRITE(numout,*) 'Number of north segments : ', nbdysegn 1282 IF(lwp) WRITE(numout,*) 'Number of south segments : ', nbdysegs 1292 IF (lwp) THEN 1293 WRITE(numout,*) ' ' 1294 WRITE(numout,*) 'bdy_ctl_seg: Check analytical segments' 1295 WRITE(numout,*) '~~~~~~~~~~~~' 1296 ! 1297 WRITE(numout,*) 'Number of east segments : ', nbdysege 1298 WRITE(numout,*) 'Number of west segments : ', nbdysegw 1299 WRITE(numout,*) 'Number of north segments : ', nbdysegn 1300 WRITE(numout,*) 'Number of south segments : ', nbdysegs 1301 IF(lflush) CALL FLUSH(numout) 1302 ENDIF 1283 1303 ! 1. Check bounds 1284 1304 !---------------- … … 1322 1342 ! 2. Look for segment crossings 1323 1343 !------------------------------ 1324 IF (lwp) WRITE(numout,*) '**Look for segments corners :' 1344 IF (lwp) THEN 1345 WRITE(numout,*) '**Look for segments corners :' 1346 IF(lflush) CALL FLUSH(numout) 1347 ENDIF 1325 1348 ! 1326 1349 itest = 0 ! corner number … … 1482 1505 ELSE 1483 1506 ! This is a corner 1484 IF(lwp) WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib) 1507 IF(lwp) THEN 1508 WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib) 1509 IF(lflush) CALL FLUSH(numout) 1510 ENDIF 1485 1511 CALL bdy_ctl_corn(npckgw(ib), icornw(ib,1)) 1486 1512 itest=itest+1 … … 1494 1520 ELSE 1495 1521 ! This is a corner 1496 IF(lwp) WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib) 1522 IF(lwp) THEN 1523 WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib) 1524 IF(lflush) CALL FLUSH(numout) 1525 ENDIF 1497 1526 CALL bdy_ctl_corn(npckgw(ib), icornw(ib,2)) 1498 1527 itest=itest+1 … … 1522 1551 ELSE 1523 1552 ! This is a corner 1524 IF(lwp) WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib) 1553 IF(lwp) THEN 1554 WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib) 1555 IF(lflush) CALL FLUSH(numout) 1556 ENDIF 1525 1557 CALL bdy_ctl_corn(npckge(ib), icorne(ib,1)) 1526 1558 itest=itest+1 … … 1534 1566 ELSE 1535 1567 ! This is a corner 1536 IF(lwp) WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib) 1568 IF(lwp) THEN 1569 WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib) 1570 IF(lflush) CALL FLUSH(numout) 1571 ENDIF 1537 1572 CALL bdy_ctl_corn(npckge(ib), icorne(ib,2)) 1538 1573 itest=itest+1 … … 1593 1628 END DO 1594 1629 ! 1595 IF ((itest==0).AND.(lwp)) WRITE(numout,*) 'NO open boundary corner found' 1630 IF ((itest==0).AND.(lwp)) THEN 1631 WRITE(numout,*) 'NO open boundary corner found' 1632 IF(lflush) CALL FLUSH(numout) 1633 ENDIF 1596 1634 ! 1597 1635 ! Other tests TBD:
Note: See TracChangeset
for help on using the changeset viewer.
