real x(3),Bout(3)
      REAL*8 TH,PH
      TH=0
      PH=0
      
      call trgInit('trg_field_map_extended.dat',tH,pH)
      do i=0,0
         do iz=-100,100
            x(1) =i
            x(2) = i
            x(3) =iz-182.5
            call GUFLD (x,Bout)
            write(*,*)x(2),IZ,Bout,
     ,           SQRT(Bout(1)**2+Bout(2)**2+Bout(3)**2)
         enddo
      enddo
      end