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