Return to mvefil.f CVS log

Up to [Research Unix] / researchv10no / cmd / view2d / cray

Annotation of researchv10no/cmd/view2d/cray/mvefil.f, revision 1.1

1.1     ! root        1:       subroutine mvefil (unit,time,nx,ny,fmin,fmax,outsid,f)
        !             2: *
        !             3: *     mvefil writes the array f(nx,ny) in a format for view2d movies.
        !             4: *     input:
        !             5: *          unit   = integer variable specifying output unit
        !             6: *          time   = real variable specifying time
        !             7: *          nx     = integer variable specifying number of x points
        !             8: *          ny     = integer variable specifying number of y points
        !             9: *          fmin   = real variable specifying minimum of f(i,j)
        !            10: *          fmax   = real variable specifying maximum of f(i,j)
        !            11: *               If fmin=fmax on input, the range will be set inside.
        !            12: *          outsid  = real variable specifying points outside region
        !            13: *                   if f(i,j) <= outsid => point is outside region
        !            14: *          f      = real array of points to be plotted
        !            15: *                   f(i,j)  i = 1,...,nx  j = 1,...,ny
        !            16: *   suggested jcl:
        !            17: * DISPOSE,DN=FT30,MF=XH,DC=ST,TEXT='|/b2/c1127/ehg/rks astro/swift file'.    
        !            18: *
        !            19: *  Each data point is represented in the packed ASCII character format
        !            20: *     by a 12 bit integer between 0 and 4095.  
        !            21: *             f(i,j)=fmin+(fmax-fmin)*icode(i,j)/4095
        !            22: *     icode is packed as two 6 bit integers ( / hi / lo / ).  Each 6 bits
        !            23: *     is represented as an printable ASCII character, 6 bits + ! (033 octal).
        !            24: *     A point which lies outside the computation region is flaged with (~~).
        !            25: *     written by R. Kent Smith,   May 1984
        !            26: *       to supply programs written by Andrew Hume and Eric Grosse
        !            27:       integer bits,shift,lbuf
        !            28:       parameter (bits=64,lbuf=128)
        !            29:       parameter (shift=2**(bits-16))
        !            30:       integer high,low,zero
        !            31:       data high,low,zero,otside/7700b,77b,2r!!,2r~~/
        !            32:       integer hi,lo,buf(lbuf)
        !            33: *
        !            34:       integer unit,nx,ny
        !            35:       real time
        !            36:       real fmin,fmax,f(1)
        !            37: *
        !            38:       if (unit.le.0) return
        !            39: *
        !            40:       nxny=nx*ny
        !            41:       if (fmin.ge.fmax) then
        !            42: *
        !            43:          fmin=1e500
        !            44:          fmax=-fmin
        !            45:          do 110 i=1,nxny
        !            46:             c=f(i)
        !            47:             if(c.gt.outsid) then
        !            48:                fmin=amin1(c,fmin)
        !            49:                fmax=amax1(c,fmax)
        !            50:             end if
        !            51:   110       continue
        !            52: *
        !            53:          end if
        !            54: *
        !            55:       write (unit,6000) time,nx,ny,fmin,fmax
        !            56:  6000 format(1x,e12.6,2i6,1x,e12.6,1x,e12.6)
        !            57: *
        !            58:       if (fmax.gt.fmin) then
        !            59:          slope=4095.0/(fmax-fmin)
        !            60:       else
        !            61:          slope=0.0
        !            62:       end if
        !            63:       nseg=(nxny-1)/lbuf
        !            64:       num=nxny-nseg*lbuf
        !            65:       lf=0
        !            66:       do 220 iseg=1,nseg+1
        !            67: *
        !            68:          do 210 i=1,num
        !            69: *
        !            70:             code=slope*(f(lf+i)-fmin)
        !            71:             index=int(code+0.5)
        !            72: *
        !            73:             lo=low.and.index
        !            74:             hi=high.and.index
        !            75:             ibuf=4*hi+lo+zero
        !            76:             buf(i)=cvmgp(ibuf,otside,index)*shift
        !            77: *
        !            78:   210       continue
        !            79: *
        !            80:          write (unit,6100) (buf(i),i=1,num)
        !            81:  6100    format(32a2)
        !            82: *
        !            83:          lf=lf+num
        !            84:          num=lbuf
        !            85:   220    continue
        !            86: *
        !            87:       return
        !            88:       end

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.