#include    <stdio.h>

main(int argc, char *argv[])
{   int mc, nc, lat, lon; double latd, latf, lond, lonf, x, y;
    opnpinc(argc, argv);
    premsg("\nLat. Long. to UTM x,y coordinates Conversion\n");
    while (1) {
        gparmi("Select DATUM/ellipsoid: WGS (1) or OldTokyo (0)  ==> ", &mc);
        if (mc == 1) mc = 200; else mc = 0;
        gparmi("Enter Long.E of central meridian  ==> ", &lon);
        nc = (lon+183)/6; 
        if ((nc>0) && (nc<=60) && ((nc*6-183) == lon))
            { cvinit(mc+nc); break; }
    }
    if (mc == 0) printf("\nutmcal:  Old"); else printf("\nutmcal:  New");
    printf(" nc = %d (Central meridian = %d)\n", nc, lon);
    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 = %7.3f  km\n",
                lat,latf,lon,lonf,x,y);
    }
    clspin();
}