ANAM: 空中磁気探査データ各種解析プログラム

1.　１つのファイルは，１組のグリッドデータまたは複数組のグリッドデータからな
　る．但し，複数組のグリッドデータからなるファイルの２組目以降のグリッドデー
　タには，利用上の制約があり，次の場合にのみ有効とする．
　　(1) 第１組目のグリッドデータがある曲面上での物理量の分布を示し，第２組の
　　　グリッドデータがその曲面の高度分布を示す．
　　(2) 複数組のグリッドデータの羅列（一般に処理の中間結果など）であり，その
　　　値の分布の図化以外には，利用を見込まない．
2.　格子点データの値に使用する単位は，磁力値なら nＴ，高度なら ｍ，磁化強度
　なら 0.01Ａ/ｍ，重力値なら mＧal とし，その他の場合もこれらから導かれるも
　のを使用することを原則とする．また，格子点データの未定義を示す数値として
　は，有効数字が 99999 となる正の値（実データ値の変域に応じて桁移動する）を
　用いることを推奨する．
3.　グリッドデータファイル上のすべてのレコード（行）の長さは，改行コードを除
　いて 80バイト以下とする．また，可能な限り改行コードを除いて 79バイト以下に
　することを推奨する．
4.　グリッドデータファイル上のデータには，2バイトコード文字（全角文字）を使用
　せず，また，改行コード以外の制御コードも含めないものとする．

各組のグリッドデータは，下記のように構成される．
a) コメント ： 第１ヘッダーの前に任意行数のコメント行をおくことができる．
        コメント行は，１行 80バイト以下で，先頭桁を '#' とする．
b) 第１ヘッダー（地域名，使用座標系情報） [固定書式] 
    FORMAT(a8, i4,4x, 2i8, 2i8)
  area: 地域名（もしくはデータ種別）を示す 8文字の文字列（先頭は '#' 以外）
  nc :  座標系番号：地図投影座標変換サブルーチンで使用する図法番号
            0 ： 国土調査法に基づく新平面直角座標（ガウスクルィーゲル図法）
         1-60 ： ＵＴＭ座標（ガウスクルィーゲル図法，ゾーン番号の 1-60 に対応）
           61 ： 北極ＵＰＳ座標（Universal Polar Stereographic Projection）
           62 ： 南極ＵＰＳ座標（   ＝ 国際極心平射図法                  ）
           65 ： UTM座標と同様であるが非標準の中央経線を指定できる
           70 ： メルカトール図法
           71 ： ランベルト等角円錐図法（１基本緯線）
           72 ： ランベルト等角円錐図法（２基本緯線）
          100 ： ランベルト等積方位図法（楕円体と表面積の等しい球からの投影）
          109 ： ランベルト等積方位図法（楕円体の長半径を半径とする球から投影）
          199 ： 緯経度の分単位の値をそのまま km 単位の北・東向き座標値とみなす
        この番号は，ベッセル楕円体(東京測地系)に対するものであり，
        GRS楕円体(WGS, 世界測地系)に対しては 200 を加算したものとする．
  ig, kg : 座標原点の緯度・経度（分単位） [1≦nc≦62 の場合は不要] 
  i1, i2 : ２つの基本緯線の緯度（分単位） [nc＝72 の場合にのみ必要] 
    　各図法とも一般に，座標原点で座標値は X＝Y＝0 であるが，
    ＵＴＭ座標では，基準点で X＝0，Y＝500,000（ｍ），
    ＵＰＳ座標では，極で X＝Y＝2,000,000（ｍ）である．
c) 第２ヘッダー（格子位置情報と未定義値，高度） [自由書式] 
    FORMAT(2i12, 2i6, 2i6, 1x, f7.1, 1x, f7.0) を標準とする．
  ixs, iys : データ格子の南西隅の座標値（北向き(X)・東向き(Y)，ｍ単位） [整数] 
  mszx,mszy: 北方向・東方向の格子間隔（ｍ単位） [整数] 
  mxn, myn : 北方向・東方向の格子点数（両端を含む） [整数] 
  vnul :     格子点データの値が未定義であることを示すために用いる数値 [実数] 
  alt :      観測高度値（ｍ単位） [実数] 
             （値が 0. のときは，観測面の高度分布が第2組のグリッドデータとして
               与えられること示し，値が 負(-1.)のときは，高度が未定義であると
               みなされる．第2組の高度分布データに対しては無意味であるが，-1. 
               の値を埋めておくものとする．）
d) 格子点データ [自由書式] 
    FORMAT((f7.1, 9(1x,f7.1))) を標準とする．
  全格子点データを，南西隅からはじめて北方向に向う順序（北端まで到達して始めて
  すぐ東側の列のデータに移る）で，列記する．
  このデータを読みとるプログラムは，
      read(10,*) ((f(i,k),i=1,mxn),k=1,myn)
  のごとくであるが，書き出すプログラムは，
      do 1 k=1,myn
        write(10,'((f7.1,9(1x,f7.1)))') (f(i,k),i=1,mxn)
    1 continue
  のように，列を移る位置で改行しておくことを推奨する．

# Areaname: Kobe-Kyoto
# Survey Date: 1995.12.07-12.27
&A-01
 2079.0222N  8116.2764E   277.8m   -45.1nT
 2079.0405N  8116.3164E   278.5m   -44.6nT
 2079.0588N  8116.3564E   279.1m   -44.4nT
        .......
 2087.3958N  8134.2559E   275.4m   -48.3nT
 2087.4158N  8134.2964E   275.4m   -53.6nT
&    C-2r
 2088.2712N  8134.3799E   279.1m   -44.9nT
 2088.2563N  8134.3384E   279.3m   -40.8nT
 2088.2407N  8134.2964E   279.6m   -40.1nT
        .......
     先頭桁が，'#' の行は，単なるコメント行であり，通常，ファイルの先頭にのみ
       置かれる．（一連の個別点データ(測線データ)の途中に，置いてはならない．）
     先頭桁が，'&' または '%' の行は測線の開始を示す．
       第2〜9桁に，測線名を記述する．それ以外の記述内容とその形式に制約は
       設けない．（開始時刻や含まれる測点数などを記述する場合が多い．）
       ファイルの末尾に，含まれる測点数が 0 の測線を表示して，データの末尾を
       明示する場合がある．
     それ以外は，個別の点データであり，緯度(分単位)・経度(分単位)・高度(m単位)
       ・残差磁力値(nT単位)の値が，
          format(f10.4, 1hN, f11.4, 1hE, f8.1, 1hm, f8.1, 2hnT)
       の形式（計42桁）で入っている．