オープンとクローズド領域の境界の検定
subroutine zsub33(ia,aa)
説明：ainte1aの結果と用いて、地球磁気圏の３次元磁力線の中から、
オープン領域とクローズド領域の間の狭い領域に存在する磁力線を選び出して
、その選択された磁力線だけ描くように検定結果をainte1に渡す。
入力変数：
      nxg=ia(1)：磁力線の出発点のx方向のデータ数
      nyg=ia(2)：磁力線の出発点のy方向のデータ数
関連するサブルーチン：
　　　ainte1a(ia,aa,f,p)
　　　ainte1(ia,aa,f,p)