2018-10 Update for AMSS3 (GSJ Open-file Report no.648)

Program dpam/pltrk.f90 :

  The line 85 of the source
    call gparmi('Select Coordinate Number (+200 for WGS) ==> ', ncc)
  should be read as
    call gparmi('Select Coordinate Number (+800 for Old-TOKYO) ==> ', ncc)

Program dpam/dmaq2asc.c :

  As there are many lines to be changed, download the updated version source
 dmaq2asc20180818.c, and replace the source with it.

Program anam/plamagc.f90 :

  a) The line 37 of the source
          nc = ncc; ncc = nc + 800
    should be read as
          nc = ncc; ncc = nc + 800;  mc = 0
  b) The line 39 of the source
          ncc = ncc - 200; nc = ncc
    should be read as
          ncc = ncc - 200; nc = ncc;  mc = 1
  c) The existing 2 lines 43-44 of the source
        nc = ncc
        if ((ncc >= 800) .and. (ncc < 1000)) nc = ncc - 800
    should be read as following 4 lines
        nc = ncc;  mc = 1
        if ((ncc >= 800) .and. (ncc < 1000)) then
          nc = ncc - 800;  mc = 0
        endif

Program utils/hgeoidcal.c :

  a) The line 16 of the source
      double latd, latf, latm, lond, lonf, lonm, gi, gk, hg;
    should be read as
      double latd, latf, latm, lond, lonf, lonm, gi, gk, fido, fkdo, hg;
  b) The line 35 of the source
        gi = latd*60.+latf; gk = lond*60.+lonf;
    should be read as
        gi = latm = latd*60.+latf; gk = lonm = lond*60.+lonf;
  c) The lines 38-39 of the source
        ido = (int)(gi/60.); latm = gi - ido*60.;
        kdo = (int)(gk/60.); lonm = gk - kdo*60.;
    should be read as
        ido = (int)(gi/60.); fido = gi - ido*60.;
        kdo = (int)(gk/60.); fkdo = gk - kdo*60.;
  d) The lines 41-42 of the source
               ido, latm, cns, kdo, lonm, cew);
        if (hgeoid(latf, lonf, &hg) < 0) {
    should be read as
               ido, fido, cns, kdo, fkdo, cew);
        if (hgeoid(latm, lonm, &hg) < 0) {