Map Projection (XYCONV)

NAME
    xyconv - Lat.Lon. to Japanese Transverse Mercator (Old Tokyo Datum) X,Y
    nxycnv - Lat.Lon. to Japanese Transverse Mercator (New JGD2000 Datum) X,Y
    utm    - Lat.Lon. to Universal Transverse Mercator X,Y
    ikconv - Japanese Transverse Mercator (Old Tokyo Datum) X,Y to Lat.Lon.
    nikcnv - Japanese Transverse Mercator (New JGD2000 Datum) X,Y to Lat.Lon.
    utmik  - Universal Transverse Mercator X,Y to Lat.Lon.
                    [i: Ido = Latitude ;  k: Keido = Longitude]
 (Note) For Universal Transverse Mercator (UTM), X=0. on the equator,
       Y=500.(km) on the central meridian, while X=Y=0. at the origin for
       Japanese Transverse Mercator (JTM).

SYNOPSIS

      call xyconv(gi, gk, fi, fk, x, y)
      call nxycnv(gi, gk, fi, fk, x, y)
      call utm(nc, fi, fk, x, y)
  gi, gk  [float] Lat. and Lon. (in minutes) of origin
  nc      [int]   Zone number of UTM coordinate (1-60) for WGS geodetic system
                  (GRS80 ellipsoid) [= New JGD2000 Datum], or
                  plus 800 (801-860) for Old Tokyo Datum (Bessel's ellipsoid)
  fi, fk  [float] Lat. and Lon. (in minutes) of the point to be calculated
  x,  y   [float] Calculated results. Northword(X) and Eastword(Y) coordinates
                  values (in km) are given.

      call ikconv(gi, gk, x, y, fi, fk)
      call nikcnv(gi, gk, x, y, fi, fk)
      call utmik(nc, x, y, fi, fk)
  gi, gk  [float] Lat. and Lon. (in minutes) of origin
  nc      [int]   Zone number of UTM coordinate (1-60) for WGS geodetic system
                  or plus 800 (801-860) for Old Tokyo Datum
  x,  y   [float] Northword(X) and Eastword(Y) coordinates (in km)
                  of the point to be calculated.
  fi, fk  [float] Calculated results. Lat. and Lon. (in minutes) are given.

NAME
    cvinit / cvdinit - Specify projection method
    cviken / cvdiken - Convert Lat.Lon. into plane rectangular coordinates
    cvenik / cvdenik - Convert plane rectangular coordinates into Lat.Lon.

    Two models of Earth's ellipsoid are implemented:
     1. GRS80 ellipsoid :    a = 6378.137 km ,    flattening = 1/298.257222
     2. Bessel's ellipsoid : a = 6377.397155 km , flattening = 1/299.152813

      Following projection numbers are available for WGS geodetic system
    (GRS80 ellipsoid) [= New JGD2000 Datum], while the numbers for Bessel's
    ellipsoid [Old Tokyo datum] are plus 800 for each projection.
         0    Japanese Transverse Mercator [JTM]
      1-60    Universal Transverse Mercator [UTM] (with zone number)
        61    Universal Polar Stereographic Projection [UPS] North Pole
        62    Universal Polar Stereographic Projection [UPS] South Pole
        65    Same as UTM but non-standard Central Meridian
        70    Mercator projection
        71    Lambert Conformal Conic projection with 1 Standard Pararrel
        72    Lambert Conformal Conic projection with 2 Standard Pararrels
       100    Lambert Azimuthal Equal Area projection from a sphere
               with the same surface area as the ellipsoid
       109    Lambert Azimuthal Equal Area projection from a sphere
               with the same eqatorial radius as the ellipsoid
       199    Latitude/Longitude in minutes are regarded as distance in km

  (Note) Although X=Y=0. at the origin in general,
       for UTM, X=0. on the equator, and Y=500.(km) on the central meridian,
       for UPS, X=Y=2000.(km) at the pole.

SYNOPSIS

      call cvinit(nc, gi, gk, sp1, sp2)
      call cvdinit(nc, dgi, dgk, dp1, dp2)
  nc        [int]     Projection number
  gi,  gk   [float]   Lat. and Lon. (in minutes) of origin (X=Y=0.)
                       (not refered in cases of UTM/UPS)
  dgi, dgk  [double]   [ditto]
  sp1, sp2  [float]   Lat. (in minutes) of two standard pararells
                       (not refered in cases other than projection 72,872)
  dp1, dp2  [double]   [ditto]

      call cviken(fi, fk, ye, xn)
      call cvdiken(dfi, dfk, dye, dxn)
  fi,  fk   [float]   Lat. and Lon. (in minutes) of the point to be calculated
  dfi, dfk  [double]   [ditto]
  ye,  xn   [float]   Calculated results. Eastword(Y) and Northword(X)
                      coordinates values (in km) are given.
  dye, dxn  [double]   [ditto]

      call cvenik(ye, xn, fi, fk)
      call cvdenik(dye, dxn, dfi, dfk)
  ye,  xn   [float]   Eastword(Y) and Northword(X) coordinates (in km)
                      of the point to be calculated.
  dye, dxn  [double]   [ditto]
  fi,  fk   [float]   Calculated results. Lat. and Lon. (in minutes) are given.
  dfi, dfk  [double]   [ditto]

 (Note)  Be aware that the sequence of coordinates values (X, Y) arguments
       in  cviken / cvenik / cvdiken / cvdenik is different from those
       in  utm / utmik / xyconv / ikconv / nxycnv / nikcnv .