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