DCL:MATH1:MAPLIB : 地図投影変換:メソッド(元サブルーチン)の説明
ボンヌ図法の投影, 逆投影を行なう.
CALL MPFBON(XLON,YLAT,X,Y)
CALL MPIBON(X,Y,XLON,YLAT)
CALL MPSBON(YLAT1)
xlon,ylat (R) 経度, 緯度. x,y (R) x, y座標. ylat1 (R) 標準緯線の緯度 
x = r sin(kλ)
y = -r cos(kλ)
   
              
   
ここで, r, k は
    
k  =   sin(θ)/r
r   =  
  θ-θ1 
   + tan(θ1)  
              
 
  
である (θは余緯度, θ1は標準緯線の余緯度).