#include    <stdio.h>

main(int argc, char *argv[])
{   int nc=0, mc, i, k, lat, lon;
    int io[24] = { 0, 1980, 1980, 2160, 1980, 2160, 2160, 2160, 2160,
                    2160, 2400, 2640, 2640, 2640, 1560, 1560, 1560,
                    1560, 2100, 2100, 2100, 1560, 1560, 1560 } ;
    int ko[24] = { 0, 7770, 7860, 7930, 8010, 8060, 8160, 8230, 8310,
                    8390, 8450, 8415, 8535, 8655, 8520, 7650, 7440,
                    7860, 8400, 8400, 8400, 7650, 7440, 7860 } ;
    double glat, glon, latd, latf, lond, lonf, x, y;
    char out[80];
    opnpinc(argc, argv);
    premsg("\nLat. Long. to XY coordinates Conversion\n");
    while ((nc<1) || (nc>23)) {
        gparmi("Select DATUM/ellipsoid: WGS (1) or OldTokyo (0)  ==> ", &mc);
        if (mc == 1) mc = 200; else mc = 0;
        gparmi("Enter coordinates number (1-23) ==> ", &nc);
    }
    i = io[nc]; k = ko[nc]; glat = (double)i; glon = (double)k;
    sprintf(out, "    origin : %d:%02d'N, %d:%02d'E", i/60,i%60,k/60,k%60);
    premsg(out);
    if (mc == 0) printf("\nxycal:  Old"); else printf("\nxycal:  New");
    printf(" XY-%d  Org.(%d:%02d'N, %d:%02d'E)\n", nc,i/60,i%60,k/60,k%60);
    cvinit(mc, glat, glon);
    while (1) {
        latd = 99999.; latf = 0.; lond = 0.; lonf = 0.;
        prompt("\n  Latitude ?  (deg. min.)\n");
        gparmd2("   < Hit RETURN to Quit > ==> ", &latd,&latf);
        if (latd == 99999.) { fputc('\n',stderr); break; }
        gparmd2("  Longitude ? (deg. min.) ==> ", &lond,&lonf);
        latf = latd*60.+latf; lonf = lond*60.+lonf;
        cviken(latf, lonf, &y, &x);
        lat = (int)(latf/60.); latf = latf - lat*60.;
        lon = (int)(lonf/60.); lonf = lonf - lon*60.;
        printf("  %4d:%7.3f\'N  %4d:%7.3f\'E  :  x = %8.3f  y =%8.3f  km\n",
                lat,latf,lon,lonf,x,y);
    }
    clspin();
}