磁力線を予め描いて判定
subroutine ainte1a(ia,aa,f,p)
説明：３次元のＭＨＤシミュレーションデータf(i1)から３成分の磁場データを用いて磁力線を追跡し、
open, closed, detached 磁力線の判定を行う
入力変数：
      nx=ia(1)：磁力線の出発点のx方向のデータ数
      ny=ia(2)：磁力線の出発点のy方向のデータ数
      mx=ia(3)：３次元データf(i1)のx方向のデータ数
      my=ia(4)：３次元データf(i1)のy方向のデータ数
      mz=ia(5)：３次元データf(i1)のz方向のデータ数
      nxz=ia(6)：地球の中心を調節するパラメータ
      mi=ia(7)：入力変数の最初の位置
      mo=ia(8)：出力変数の最初の位置
      th0=aa(1)*pi/180.0：磁力線の出発点の最低緯度
      aru=aa(2)：内側境界の遷移を特徴づけるパラメータ
      ar1=aa(3)：内側境界の遷移領域の外端半径
      arb=aa(4)：内側境界の遷移領域の内端半径
      xxl=aa(5)：x方向のシミュレーションボックスの長さ
      yyl=aa(6)：y方向のシミュレーションボックスの長さ
      zzl=aa(7)：z方向のシミュレーションボックスの長さ
      b0=aa(8)：ダイポール磁場の大きさ（通常 1.0）
      gx0=aa(13)：グラフを描く原点のx座標
      gy0=aa(14)：グラフを描く原点のy座標
      gxl=aa(15)：グラフを描くx方向の長さ
      gyl=aa(16)：グラフを描くy方向の長さ
      gth=aa(17)*pi/180.0：グラフの座標の回転角度
      gxmi=aa(18)：グラフを描くx方向の最小値
      gxma=aa(19)：グラフを描くx方向の最大値
      ep1=aa(20)：判定の最小閾値
使用している主なサブルーチン：
　　　quant1(x,y,z,hx,hy,hz,1,mx,my,mz,nxp,aa,f,q)：
      　３次元格子点のＭＨＤデータf(i1)から任意の場所(x,y,z)
      　のＭＨＤデータq(8)を計算