525 |
vv(i)=v(i)/100. !inter_B.f works in meters |
vv(i)=v(i)/100. !inter_B.f works in meters |
526 |
enddo |
enddo |
527 |
c inter_B: coordinates in m, B field in Tesla |
c inter_B: coordinates in m, B field in Tesla |
528 |
|
c print*,'GUFLD: v ',v |
529 |
call inter_B(vv(1),vv(2),vv(3),ff) |
call inter_B(vv(1),vv(2),vv(3),ff) |
530 |
do i=1,3 !change back the field in kGauss |
do i=1,3 !change back the field in kGauss |
531 |
f(i)=ff(i)*10. |
f(i)=ff(i)*10. |
532 |
enddo |
enddo |
533 |
|
c print*,'GUFLD: b ',f |
534 |
|
c print*,vv,ff |
535 |
|
|
536 |
return |
return |
537 |
end |
end |