mflxd_cylind

    円柱型永久磁石が生成する磁束密度Bを計算する倍精度サブルーチン．但し，永久磁石の中心対称軸をz軸とする円柱座標系(r, θ, z)を用い，磁化ベクトルはM=Mezを満たすと仮定する．
     
     
  • サブルーチン名： mflxd_cylind

    •  
  • 呼び出し方：call mflxd_cylind(n_gl,gzi,weit,rpm,zpm,r,z,br,bz)

    •  
  • 引　数：
    • 引数
    意　味
    入出力の区分
    備　考
    n_gl Gauss-Legendre積分の積分点数
    入力
    		 整数型
    gzi(n_gl) Gauss-Legendre積分の積分点の座標を格納した配列
    入力
    		 倍精度実数型
    weit(n_gl) Gauss-Legendre積分の重みを格納した配列
    入力
    		 倍精度実数型
    rpm 永久磁石の半径
    入力
    		 倍精度実数型
    zpm(2) 永久磁石のz座標の下限と上限．zpm(1) < zpm(2)を満たす必要あり．
    入力
    		 倍精度実数型
    r,z 観測点の(r,z)座標
    入力
    		 倍精度実数型
    br, bz 磁束密度のr, z成分
    出力
    		 倍精度実数型
  • 注　意
    1. 本ルーチンはB/(μ0M)を出力する．
    2. 本ルーチンを呼び出す前には，Gauss-Legendre積分の積分点と重みを計算するために，予めdawg2を呼び出しておく必要がある．