vtf-logo

fsi/sfc-amroc/TubeCJBurnFlaps/src/lset3.f

c
c     Copyright (C) 2003-2007 California Institute of Technology
c     Ralf Deiterding, ralf@amroc.net
c
c     Set sign behind entry
c
c     =====================================================
      subroutine ls(maxmx,maxmy,maxmz,mbc,mx,my,mz,x,y,z,
     &     dx,dy,dz,phi,t)
c     =====================================================
      implicit double precision (a-h,o-z)
c     
      include  "cuser.i"
c     
      dimension phi(1-mbc:maxmx+mbc, 1-mbc:maxmy+mbc, 
     &     1-mbc:maxmz+mbc)
      dimension x(1-mbc:maxmx+mbc),y(1-mbc:maxmy+mbc),
     &     z(1-mbc:maxmz+mbc)
c     
      do 60 k = 1-mbc, mz+mbc
         if (z(k).gt.zt0) goto 60
         do 50 i = 1-mbc, mx+mbc
            do 50 j = 1-mbc, my+mbc
               if (dsqrt(x(i)**2+y(j)**2).gt.rf)
     &              phi(i,j,k) = -phi(i,j,k)
 50      continue
 60   continue
c     
      return
      end
c