| 521 |
real*8 vv(3),ff(3) !inter_B.f works in double precision |
real*8 vv(3),ff(3) !inter_B.f works in double precision |
| 522 |
|
|
| 523 |
|
|
| 524 |
do i=1,3 |
do i=1,3 |
| 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 |
call inter_B(vv(1),vv(2),vv(3),ff) |
c$$$ print*,'GUFLD: v ',v |
| 529 |
do i=1,3 !change back the field in kGauss |
call inter_B(vv(1),vv(2),vv(3),ff) |
| 530 |
f(i)=ff(i)*10. |
do i=1,3 !change back the field in kGauss |
| 531 |
enddo |
f(i)=ff(i)*10. |
| 532 |
|
enddo |
| 533 |
|
c$$$ print*,'GUFLD: b ',f |
| 534 |
|
|
| 535 |
return |
return |
| 536 |
end |
end |
| 537 |
|
|