Symbol Description Uses
jp 454



File Line number Procedure Code
./pkg/exch2/w2_print_e2setup.F 171 W2_PRINT_E2SETUP
      jp = 0
./pkg/exch2/w2_print_e2setup.F 179 W2_PRINT_E2SETUP
         IF ( js .NE. jp )
./pkg/exch2/w2_print_e2setup.F 183 W2_PRINT_E2SETUP
         jp = js
./pkg/exch2/w2_print_e2setup.F 42 W2_PRINT_E2SETUP
      INTEGER i, j, js, jp, jt, ii, is, it, ns, nt, k, tx, ty
./pkg/exch2/w2_print_e2setup.F 58 W2_PRINT_E2SETUP
      jp = 0
./pkg/exch2/w2_print_e2setup.F 68 W2_PRINT_E2SETUP
         IF ( js.NE.jp ) THEN
./pkg/exch2/w2_print_e2setup.F 69 W2_PRINT_E2SETUP
           IF ( jp.NE.0 ) THEN
./pkg/exch2/w2_print_e2setup.F 91 W2_PRINT_E2SETUP
           jp = js
./pkg/flt/flt_exch2.F 183 FLT_EXCH2
             DO jp = npNew+1,npart_tile(bi,bj)
./pkg/flt/flt_exch2.F 184 FLT_EXCH2
              IF ( npart(jp,bi,bj).NE.0. _d 0 ) THEN
./pkg/flt/flt_exch2.F 188 FLT_EXCH2
                npart (ip,bi,bj) =  npart (jp,bi,bj)
./pkg/flt/flt_exch2.F 189 FLT_EXCH2
                tstart(ip,bi,bj) =  tstart(jp,bi,bj)
./pkg/flt/flt_exch2.F 190 FLT_EXCH2
                ipart (ip,bi,bj) =  ipart (jp,bi,bj)
./pkg/flt/flt_exch2.F 191 FLT_EXCH2
                jpart (ip,bi,bj) =  jpart (jp,bi,bj)
./pkg/flt/flt_exch2.F 192 FLT_EXCH2
                kpart (ip,bi,bj) =  kpart (jp,bi,bj)
./pkg/flt/flt_exch2.F 193 FLT_EXCH2
                kfloat(ip,bi,bj) =  kfloat(jp,bi,bj)
./pkg/flt/flt_exch2.F 194 FLT_EXCH2
                iup   (ip,bi,bj) =  iup   (jp,bi,bj)
./pkg/flt/flt_exch2.F 195 FLT_EXCH2
                itop  (ip,bi,bj) =  itop  (jp,bi,bj)
./pkg/flt/flt_exch2.F 196 FLT_EXCH2
                tend  (ip,bi,bj) =  tend  (jp,bi,bj)
./pkg/flt/flt_exch2.F 312 FLT_EXCH2
             DO jp = npNew+1,npart_tile(bi,bj)
./pkg/flt/flt_exch2.F 313 FLT_EXCH2
              IF ( npart(jp,bi,bj).NE.0. _d 0 ) THEN
./pkg/flt/flt_exch2.F 317 FLT_EXCH2
                npart (ip,bi,bj) =  npart (jp,bi,bj)
./pkg/flt/flt_exch2.F 318 FLT_EXCH2
                tstart(ip,bi,bj) =  tstart(jp,bi,bj)
./pkg/flt/flt_exch2.F 319 FLT_EXCH2
                ipart (ip,bi,bj) =  ipart (jp,bi,bj)
./pkg/flt/flt_exch2.F 320 FLT_EXCH2
                jpart (ip,bi,bj) =  jpart (jp,bi,bj)
./pkg/flt/flt_exch2.F 321 FLT_EXCH2
                kpart (ip,bi,bj) =  kpart (jp,bi,bj)
./pkg/flt/flt_exch2.F 322 FLT_EXCH2
                kfloat(ip,bi,bj) =  kfloat(jp,bi,bj)
./pkg/flt/flt_exch2.F 323 FLT_EXCH2
                iup   (ip,bi,bj) =  iup   (jp,bi,bj)
./pkg/flt/flt_exch2.F 324 FLT_EXCH2
                itop  (ip,bi,bj) =  itop  (jp,bi,bj)
./pkg/flt/flt_exch2.F 325 FLT_EXCH2
                tend  (ip,bi,bj) =  tend  (jp,bi,bj)
./pkg/flt/flt_exch2.F 379 FLT_EXCH2
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exch2.F 380 FLT_EXCH2
               npart( jp,bi,bj) = fltbuf_recv(ic+1,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 381 FLT_EXCH2
               tstart(jp,bi,bj) = fltbuf_recv(ic+2,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 382 FLT_EXCH2
               ipart( jp,bi,bj) = fltbuf_recv(ic+3,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 383 FLT_EXCH2
               jpart( jp,bi,bj) = fltbuf_recv(ic+4,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 384 FLT_EXCH2
               kpart( jp,bi,bj) = fltbuf_recv(ic+5,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 385 FLT_EXCH2
               kfloat(jp,bi,bj) = fltbuf_recv(ic+6,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 386 FLT_EXCH2
               iup(   jp,bi,bj) = fltbuf_recv(ic+7,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 387 FLT_EXCH2
               itop(  jp,bi,bj) = fltbuf_recv(ic+8,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 388 FLT_EXCH2
               tend(  jp,bi,bj) = fltbuf_recv(ic+9,bi,bj,W2_EAST)
./pkg/flt/flt_exch2.F 408 FLT_EXCH2
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exch2.F 409 FLT_EXCH2
               npart( jp,bi,bj) = fltbuf_recv(ic+1,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 410 FLT_EXCH2
               tstart(jp,bi,bj) = fltbuf_recv(ic+2,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 411 FLT_EXCH2
               ipart( jp,bi,bj) = fltbuf_recv(ic+3,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 412 FLT_EXCH2
               jpart( jp,bi,bj) = fltbuf_recv(ic+4,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 413 FLT_EXCH2
               kpart( jp,bi,bj) = fltbuf_recv(ic+5,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 414 FLT_EXCH2
               kfloat(jp,bi,bj) = fltbuf_recv(ic+6,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 415 FLT_EXCH2
               iup(   jp,bi,bj) = fltbuf_recv(ic+7,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 416 FLT_EXCH2
               itop(  jp,bi,bj) = fltbuf_recv(ic+8,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 417 FLT_EXCH2
               tend(  jp,bi,bj) = fltbuf_recv(ic+9,bi,bj,W2_WEST)
./pkg/flt/flt_exch2.F 437 FLT_EXCH2
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exch2.F 438 FLT_EXCH2
               npart( jp,bi,bj) = fltbuf_recv(ic+1,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 439 FLT_EXCH2
               tstart(jp,bi,bj) = fltbuf_recv(ic+2,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 43 FLT_EXCH2
      INTEGER ip, jp, jl, npNew
./pkg/flt/flt_exch2.F 440 FLT_EXCH2
               ipart( jp,bi,bj) = fltbuf_recv(ic+3,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 441 FLT_EXCH2
               jpart( jp,bi,bj) = fltbuf_recv(ic+4,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 442 FLT_EXCH2
               kpart( jp,bi,bj) = fltbuf_recv(ic+5,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 443 FLT_EXCH2
               kfloat(jp,bi,bj) = fltbuf_recv(ic+6,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 444 FLT_EXCH2
               iup(   jp,bi,bj) = fltbuf_recv(ic+7,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 445 FLT_EXCH2
               itop(  jp,bi,bj) = fltbuf_recv(ic+8,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 446 FLT_EXCH2
               tend(  jp,bi,bj) = fltbuf_recv(ic+9,bi,bj,W2_NORTH)
./pkg/flt/flt_exch2.F 466 FLT_EXCH2
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exch2.F 467 FLT_EXCH2
               npart( jp,bi,bj) = fltbuf_recv(ic+1,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 468 FLT_EXCH2
               tstart(jp,bi,bj) = fltbuf_recv(ic+2,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 469 FLT_EXCH2
               ipart( jp,bi,bj) = fltbuf_recv(ic+3,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 470 FLT_EXCH2
               jpart( jp,bi,bj) = fltbuf_recv(ic+4,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 471 FLT_EXCH2
               kpart( jp,bi,bj) = fltbuf_recv(ic+5,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 472 FLT_EXCH2
               kfloat(jp,bi,bj) = fltbuf_recv(ic+6,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 473 FLT_EXCH2
               iup(   jp,bi,bj) = fltbuf_recv(ic+7,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 474 FLT_EXCH2
               itop(  jp,bi,bj) = fltbuf_recv(ic+8,bi,bj,W2_SOUTH)
./pkg/flt/flt_exch2.F 475 FLT_EXCH2
               tend(  jp,bi,bj) = fltbuf_recv(ic+9,bi,bj,W2_SOUTH)
./pkg/flt/flt_exchg.F 187 FLT_EXCHG
             DO jp = npNew+1,npart_tile(bi,bj)
./pkg/flt/flt_exchg.F 188 FLT_EXCHG
              IF ( npart(jp,bi,bj).NE.0. _d 0 ) THEN
./pkg/flt/flt_exchg.F 192 FLT_EXCHG
                npart (ip,bi,bj) =  npart (jp,bi,bj)
./pkg/flt/flt_exchg.F 193 FLT_EXCHG
                tstart(ip,bi,bj) =  tstart(jp,bi,bj)
./pkg/flt/flt_exchg.F 194 FLT_EXCHG
                ipart (ip,bi,bj) =  ipart (jp,bi,bj)
./pkg/flt/flt_exchg.F 195 FLT_EXCHG
                jpart (ip,bi,bj) =  jpart (jp,bi,bj)
./pkg/flt/flt_exchg.F 196 FLT_EXCHG
                kpart (ip,bi,bj) =  kpart (jp,bi,bj)
./pkg/flt/flt_exchg.F 197 FLT_EXCHG
                kfloat(ip,bi,bj) =  kfloat(jp,bi,bj)
./pkg/flt/flt_exchg.F 198 FLT_EXCHG
                iup   (ip,bi,bj) =  iup   (jp,bi,bj)
./pkg/flt/flt_exchg.F 199 FLT_EXCHG
                itop  (ip,bi,bj) =  itop  (jp,bi,bj)
./pkg/flt/flt_exchg.F 200 FLT_EXCHG
                tend  (ip,bi,bj) =  tend  (jp,bi,bj)
./pkg/flt/flt_exchg.F 240 FLT_EXCHG
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exchg.F 241 FLT_EXCHG
               npart( jp,bi,bj) = fltbuf_recvE(ic+1,bi,bj)
./pkg/flt/flt_exchg.F 242 FLT_EXCHG
               tstart(jp,bi,bj) = fltbuf_recvE(ic+2,bi,bj)
./pkg/flt/flt_exchg.F 243 FLT_EXCHG
               ipart( jp,bi,bj) = fltbuf_recvE(ic+3,bi,bj)
./pkg/flt/flt_exchg.F 244 FLT_EXCHG
               jpart( jp,bi,bj) = fltbuf_recvE(ic+4,bi,bj)
./pkg/flt/flt_exchg.F 245 FLT_EXCHG
               kpart( jp,bi,bj) = fltbuf_recvE(ic+5,bi,bj)
./pkg/flt/flt_exchg.F 246 FLT_EXCHG
               kfloat(jp,bi,bj) = fltbuf_recvE(ic+6,bi,bj)
./pkg/flt/flt_exchg.F 247 FLT_EXCHG
               iup(   jp,bi,bj) = fltbuf_recvE(ic+7,bi,bj)
./pkg/flt/flt_exchg.F 248 FLT_EXCHG
               itop(  jp,bi,bj) = fltbuf_recvE(ic+8,bi,bj)
./pkg/flt/flt_exchg.F 249 FLT_EXCHG
               tend(  jp,bi,bj) = fltbuf_recvE(ic+9,bi,bj)
./pkg/flt/flt_exchg.F 269 FLT_EXCHG
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exchg.F 270 FLT_EXCHG
               npart( jp,bi,bj) = fltbuf_recvW(ic+1,bi,bj)
./pkg/flt/flt_exchg.F 271 FLT_EXCHG
               tstart(jp,bi,bj) = fltbuf_recvW(ic+2,bi,bj)
./pkg/flt/flt_exchg.F 272 FLT_EXCHG
               ipart( jp,bi,bj) = fltbuf_recvW(ic+3,bi,bj)
./pkg/flt/flt_exchg.F 273 FLT_EXCHG
               jpart( jp,bi,bj) = fltbuf_recvW(ic+4,bi,bj)
./pkg/flt/flt_exchg.F 274 FLT_EXCHG
               kpart( jp,bi,bj) = fltbuf_recvW(ic+5,bi,bj)
./pkg/flt/flt_exchg.F 275 FLT_EXCHG
               kfloat(jp,bi,bj) = fltbuf_recvW(ic+6,bi,bj)
./pkg/flt/flt_exchg.F 276 FLT_EXCHG
               iup(   jp,bi,bj) = fltbuf_recvW(ic+7,bi,bj)
./pkg/flt/flt_exchg.F 277 FLT_EXCHG
               itop(  jp,bi,bj) = fltbuf_recvW(ic+8,bi,bj)
./pkg/flt/flt_exchg.F 278 FLT_EXCHG
               tend(  jp,bi,bj) = fltbuf_recvW(ic+9,bi,bj)
./pkg/flt/flt_exchg.F 430 FLT_EXCHG
             DO jp = npNew+1,npart_tile(bi,bj)
./pkg/flt/flt_exchg.F 431 FLT_EXCHG
              IF ( npart(jp,bi,bj).NE.0. _d 0 ) THEN
./pkg/flt/flt_exchg.F 435 FLT_EXCHG
                npart (ip,bi,bj) =  npart (jp,bi,bj)
./pkg/flt/flt_exchg.F 436 FLT_EXCHG
                tstart(ip,bi,bj) =  tstart(jp,bi,bj)
./pkg/flt/flt_exchg.F 437 FLT_EXCHG
                ipart (ip,bi,bj) =  ipart (jp,bi,bj)
./pkg/flt/flt_exchg.F 438 FLT_EXCHG
                jpart (ip,bi,bj) =  jpart (jp,bi,bj)
./pkg/flt/flt_exchg.F 439 FLT_EXCHG
                kpart (ip,bi,bj) =  kpart (jp,bi,bj)
./pkg/flt/flt_exchg.F 440 FLT_EXCHG
                kfloat(ip,bi,bj) =  kfloat(jp,bi,bj)
./pkg/flt/flt_exchg.F 441 FLT_EXCHG
                iup   (ip,bi,bj) =  iup   (jp,bi,bj)
./pkg/flt/flt_exchg.F 442 FLT_EXCHG
                itop  (ip,bi,bj) =  itop  (jp,bi,bj)
./pkg/flt/flt_exchg.F 443 FLT_EXCHG
                tend  (ip,bi,bj) =  tend  (jp,bi,bj)
./pkg/flt/flt_exchg.F 47 FLT_EXCHG
      INTEGER ip, jp, jl, m, npNew
./pkg/flt/flt_exchg.F 483 FLT_EXCHG
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exchg.F 484 FLT_EXCHG
               npart( jp,bi,bj) = fltbuf_recvE(ic+1,bi,bj)
./pkg/flt/flt_exchg.F 485 FLT_EXCHG
               tstart(jp,bi,bj) = fltbuf_recvE(ic+2,bi,bj)
./pkg/flt/flt_exchg.F 486 FLT_EXCHG
               ipart( jp,bi,bj) = fltbuf_recvE(ic+3,bi,bj)
./pkg/flt/flt_exchg.F 487 FLT_EXCHG
               jpart( jp,bi,bj) = fltbuf_recvE(ic+4,bi,bj)
./pkg/flt/flt_exchg.F 488 FLT_EXCHG
               kpart( jp,bi,bj) = fltbuf_recvE(ic+5,bi,bj)
./pkg/flt/flt_exchg.F 489 FLT_EXCHG
               kfloat(jp,bi,bj) = fltbuf_recvE(ic+6,bi,bj)
./pkg/flt/flt_exchg.F 490 FLT_EXCHG
               iup(   jp,bi,bj) = fltbuf_recvE(ic+7,bi,bj)
./pkg/flt/flt_exchg.F 491 FLT_EXCHG
               itop(  jp,bi,bj) = fltbuf_recvE(ic+8,bi,bj)
./pkg/flt/flt_exchg.F 492 FLT_EXCHG
               tend(  jp,bi,bj) = fltbuf_recvE(ic+9,bi,bj)
./pkg/flt/flt_exchg.F 512 FLT_EXCHG
               jp = npart_tile(bi,bj)
./pkg/flt/flt_exchg.F 513 FLT_EXCHG
               npart( jp,bi,bj) = fltbuf_recvW(ic+1,bi,bj)
./pkg/flt/flt_exchg.F 514 FLT_EXCHG
               tstart(jp,bi,bj) = fltbuf_recvW(ic+2,bi,bj)
./pkg/flt/flt_exchg.F 515 FLT_EXCHG
               ipart( jp,bi,bj) = fltbuf_recvW(ic+3,bi,bj)
./pkg/flt/flt_exchg.F 516 FLT_EXCHG
               jpart( jp,bi,bj) = fltbuf_recvW(ic+4,bi,bj)
./pkg/flt/flt_exchg.F 517 FLT_EXCHG
               kpart( jp,bi,bj) = fltbuf_recvW(ic+5,bi,bj)
./pkg/flt/flt_exchg.F 518 FLT_EXCHG
               kfloat(jp,bi,bj) = fltbuf_recvW(ic+6,bi,bj)
./pkg/flt/flt_exchg.F 519 FLT_EXCHG
               iup(   jp,bi,bj) = fltbuf_recvW(ic+7,bi,bj)
./pkg/flt/flt_exchg.F 520 FLT_EXCHG
               itop(  jp,bi,bj) = fltbuf_recvW(ic+8,bi,bj)
./pkg/flt/flt_exchg.F 521 FLT_EXCHG
               tend(  jp,bi,bj) = fltbuf_recvW(ic+9,bi,bj)
./pkg/mdsio/mdsio_gl.F 1011 MDSREADFIELD_2D_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl.F 1019 MDSREADFIELD_2D_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl.F 1057 MDSREADFIELD_2D_GL
            jG=bj+(jp-1)*nsy
./pkg/mdsio/mdsio_gl.F 1106 MDSREADFIELD_2D_GL
        arr_gl(ii,bi,ip,j,bj,jp,k)=arr(ii,j,k,bi,bj)
./pkg/mdsio/mdsio_gl.F 1172 MDSREADFIELD_2D_GL
        DO jp=1,nPy
./pkg/mdsio/mdsio_gl.F 1177 MDSREADFIELD_2D_GL
             JJ=((jp-1)*nSy+(bj-1))*sNy+J
./pkg/mdsio/mdsio_gl.F 1180 MDSREADFIELD_2D_GL
              arr_gl(i,bi,ip,j,bj,jp,k) = global(II,JJ)
./pkg/mdsio/mdsio_gl.F 1280 MDSWRITEFIELD_2D_GL
      integer ip,jp,iG,jG,irec,bi,bj,i,j,k,dUnit,IL
./pkg/mdsio/mdsio_gl.F 1377 MDSWRITEFIELD_2D_GL
        do jp=1,nPy
./pkg/mdsio/mdsio_gl.F 1382 MDSWRITEFIELD_2D_GL
           JJ=((jp-1)*nSy+(bj-1))*sNy+J
./pkg/mdsio/mdsio_gl.F 1385 MDSWRITEFIELD_2D_GL
            global(II,JJ) = arr_gl(i,bi,ip,j,bj,jp,k)
./pkg/mdsio/mdsio_gl.F 1458 MDSWRITEFIELD_2D_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl.F 1465 MDSWRITEFIELD_2D_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl.F 1483 MDSWRITEFIELD_2D_GL
                arr(i,j,k,bi,bj)=arr_gl(i,bi,ip,j,bj,jp,k)
./pkg/mdsio/mdsio_gl.F 1545 MDSWRITEFIELD_2D_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl.F 1552 MDSWRITEFIELD_2D_GL
         dimList(2,2)=((jp-1)*nSy+(bj-1))*sNy+1
./pkg/mdsio/mdsio_gl.F 1553 MDSWRITEFIELD_2D_GL
         dimList(3,2)=((jp-1)*nSy+bj)*sNy
./pkg/mdsio/mdsio_gl.F 224 MDSREADFIELD_3D_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl.F 232 MDSREADFIELD_3D_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl.F 270 MDSREADFIELD_3D_GL
            jG=bj+(jp-1)*nsy
./pkg/mdsio/mdsio_gl.F 319 MDSREADFIELD_3D_GL
        arr_gl(ii,bi,ip,j,bj,jp,k)=arr(ii,j,k,bi,bj)
./pkg/mdsio/mdsio_gl.F 385 MDSREADFIELD_3D_GL
        DO jp=1,nPy
./pkg/mdsio/mdsio_gl.F 390 MDSREADFIELD_3D_GL
             JJ=((jp-1)*nSy+(bj-1))*sNy+J
./pkg/mdsio/mdsio_gl.F 393 MDSREADFIELD_3D_GL
              arr_gl(i,bi,ip,j,bj,jp,k) = global(II,JJ)
./pkg/mdsio/mdsio_gl.F 492 MDSWRITEFIELD_3D_GL
      integer ip,jp,iG,jG,irec,bi,bj,i,j,k,dUnit,IL
./pkg/mdsio/mdsio_gl.F 588 MDSWRITEFIELD_3D_GL
        do jp=1,nPy
./pkg/mdsio/mdsio_gl.F 593 MDSWRITEFIELD_3D_GL
           JJ=((jp-1)*nSy+(bj-1))*sNy+J
./pkg/mdsio/mdsio_gl.F 596 MDSWRITEFIELD_3D_GL
            global(II,JJ) = arr_gl(i,bi,ip,j,bj,jp,k)
./pkg/mdsio/mdsio_gl.F 669 MDSWRITEFIELD_3D_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl.F 676 MDSWRITEFIELD_3D_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl.F 694 MDSWRITEFIELD_3D_GL
                arr(i,j,k,bi,bj)=arr_gl(i,bi,ip,j,bj,jp,k)
./pkg/mdsio/mdsio_gl.F 73 MDSREADFIELD_3D_GL
      integer ip,jp,iG,jG,irec,bi,bj,i,j,k,dUnit,IL
./pkg/mdsio/mdsio_gl.F 756 MDSWRITEFIELD_3D_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl.F 763 MDSWRITEFIELD_3D_GL
         dimList(2,2)=((jp-1)*nSy+(bj-1))*sNy+1
./pkg/mdsio/mdsio_gl.F 764 MDSWRITEFIELD_3D_GL
         dimList(3,2)=((jp-1)*nSy+bj)*sNy
./pkg/mdsio/mdsio_gl.F 860 MDSREADFIELD_2D_GL
      integer ip,jp,iG,jG,irec,bi,bj,i,j,k,dUnit,IL
./pkg/mdsio/mdsio_gl_slice.F 138 MDSREADFIELD_XZ_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl_slice.F 146 MDSREADFIELD_XZ_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl_slice.F 222 MDSREADFIELD_XZ_GL
        arr_gl(ii,bi,ip,bj,jp,k)=arr(ii,k,bi,bj)
./pkg/mdsio/mdsio_gl_slice.F 313 MDSREADFIELD_YZ_GL
      integer ip,jp,iG,jG,irec,bi,bj,jj,k,dUnit,IL
./pkg/mdsio/mdsio_gl_slice.F 378 MDSREADFIELD_YZ_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl_slice.F 386 MDSREADFIELD_YZ_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl_slice.F 462 MDSREADFIELD_YZ_GL
        arr_gl(bi,ip,jj,bj,jp,k)=arr(jj,k,bi,bj)
./pkg/mdsio/mdsio_gl_slice.F 567 MDSWRITEFIELD_XZ_GL
      integer ip,jp,iG,jG,irec,bi,bj,ii,k,dUnit,IL
./pkg/mdsio/mdsio_gl_slice.F 615 MDSWRITEFIELD_XZ_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl_slice.F 622 MDSWRITEFIELD_XZ_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl_slice.F 639 MDSWRITEFIELD_XZ_GL
               arr(ii,k,bi,bj)=arr_gl(ii,bi,ip,bj,jp,k)
./pkg/mdsio/mdsio_gl_slice.F 699 MDSWRITEFIELD_XZ_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl_slice.F 706 MDSWRITEFIELD_XZ_GL
         dimList(2,2)=(jp-1)*nSy+bj
./pkg/mdsio/mdsio_gl_slice.F 707 MDSWRITEFIELD_XZ_GL
         dimList(3,2)= jp*nSy+bj
./pkg/mdsio/mdsio_gl_slice.F 72 MDSREADFIELD_XZ_GL
      integer ip,jp,iG,jG,irec,bi,bj,ii,k,dUnit,IL
./pkg/mdsio/mdsio_gl_slice.F 808 MDSWRITEFIELD_YZ_GL
      integer ip,jp,iG,jG,irec,bi,bj,jj,k,dUnit,IL
./pkg/mdsio/mdsio_gl_slice.F 856 MDSWRITEFIELD_YZ_GL
      do jp=1,nPy
./pkg/mdsio/mdsio_gl_slice.F 863 MDSWRITEFIELD_YZ_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl_slice.F 880 MDSWRITEFIELD_YZ_GL
               arr(jj,k,bi,bj)=arr_gl(bi,ip,jj,bj,jp,k)
./pkg/mdsio/mdsio_gl_slice.F 940 MDSWRITEFIELD_YZ_GL
         jG=bj+(jp-1)*nsy ! Kludge until unstructered tiles
./pkg/mdsio/mdsio_gl_slice.F 947 MDSWRITEFIELD_YZ_GL
         dimList(2,2)=((jp-1)*nSy+(bj-1))*sNy+1
./pkg/mdsio/mdsio_gl_slice.F 948 MDSWRITEFIELD_YZ_GL
         dimList(3,2)=((jp-1)*nSy+bj)*sNy
./model/src/solve_uv_tridiago.F 349 SOLVE_UV_TRIDIAGO
            jp = j+1
./model/src/solve_uv_tridiago.F 352 SOLVE_UV_TRIDIAGO
     &                       - gamV(i,j,bi,bj)*yy_V(i,jp,bi,bj)
./model/src/solve_uv_tridiago.F 354 SOLVE_UV_TRIDIAGO
     &                       - gamV(i,j,bi,bj)*alpV(i,jp,bi,bj)
./model/src/solve_uv_tridiago.F 355 SOLVE_UV_TRIDIAGO
            gamV(i,j,bi,bj) = -gamV(i,j,bi,bj)*gamV(i,jp,bi,bj)
./model/src/solve_uv_tridiago.F 80 SOLVE_UV_TRIDIAGO
      INTEGER jj, jm, jp
./pkg/ctrl/ctrl_set_globfld_xy.F 46 CTRL_SET_GLOBFLD_XY
      integer ip,jp
./pkg/ctrl/ctrl_set_globfld_xyz.F 41 CTRL_SET_GLOBFLD_XYZ
      integer ip,jp
./pkg/ctrl/ctrl_set_globfld_xz.F 46 CTRL_SET_GLOBFLD_XZ
      integer ip,jp
./pkg/ctrl/ctrl_set_globfld_yz.F 46 CTRL_SET_GLOBFLD_YZ
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_xy.F 101 CTRL_SET_PACK_XY
             globfld3d  (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xy.F 102 CTRL_SET_PACK_XY
             globmsk    (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xy.F 103 CTRL_SET_PACK_XY
             globfldtmp2(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xy.F 104 CTRL_SET_PACK_XY
             globfldtmp3(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xy.F 153 CTRL_SET_PACK_XY
        do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xy.F 159 CTRL_SET_PACK_XY
              weightfld2d(i,bi,ip,j,bj,jp) = 1. _d 0
./pkg/ctrl/ctrl_set_pack_xy.F 186 CTRL_SET_PACK_XY
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xy.F 192 CTRL_SET_PACK_XY
                  if (globmsk(i,bi,ip,j,bj,jp,1)  .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_xy.F 195 CTRL_SET_PACK_XY
                     globfldtmp3(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_pack_xy.F 196 CTRL_SET_PACK_XY
     &                    globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 202 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k) *
./pkg/ctrl/ctrl_set_pack_xy.F 203 CTRL_SET_PACK_XY
     &                       sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_pack_xy.F 207 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k) /
./pkg/ctrl/ctrl_set_pack_xy.F 208 CTRL_SET_PACK_XY
     &                       sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_pack_xy.F 214 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 218 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 225 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 229 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 234 CTRL_SET_PACK_XY
                     globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_xy.F 280 CTRL_SET_PACK_XY
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xy.F 286 CTRL_SET_PACK_XY
                  if (globmsk(i,bi,ip,j,bj,jp,k)  .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_xy.F 289 CTRL_SET_PACK_XY
                     globfldtmp3(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_pack_xy.F 290 CTRL_SET_PACK_XY
     &                    globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 296 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k) *
./pkg/ctrl/ctrl_set_pack_xy.F 297 CTRL_SET_PACK_XY
     &                       sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_pack_xy.F 301 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k) /
./pkg/ctrl/ctrl_set_pack_xy.F 302 CTRL_SET_PACK_XY
     &                       sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_pack_xy.F 308 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 312 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 319 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 323 CTRL_SET_PACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xy.F 328 CTRL_SET_PACK_XY
                     globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_xy.F 372 CTRL_SET_PACK_XY
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_xy.F 48 CTRL_SET_PACK_XY
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_xy.F 95 CTRL_SET_PACK_XY
       do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xyz.F 104 CTRL_SET_PACK_XYZ
       do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xyz.F 110 CTRL_SET_PACK_XYZ
             globfld3d  (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xyz.F 111 CTRL_SET_PACK_XYZ
             globmsk    (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xyz.F 112 CTRL_SET_PACK_XYZ
             globfldtmp2(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xyz.F 113 CTRL_SET_PACK_XYZ
             globfldtmp3(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xyz.F 163 CTRL_SET_PACK_XYZ
        do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xyz.F 169 CTRL_SET_PACK_XYZ
              weightfld3d(i,bi,ip,j,bj,jp,k) = 1. _d 0
./pkg/ctrl/ctrl_set_pack_xyz.F 202 CTRL_SET_PACK_XYZ
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xyz.F 208 CTRL_SET_PACK_XYZ
                  if (globmsk(i,bi,ip,j,bj,jp,k) .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_xyz.F 211 CTRL_SET_PACK_XYZ
                     globfldtmp3(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_pack_xyz.F 212 CTRL_SET_PACK_XYZ
     &                    globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xyz.F 218 CTRL_SET_PACK_XYZ
     &                       * globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xyz.F 220 CTRL_SET_PACK_XYZ
     &                       * sqrt(weightfld3d(i,bi,ip,j,bj,jp,k))
./pkg/ctrl/ctrl_set_pack_xyz.F 226 CTRL_SET_PACK_XYZ
     &                       * globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xyz.F 228 CTRL_SET_PACK_XYZ
     &                       / sqrt(weightfld3d(i,bi,ip,j,bj,jp,k))
./pkg/ctrl/ctrl_set_pack_xyz.F 234 CTRL_SET_PACK_XYZ
                     globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_xyz.F 237 CTRL_SET_PACK_XYZ
                     cbuff(cbuffindex) = globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xyz.F 240 CTRL_SET_PACK_XYZ
                     cbuff(cbuffindex) = globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xyz.F 284 CTRL_SET_PACK_XYZ
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_xyz.F 52 CTRL_SET_PACK_XYZ
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_xz.F 103 CTRL_SET_PACK_XZ
                        globfldxz  (i,bi,ip,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xz.F 104 CTRL_SET_PACK_XZ
                        globfldtmp2(i,bi,ip,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_pack_xz.F 105 CTRL_SET_PACK_XZ
                        globfldtmp3(i,bi,ip,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_pack_xz.F 107 CTRL_SET_PACK_XZ
                           globmskxz(i,bi,ip,bj,jp,k,iobcs) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xz.F 117 CTRL_SET_PACK_XZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xz.F 123 CTRL_SET_PACK_XZ
                           globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_xz.F 196 CTRL_SET_PACK_XZ
             do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xz.F 203 CTRL_SET_PACK_XZ
                  if (globmskxz(i,bi,ip,bj,jp,k,iobcs)  .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_xz.F 206 CTRL_SET_PACK_XZ
                     globfldtmp3(i,bi,ip,bj,jp) =
./pkg/ctrl/ctrl_set_pack_xz.F 207 CTRL_SET_PACK_XZ
     &                    globfld3d(i,bi,ip,jj,bj,jp,kk)
./pkg/ctrl/ctrl_set_pack_xz.F 212 CTRL_SET_PACK_XZ
     &                       globfld3d(i,bi,ip,jj,bj,jp,kk) *
./pkg/ctrl/ctrl_set_pack_xz.F 214 CTRL_SET_PACK_XZ
     &                       sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_xz.F 220 CTRL_SET_PACK_XZ
     &                       globfld3d(i,bi,ip,jj,bj,jp,kk) /
./pkg/ctrl/ctrl_set_pack_xz.F 222 CTRL_SET_PACK_XZ
     &                       sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_xz.F 228 CTRL_SET_PACK_XZ
                     globfldtmp2(i,bi,ip,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_xz.F 231 CTRL_SET_PACK_XZ
                     cbuff(cbuffindex) = globfld3d(i,bi,ip,jj,bj,jp,kk)
./pkg/ctrl/ctrl_set_pack_xz.F 271 CTRL_SET_PACK_XZ
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_xz.F 276 CTRL_SET_PACK_XZ
                 if (globmskxz(i,bi,ip,bj,jp,k,iobcs)  .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_xz.F 279 CTRL_SET_PACK_XZ
                     globfldtmp3(i,bi,ip,bj,jp) =
./pkg/ctrl/ctrl_set_pack_xz.F 280 CTRL_SET_PACK_XZ
     &                    globfldxz(i,bi,ip,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xz.F 285 CTRL_SET_PACK_XZ
     &                       globfldxz(i,bi,ip,bj,jp,k) *
./pkg/ctrl/ctrl_set_pack_xz.F 287 CTRL_SET_PACK_XZ
     &                       sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_xz.F 293 CTRL_SET_PACK_XZ
     &                       globfldxz(i,bi,ip,bj,jp,k) /
./pkg/ctrl/ctrl_set_pack_xz.F 295 CTRL_SET_PACK_XZ
     &                       sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_xz.F 301 CTRL_SET_PACK_XZ
                     globfldtmp2(i,bi,ip,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_xz.F 304 CTRL_SET_PACK_XZ
                     cbuff(cbuffindex) = globfldxz(i,bi,ip,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_xz.F 54 CTRL_SET_PACK_XZ
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_xz.F 98 CTRL_SET_PACK_XZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_yz.F 103 CTRL_SET_PACK_YZ
                        globfldyz  (bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_yz.F 104 CTRL_SET_PACK_YZ
                        globfldtmp2(bi,ip,j,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_pack_yz.F 105 CTRL_SET_PACK_YZ
                        globfldtmp3(bi,ip,j,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_pack_yz.F 107 CTRL_SET_PACK_YZ
                           globmskyz(bi,ip,j,bj,jp,k,iobcs) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_yz.F 117 CTRL_SET_PACK_YZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_yz.F 123 CTRL_SET_PACK_YZ
                           globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_pack_yz.F 199 CTRL_SET_PACK_YZ
             do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_yz.F 207 CTRL_SET_PACK_YZ
                  if (globmskyz(bi,ip,j,bj,jp,k,iobcs)  .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_yz.F 210 CTRL_SET_PACK_YZ
                     globfldtmp3(bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_pack_yz.F 211 CTRL_SET_PACK_YZ
     &                    globfld3d(ii,bi,ip,jj,bj,jp,kk)
./pkg/ctrl/ctrl_set_pack_yz.F 216 CTRL_SET_PACK_YZ
     &                       globfld3d(ii,bi,ip,jj,bj,jp,kk) *
./pkg/ctrl/ctrl_set_pack_yz.F 218 CTRL_SET_PACK_YZ
     &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_yz.F 224 CTRL_SET_PACK_YZ
     &                       globfld3d(ii,bi,ip,jj,bj,jp,kk) /
./pkg/ctrl/ctrl_set_pack_yz.F 226 CTRL_SET_PACK_YZ
     &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_yz.F 232 CTRL_SET_PACK_YZ
                     globfldtmp2(bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_yz.F 235 CTRL_SET_PACK_YZ
                     cbuff(cbuffindex) = globfld3d(ii,bi,ip,jj,bj,jp,kk)
./pkg/ctrl/ctrl_set_pack_yz.F 275 CTRL_SET_PACK_YZ
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_pack_yz.F 280 CTRL_SET_PACK_YZ
                 if (globmskyz(bi,ip,j,bj,jp,k,iobcs)  .ne. 0. ) then
./pkg/ctrl/ctrl_set_pack_yz.F 283 CTRL_SET_PACK_YZ
                     globfldtmp3(bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_pack_yz.F 284 CTRL_SET_PACK_YZ
     &                    globfldyz(bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_yz.F 289 CTRL_SET_PACK_YZ
     &                       globfldyz(bi,ip,j,bj,jp,k) *
./pkg/ctrl/ctrl_set_pack_yz.F 291 CTRL_SET_PACK_YZ
     &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_yz.F 297 CTRL_SET_PACK_YZ
     &                       globfldyz(bi,ip,j,bj,jp,k) /
./pkg/ctrl/ctrl_set_pack_yz.F 299 CTRL_SET_PACK_YZ
     &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_pack_yz.F 305 CTRL_SET_PACK_YZ
                     globfldtmp2(bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_pack_yz.F 308 CTRL_SET_PACK_YZ
                     cbuff(cbuffindex) = globfldyz(bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_pack_yz.F 54 CTRL_SET_PACK_YZ
      integer ip,jp
./pkg/ctrl/ctrl_set_pack_yz.F 98 CTRL_SET_PACK_YZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xy.F 102 CTRL_SET_UNPACK_XY
             globfld3d  (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 103 CTRL_SET_UNPACK_XY
             globmsk    (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 104 CTRL_SET_UNPACK_XY
             globfldtmp2(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 105 CTRL_SET_UNPACK_XY
             globfldtmp3(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 154 CTRL_SET_UNPACK_XY
        do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xy.F 160 CTRL_SET_UNPACK_XY
              weightfld2d(i,bi,ip,j,bj,jp) = 1. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 210 CTRL_SET_UNPACK_XY
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xy.F 216 CTRL_SET_UNPACK_XY
                  if ( globmsk(i,bi,ip,j,bj,jp,1) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_xy.F 218 CTRL_SET_UNPACK_XY
                     globfld3d(i,bi,ip,j,bj,jp,k) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xy.F 220 CTRL_SET_UNPACK_XY
                     globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xy.F 224 CTRL_SET_UNPACK_XY
                     globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 227 CTRL_SET_UNPACK_XY
                     globfldtmp2(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xy.F 233 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 234 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 236 CTRL_SET_UNPACK_XY
     &                       / sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_unpack_xy.F 240 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 241 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 243 CTRL_SET_UNPACK_XY
     &                       * sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_unpack_xy.F 249 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 250 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 253 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 254 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 259 CTRL_SET_UNPACK_XY
                     globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 262 CTRL_SET_UNPACK_XY
                  globfldtmp3(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xy.F 263 CTRL_SET_UNPACK_XY
     &                 globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 320 CTRL_SET_UNPACK_XY
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xy.F 326 CTRL_SET_UNPACK_XY
                  if ( globmsk(i,bi,ip,j,bj,jp,k) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_xy.F 328 CTRL_SET_UNPACK_XY
                     globfld3d(i,bi,ip,j,bj,jp,k) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xy.F 330 CTRL_SET_UNPACK_XY
                     globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xy.F 334 CTRL_SET_UNPACK_XY
                     globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 337 CTRL_SET_UNPACK_XY
                     globfldtmp2(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xy.F 343 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 344 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 346 CTRL_SET_UNPACK_XY
     &                       / sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_unpack_xy.F 350 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 351 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 353 CTRL_SET_UNPACK_XY
     &                       * sqrt(weightfld2d(i,bi,ip,j,bj,jp))
./pkg/ctrl/ctrl_set_unpack_xy.F 359 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 360 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 363 CTRL_SET_UNPACK_XY
                        globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xy.F 364 CTRL_SET_UNPACK_XY
     &                       globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 369 CTRL_SET_UNPACK_XY
                     globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xy.F 372 CTRL_SET_UNPACK_XY
                  globfldtmp3(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xy.F 373 CTRL_SET_UNPACK_XY
     &                 globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xy.F 406 CTRL_SET_UNPACK_XY
      integer ip,jp
./pkg/ctrl/ctrl_set_unpack_xy.F 49 CTRL_SET_UNPACK_XY
      integer ip,jp
./pkg/ctrl/ctrl_set_unpack_xy.F 96 CTRL_SET_UNPACK_XY
       do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xyz.F 106 CTRL_SET_UNPACK_XYZ
       do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xyz.F 112 CTRL_SET_UNPACK_XYZ
             globfld3d  (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xyz.F 113 CTRL_SET_UNPACK_XYZ
             globmsk    (i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xyz.F 114 CTRL_SET_UNPACK_XYZ
             globfldtmp2(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xyz.F 115 CTRL_SET_UNPACK_XYZ
             globfldtmp3(i,bi,ip,j,bj,jp)   = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xyz.F 172 CTRL_SET_UNPACK_XYZ
        do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xyz.F 178 CTRL_SET_UNPACK_XYZ
              weightfld3d(i,bi,ip,j,bj,jp,k) = 1. _d 0
./pkg/ctrl/ctrl_set_unpack_xyz.F 233 CTRL_SET_UNPACK_XYZ
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xyz.F 239 CTRL_SET_UNPACK_XYZ
                  if ( globmsk(i,bi,ip,j,bj,jp,k) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_xyz.F 241 CTRL_SET_UNPACK_XYZ
                     globfld3d(i,bi,ip,j,bj,jp,k) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xyz.F 243 CTRL_SET_UNPACK_XYZ
                     globfldtmp2(i,bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xyz.F 247 CTRL_SET_UNPACK_XYZ
                     globfld3d(i,bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xyz.F 250 CTRL_SET_UNPACK_XYZ
                     globfldtmp2(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xyz.F 257 CTRL_SET_UNPACK_XYZ
                        globfld3d(i,bi,ip,j,bj,jp,k) = delZnorm
./pkg/ctrl/ctrl_set_unpack_xyz.F 258 CTRL_SET_UNPACK_XYZ
     &                       * globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xyz.F 260 CTRL_SET_UNPACK_XYZ
     &                       / sqrt(weightfld3d(i,bi,ip,j,bj,jp,k))
./pkg/ctrl/ctrl_set_unpack_xyz.F 265 CTRL_SET_UNPACK_XYZ
                        globfld3d(i,bi,ip,j,bj,jp,k) = 1/delZnorm
./pkg/ctrl/ctrl_set_unpack_xyz.F 266 CTRL_SET_UNPACK_XYZ
     &                       * globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xyz.F 268 CTRL_SET_UNPACK_XYZ
     &                       * sqrt(weightfld3d(i,bi,ip,j,bj,jp,k))
./pkg/ctrl/ctrl_set_unpack_xyz.F 276 CTRL_SET_UNPACK_XYZ
                     globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xyz.F 279 CTRL_SET_UNPACK_XYZ
                  globfldtmp3(i,bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xyz.F 280 CTRL_SET_UNPACK_XYZ
     &                 globfld3d(i,bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xyz.F 313 CTRL_SET_UNPACK_XYZ
      integer ip,jp
./pkg/ctrl/ctrl_set_unpack_xyz.F 52 CTRL_SET_UNPACK_XYZ
      integer ip,jp
./pkg/ctrl/ctrl_set_unpack_xz.F 103 CTRL_SET_UNPACK_XZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xz.F 108 CTRL_SET_UNPACK_XZ
                        globfldxz  (i,bi,ip,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xz.F 109 CTRL_SET_UNPACK_XZ
                        globfldtmp2(i,bi,ip,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_unpack_xz.F 110 CTRL_SET_UNPACK_XZ
                        globfldtmp3(i,bi,ip,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_unpack_xz.F 112 CTRL_SET_UNPACK_XZ
                           globmskxz(i,bi,ip,bj,jp,k,iobcs) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xz.F 122 CTRL_SET_UNPACK_XZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xz.F 128 CTRL_SET_UNPACK_XZ
                           globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xz.F 225 CTRL_SET_UNPACK_XZ
               do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xz.F 230 CTRL_SET_UNPACK_XZ
                    if ( globmskxz(i,bi,ip,bj,jp,k,iobcs) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_xz.F 232 CTRL_SET_UNPACK_XZ
                       globfld3d(i,bi,ip,jj,bj,jp,kk) =
./pkg/ctrl/ctrl_set_unpack_xz.F 235 CTRL_SET_UNPACK_XZ
                       globfldtmp2(i,bi,ip,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xz.F 239 CTRL_SET_UNPACK_XZ
                       globfld3d(i,bi,ip,jj,bj,jp,kk) =
./pkg/ctrl/ctrl_set_unpack_xz.F 240 CTRL_SET_UNPACK_XZ
     &                      globfld3d(i,bi,ip,jj,bj,jp,kk)/
./pkg/ctrl/ctrl_set_unpack_xz.F 242 CTRL_SET_UNPACK_XZ
     &                      sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_unpack_xz.F 248 CTRL_SET_UNPACK_XZ
                       globfld3d(i,bi,ip,jj,bj,jp,kk) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xz.F 251 CTRL_SET_UNPACK_XZ
                    globfldtmp3(i,bi,ip,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xz.F 252 CTRL_SET_UNPACK_XZ
     &                   globfld3d(i,bi,ip,jj,bj,jp,kk)
./pkg/ctrl/ctrl_set_unpack_xz.F 310 CTRL_SET_UNPACK_XZ
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_xz.F 315 CTRL_SET_UNPACK_XZ
                  if ( globmskxz(i,bi,ip,bj,jp,k,iobcs) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_xz.F 317 CTRL_SET_UNPACK_XZ
                     globfldxz(i,bi,ip,bj,jp,k) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xz.F 319 CTRL_SET_UNPACK_XZ
                     globfldtmp2(i,bi,ip,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_xz.F 322 CTRL_SET_UNPACK_XZ
                     globfldxz(i,bi,ip,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_xz.F 323 CTRL_SET_UNPACK_XZ
     &                    globfldxz(i,bi,ip,bj,jp,k)/
./pkg/ctrl/ctrl_set_unpack_xz.F 325 CTRL_SET_UNPACK_XZ
     &                    sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_unpack_xz.F 331 CTRL_SET_UNPACK_XZ
                     globfldxz(i,bi,ip,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_xz.F 334 CTRL_SET_UNPACK_XZ
                  globfldtmp3(i,bi,ip,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_xz.F 335 CTRL_SET_UNPACK_XZ
     &                 globfldxz(i,bi,ip,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_xz.F 55 CTRL_SET_UNPACK_XZ
      integer ip,jp
./pkg/ctrl/ctrl_set_unpack_yz.F 103 CTRL_SET_UNPACK_YZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_yz.F 108 CTRL_SET_UNPACK_YZ
                        globfldyz  (bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_yz.F 109 CTRL_SET_UNPACK_YZ
                        globfldtmp2(bi,ip,j,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_unpack_yz.F 110 CTRL_SET_UNPACK_YZ
                        globfldtmp3(bi,ip,j,bj,jp)   = 0.
./pkg/ctrl/ctrl_set_unpack_yz.F 112 CTRL_SET_UNPACK_YZ
                           globmskyz(bi,ip,j,bj,jp,k,iobcs) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_yz.F 122 CTRL_SET_UNPACK_YZ
         do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_yz.F 128 CTRL_SET_UNPACK_YZ
                           globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_yz.F 222 CTRL_SET_UNPACK_YZ
               do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_yz.F 230 CTRL_SET_UNPACK_YZ
                    if ( globmskyz(bi,ip,j,bj,jp,k,iobcs) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_yz.F 232 CTRL_SET_UNPACK_YZ
                       globfld3d(ii,bi,ip,jj,bj,jp,kk) =
./pkg/ctrl/ctrl_set_unpack_yz.F 235 CTRL_SET_UNPACK_YZ
                       globfldtmp2(bi,ip,jj,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_yz.F 239 CTRL_SET_UNPACK_YZ
                       globfld3d(ii,bi,ip,jj,bj,jp,kk) =
./pkg/ctrl/ctrl_set_unpack_yz.F 240 CTRL_SET_UNPACK_YZ
     &                      globfld3d(ii,bi,ip,jj,bj,jp,kk)/
./pkg/ctrl/ctrl_set_unpack_yz.F 242 CTRL_SET_UNPACK_YZ
     &                      sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_unpack_yz.F 248 CTRL_SET_UNPACK_YZ
                       globfld3d(ii,bi,ip,jj,bj,jp,kk) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_yz.F 251 CTRL_SET_UNPACK_YZ
                    globfldtmp3(bi,ip,jj,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_yz.F 252 CTRL_SET_UNPACK_YZ
     &                   globfld3d(ii,bi,ip,jj,bj,jp,kk)
./pkg/ctrl/ctrl_set_unpack_yz.F 309 CTRL_SET_UNPACK_YZ
            do jp = 1,nPy
./pkg/ctrl/ctrl_set_unpack_yz.F 314 CTRL_SET_UNPACK_YZ
                  if ( globmskyz(bi,ip,j,bj,jp,k,iobcs) .ne. 0. ) then
./pkg/ctrl/ctrl_set_unpack_yz.F 316 CTRL_SET_UNPACK_YZ
                     globfldyz(bi,ip,j,bj,jp,k) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_yz.F 318 CTRL_SET_UNPACK_YZ
                     globfldtmp2(bi,ip,j,bj,jp) = cbuff(cbuffindex)
./pkg/ctrl/ctrl_set_unpack_yz.F 321 CTRL_SET_UNPACK_YZ
                     globfldyz(bi,ip,j,bj,jp,k) =
./pkg/ctrl/ctrl_set_unpack_yz.F 322 CTRL_SET_UNPACK_YZ
     &                    globfldyz(bi,ip,j,bj,jp,k)/
./pkg/ctrl/ctrl_set_unpack_yz.F 324 CTRL_SET_UNPACK_YZ
     &                    sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
./pkg/ctrl/ctrl_set_unpack_yz.F 330 CTRL_SET_UNPACK_YZ
                     globfldyz(bi,ip,j,bj,jp,k) = 0. _d 0
./pkg/ctrl/ctrl_set_unpack_yz.F 333 CTRL_SET_UNPACK_YZ
                  globfldtmp3(bi,ip,j,bj,jp) =
./pkg/ctrl/ctrl_set_unpack_yz.F 334 CTRL_SET_UNPACK_YZ
     &                 globfldyz(bi,ip,j,bj,jp,k)
./pkg/ctrl/ctrl_set_unpack_yz.F 55 CTRL_SET_UNPACK_YZ
      integer ip,jp
./lsopt/hessupd.F 47 HESSUPD
      integer jfin, i, j, jp
./lsopt/hessupd.F 70 HESSUPD
         jp = j
./lsopt/hessupd.F 71 HESSUPD
         if (jp.gt.mupd) jp = jp-mupd
./lsopt/hessupd.F 72 HESSUPD
         call DOSTORE( nn, xdiff, .false., 2*jp+3 )
./lsopt/hessupd.F 74 HESSUPD
         call DOSTORE( nn, xdiff, .false., 2*jp+2 )
./lsopt/hessupd.F 75 HESSUPD
         alpha(jp) = r
./lsopt/hessupd.F 95 HESSUPD
         jp = j
./lsopt/hessupd.F 96 HESSUPD
         if (jp .gt. mupd) jp = jp-mupd
./lsopt/hessupd.F 97 HESSUPD
         call DOSTORE( nn, xdiff, .false., 2*jp+2 )
./lsopt/hessupd.F 98 HESSUPD
         r = alpha(jp) - DDOT( nn, dd,1 , xdiff, 1 )
./lsopt/hessupd.F 99 HESSUPD
         call DOSTORE( nn, xdiff, .false., 2*jp+3 )