| 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 |