ihi = ilnblnk(data_errfile)
open(gwunit, file=data_errfile(ilo:ihi),
integer ilo,ihi
ihi = ilnblnk(geoid_covariancefile)
& geoid_covariancefile(ilo:ihi)
open(egmunit, file = geoid_covariancefile(ilo:ihi),
I data_errfile(ilo:ihi),