/* amedas_sub.f -- translated by f2c (version 20100827).
   You must link the resulting object file with libf2c:
      on Microsoft Windows system, link with libf2c.lib;
      on Linux or Unix systems, link with .../path/to/libf2c.a -lm
      or, if you install libf2c.a in a standard place, with -lf2c -lm
      -- in that order, at the end of the command line, as in
              cc *.o -lf2c -lm
      Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,

              http://www.netlib.org/f2c/libf2c.zip
*/

#include <stdlib.h>
#include "f2c.h"

int amdhr (int ihx, int idx, int nx, int *ista, int iyy, int imm,
	   int imiss, float rmiss, int *itype, int *ir, int *id, int *iw,
	   float *s, float *t){
  int i;
  long ihx_long, idx_long, nx_long, iyy_long, imm_long, imiss_long;
  long *ista_long, *itype_long, *ir_long, *id_long, *iw_long;
  ista_long = malloc(nx*sizeof(long));
  itype_long = malloc(idx*nx*sizeof(long));
  ir_long = malloc(ihx*idx*nx*sizeof(long));
  id_long = malloc(ihx*idx*nx*sizeof(long));
  iw_long = malloc(ihx*idx*nx*sizeof(long));
  ihx_long = ihx;
  idx_long = idx;
  nx_long = nx;
  iyy_long = iyy;
  imm_long = imm;
  imiss_long = imiss;
  for (i=0; i<=nx-1; i++){
    ista_long[i] = ista[i];
  }
  amdhr_(&ihx_long, &idx_long, &nx_long, ista_long, &iyy_long, &imm_long,
	 &imiss_long, &rmiss, itype_long,
	 ir_long, id_long, iw_long, s, t);
  for (i=0; i<=idx*nx-1; i++){
    itype[i] = itype_long[i];
  }
  for (i=0; i<=ihx*idx*nx-1; i++){
    ir[i] = ir_long[i];
    id[i] = id_long[i];
    iw[i] = iw_long[i];
  }
  free (ista_long);
  free (itype_long);
  free (ir_long);
  free (id_long);
  free (iw_long);
  return 0;
}

int amdid (int nx, int *ista, float *xx, float *yy, float *hh){
  int i;
  long nx_long;
  long *ista_long;
  ista_long = malloc(nx*sizeof(long));
  nx_long = nx;
  amdid_(&nx_long, ista_long, xx, yy, hh);
  for (i=0; i<=nx-1; i++){
    ista[i] = ista_long[i];
  }
  free (ista_long);
  return 0;
}

/* Table of constant values */

static integer c__1 = 1;
static integer c__2 = 2;
static integer c__1000 = 1000;
static integer c__3 = 3;
static integer c__4 = 4;
static integer c__5 = 5;
static integer c__7 = 7;
static integer c__6 = 6;
static integer c__9 = 9;
static integer c__11 = 11;
static real c_b137 = .1f;
static integer c__8 = 8;
static integer c__15 = 15;
static integer c__16 = 16;
static integer c__13 = 13;

/* Subroutine */ int amdhr_(integer *ihx, integer *idx, integer *nx, integer *
	ista, integer *iyy, integer *imm, integer *imiss, real *rmiss, 
	integer *itype, integer *ir, integer *id, integer *iw, real *s, real *
	t)
{
    /* System generated locals */
    integer itype_dim1, itype_offset, ir_dim1, ir_dim2, ir_offset, id_dim1, 
	    id_dim2, id_offset, iw_dim1, iw_dim2, iw_offset, s_dim1, s_dim2, 
	    s_offset, t_dim1, t_dim2, t_offset;

    /* Local variables */
    extern /* Subroutine */ int amdhr1_(integer *, integer *, integer *, 
	    integer *, integer *, integer *, integer *, real *, integer *, 
	    integer *, integer *, integer *, real *, real *), amdhr2_(integer 
	    *, integer *, integer *, integer *, integer *, integer *, integer 
	    *, real *, integer *, integer *, integer *, integer *, real *, 
	    real *), amdhr3_(integer *, integer *, integer *, integer *, 
	    integer *, integer *, integer *, real *, integer *, integer *, 
	    integer *, integer *, real *, real *), amdhr4_(integer *, integer 
	    *, integer *, integer *, integer *, integer *, integer *, real *, 
	    integer *, integer *, integer *, integer *, real *, real *);

    /* Parameter adjustments */
    t_dim1 = *ihx;
    t_dim2 = *idx;
    t_offset = 1 + t_dim1 * (1 + t_dim2);
    t -= t_offset;
    s_dim1 = *ihx;
    s_dim2 = *idx;
    s_offset = 1 + s_dim1 * (1 + s_dim2);
    s -= s_offset;
    iw_dim1 = *ihx;
    iw_dim2 = *idx;
    iw_offset = 1 + iw_dim1 * (1 + iw_dim2);
    iw -= iw_offset;
    id_dim1 = *ihx;
    id_dim2 = *idx;
    id_offset = 1 + id_dim1 * (1 + id_dim2);
    id -= id_offset;
    ir_dim1 = *ihx;
    ir_dim2 = *idx;
    ir_offset = 1 + ir_dim1 * (1 + ir_dim2);
    ir -= ir_offset;
    itype_dim1 = *idx;
    itype_offset = 1 + itype_dim1;
    itype -= itype_offset;
    --ista;

    /* Function Body */
    if (*iyy >= 1975 && *iyy <= 2000) {
	amdhr1_(ihx, idx, nx, &ista[1], iyy, imm, imiss, rmiss, &itype[
		itype_offset], &ir[ir_offset], &id[id_offset], &iw[iw_offset],
		 &s[s_offset], &t[t_offset]);
    }
    if (*iyy >= 2001 && *iyy <= 2003) {
	amdhr2_(ihx, idx, nx, &ista[1], iyy, imm, imiss, rmiss, &itype[
		itype_offset], &ir[ir_offset], &id[id_offset], &iw[iw_offset],
		 &s[s_offset], &t[t_offset]);
    }
    if (*iyy >= 2004 && *iyy <= 2009) {
	amdhr3_(ihx, idx, nx, &ista[1], iyy, imm, imiss, rmiss, &itype[
		itype_offset], &ir[ir_offset], &id[id_offset], &iw[iw_offset],
		 &s[s_offset], &t[t_offset]);
    }
/*      IF ((IYY.GE.2010).AND.(IYY.LE.2010)) THEN */
    if (*iyy >= 2010) {
	amdhr4_(ihx, idx, nx, &ista[1], iyy, imm, imiss, rmiss, &itype[
		itype_offset], &ir[ir_offset], &id[id_offset], &iw[iw_offset],
		 &s[s_offset], &t[t_offset]);
    }
    return 0;
} /* amdhr_ */

/* ======================================== */
/* Subroutine */ int amdhr1_(integer *ihx, integer *idx, integer *nx, integer 
	*ista, integer *iyy, integer *imm, integer *imiss, real *rmiss, 
	integer *itype, integer *ir, integer *id, integer *iw, real *s, real *
	t)
{
    /* System generated locals */
    shortint s__1;
    integer itype_dim1, itype_offset, ir_dim1, ir_dim2, ir_offset, id_dim1, 
	    id_dim2, id_offset, iw_dim1, iw_dim2, iw_offset, s_dim1, s_dim2, 
	    s_offset, t_dim1, t_dim2, t_offset, i__1, i__2, i__3;
    char ch__1[80];
    olist o__1;
    cllist cl__1;

    /* Builtin functions */
    integer f_open(olist *), s_rdue(cilist *), do_uio(integer *, char *, 
	    ftnlen), e_rdue(void), f_clos(cllist *);

    /* Local variables */
    static integer i__, n, nn;
    static shortint ir2[24], is2[24], it2[24], iw2[24];
    static integer idd, ihh;
    static shortint idd2, imm2;
    extern shortint idv2_(integer *, shortint *), inv2_(shortint *);
    extern integer inv4_(integer *);
    static shortint iyy2;
    extern /* Character */ VOID amdf1_(char *, ftnlen, integer *, integer *);
    static shortint irec2;
    static integer ista1, ista4;
    static shortint itype2;

    /* Fortran I/O blocks */
    static cilist io___5 = { 1, 10, 0, 0, 0 };


    /* Parameter adjustments */
    t_dim1 = *ihx;
    t_dim2 = *idx;
    t_offset = 1 + t_dim1 * (1 + t_dim2);
    t -= t_offset;
    s_dim1 = *ihx;
    s_dim2 = *idx;
    s_offset = 1 + s_dim1 * (1 + s_dim2);
    s -= s_offset;
    iw_dim1 = *ihx;
    iw_dim2 = *idx;
    iw_offset = 1 + iw_dim1 * (1 + iw_dim2);
    iw -= iw_offset;
    id_dim1 = *ihx;
    id_dim2 = *idx;
    id_offset = 1 + id_dim1 * (1 + id_dim2);
    id -= id_offset;
    ir_dim1 = *ihx;
    ir_dim2 = *idx;
    ir_offset = 1 + ir_dim1 * (1 + ir_dim2);
    ir -= ir_offset;
    itype_dim1 = *idx;
    itype_offset = 1 + itype_dim1;
    itype -= itype_offset;
    --ista;

    /* Function Body */
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	i__2 = *idx;
	for (idd = 1; idd <= i__2; ++idd) {
	    itype[idd + n * itype_dim1] = 0;
	    i__3 = *ihx;
	    for (ihh = 1; ihh <= i__3; ++ihh) {
		ir[ihh + (idd + n * ir_dim2) * ir_dim1] = *imiss;
		id[ihh + (idd + n * id_dim2) * id_dim1] = *imiss;
		iw[ihh + (idd + n * iw_dim2) * iw_dim1] = *imiss;
		s[ihh + (idd + n * s_dim2) * s_dim1] = *rmiss;
		t[ihh + (idd + n * t_dim2) * t_dim1] = *rmiss;
/* L13: */
	    }
/* L12: */
	}
/* L11: */
    }
    o__1.oerr = 1;
    o__1.ounit = 10;
    o__1.ofnmlen = 80;
    amdf1_(ch__1, (ftnlen)80, iyy, imm);
    o__1.ofnm = ch__1;
    o__1.orl = 256;
    o__1.osta = "OLD";
    o__1.oacc = "DIRECT";
    o__1.ofm = "UNFORMATTED";
    o__1.oblnk = 0;
    i__1 = f_open(&o__1);
    if (i__1 != 0) {
	goto L29;
    }
    i__ = 1;
L21:
    io___5.cirec = i__;
    i__1 = s_rdue(&io___5);
    if (i__1 != 0) {
	goto L29;
    }
    i__1 = do_uio(&c__1, (char *)&irec2, (ftnlen)sizeof(shortint));
    if (i__1 != 0) {
	goto L29;
    }
    i__1 = do_uio(&c__1, (char *)&ista4, (ftnlen)sizeof(integer));
    if (i__1 != 0) {
	goto L29;
    }
    i__1 = do_uio(&c__1, (char *)&itype2, (ftnlen)sizeof(shortint));
    if (i__1 != 0) {
	goto L29;
    }
    i__1 = do_uio(&c__1, (char *)&iyy2, (ftnlen)sizeof(shortint));
    if (i__1 != 0) {
	goto L29;
    }
    i__1 = do_uio(&c__1, (char *)&imm2, (ftnlen)sizeof(shortint));
    if (i__1 != 0) {
	goto L29;
    }
    i__1 = do_uio(&c__1, (char *)&idd2, (ftnlen)sizeof(shortint));
    if (i__1 != 0) {
	goto L29;
    }
    for (ihh = 1; ihh <= 24; ++ihh) {
	i__1 = do_uio(&c__1, (char *)&ir2[ihh - 1], (ftnlen)sizeof(shortint));
	if (i__1 != 0) {
	    goto L29;
	}
	i__1 = do_uio(&c__1, (char *)&iw2[ihh - 1], (ftnlen)sizeof(shortint));
	if (i__1 != 0) {
	    goto L29;
	}
	i__1 = do_uio(&c__1, (char *)&is2[ihh - 1], (ftnlen)sizeof(shortint));
	if (i__1 != 0) {
	    goto L29;
	}
	i__1 = do_uio(&c__1, (char *)&it2[ihh - 1], (ftnlen)sizeof(shortint));
	if (i__1 != 0) {
	    goto L29;
	}
    }
    i__1 = e_rdue();
    if (i__1 != 0) {
	goto L29;
    }
    ista1 = inv4_(&ista4);
    n = 0;
    i__1 = *nx;
    for (nn = 1; nn <= i__1; ++nn) {
	if (ista1 == ista[nn]) {
	    n = nn;
	    goto L39;
	}
/* L31: */
    }
L39:
    if (n != 0) {
	idd = inv2_(&idd2);
	itype[idd + n * itype_dim1] = inv2_(&itype2);
	if (itype[idd + n * itype_dim1] == 1 || itype[idd + n * itype_dim1] ==
		 4) {
	    i__1 = *ihx;
	    for (ihh = 1; ihh <= i__1; ++ihh) {
		if (inv2_(&ir2[ihh - 1]) != 32767) {
		    ir[ihh + (idd + n * ir_dim2) * ir_dim1] = inv2_(&ir2[ihh 
			    - 1]);
		}
/* L41: */
	    }
	}
	if (itype[idd + n * itype_dim1] == 4) {
	    i__1 = *ihx;
	    for (ihh = 1; ihh <= i__1; ++ihh) {
		if (inv2_(&iw2[ihh - 1]) != -1) {
		    s__1 = inv2_(&iw2[ihh - 1]);
		    id[ihh + (idd + n * id_dim2) * id_dim1] = idv2_(&c__2, &
			    s__1);
		    s__1 = inv2_(&iw2[ihh - 1]);
		    iw[ihh + (idd + n * iw_dim2) * iw_dim1] = idv2_(&c__1, &
			    s__1);
		}
		if (inv2_(&is2[ihh - 1]) != 32767) {
		    s[ihh + (idd + n * s_dim2) * s_dim1] = (real) inv2_(&is2[
			    ihh - 1]) * .1f;
		}
		if (inv2_(&it2[ihh - 1]) != 32767) {
		    t[ihh + (idd + n * t_dim2) * t_dim1] = (real) inv2_(&it2[
			    ihh - 1]) * .1f;
		}
/* L42: */
	    }
	}
    }
    ++i__;
    goto L21;
L29:
    cl__1.cerr = 0;
    cl__1.cunit = 10;
    cl__1.csta = 0;
    f_clos(&cl__1);
    return 0;
} /* amdhr1_ */

/* Subroutine */ int amdhr2_(integer *ihx, integer *idx, integer *nx, integer 
	*ista, integer *iyy, integer *imm, integer *imiss, real *rmiss, 
	integer *itype, integer *ir, integer *id, integer *iw, real *s, real *
	t)
{
    /* System generated locals */
    integer itype_dim1, itype_offset, ir_dim1, ir_dim2, ir_offset, id_dim1, 
	    id_dim2, id_offset, iw_dim1, iw_dim2, iw_offset, s_dim1, s_dim2, 
	    s_offset, t_dim1, t_dim2, t_offset, i__1, i__2, i__3;
    char ch__1[80];
    olist o__1;
    cllist cl__1;

    /* Builtin functions */
    integer f_open(olist *), s_rsfe(cilist *), e_rsfe(void), do_fio(integer *,
	     char *, ftnlen), i_nint(real *), f_clos(cllist *), s_wsle(cilist 
	    *), do_lio(integer *, integer *, char *, ftnlen), e_wsle(void);
    /* Subroutine */ int s_stop(char *, ftnlen);

    /* Local variables */
    static integer n;
    static char a1[1*1000];
    static integer ii, nd1;
    static real rd1;
    static integer nr1, ns1, nt1;
    static real rr1;
    static integer nw1;
    static real rs1, rt1, rw1;
    static integer idd, ndd, ihh;
    static real rdd;
    static integer nhh;
    static real rhh;
    extern /* Subroutine */ int amdcr_(integer *, char *, integer *, real *, 
	    integer *, ftnlen);
    extern /* Character */ VOID amdf2h_(char *, ftnlen, integer *, integer *, 
	    integer *);

    /* Fortran I/O blocks */
    static cilist io___21 = { 1, 10, 1, "()", 0 };
    static cilist io___22 = { 1, 10, 1, "()", 0 };
    static cilist io___23 = { 1, 10, 1, "()", 0 };
    static cilist io___24 = { 1, 10, 1, "()", 0 };
    static cilist io___25 = { 1, 10, 1, "(1000A1)", 0 };
    static cilist io___42 = { 0, 6, 0, 0, 0 };


    /* Parameter adjustments */
    t_dim1 = *ihx;
    t_dim2 = *idx;
    t_offset = 1 + t_dim1 * (1 + t_dim2);
    t -= t_offset;
    s_dim1 = *ihx;
    s_dim2 = *idx;
    s_offset = 1 + s_dim1 * (1 + s_dim2);
    s -= s_offset;
    iw_dim1 = *ihx;
    iw_dim2 = *idx;
    iw_offset = 1 + iw_dim1 * (1 + iw_dim2);
    iw -= iw_offset;
    id_dim1 = *ihx;
    id_dim2 = *idx;
    id_offset = 1 + id_dim1 * (1 + id_dim2);
    id -= id_offset;
    ir_dim1 = *ihx;
    ir_dim2 = *idx;
    ir_offset = 1 + ir_dim1 * (1 + ir_dim2);
    ir -= ir_offset;
    itype_dim1 = *idx;
    itype_offset = 1 + itype_dim1;
    itype -= itype_offset;
    --ista;

    /* Function Body */
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	i__2 = *idx;
	for (idd = 1; idd <= i__2; ++idd) {
	    itype[idd + n * itype_dim1] = 0;
	    i__3 = *ihx;
	    for (ihh = 1; ihh <= i__3; ++ihh) {
		ir[ihh + (idd + n * ir_dim2) * ir_dim1] = *imiss;
		id[ihh + (idd + n * id_dim2) * id_dim1] = *imiss;
		iw[ihh + (idd + n * iw_dim2) * iw_dim1] = *imiss;
		s[ihh + (idd + n * s_dim2) * s_dim1] = *rmiss;
		t[ihh + (idd + n * t_dim2) * t_dim1] = *rmiss;
/* L13: */
	    }
/* L12: */
	}
/* L11: */
    }
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	o__1.oerr = 1;
	o__1.ounit = 10;
	o__1.ofnmlen = 80;
	amdf2h_(ch__1, (ftnlen)80, iyy, imm, &ista[n]);
	o__1.ofnm = ch__1;
	o__1.orl = 0;
	o__1.osta = "OLD";
	o__1.oacc = 0;
	o__1.ofm = "FORMATTED";
	o__1.oblnk = 0;
	i__2 = f_open(&o__1);
	if (i__2 != 0) {
	    goto L49;
	}
	i__2 = s_rsfe(&io___21);
	if (i__2 != 0) {
	    goto L100001;
	}
	i__2 = e_rsfe();
L100001:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	i__2 = s_rsfe(&io___22);
	if (i__2 != 0) {
	    goto L100002;
	}
	i__2 = e_rsfe();
L100002:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	i__2 = s_rsfe(&io___23);
	if (i__2 != 0) {
	    goto L100003;
	}
	i__2 = e_rsfe();
L100003:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	i__2 = s_rsfe(&io___24);
	if (i__2 != 0) {
	    goto L100004;
	}
	i__2 = e_rsfe();
L100004:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
L41:
	i__2 = s_rsfe(&io___25);
	if (i__2 != 0) {
	    goto L100005;
	}
	for (ii = 1; ii <= 1000; ++ii) {
	    i__2 = do_fio(&c__1, a1 + (ii - 1), (ftnlen)1);
	    if (i__2 != 0) {
		goto L100005;
	    }
	}
	i__2 = e_rsfe();
L100005:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	amdcr_(&c__1000, a1, &c__1, &rdd, &ndd, (ftnlen)1);
	idd = i_nint(&rdd);
	amdcr_(&c__1000, a1, &c__2, &rhh, &nhh, (ftnlen)1);
	ihh = i_nint(&rhh);
	itype[idd + n * itype_dim1] = 4;
	if (idd < 1 || idd > *idx) {
	    goto L41;
	}
	if (ihh < 1 || ihh > *ihx) {
	    goto L41;
	}
	amdcr_(&c__1000, a1, &c__3, &rr1, &nr1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__4, &rd1, &nd1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__5, &rw1, &nw1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__7, &rs1, &ns1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__6, &rt1, &nt1, (ftnlen)1);
	if (nr1 >= 1) {
	    ir[ihh + (idd + n * ir_dim2) * ir_dim1] = i_nint(&rr1);
	}
	if (nd1 >= 1) {
	    id[ihh + (idd + n * id_dim2) * id_dim1] = i_nint(&rd1);
	}
	if (nw1 >= 1) {
	    iw[ihh + (idd + n * iw_dim2) * iw_dim1] = i_nint(&rw1);
	}
	if (ns1 >= 1) {
	    s[ihh + (idd + n * s_dim2) * s_dim1] = rs1;
	}
	if (nt1 >= 1) {
	    t[ihh + (idd + n * t_dim2) * t_dim1] = rt1;
	}
	goto L41;
L49:
	cl__1.cerr = 0;
	cl__1.cunit = 10;
	cl__1.csta = 0;
	f_clos(&cl__1);
/* L21: */
    }
    return 0;
L99:
    cl__1.cerr = 0;
    cl__1.cunit = 10;
    cl__1.csta = 0;
    f_clos(&cl__1);
    s_wsle(&io___42);
    do_lio(&c__9, &c__1, "Can not read file.", (ftnlen)18);
    e_wsle();
    s_stop("", (ftnlen)0);
    return 0;
} /* amdhr2_ */

/* Subroutine */ int amdhr3_(integer *ihx, integer *idx, integer *nx, integer 
	*ista, integer *iyy, integer *imm, integer *imiss, real *rmiss, 
	integer *itype, integer *ir, integer *id, integer *iw, real *s, real *
	t)
{
    /* System generated locals */
    integer itype_dim1, itype_offset, ir_dim1, ir_dim2, ir_offset, id_dim1, 
	    id_dim2, id_offset, iw_dim1, iw_dim2, iw_offset, s_dim1, s_dim2, 
	    s_offset, t_dim1, t_dim2, t_offset, i__1, i__2, i__3;
    char ch__1[80];
    olist o__1;
    cllist cl__1;

    /* Builtin functions */
    integer f_open(olist *), s_rsfe(cilist *), e_rsfe(void), do_fio(integer *,
	     char *, ftnlen), i_nint(real *), f_clos(cllist *), s_wsle(cilist 
	    *), do_lio(integer *, integer *, char *, ftnlen), e_wsle(void);
    /* Subroutine */ int s_stop(char *, ftnlen);

    /* Local variables */
    static integer n;
    static char a1[1*1000];
    static integer ii, nd1;
    static real rd1;
    static integer nr1, ns1, nt1;
    static real rr1;
    static integer nw1;
    static real rs1, rt1, rw1;
    static integer idd, ndd, ihh;
    static real rdd;
    static integer nhh;
    static real rhh;
    extern /* Subroutine */ int amdcr_(integer *, char *, integer *, real *, 
	    integer *, ftnlen);
    extern /* Character */ VOID amdf3h_(char *, ftnlen, integer *, integer *, 
	    integer *);

    /* Fortran I/O blocks */
    static cilist io___46 = { 1, 10, 1, "()", 0 };
    static cilist io___47 = { 1, 10, 1, "()", 0 };
    static cilist io___48 = { 1, 10, 1, "()", 0 };
    static cilist io___49 = { 1, 10, 1, "()", 0 };
    static cilist io___50 = { 1, 10, 1, "(1000A1)", 0 };
    static cilist io___67 = { 0, 6, 0, 0, 0 };


    /* Parameter adjustments */
    t_dim1 = *ihx;
    t_dim2 = *idx;
    t_offset = 1 + t_dim1 * (1 + t_dim2);
    t -= t_offset;
    s_dim1 = *ihx;
    s_dim2 = *idx;
    s_offset = 1 + s_dim1 * (1 + s_dim2);
    s -= s_offset;
    iw_dim1 = *ihx;
    iw_dim2 = *idx;
    iw_offset = 1 + iw_dim1 * (1 + iw_dim2);
    iw -= iw_offset;
    id_dim1 = *ihx;
    id_dim2 = *idx;
    id_offset = 1 + id_dim1 * (1 + id_dim2);
    id -= id_offset;
    ir_dim1 = *ihx;
    ir_dim2 = *idx;
    ir_offset = 1 + ir_dim1 * (1 + ir_dim2);
    ir -= ir_offset;
    itype_dim1 = *idx;
    itype_offset = 1 + itype_dim1;
    itype -= itype_offset;
    --ista;

    /* Function Body */
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	i__2 = *idx;
	for (idd = 1; idd <= i__2; ++idd) {
	    itype[idd + n * itype_dim1] = 0;
	    i__3 = *ihx;
	    for (ihh = 1; ihh <= i__3; ++ihh) {
		ir[ihh + (idd + n * ir_dim2) * ir_dim1] = *imiss;
		id[ihh + (idd + n * id_dim2) * id_dim1] = *imiss;
		iw[ihh + (idd + n * iw_dim2) * iw_dim1] = *imiss;
		s[ihh + (idd + n * s_dim2) * s_dim1] = *rmiss;
		t[ihh + (idd + n * t_dim2) * t_dim1] = *rmiss;
/* L13: */
	    }
/* L12: */
	}
/* L11: */
    }
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	o__1.oerr = 1;
	o__1.ounit = 10;
	o__1.ofnmlen = 80;
	amdf3h_(ch__1, (ftnlen)80, iyy, imm, &ista[n]);
	o__1.ofnm = ch__1;
	o__1.orl = 0;
	o__1.osta = "OLD";
	o__1.oacc = 0;
	o__1.ofm = "FORMATTED";
	o__1.oblnk = 0;
	i__2 = f_open(&o__1);
	if (i__2 != 0) {
	    goto L49;
	}
	i__2 = s_rsfe(&io___46);
	if (i__2 != 0) {
	    goto L100006;
	}
	i__2 = e_rsfe();
L100006:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	i__2 = s_rsfe(&io___47);
	if (i__2 != 0) {
	    goto L100007;
	}
	i__2 = e_rsfe();
L100007:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	i__2 = s_rsfe(&io___48);
	if (i__2 != 0) {
	    goto L100008;
	}
	i__2 = e_rsfe();
L100008:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	i__2 = s_rsfe(&io___49);
	if (i__2 != 0) {
	    goto L100009;
	}
	i__2 = e_rsfe();
L100009:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
L41:
	i__2 = s_rsfe(&io___50);
	if (i__2 != 0) {
	    goto L100010;
	}
	for (ii = 1; ii <= 1000; ++ii) {
	    i__2 = do_fio(&c__1, a1 + (ii - 1), (ftnlen)1);
	    if (i__2 != 0) {
		goto L100010;
	    }
	}
	i__2 = e_rsfe();
L100010:
	if (i__2 < 0) {
	    goto L49;
	}
	if (i__2 > 0) {
	    goto L99;
	}
	amdcr_(&c__1000, a1, &c__1, &rdd, &ndd, (ftnlen)1);
	idd = i_nint(&rdd);
	amdcr_(&c__1000, a1, &c__2, &rhh, &nhh, (ftnlen)1);
	ihh = i_nint(&rhh);
	itype[idd + n * itype_dim1] = 4;
	if (idd < 1 || idd > *idx) {
	    goto L41;
	}
	if (ihh < 1 || ihh > *ihx) {
	    goto L41;
	}
	amdcr_(&c__1000, a1, &c__3, &rr1, &nr1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__5, &rd1, &nd1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__7, &rw1, &nw1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__11, &rs1, &ns1, (ftnlen)1);
	amdcr_(&c__1000, a1, &c__9, &rt1, &nt1, (ftnlen)1);
	if (nr1 >= 1) {
	    ir[ihh + (idd + n * ir_dim2) * ir_dim1] = i_nint(&rr1);
	}
	if (nd1 >= 1) {
	    id[ihh + (idd + n * id_dim2) * id_dim1] = i_nint(&rd1);
	}
	if (nw1 >= 1) {
	    iw[ihh + (idd + n * iw_dim2) * iw_dim1] = i_nint(&rw1);
	}
	if (ns1 >= 1) {
	    s[ihh + (idd + n * s_dim2) * s_dim1] = rs1;
	}
	if (nt1 >= 1) {
	    t[ihh + (idd + n * t_dim2) * t_dim1] = rt1;
	}
	goto L41;
L49:
	cl__1.cerr = 0;
	cl__1.cunit = 10;
	cl__1.csta = 0;
	f_clos(&cl__1);
/* L21: */
    }
    return 0;
L99:
    cl__1.cerr = 0;
    cl__1.cunit = 10;
    cl__1.csta = 0;
    f_clos(&cl__1);
    s_wsle(&io___67);
    do_lio(&c__9, &c__1, "Can not read file.", (ftnlen)18);
    e_wsle();
    s_stop("", (ftnlen)0);
    return 0;
} /* amdhr3_ */

/* Subroutine */ int amdhr4_(integer *ihx, integer *idx, integer *nx, integer 
	*ista, integer *iyy, integer *imm, integer *imiss, real *rmiss, 
	integer *itype, integer *ir, integer *id, integer *iw, real *s, real *
	t)
{
    /* System generated locals */
    integer itype_dim1, itype_offset, ir_dim1, ir_dim2, ir_offset, id_dim1, 
	    id_dim2, id_offset, iw_dim1, iw_dim2, iw_offset, s_dim1, s_dim2, 
	    s_offset, t_dim1, t_dim2, t_offset, i__1, i__2, i__3, i__4;
    real r__1;
    char ch__1[80];
    olist o__1;
    cllist cl__1;

    /* Builtin functions */
    integer f_open(olist *), s_rdue(cilist *), do_uio(integer *, char *, 
	    ftnlen), e_rdue(void), i_nint(real *), f_clos(cllist *), s_wsle(
	    cilist *), do_lio(integer *, integer *, char *, ftnlen), e_wsle(
	    void);
    /* Subroutine */ int s_stop(char *, ftnlen);

    /* Local variables */
    static integer i__, n;
    static char a1[1*262];
    static integer ii, id1, ir1, is1, it1, iw1, idd, ihh;
    extern /* Character */ VOID amdf4h_(char *, ftnlen, integer *, integer *, 
	    integer *);

    /* Fortran I/O blocks */
    static cilist io___72 = { 1, 10, 0, 0, 0 };
    static cilist io___80 = { 0, 6, 0, 0, 0 };


    /* Parameter adjustments */
    t_dim1 = *ihx;
    t_dim2 = *idx;
    t_offset = 1 + t_dim1 * (1 + t_dim2);
    t -= t_offset;
    s_dim1 = *ihx;
    s_dim2 = *idx;
    s_offset = 1 + s_dim1 * (1 + s_dim2);
    s -= s_offset;
    iw_dim1 = *ihx;
    iw_dim2 = *idx;
    iw_offset = 1 + iw_dim1 * (1 + iw_dim2);
    iw -= iw_offset;
    id_dim1 = *ihx;
    id_dim2 = *idx;
    id_offset = 1 + id_dim1 * (1 + id_dim2);
    id -= id_offset;
    ir_dim1 = *ihx;
    ir_dim2 = *idx;
    ir_offset = 1 + ir_dim1 * (1 + ir_dim2);
    ir -= ir_offset;
    itype_dim1 = *idx;
    itype_offset = 1 + itype_dim1;
    itype -= itype_offset;
    --ista;

    /* Function Body */
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	i__2 = *idx;
	for (idd = 1; idd <= i__2; ++idd) {
	    itype[idd + n * itype_dim1] = 0;
	    i__3 = *ihx;
	    for (ihh = 1; ihh <= i__3; ++ihh) {
		ir[ihh + (idd + n * ir_dim2) * ir_dim1] = *imiss;
		id[ihh + (idd + n * id_dim2) * id_dim1] = *imiss;
		iw[ihh + (idd + n * iw_dim2) * iw_dim1] = *imiss;
		s[ihh + (idd + n * s_dim2) * s_dim1] = *rmiss;
		t[ihh + (idd + n * t_dim2) * t_dim1] = *rmiss;
/* L13: */
	    }
/* L12: */
	}
/* L11: */
    }
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
	o__1.oerr = 1;
	o__1.ounit = 10;
	o__1.ofnmlen = 80;
	amdf4h_(ch__1, (ftnlen)80, iyy, imm, &ista[n]);
	o__1.ofnm = ch__1;
	o__1.orl = 262;
	o__1.osta = "OLD";
	o__1.oacc = "DIRECT";
	o__1.ofm = "UNFORMATTED";
	o__1.oblnk = 0;
	i__2 = f_open(&o__1);
	if (i__2 != 0) {
	    goto L49;
	}
	i__2 = *idx;
	for (idd = 1; idd <= i__2; ++idd) {
	    i__3 = *ihx;
	    for (ihh = 1; ihh <= i__3; ++ihh) {
		i__ = ihh * 6 + (idd - 1) * 144;
		io___72.cirec = i__;
		i__4 = s_rdue(&io___72);
		if (i__4 != 0) {
		    goto L49;
		}
		for (ii = 1; ii <= 262; ++ii) {
		    i__4 = do_uio(&c__1, a1 + (ii - 1), (ftnlen)1);
		    if (i__4 != 0) {
			goto L49;
		    }
		}
		i__4 = e_rdue();
		if (i__4 != 0) {
		    goto L49;
		}
		itype[idd + n * itype_dim1] = 4;
		if (*(unsigned char *)&a1[66] <= 11) {
		    ir1 = (*(unsigned char *)&a1[65] << 24) + (*(unsigned 
			    char *)&a1[64] << 16) + (*(unsigned char *)&a1[63]
			     << 8) + *(unsigned char *)&a1[62];
		    r__1 = (real) ir1 * .1f;
		    ir[ihh + (idd + n * ir_dim2) * ir_dim1] = i_nint(&r__1);
		}
		if (*(unsigned char *)&a1[130] <= 11) {
		    id1 = (*(unsigned char *)&a1[129] << 24) + (*(unsigned 
			    char *)&a1[128] << 16) + (*(unsigned char *)&a1[
			    127] << 8) + *(unsigned char *)&a1[126];
		    id[ihh + (idd + n * id_dim2) * id_dim1] = id1;
		}
		if (*(unsigned char *)&a1[158] <= 11) {
		    iw1 = (*(unsigned char *)&a1[157] << 24) + (*(unsigned 
			    char *)&a1[156] << 16) + (*(unsigned char *)&a1[
			    155] << 8) + *(unsigned char *)&a1[154];
		    r__1 = (real) iw1 * .1f;
		    iw[ihh + (idd + n * iw_dim2) * iw_dim1] = i_nint(&r__1);
		}
		if (*(unsigned char *)&a1[190] <= 11) {
		    it1 = (*(unsigned char *)&a1[189] << 24) + (*(unsigned 
			    char *)&a1[188] << 16) + (*(unsigned char *)&a1[
			    187] << 8) + *(unsigned char *)&a1[186];
		    t[ihh + (idd + n * t_dim2) * t_dim1] = (real) it1 * .1f;
		}
		if (*(unsigned char *)&a1[242] <= 11) {
		    is1 = (*(unsigned char *)&a1[241] << 24) + (*(unsigned 
			    char *)&a1[240] << 16) + (*(unsigned char *)&a1[
			    239] << 8) + *(unsigned char *)&a1[238];
		    s[ihh + (idd + n * s_dim2) * s_dim1] = (real) is1 * .1f;
		}
/* L32: */
	    }
/* L31: */
	}
L49:
	cl__1.cerr = 0;
	cl__1.cunit = 10;
	cl__1.csta = 0;
	f_clos(&cl__1);
/* L21: */
    }
    return 0;
/* L99: */
    cl__1.cerr = 0;
    cl__1.cunit = 10;
    cl__1.csta = 0;
    f_clos(&cl__1);
    s_wsle(&io___80);
    do_lio(&c__9, &c__1, "Can not read file.", (ftnlen)18);
    e_wsle();
    s_stop("", (ftnlen)0);
    return 0;
} /* amdhr4_ */

/* Subroutine */ int amdcr_(integer *ix, char *a, integer *n, real *g, 
	integer *ng, ftnlen a_len)
{
    /* System generated locals */
    integer i__1;

    /* Builtin functions */
    double pow_ri(real *, integer *);

    /* Local variables */
    static integer i__, i1, i2, iadd, ifig, icrd, isign, icheck, icount, 
	    ipoint;

    /* Parameter adjustments */
    --a;

    /* Function Body */
    icrd = 0;
    isign = 1;
    icheck = 0;
    ipoint = 0;
    ifig = 0;
    i1 = 0;
    i2 = 0;
    icount = 0;
    *g = 0.f;
    *ng = 0;
    i__1 = *ix;
    for (i__ = 1; i__ <= i__1; ++i__) {
	if (icount == *n - 1) {
	    i1 = i__;
	    goto L12;
	}
	if (*(unsigned char *)&a[i__] == ',') {
	    ++icount;
	}
/* L11: */
    }
L12:
    icount = 0;
    i__1 = *ix;
    for (i__ = 1; i__ <= i__1; ++i__) {
	if (icount == *n) {
	    i2 = i__ - 2;
	    goto L14;
	}
	if (*(unsigned char *)&a[i__] == ',') {
	    ++icount;
	}
/* L13: */
    }
L14:
    if (i1 > i2) {
	return 0;
    }
    i__1 = i2;
    for (i__ = i1; i__ <= i__1; ++i__) {
	iadd = -1;
	if (*(unsigned char *)&a[i__] >= 48 && *(unsigned char *)&a[i__] <= 
		57) {
	    icheck = 1;
	    iadd = *(unsigned char *)&a[i__] - 48;
	    icrd = icrd * 10 + iadd;
	    if (ipoint == 1) {
		++ifig;
	    }
	}
	if (icheck == 0 && *(unsigned char *)&a[i__] == ' ') {
	    iadd = 0;
	}
	if (icheck == 0 && *(unsigned char *)&a[i__] == '-') {
	    icheck = 1;
	    isign = -1;
	    iadd = 0;
	}
	if (ipoint == 0 && *(unsigned char *)&a[i__] == '.') {
	    icheck = 1;
	    ipoint = 1;
	    iadd = 0;
	}
	if (iadd < 0) {
	    return 0;
	}
/* L21: */
    }
    *g = (real) isign * (real) icrd * pow_ri(&c_b137, &ifig);
    *ng = 1;
    return 0;
} /* amdcr_ */

/* Character */ VOID amdf1_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im)
{
    /* System generated locals */
    address a__1[4], a__2[2], a__3[6], a__4[8];
    integer i__1[4], i__2[2], i__3[6], i__4[8];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cm[2], cy[4];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__1[2] = 1, a__1[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__1[3] = 1, a__1[3] = ch__5;
    s_cat(cy, a__1, i__1, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    s_cat(cm, a__2, i__2, &c__2, (ftnlen)2);
    if (*iy >= 1976 && *iy <= 1978) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1976/AMEDAS/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/AMD";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy >= 1979 && *iy <= 1982) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1979/AMEDAS/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/AMD";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy >= 1983 && *iy <= 1986) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1983/AMEDAS/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/AMD";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy >= 1987 && *iy <= 1990) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1987/AMEDAS/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/AMD";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy >= 1991 && *iy <= 1994) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1991/AMEDAS/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/AMD";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy >= 1995 && *iy <= 1997) {
/* Writing concatenation */
	i__4[0] = 2, a__4[0] = "CD";
	i__4[1] = 4, a__4[1] = cy;
	i__4[2] = 8, a__4[2] = "/AMEDAS/";
	i__4[3] = 4, a__4[3] = cy;
	i__4[4] = 4, a__4[4] = "/AMD";
	i__4[5] = 4, a__4[5] = cy;
	i__4[6] = 1, a__4[6] = ".";
	i__4[7] = 2, a__4[7] = cm;
	s_cat(ret_val, a__4, i__4, &c__8, (ftnlen)80);
    }
    if (*iy == 1998) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1998/amedas/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/amd";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy == 1999) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD1999/amedas/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/Amd";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (*iy == 2000) {
/* Writing concatenation */
	i__3[0] = 14, a__3[0] = "CD2000/amedas/";
	i__3[1] = 4, a__3[1] = cy;
	i__3[2] = 4, a__3[2] = "/amd";
	i__3[3] = 4, a__3[3] = cy;
	i__3[4] = 1, a__3[4] = ".";
	i__3[5] = 2, a__3[5] = cm;
	s_cat(ret_val, a__3, i__3, &c__6, (ftnlen)80);
    }
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__2[0] = 7, a__2[0] = "AMeDAS/";
	i__2[1] = 80, a__2[1] = ret_val;
	s_cat(ret_val, a__2, i__2, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf1_ */

/* Character */ VOID amdf2d_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im, integer *ista)
{
    /* System generated locals */
    address a__1[2], a__2[4], a__3[5], a__4[15];
    integer i__1[2], i__2[4], i__3[5], i__4[15];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1], ch__6[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cm[2], cy2[2], cy4[4], csta2[2], csta5[5];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cy2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__2[2] = 1, a__2[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__2[3] = 1, a__2[3] = ch__5;
    s_cat(cy4, a__2, i__2, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cm, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(csta2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__3[0] = 1, a__3[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__3[1] = 1, a__3[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *ista / 100 % 10 + 48;
    i__3[2] = 1, a__3[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *ista / 10 % 10 + 48;
    i__3[3] = 1, a__3[3] = ch__5;
    *(unsigned char *)&ch__6[0] = *ista % 10 + 48;
    i__3[4] = 1, a__3[4] = ch__6;
    s_cat(csta5, a__3, i__3, &c__5, (ftnlen)5);
/* Writing concatenation */
    i__4[0] = 2, a__4[0] = "CD";
    i__4[1] = 4, a__4[1] = cy4;
    i__4[2] = 7, a__4[2] = "/amedas";
    i__4[3] = 1, a__4[3] = "/";
    i__4[4] = 4, a__4[4] = cy4;
    i__4[5] = 2, a__4[5] = cm;
    i__4[6] = 1, a__4[6] = "/";
    i__4[7] = 4, a__4[7] = "area";
    i__4[8] = 2, a__4[8] = csta2;
    i__4[9] = 2, a__4[9] = "/d";
    i__4[10] = 5, a__4[10] = csta5;
    i__4[11] = 1, a__4[11] = "_";
    i__4[12] = 4, a__4[12] = cy4;
    i__4[13] = 2, a__4[13] = cm;
    i__4[14] = 4, a__4[14] = ".csv";
    s_cat(ret_val, a__4, i__4, &c__15, (ftnlen)80);
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__1[0] = 7, a__1[0] = "AMeDAS/";
	i__1[1] = 80, a__1[1] = ret_val;
	s_cat(ret_val, a__1, i__1, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf2d_ */

/* Character */ VOID amdf2h_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im, integer *ista)
{
    /* System generated locals */
    address a__1[2], a__2[4], a__3[5], a__4[15];
    integer i__1[2], i__2[4], i__3[5], i__4[15];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1], ch__6[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cm[2], cy2[2], cy4[4], csta2[2], csta5[5];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cy2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__2[2] = 1, a__2[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__2[3] = 1, a__2[3] = ch__5;
    s_cat(cy4, a__2, i__2, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cm, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(csta2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__3[0] = 1, a__3[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__3[1] = 1, a__3[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *ista / 100 % 10 + 48;
    i__3[2] = 1, a__3[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *ista / 10 % 10 + 48;
    i__3[3] = 1, a__3[3] = ch__5;
    *(unsigned char *)&ch__6[0] = *ista % 10 + 48;
    i__3[4] = 1, a__3[4] = ch__6;
    s_cat(csta5, a__3, i__3, &c__5, (ftnlen)5);
/* Writing concatenation */
    i__4[0] = 2, a__4[0] = "CD";
    i__4[1] = 4, a__4[1] = cy4;
    i__4[2] = 7, a__4[2] = "/amedas";
    i__4[3] = 1, a__4[3] = "/";
    i__4[4] = 4, a__4[4] = cy4;
    i__4[5] = 2, a__4[5] = cm;
    i__4[6] = 1, a__4[6] = "/";
    i__4[7] = 4, a__4[7] = "area";
    i__4[8] = 2, a__4[8] = csta2;
    i__4[9] = 2, a__4[9] = "/h";
    i__4[10] = 5, a__4[10] = csta5;
    i__4[11] = 1, a__4[11] = "_";
    i__4[12] = 4, a__4[12] = cy4;
    i__4[13] = 2, a__4[13] = cm;
    i__4[14] = 4, a__4[14] = ".csv";
    s_cat(ret_val, a__4, i__4, &c__15, (ftnlen)80);
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__1[0] = 7, a__1[0] = "AMeDAS/";
	i__1[1] = 80, a__1[1] = ret_val;
	s_cat(ret_val, a__1, i__1, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf2h_ */

/* Character */ VOID amdf3d_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im, integer *ista)
{
    /* System generated locals */
    address a__1[2], a__2[4], a__3[5], a__4[16], a__5[15];
    integer i__1[2], i__2[4], i__3[5], i__4[16], i__5[15];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1], ch__6[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cd[1], cm[2], cy2[2], cy4[4], csta2[2], csta5[5];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cy2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__2[2] = 1, a__2[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__2[3] = 1, a__2[3] = ch__5;
    s_cat(cy4, a__2, i__2, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cm, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(csta2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__3[0] = 1, a__3[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__3[1] = 1, a__3[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *ista / 100 % 10 + 48;
    i__3[2] = 1, a__3[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *ista / 10 % 10 + 48;
    i__3[3] = 1, a__3[3] = ch__5;
    *(unsigned char *)&ch__6[0] = *ista % 10 + 48;
    i__3[4] = 1, a__3[4] = ch__6;
    s_cat(csta5, a__3, i__3, &c__5, (ftnlen)5);
    if (*iy <= 2008) {
	if (*im <= 6) {
	    *(unsigned char *)cd = '1';
	}
	if (*im > 6) {
	    *(unsigned char *)cd = '2';
	}
	if (*iy >= 2008) {
	    *(unsigned char *)cd = '1';
	}
/* Writing concatenation */
	i__4[0] = 2, a__4[0] = "CD";
	i__4[1] = 4, a__4[1] = cy4;
	i__4[2] = 1, a__4[2] = cd;
	i__4[3] = 7, a__4[3] = "/amedas";
	i__4[4] = 7, a__4[4] = "/daily/";
	i__4[5] = 4, a__4[5] = cy4;
	i__4[6] = 2, a__4[6] = cm;
	i__4[7] = 1, a__4[7] = "/";
	i__4[8] = 4, a__4[8] = "area";
	i__4[9] = 2, a__4[9] = csta2;
	i__4[10] = 3, a__4[10] = "/ad";
	i__4[11] = 5, a__4[11] = csta5;
	i__4[12] = 1, a__4[12] = "_";
	i__4[13] = 4, a__4[13] = cy4;
	i__4[14] = 2, a__4[14] = cm;
	i__4[15] = 4, a__4[15] = ".csv";
	s_cat(ret_val, a__4, i__4, &c__16, (ftnlen)80);
    } else {
/* Writing concatenation */
	i__5[0] = 2, a__5[0] = "CD";
	i__5[1] = 4, a__5[1] = cy4;
	i__5[2] = 7, a__5[2] = "/amedas";
	i__5[3] = 7, a__5[3] = "/daily/";
	i__5[4] = 4, a__5[4] = cy4;
	i__5[5] = 2, a__5[5] = cm;
	i__5[6] = 1, a__5[6] = "/";
	i__5[7] = 4, a__5[7] = "area";
	i__5[8] = 2, a__5[8] = csta2;
	i__5[9] = 3, a__5[9] = "/ad";
	i__5[10] = 5, a__5[10] = csta5;
	i__5[11] = 1, a__5[11] = "_";
	i__5[12] = 4, a__5[12] = cy4;
	i__5[13] = 2, a__5[13] = cm;
	i__5[14] = 4, a__5[14] = ".csv";
	s_cat(ret_val, a__5, i__5, &c__15, (ftnlen)80);
    }
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__1[0] = 7, a__1[0] = "AMeDAS/";
	i__1[1] = 80, a__1[1] = ret_val;
	s_cat(ret_val, a__1, i__1, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf3d_ */

/* Character */ VOID amdf3h_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im, integer *ista)
{
    /* System generated locals */
    address a__1[2], a__2[4], a__3[5], a__4[16], a__5[15];
    integer i__1[2], i__2[4], i__3[5], i__4[16], i__5[15];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1], ch__6[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cd[1], cm[2], cy2[2], cy4[4], csta2[2], csta5[5];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cy2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__2[2] = 1, a__2[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__2[3] = 1, a__2[3] = ch__5;
    s_cat(cy4, a__2, i__2, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cm, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(csta2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__3[0] = 1, a__3[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__3[1] = 1, a__3[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *ista / 100 % 10 + 48;
    i__3[2] = 1, a__3[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *ista / 10 % 10 + 48;
    i__3[3] = 1, a__3[3] = ch__5;
    *(unsigned char *)&ch__6[0] = *ista % 10 + 48;
    i__3[4] = 1, a__3[4] = ch__6;
    s_cat(csta5, a__3, i__3, &c__5, (ftnlen)5);
    if (*iy <= 2008) {
	if (*im <= 6) {
	    *(unsigned char *)cd = '1';
	}
	if (*im > 6) {
	    *(unsigned char *)cd = '2';
	}
	if (*iy >= 2008) {
	    *(unsigned char *)cd = '1';
	}
/* Writing concatenation */
	i__4[0] = 2, a__4[0] = "CD";
	i__4[1] = 4, a__4[1] = cy4;
	i__4[2] = 1, a__4[2] = cd;
	i__4[3] = 7, a__4[3] = "/amedas";
	i__4[4] = 8, a__4[4] = "/hourly/";
	i__4[5] = 4, a__4[5] = cy4;
	i__4[6] = 2, a__4[6] = cm;
	i__4[7] = 1, a__4[7] = "/";
	i__4[8] = 4, a__4[8] = "area";
	i__4[9] = 2, a__4[9] = csta2;
	i__4[10] = 3, a__4[10] = "/ah";
	i__4[11] = 5, a__4[11] = csta5;
	i__4[12] = 1, a__4[12] = "_";
	i__4[13] = 4, a__4[13] = cy4;
	i__4[14] = 2, a__4[14] = cm;
	i__4[15] = 4, a__4[15] = ".csv";
	s_cat(ret_val, a__4, i__4, &c__16, (ftnlen)80);
    } else {
/* Writing concatenation */
	i__5[0] = 3, a__5[0] = "DVD";
	i__5[1] = 4, a__5[1] = cy4;
	i__5[2] = 7, a__5[2] = "/amedas";
	i__5[3] = 8, a__5[3] = "/hourly/";
	i__5[4] = 4, a__5[4] = cy4;
	i__5[5] = 2, a__5[5] = cm;
	i__5[6] = 1, a__5[6] = "/";
	i__5[7] = 4, a__5[7] = "area";
	i__5[8] = 2, a__5[8] = csta2;
	i__5[9] = 3, a__5[9] = "/ah";
	i__5[10] = 5, a__5[10] = csta5;
	i__5[11] = 1, a__5[11] = "_";
	i__5[12] = 4, a__5[12] = cy4;
	i__5[13] = 2, a__5[13] = cm;
	i__5[14] = 4, a__5[14] = ".csv";
	s_cat(ret_val, a__5, i__5, &c__15, (ftnlen)80);
    }
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__1[0] = 7, a__1[0] = "AMeDAS/";
	i__1[1] = 80, a__1[1] = ret_val;
	s_cat(ret_val, a__1, i__1, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf3h_ */

/* Character */ VOID amdf4d_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im, integer *ista)
{
    /* System generated locals */
    address a__1[2], a__2[4], a__3[5], a__4[13];
    integer i__1[2], i__2[4], i__3[5], i__4[13];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1], ch__6[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cm[2], cy2[2], cy4[4], csta5[5];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cy2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__2[2] = 1, a__2[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__2[3] = 1, a__2[3] = ch__5;
    s_cat(cy4, a__2, i__2, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cm, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__3[0] = 1, a__3[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__3[1] = 1, a__3[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *ista / 100 % 10 + 48;
    i__3[2] = 1, a__3[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *ista / 10 % 10 + 48;
    i__3[3] = 1, a__3[3] = ch__5;
    *(unsigned char *)&ch__6[0] = *ista % 10 + 48;
    i__3[4] = 1, a__3[4] = ch__6;
    s_cat(csta5, a__3, i__3, &c__5, (ftnlen)5);
/* Writing concatenation */
    i__4[0] = 3, a__4[0] = "DVD";
    i__4[1] = 4, a__4[1] = cy4;
    i__4[2] = 7, a__4[2] = "/amedas";
    i__4[3] = 7, a__4[3] = "/daily/";
    i__4[4] = 4, a__4[4] = cy4;
    i__4[5] = 1, a__4[5] = "/";
    i__4[6] = 2, a__4[6] = cm;
    i__4[7] = 1, a__4[7] = "/";
    i__4[8] = 6, a__4[8] = "amd_d_";
    i__4[9] = 4, a__4[9] = cy4;
    i__4[10] = 2, a__4[10] = cm;
    i__4[11] = 1, a__4[11] = ".";
    i__4[12] = 5, a__4[12] = csta5;
    s_cat(ret_val, a__4, i__4, &c__13, (ftnlen)80);
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__1[0] = 7, a__1[0] = "AMeDAS/";
	i__1[1] = 80, a__1[1] = ret_val;
	s_cat(ret_val, a__1, i__1, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf4d_ */

/* Character */ VOID amdf4h_(char *ret_val, ftnlen ret_val_len, integer *iy, 
	integer *im, integer *ista)
{
    /* System generated locals */
    address a__1[2], a__2[4], a__3[5], a__4[13];
    integer i__1[2], i__2[4], i__3[5], i__4[13];
    char ch__2[1], ch__3[1], ch__4[1], ch__5[1], ch__6[1];

    /* Builtin functions */
    /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen), s_cat(char *,
	     char **, integer *, integer *, ftnlen);
    integer s_cmp(char *, char *, ftnlen, ftnlen);

    /* Local variables */
    static char cm[2], cy2[2], cy4[4], csta5[5];

    s_copy(ret_val, " ", (ftnlen)80, (ftnlen)1);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cy2, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *iy / 1000 % 10 + 48;
    i__2[0] = 1, a__2[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *iy / 100 % 10 + 48;
    i__2[1] = 1, a__2[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *iy / 10 % 10 + 48;
    i__2[2] = 1, a__2[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *iy % 10 + 48;
    i__2[3] = 1, a__2[3] = ch__5;
    s_cat(cy4, a__2, i__2, &c__4, (ftnlen)4);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *im / 10 % 10 + 48;
    i__1[0] = 1, a__1[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *im % 10 + 48;
    i__1[1] = 1, a__1[1] = ch__3;
    s_cat(cm, a__1, i__1, &c__2, (ftnlen)2);
/* Writing concatenation */
    *(unsigned char *)&ch__2[0] = *ista / 10000 % 10 + 48;
    i__3[0] = 1, a__3[0] = ch__2;
    *(unsigned char *)&ch__3[0] = *ista / 1000 % 10 + 48;
    i__3[1] = 1, a__3[1] = ch__3;
    *(unsigned char *)&ch__4[0] = *ista / 100 % 10 + 48;
    i__3[2] = 1, a__3[2] = ch__4;
    *(unsigned char *)&ch__5[0] = *ista / 10 % 10 + 48;
    i__3[3] = 1, a__3[3] = ch__5;
    *(unsigned char *)&ch__6[0] = *ista % 10 + 48;
    i__3[4] = 1, a__3[4] = ch__6;
    s_cat(csta5, a__3, i__3, &c__5, (ftnlen)5);
/* Writing concatenation */
    i__4[0] = 3, a__4[0] = "DVD";
    i__4[1] = 4, a__4[1] = cy4;
    i__4[2] = 7, a__4[2] = "/amedas";
    i__4[3] = 9, a__4[3] = "/10min_h/";
    i__4[4] = 4, a__4[4] = cy4;
    i__4[5] = 1, a__4[5] = "/";
    i__4[6] = 2, a__4[6] = cm;
    i__4[7] = 1, a__4[7] = "/";
    i__4[8] = 11, a__4[8] = "amd_10minh_";
    i__4[9] = 4, a__4[9] = cy4;
    i__4[10] = 2, a__4[10] = cm;
    i__4[11] = 1, a__4[11] = ".";
    i__4[12] = 5, a__4[12] = csta5;
    s_cat(ret_val, a__4, i__4, &c__13, (ftnlen)80);
    if (s_cmp(ret_val, " ", (ftnlen)80, (ftnlen)1) != 0) {
/* Writing concatenation */
	i__1[0] = 7, a__1[0] = "AMeDAS/";
	i__1[1] = 80, a__1[1] = ret_val;
	s_cat(ret_val, a__1, i__1, &c__2, (ftnlen)80);
    }
    return ;
} /* amdf4h_ */

shortint inv2_(shortint *i__)
{
    /* System generated locals */
    shortint ret_val;

    /* Local variables */
    static shortint i1, i2;
    static integer n1;

    if (*i__ >= 0) {
	n1 = *i__;
    } else {
	n1 = *i__ + 65536;
    }
    i1 = (shortint) (n1 / 256);
    i2 = (shortint) (n1 % 256);
/* -- For big-endian format -- */
/*      INV2 = 256*I2 + I1 */
/* -- For little-endian format -- */
    ret_val = (shortint) ((i1 << 8) + i2);
    return ret_val;
} /* inv2_ */

integer inv4_(integer *i__)
{
    /* System generated locals */
    integer ret_val;

    /* Local variables */
    static integer i1, i2, i3, i4, n0, n1;

    n1 = *i__;
    n0 = 0;
    if (*i__ < 0) {
	n1 += 1073741824;
	n1 += 1073741824;
	n0 = 1;
    }
    i1 = n1 / 16777216;
    i2 = n1 % 16777216 / 65536;
    i3 = n1 % 65536 / 256;
    i4 = n1 % 256;
    i1 += n0 << 7;
/* -- For big-endian format -- */
/*      INV4 = (256**3)*I4 + (256**2)*I3 + 256*I2 + I1 */
/* -- For little-endian format -- */
    ret_val = (i1 << 24) + (i2 << 16) + (i3 << 8) + i4;
    return ret_val;
} /* inv4_ */

shortint idv2_(integer *m, shortint *i__)
{
    /* System generated locals */
    shortint ret_val;

    /* Local variables */
    static integer n;
    static shortint i1, i2;

    n = *i__;
    i1 = (shortint) (n / 256);
    i2 = (shortint) (n % 256);
    if (*m == 1) {
	ret_val = i1;
    }
    if (*m == 2) {
	ret_val = i2;
    }
    return ret_val;
} /* idv2_ */

/* ======================================== */
/* Subroutine */ int amdid_(integer *nx, integer *ista, real *xx, real *yy, 
	real *hh)
{
    /* Format strings */
    static char fmt_101[] = "(1x,i5,1x,a24,1x,i2,1x,f4.1,1x,i3,1x,f4.1,1x,i4)"
	    ;

    /* System generated locals */
    integer i__1, i__2;
    olist o__1;
    cllist cl__1;

    /* Builtin functions */
    integer f_open(olist *), s_rsfe(cilist *), do_fio(integer *, char *, 
	    ftnlen), e_rsfe(void), f_clos(cllist *);

    /* Local variables */
    static integer n, ih, ix, iy;
    static real rx, ry;
    static char csta__[24];

    /* Fortran I/O blocks */
    static cilist io___136 = { 0, 10, 1, fmt_101, 0 };


/*      CHARACTER*24 CSTA(NX) */
    /* Parameter adjustments */
    --hh;
    --yy;
    --xx;
    --ista;

    /* Function Body */
    o__1.oerr = 0;
    o__1.ounit = 10;
    o__1.ofnmlen = 9;
    o__1.ofnm = "index.txt";
    o__1.orl = 0;
    o__1.osta = "OLD";
    o__1.oacc = 0;
    o__1.ofm = "FORMATTED";
    o__1.oblnk = 0;
    f_open(&o__1);
    i__1 = *nx;
    for (n = 1; n <= i__1; ++n) {
/* L101: */
	i__2 = s_rsfe(&io___136);
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, (char *)&ista[n], (ftnlen)sizeof(integer));
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, csta__, (ftnlen)24);
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, (char *)&iy, (ftnlen)sizeof(integer));
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, (char *)&ry, (ftnlen)sizeof(real));
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, (char *)&ix, (ftnlen)sizeof(integer));
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, (char *)&rx, (ftnlen)sizeof(real));
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = do_fio(&c__1, (char *)&ih, (ftnlen)sizeof(integer));
	if (i__2 != 0) {
	    goto L19;
	}
	i__2 = e_rsfe();
	if (i__2 != 0) {
	    goto L19;
	}
	yy[n] = (real) iy + ry / 60.f;
	xx[n] = (real) ix + rx / 60.f;
	hh[n] = (real) ih;
/* L11: */
    }
L19:
    cl__1.cerr = 0;
    cl__1.cunit = 10;
    cl__1.csta = 0;
    f_clos(&cl__1);
    return 0;
} /* amdid_ */