Symbol Description Uses
jp1 190


File Line number Procedure Code
./pkg/kpp/kpp_calc.F 202 KPP_CALC 
      integer i, j, k, kp1, km1, im1, ip1, jm1, jp1
./pkg/kpp/kpp_calc.F 427 KPP_CALC 
        jp1 = j + 1
./pkg/kpp/kpp_calc.F 438 KPP_CALC 
     &        (vVel(i,  jp1,k,bi,bj)-vVel(i,  jp1,kp1,bi,bj)) *
./pkg/kpp/kpp_calc.F 439 KPP_CALC 
     &        (vVel(i,  jp1,k,bi,bj)-vVel(i,  jp1,kp1,bi,bj)) )
./pkg/kpp/kpp_calc.F 446 KPP_CALC 
     &        (uVel(i,  jp1,k,bi,bj)-uVel(i,  jp1,kp1,bi,bj)) *
./pkg/kpp/kpp_calc.F 447 KPP_CALC 
     &        (uVel(i,  jp1,k,bi,bj)-uVel(i,  jp1,kp1,bi,bj)) +
./pkg/kpp/kpp_calc.F 448 KPP_CALC 
     &        (uVel(ip1,jp1,k,bi,bj)-uVel(ip1,jp1,kp1,bi,bj)) *
./pkg/kpp/kpp_calc.F 449 KPP_CALC 
     &        (uVel(ip1,jp1,k,bi,bj)-uVel(ip1,jp1,kp1,bi,bj)) +
./pkg/kpp/kpp_calc.F 452 KPP_CALC 
     &        (vVel(im1,jp1,k,bi,bj)-vVel(im1,jp1,kp1,bi,bj)) *
./pkg/kpp/kpp_calc.F 453 KPP_CALC 
     &        (vVel(im1,jp1,k,bi,bj)-vVel(im1,jp1,kp1,bi,bj)) +
./pkg/kpp/kpp_calc.F 456 KPP_CALC 
     &        (vVel(ip1,jp1,k,bi,bj)-vVel(ip1,jp1,kp1,bi,bj)) *
./pkg/kpp/kpp_calc.F 457 KPP_CALC 
     &        (vVel(ip1,jp1,k,bi,bj)-vVel(ip1,jp1,kp1,bi,bj)) )
./pkg/kpp/kpp_forcing_surf.F 132 KPP_FORCING_SURF 
      integer i, j, k, im1, ip1, jm1, jp1
./pkg/kpp/kpp_forcing_surf.F 171 KPP_FORCING_SURF 
       jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 177 KPP_FORCING_SURF 
     &   (surfForcV(i,j,bi,bj) + surfForcV(i,jp1,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 178 KPP_FORCING_SURF 
     &   (surfForcV(i,j,bi,bj) + surfForcV(i,jp1,bi,bj))
./pkg/kpp/kpp_forcing_surf.F 185 KPP_FORCING_SURF 
       jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 200 KPP_FORCING_SURF 
       jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 217 KPP_FORCING_SURF 
            jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 270 KPP_FORCING_SURF 
       jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 308 KPP_FORCING_SURF 
     &       (vVel(i,  jp1,1,bi,bj)-vVel(i,  jp1,2,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 309 KPP_FORCING_SURF 
     &       (vVel(i,  jp1,1,bi,bj)-vVel(i,  jp1,2,bi,bj)) )
./pkg/kpp/kpp_forcing_surf.F 328 KPP_FORCING_SURF 
        vRef(i,j) = p5 * ( vVel(i,j,1,bi,bj) + vVel(i,jp1,1,bi,bj) )
./pkg/kpp/kpp_forcing_surf.F 333 KPP_FORCING_SURF 
     &        surfForcV(i,jp1,bi,bj) ) * p5 *recip_drF(1)
./pkg/kpp/kpp_forcing_surf.F 354 KPP_FORCING_SURF 
        jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 365 KPP_FORCING_SURF 
     $        (vRef(i,j) - vVel(i,  jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 366 KPP_FORCING_SURF 
     $        (vRef(i,j) - vVel(i,  jp1,k,bi,bj)) )
./pkg/kpp/kpp_forcing_surf.F 373 KPP_FORCING_SURF 
     $        (uRef(i,j) - uVel(i,  jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 374 KPP_FORCING_SURF 
     $        (uRef(i,j) - uVel(i,  jp1,k,bi,bj)) +
./pkg/kpp/kpp_forcing_surf.F 375 KPP_FORCING_SURF 
     $        (uRef(i,j) - uVel(ip1,jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 376 KPP_FORCING_SURF 
     $        (uRef(i,j) - uVel(ip1,jp1,k,bi,bj)) +
./pkg/kpp/kpp_forcing_surf.F 379 KPP_FORCING_SURF 
     $        (vRef(i,j) - vVel(im1,jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 380 KPP_FORCING_SURF 
     $        (vRef(i,j) - vVel(im1,jp1,k,bi,bj)) +
./pkg/kpp/kpp_forcing_surf.F 383 KPP_FORCING_SURF 
     $        (vRef(i,j) - vVel(ip1,jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 384 KPP_FORCING_SURF 
     $        (vRef(i,j) - vVel(ip1,jp1,k,bi,bj)) )
./pkg/kpp/kpp_forcing_surf.F 395 KPP_FORCING_SURF 
        jp1 = j + 1
./pkg/kpp/kpp_forcing_surf.F 406 KPP_FORCING_SURF 
     $        (vVel(i,  jp1,1,bi,bj)-vVel(i,  jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 407 KPP_FORCING_SURF 
     $        (vVel(i,  jp1,1,bi,bj)-vVel(i,  jp1,k,bi,bj)) )
./pkg/kpp/kpp_forcing_surf.F 414 KPP_FORCING_SURF 
     $        (uVel(i,  jp1,1,bi,bj)-uVel(i,  jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 415 KPP_FORCING_SURF 
     $        (uVel(i,  jp1,1,bi,bj)-uVel(i,  jp1,k,bi,bj)) +
./pkg/kpp/kpp_forcing_surf.F 416 KPP_FORCING_SURF 
     $        (uVel(ip1,jp1,1,bi,bj)-uVel(ip1,jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 417 KPP_FORCING_SURF 
     $        (uVel(ip1,jp1,1,bi,bj)-uVel(ip1,jp1,k,bi,bj)) +
./pkg/kpp/kpp_forcing_surf.F 420 KPP_FORCING_SURF 
     $        (vVel(im1,jp1,1,bi,bj)-vVel(im1,jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 421 KPP_FORCING_SURF 
     $        (vVel(im1,jp1,1,bi,bj)-vVel(im1,jp1,k,bi,bj)) +
./pkg/kpp/kpp_forcing_surf.F 424 KPP_FORCING_SURF 
     $        (vVel(ip1,jp1,1,bi,bj)-vVel(ip1,jp1,k,bi,bj)) *
./pkg/kpp/kpp_forcing_surf.F 425 KPP_FORCING_SURF 
     $        (vVel(ip1,jp1,1,bi,bj)-vVel(ip1,jp1,k,bi,bj)) )
./pkg/kpp/kpp_routines.F 1094 SMOOTH_HORIZ 
      integer i, j, im1, ip1, jm1, jp1
./pkg/kpp/kpp_routines.F 1106 SMOOTH_HORIZ 
         jp1 = j+1
./pkg/kpp/kpp_routines.F 1115 SMOOTH_HORIZ 
     &                     maskC(i  ,jp1,k,bi,bj) ) +
./pkg/kpp/kpp_routines.F 1117 SMOOTH_HORIZ 
     &                     maskC(im1,jp1,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1119 SMOOTH_HORIZ 
     &                     maskC(ip1,jp1,k,bi,bj) )
./pkg/kpp/kpp_routines.F 1126 SMOOTH_HORIZ 
     &                     fld(i  ,jp1)*maskC(i  ,jp1,k,bi,bj))+
./pkg/kpp/kpp_routines.F 1128 SMOOTH_HORIZ 
     &                     fld(im1,jp1)*maskC(im1,jp1,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1130 SMOOTH_HORIZ 
     &                     fld(ip1,jp1)*maskC(ip1,jp1,k,bi,bj)))
./verification/advect_cs/code/ini_vel.F 37 INI_VEL 
      INTEGER ip1,jp1
./verification/advect_cs/code/ini_vel.F 56 INI_VEL 
          jp1=MIN(j+1,sNy+Oly)
./verification/advect_cs/code/ini_vel.F 60 INI_VEL 
     &      +(psi(i,j,bi,bj)-psi(i,jp1,bi,bj))*recip_dyG(i,j,bi,bj)
./verification/solid-body.cs-32x32x1/code/ini_vel.F 35 INI_VEL 
      INTEGER ip1,jp1
./verification/solid-body.cs-32x32x1/code/ini_vel.F 53 INI_VEL 
          jp1=MIN(j+1,sNy+Oly)
./verification/solid-body.cs-32x32x1/code/ini_vel.F 57 INI_VEL 
     &      +(psi(i,j,bi,bj)-psi(i,jp1,bi,bj))*recip_dyg(i,j,bi,bj)
./pkg/ctrl/ctrl_getobcsn.F 103 CTRL_GETOBCSN 
      jp1  = 0
./pkg/ctrl/ctrl_getobcsn.F 167 CTRL_GETOBCSN 
     &                maskS(i,j+jp1,k,bi,bj) * delR(k) + vtop
./pkg/ctrl/ctrl_getobcsn.F 169 CTRL_GETOBCSN 
                      if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_getobcsn.F 298 CTRL_GETOBCSN 
     &                maskS(i,j+jp1,k,bi,bj) * delR(k) + vtop
./pkg/ctrl/ctrl_getobcsn.F 300 CTRL_GETOBCSN 
                      if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_getobcsn.F 377 CTRL_GETOBCSN 
     &                         *maskS(i,j+jp1,k,bi,bj)
./pkg/ctrl/ctrl_getobcsn.F 384 CTRL_GETOBCSN 
     &                         *maskS(i,j+jp1,k,bi,bj)
./pkg/ctrl/ctrl_getobcsn.F 398 CTRL_GETOBCSN 
     &                         *maskS(i,j+jp1,k,bi,bj)
./pkg/ctrl/ctrl_getobcsn.F 62 CTRL_GETOBCSN 
      integer jp1
./pkg/ctrl/ctrl_getobcss.F 100 CTRL_GETOBCSS 
      jp1  = 1
./pkg/ctrl/ctrl_getobcss.F 156 CTRL_GETOBCSS 
     &                maskS(i,j+jp1,k,bi,bj) * delR(k) + vtop
./pkg/ctrl/ctrl_getobcss.F 158 CTRL_GETOBCSS 
                      if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_getobcss.F 275 CTRL_GETOBCSS 
     &                maskS(i,j+jp1,k,bi,bj) * delR(k) + vtop
./pkg/ctrl/ctrl_getobcss.F 277 CTRL_GETOBCSS 
                      if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_getobcss.F 352 CTRL_GETOBCSS 
     &                      *maskS(i,j+jp1,k,bi,bj)
./pkg/ctrl/ctrl_getobcss.F 358 CTRL_GETOBCSS 
     &                      *maskS(i,j+jp1,k,bi,bj)
./pkg/ctrl/ctrl_getobcss.F 370 CTRL_GETOBCSS 
     &                      *maskS(i,j+jp1,k,bi,bj)
./pkg/ctrl/ctrl_getobcss.F 70 CTRL_GETOBCSS 
      integer jp1
./pkg/ctrl/ctrl_mask_set_xz.F 107 CTRL_MASK_SET_XZ 
                    if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_mask_set_xz.F 36 CTRL_MASK_SET_XZ 
      integer jp1
./pkg/ctrl/ctrl_mask_set_xz.F 8 CTRL_MASK_SET_XZ 
     &     jp1, OB_J, nwetobcs, ymaskobcs, mythid )
./pkg/ctrl/ctrl_obcsbal.F 66 CTRL_OBCSBAL 
      integer jp1
./pkg/ctrl/ctrl_obcsvol.F 168 CTRL_OBCSVOL 
        jp1 = 0
./pkg/ctrl/ctrl_obcsvol.F 185 CTRL_OBCSVOL 
                if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_obcsvol.F 190 CTRL_OBCSVOL 
                    tmpflux = tmpflux - tmpx* delR(k)*dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 191 CTRL_OBCSVOL 
                    tmparea = tmparea + delR(k) * dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 204 CTRL_OBCSVOL 
        jp1 = 1
./pkg/ctrl/ctrl_obcsvol.F 221 CTRL_OBCSVOL 
                if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_obcsvol.F 226 CTRL_OBCSVOL 
                    tmpflux = tmpflux + tmpx* delR(k)*dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 227 CTRL_OBCSVOL 
                    tmparea = tmparea + delR(k) * dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 327 CTRL_OBCSVOL 
        jp1 = 0
./pkg/ctrl/ctrl_obcsvol.F 343 CTRL_OBCSVOL 
                if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_obcsvol.F 348 CTRL_OBCSVOL 
                    tmpflux = tmpflux - tmpx* delR(k)*dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 349 CTRL_OBCSVOL 
                    tmparea = tmparea + delR(k) * dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 362 CTRL_OBCSVOL 
        jp1 = 1
./pkg/ctrl/ctrl_obcsvol.F 379 CTRL_OBCSVOL 
                if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_obcsvol.F 384 CTRL_OBCSVOL 
                    tmpflux = tmpflux + tmpx* delR(k)*dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 385 CTRL_OBCSVOL 
                    tmparea = tmparea + delR(k) * dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_obcsvol.F 63 CTRL_OBCSVOL 
      integer jp1
./pkg/ctrl/ctrl_smooth.F 28 CTRL_SMOOTH 
      integer i, j, im1, ip1, jm1, jp1
./pkg/ctrl/ctrl_smooth.F 46 CTRL_SMOOTH 
         jp1 = j+1
./pkg/ctrl/ctrl_smooth.F 55 CTRL_SMOOTH 
     &                     mask(i  ,jp1,k,bi,bj) ) +
./pkg/ctrl/ctrl_smooth.F 57 CTRL_SMOOTH 
     &                     mask(im1,jp1,k,bi,bj)   +
./pkg/ctrl/ctrl_smooth.F 59 CTRL_SMOOTH 
     &                     mask(ip1,jp1,k,bi,bj) )
./pkg/ctrl/ctrl_smooth.F 66 CTRL_SMOOTH 
     &                     fld(i  ,jp1,bi,bj)*mask(i  ,jp1,k,bi,bj))+
./pkg/ctrl/ctrl_smooth.F 68 CTRL_SMOOTH 
     &                     fld(im1,jp1,bi,bj)*mask(im1,jp1,k,bi,bj) +
./pkg/ctrl/ctrl_smooth.F 70 CTRL_SMOOTH 
     &                     fld(ip1,jp1,bi,bj)*mask(ip1,jp1,k,bi,bj)))
./pkg/ctrl/ctrl_volflux.F 168 CTRL_VOLFLUX 
        jp1 = 0
./pkg/ctrl/ctrl_volflux.F 187 CTRL_VOLFLUX 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_volflux.F 198 CTRL_VOLFLUX 
                        tmparea = tmparea + delR(k) * dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_volflux.F 202 CTRL_VOLFLUX 
                    tmpflux = tmpflux -tmpx*delR(k)*dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_volflux.F 216 CTRL_VOLFLUX 
        jp1 = 1
./pkg/ctrl/ctrl_volflux.F 235 CTRL_VOLFLUX 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_volflux.F 247 CTRL_VOLFLUX 
                        tmparea = tmparea + delR(k) * dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_volflux.F 252 CTRL_VOLFLUX 
                      tmpflux = tmpflux +tmpx*delR(k)*dxg(i,j+jp1,bi,bj)
./pkg/ctrl/ctrl_volflux.F 69 CTRL_VOLFLUX 
      integer jp1
./pkg/ecco/cost_obcs_ageos.F 217 COST_OBCS_AGEOS 
            jp1 = 0
./pkg/ecco/cost_obcs_ageos.F 228 COST_OBCS_AGEOS 
     &       hfacC(i,j+jp1,k,bi,bj)*hfacC(i+1,j+jp1,k,bi,bj) *
./pkg/ecco/cost_obcs_ageos.F 229 COST_OBCS_AGEOS 
     &       hfacC(i-1,j+jp1,k,bi,bj)*hfacC(i,j+jp1,k+1,bi,bj)*
./pkg/ecco/cost_obcs_ageos.F 230 COST_OBCS_AGEOS 
     &       hfacC(i-1,j+jp1,k+1,bi,bj)*hfacC(i+1,j+jp1,k+1,bi,bj)*
./pkg/ecco/cost_obcs_ageos.F 231 COST_OBCS_AGEOS 
     &       hfacS(i,j+jp1,k,bi,bj)*hfacS(i,j+jp1,k+1,bi,bj)
./pkg/ecco/cost_obcs_ageos.F 253 COST_OBCS_AGEOS 
     &            (rholoc(i-1,j+jp1,bi,bj)-rholoc(i+1,j+jp1,bi,bj))
./pkg/ecco/cost_obcs_ageos.F 254 COST_OBCS_AGEOS 
     &            /(2.*dxc(i,j+jp1,bi,bj))
./pkg/ecco/cost_obcs_ageos.F 264 COST_OBCS_AGEOS 
                  xzdvel1(i,k,bi,bj) = vbar(i,j+jp1,k  ,bi,bj)
./pkg/ecco/cost_obcs_ageos.F 265 COST_OBCS_AGEOS 
     &                               - vbar(i,j+jp1,k+1,bi,bj)
./pkg/ecco/cost_obcs_ageos.F 284 COST_OBCS_AGEOS 
            jp1 = 1
./pkg/ecco/cost_obcs_ageos.F 295 COST_OBCS_AGEOS 
     &       hfacC(i,j+jp1,k,bi,bj)*hfacC(i+1,j+jp1,k,bi,bj) *
./pkg/ecco/cost_obcs_ageos.F 296 COST_OBCS_AGEOS 
     &       hfacC(i-1,j+jp1,k,bi,bj)*hfacC(i,j+jp1,k+1,bi,bj)*
./pkg/ecco/cost_obcs_ageos.F 297 COST_OBCS_AGEOS 
     &       hfacC(i-1,j+jp1,k+1,bi,bj)*hfacC(i+1,j+jp1,k+1,bi,bj)*
./pkg/ecco/cost_obcs_ageos.F 298 COST_OBCS_AGEOS 
     &       hfacS(i,j+jp1,k,bi,bj)*hfacS(i,j+jp1,k+1,bi,bj)
./pkg/ecco/cost_obcs_ageos.F 318 COST_OBCS_AGEOS 
     &            (rholoc(i-1,j+jp1,bi,bj)-rholoc(i+1,j+jp1,bi,bj))
./pkg/ecco/cost_obcs_ageos.F 319 COST_OBCS_AGEOS 
     &            /(2.*dxc(i,j+jp1,bi,bj))
./pkg/ecco/cost_obcs_ageos.F 328 COST_OBCS_AGEOS 
                 xzdvel1(i,k,bi,bj) = vbar(i,j+jp1,k  ,bi,bj)
./pkg/ecco/cost_obcs_ageos.F 329 COST_OBCS_AGEOS 
     &                               - vbar(i,j+jp1,k+1,bi,bj)
./pkg/ecco/cost_obcs_ageos.F 69 COST_OBCS_AGEOS 
      integer ip1, jp1
./pkg/ecco/cost_obcsn.F 101 COST_OBCSN 
      jp1 = 0
./pkg/ecco/cost_obcsn.F 130 COST_OBCSN 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcsn.F 167 COST_OBCSN 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcsn.F 240 COST_OBCSN 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcsn.F 61 COST_OBCSN 
      integer jp1
./pkg/ecco/cost_obcss.F 101 COST_OBCSS 
      jp1 = 1
./pkg/ecco/cost_obcss.F 130 COST_OBCSS 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcss.F 167 COST_OBCSS 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcss.F 240 COST_OBCSS 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcss.F 61 COST_OBCSS 
      integer jp1
./pkg/ecco/cost_obcsvol.F 182 COST_OBCSVOL 
        jp1 = 0
./pkg/ecco/cost_obcsvol.F 209 COST_OBCSVOL 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcsvol.F 211 COST_OBCSVOL 
                    fctile = fctile - tmpx* delZ(k) *dxg(i,j+jp1,bi,bj)
./pkg/ecco/cost_obcsvol.F 225 COST_OBCSVOL 
        jp1 = 1
./pkg/ecco/cost_obcsvol.F 251 COST_OBCSVOL 
                  if (maskS(i,j+jp1,k,bi,bj) .ne. 0.) then
./pkg/ecco/cost_obcsvol.F 253 COST_OBCSVOL 
                    fctile = fctile + tmpx* delZ(k) *dxg(i,j+jp1,bi,bj)
./pkg/ecco/cost_obcsvol.F 66 COST_OBCSVOL 
      integer jp1