result = read3fromopa(unit, params, num)
result = readoldopadistcoast(filename, jpiglo, jpjglo, jpk, IBLOC=IBLOC, JPBYT=JPBYT, NUMREC=NUMREC)