/[MITgcm]/MITgcm/pkg/flt/flt_traj.F
ViewVC logotype

Diff of /MITgcm/pkg/flt/flt_traj.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by jmc, Thu Feb 5 01:01:48 2009 UTC revision 1.7 by jmc, Tue Feb 10 21:30:21 2009 UTC
# Line 31  C     == routine arguments == Line 31  C     == routine arguments ==
31    
32  C     === Functions ==  C     === Functions ==
33        INTEGER ILNBLNK        INTEGER ILNBLNK
       _RL global2local_i  
       _RL global2local_j  
34        _RL FLT_MAP_K2R        _RL FLT_MAP_K2R
35        EXTERNAL ILNBLNK        EXTERNAL ILNBLNK
36        EXTERNAL FLT_MAP_K2R        EXTERNAL FLT_MAP_K2R
# Line 41  C     == local variables == Line 39  C     == local variables ==
39        INTEGER bi, bj, imax        INTEGER bi, bj, imax
40        PARAMETER (imax=13)        PARAMETER (imax=13)
41        INTEGER ip, kp, ii        INTEGER ip, kp, ii
42        _RL ix, jy, i0x, j0y, zz        _RL ix, jy, i0x, j0y, xx, yy, zz
43        _RL uu, vv, tt, ss, pp        _RL uu, vv, tt, ss, pp
44    
45        INTEGER ioUnit, irecord        INTEGER ioUnit, irecord
# Line 107  C total number of floats Line 105  C total number of floats
105           j0y = DFLOAT( myYGlobalLo-1 + (bj-1)*sNy )           j0y = DFLOAT( myYGlobalLo-1 + (bj-1)*sNy )
106           DO ip=1,npart_tile(bi,bj)           DO ip=1,npart_tile(bi,bj)
107    
108              ix=global2local_i(xpart(ip,bi,bj),bi,bj,myThid)              ix = ipart(ip,bi,bj)
109              jy=global2local_j(ypart(ip,bi,bj),bi,bj,myThid)              jy = jpart(ip,bi,bj)
110              zz = FLT_MAP_K2R( kpart(ip,bi,bj),bi,bj,myThid)              CALL FLT_MAP_IJLOCAL2XY( xx, yy,
111         I                               ix, jy, bi,bj, myThid )
112                zz = FLT_MAP_K2R( kpart(ip,bi,bj),bi,bj,myThid )
113              kp = NINT(kpart(ip,bi,bj))              kp = NINT(kpart(ip,bi,bj))
114              tmp(1)  = npart(ip,bi,bj)              tmp(1)  = npart(ip,bi,bj)
115              tmp(2)  = myTime              tmp(2)  = myTime
116              tmp(3)  = xpart(ip,bi,bj)              tmp(3)  = xx
117              tmp(4)  = ypart(ip,bi,bj)              tmp(4)  = yy
118              tmp(5)  = zz              tmp(5)  = zz
119              tmp(6)  = ix + i0x              tmp(6)  = ix + i0x
120              tmp(7)  = jy + j0y              tmp(7)  = jy + j0y
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22