Changeset 6817
- Timestamp:
- 2016-07-22T02:21:07+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/nemo_v3_6_STABLE_copy-NIWA-N96-ORCA1/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90
r6237 r6817 76 76 INTEGER :: ib_bdy, ii, ij, ik, igrd, ib, ir, iseg ! dummy loop indices 77 77 INTEGER :: icount, icountr, ibr_max, ilen1, ibm1 ! local integers 78 INTEGER :: iw e, ies, iso, ino, inum, id_dummy ! - -78 INTEGER :: iw, ie, is, in, inum, id_dummy ! - - 79 79 INTEGER :: igrd_start, igrd_end, jpbdta ! - - 80 80 INTEGER :: jpbdtau, jpbdtas ! - - … … 777 777 ! is = mjg(1) + 1 ! if monotasking and no zoom, is=2 778 778 ! in = mjg(1) + nlcj-1 - 1 ! if monotasking and no zoom, in=jpjm1 779 iw e= mig(1) - jpizoom + 2 ! if monotasking and no zoom, iw=2780 ie s= mig(1) + nlci - jpizoom - 1 ! if monotasking and no zoom, ie=jpim1781 is o= mjg(1) - jpjzoom + 2 ! if monotasking and no zoom, is=2782 in o= mjg(1) + nlcj - jpjzoom - 1 ! if monotasking and no zoom, in=jpjm1779 iw = mig(1) - jpizoom + 2 ! if monotasking and no zoom, iw=2 780 ie = mig(1) + nlci - jpizoom - 1 ! if monotasking and no zoom, ie=jpim1 781 is = mjg(1) - jpjzoom + 2 ! if monotasking and no zoom, is=2 782 in = mjg(1) + nlcj - jpjzoom - 1 ! if monotasking and no zoom, in=jpjm1 783 783 784 784 ALLOCATE( nbondi_bdy(nb_bdy)) … … 853 853 ENDIF 854 854 ! check if point is in local domain 855 IF( nbidta(ib,igrd,ib_bdy) >= iw e .AND. nbidta(ib,igrd,ib_bdy) <= ies.AND. &856 & nbjdta(ib,igrd,ib_bdy) >= is o .AND. nbjdta(ib,igrd,ib_bdy) <= ino) THEN855 IF( nbidta(ib,igrd,ib_bdy) >= iw .AND. nbidta(ib,igrd,ib_bdy) <= ie .AND. & 856 & nbjdta(ib,igrd,ib_bdy) >= is .AND. nbjdta(ib,igrd,ib_bdy) <= in ) THEN 857 857 ! 858 858 icount = icount + 1 … … 897 897 DO ib = 1, nblendta(igrd,ib_bdy) 898 898 ! check if point is in local domain and equals ir 899 IF( nbidta(ib,igrd,ib_bdy) >= iw e .AND. nbidta(ib,igrd,ib_bdy) <= ies.AND. &900 & nbjdta(ib,igrd,ib_bdy) >= is o .AND. nbjdta(ib,igrd,ib_bdy) <= ino.AND. &899 IF( nbidta(ib,igrd,ib_bdy) >= iw .AND. nbidta(ib,igrd,ib_bdy) <= ie .AND. & 900 & nbjdta(ib,igrd,ib_bdy) >= is .AND. nbjdta(ib,igrd,ib_bdy) <= in .AND. & 901 901 & nbrdta(ib,igrd,ib_bdy) == ir ) THEN 902 902 ! … … 1595 1595 ELSE 1596 1596 ! This is a corner 1597 IF(lwp)WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib)1597 WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib) 1598 1598 CALL bdy_ctl_corn(npckgw(ib), icornw(ib,1)) 1599 1599 itest=itest+1 … … 1609 1609 ELSE 1610 1610 ! This is a corner 1611 IF(lwp)WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib)1611 WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib) 1612 1612 CALL bdy_ctl_corn(npckgw(ib), icornw(ib,2)) 1613 1613 itest=itest+1 … … 1639 1639 ELSE 1640 1640 ! This is a corner 1641 IF(lwp)WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib)1641 WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib) 1642 1642 CALL bdy_ctl_corn(npckge(ib), icorne(ib,1)) 1643 1643 itest=itest+1 … … 1653 1653 ELSE 1654 1654 ! This is a corner 1655 IF(lwp)WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib)1655 WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib) 1656 1656 CALL bdy_ctl_corn(npckge(ib), icorne(ib,2)) 1657 1657 itest=itest+1
Note: See TracChangeset
for help on using the changeset viewer.