ANAM: 磁気異常データ解析処理プログラム
rev.2022-03-02

English

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

各組のグリッドデータは，下記のように構成される．
a) コメント ： 第１ヘッダーの前に任意行数のコメント行をおくことができる．
        コメント行は，１行 80バイト以下で，先頭桁を '#' とする．
b) 第１ヘッダー（地域名，使用座標系情報） [固定書式] 
    FORMAT(a8, 4x,i4, 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, 世界測地系)に対するものであり，
        ベッセル楕円体(旧東京測地系)に対しては 800 を加算したものとする．
  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 k=1,myn
        write(10,'((f7.1,9(1x,f7.1)))') (f(i,k),i=1,mxn)
      enddo
  のように，列を移る位置で改行しておくことを推奨する．

# Areaname: Kobe-Kyoto
# Survey Date: 1995.12.07-12.27
&A-01
  2079.02221N  8116.27649E   277.87m   -45.15nT
  2079.04052N  8116.31640E   278.58m   -44.66nT
  2079.05883N  8116.35641E   279.19m   -44.47nT
        .......
  2087.39584N  8134.25592E   275.40m   -48.38nT
  2087.41585N  8134.29643E   275.41m   -53.69nT
&    C-2r
  2088.27126N  8134.37994E   279.12m   -44.90nT
  2088.25637N  8134.33845E   279.33m   -40.81nT
  2088.24078N  8134.29646E   279.64m   -40.12nT
        .......

     先頭桁が，'#' の行は，単なるコメント行であり，通常，ファイルの先頭にのみ
       置かれる．（一連の個別点データ(測線データ)の途中に，置いてはならない．）
     先頭桁が，'&' または '%' の行は測線の開始を示す．
       第2～9桁に，測線名を記述する．それ以外の記述内容とその形式に制約は
       設けない．（開始時刻や含まれる測点数などを記述する場合が多い．）
       ファイルの末尾に，含まれる測点数が 0 の測線を表示して，データの末尾を
       明示する場合がある．
     それ以外は，個別の点データであり，緯度(分単位)・経度(分単位)・高度(m単位)
       ・残差磁力値(nT単位)の値が，
          format(1x, f11.5, 1hN,1x, f11.5, 1hE,1x, f8.2, 1hm,1x, f8.2, 2hnT)
       の形式（計47桁）で入っている．
     北緯を意味する 'N'，東経を意味する 'E' の文字を含むが，南半球・西半球の
       場合には 緯度値・経度値をマイナス値で表現する．
     入力データの場合は数値部分の桁数は任意である．