- Timestamp:
- 2015-11-20T09:39:06+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90
r5626 r5901 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 , ie, is, in, inum, id_dummy! - -78 INTEGER :: iwe, ies, iso, ino, inum, id_dummy ! - - 79 79 INTEGER :: igrd_start, igrd_end, jpbdta ! - - 80 80 INTEGER :: jpbdtau, jpbdtas ! - - … … 801 801 ! is = mjg(1) + 1 ! if monotasking and no zoom, is=2 802 802 ! in = mjg(1) + nlcj-1 - 1 ! if monotasking and no zoom, in=jpjm1 803 iw = mig(1) - jpizoom + 2 ! if monotasking and no zoom, iw=2804 ie = mig(1) + nlci - jpizoom - 1 ! if monotasking and no zoom, ie=jpim1805 is = mjg(1) - jpjzoom + 2 ! if monotasking and no zoom, is=2806 in = mjg(1) + nlcj - jpjzoom - 1 ! if monotasking and no zoom, in=jpjm1803 iwe = mig(1) - jpizoom + 2 ! if monotasking and no zoom, iw=2 804 ies = mig(1) + nlci - jpizoom - 1 ! if monotasking and no zoom, ie=jpim1 805 iso = mjg(1) - jpjzoom + 2 ! if monotasking and no zoom, is=2 806 ino = mjg(1) + nlcj - jpjzoom - 1 ! if monotasking and no zoom, in=jpjm1 807 807 808 808 ALLOCATE( nbondi_bdy(nb_bdy)) … … 877 877 ENDIF 878 878 ! check if point is in local domain 879 IF( nbidta(ib,igrd,ib_bdy) >= iw .AND. nbidta(ib,igrd,ib_bdy) <= ie.AND. &880 & nbjdta(ib,igrd,ib_bdy) >= is .AND. nbjdta(ib,igrd,ib_bdy) <= in) THEN879 IF( nbidta(ib,igrd,ib_bdy) >= iwe .AND. nbidta(ib,igrd,ib_bdy) <= ies .AND. & 880 & nbjdta(ib,igrd,ib_bdy) >= iso .AND. nbjdta(ib,igrd,ib_bdy) <= ino ) THEN 881 881 ! 882 882 icount = icount + 1 … … 914 914 com_south_b = 0 915 915 com_north_b = 0 916 916 917 DO igrd = 1, jpbgrd 917 918 icount = 0 … … 920 921 DO ib = 1, nblendta(igrd,ib_bdy) 921 922 ! check if point is in local domain and equals ir 922 IF( nbidta(ib,igrd,ib_bdy) >= iw .AND. nbidta(ib,igrd,ib_bdy) <= ie.AND. &923 & nbjdta(ib,igrd,ib_bdy) >= is .AND. nbjdta(ib,igrd,ib_bdy) <= in.AND. &923 IF( nbidta(ib,igrd,ib_bdy) >= iwe .AND. nbidta(ib,igrd,ib_bdy) <= ies .AND. & 924 & nbjdta(ib,igrd,ib_bdy) >= iso .AND. nbjdta(ib,igrd,ib_bdy) <= ino .AND. & 924 925 & nbrdta(ib,igrd,ib_bdy) == ir ) THEN 925 926 ! … … 1618 1619 ELSE 1619 1620 ! This is a corner 1620 WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib)1621 IF(lwp) WRITE(numout,*) 'Found a South-West corner at (i,j): ', jpiwob(ib), jpjwdt(ib) 1621 1622 CALL bdy_ctl_corn(npckgw(ib), icornw(ib,1)) 1622 1623 itest=itest+1 … … 1632 1633 ELSE 1633 1634 ! This is a corner 1634 WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib)1635 IF(lwp) WRITE(numout,*) 'Found a North-West corner at (i,j): ', jpiwob(ib), jpjwft(ib) 1635 1636 CALL bdy_ctl_corn(npckgw(ib), icornw(ib,2)) 1636 1637 itest=itest+1 … … 1662 1663 ELSE 1663 1664 ! This is a corner 1664 WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib)1665 IF(lwp) WRITE(numout,*) 'Found a South-East corner at (i,j): ', jpieob(ib)+1, jpjedt(ib) 1665 1666 CALL bdy_ctl_corn(npckge(ib), icorne(ib,1)) 1666 1667 itest=itest+1 … … 1676 1677 ELSE 1677 1678 ! This is a corner 1678 WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib)1679 IF(lwp) WRITE(numout,*) 'Found a North-East corner at (i,j): ', jpieob(ib)+1, jpjeft(ib) 1679 1680 CALL bdy_ctl_corn(npckge(ib), icorne(ib,2)) 1680 1681 itest=itest+1
Note: See TracChangeset
for help on using the changeset viewer.