noz_p=10 &  npr_p=2 &  nrs_p=8 &  ntc_p=11 &  np_p=5 &  nq_p=5
eryth_noz=0L & eryth_npr=0L & eryth_nrs=0L & eryth_ntc=0L 
eryth_np=0L  & eryth_nq=0L

eryth_oz= DBLARR(noz_p)
eryth_pr= DBLARR(npr_p)
eryth_rs= DBLARR(nrs_p)
eryth_tc= DBLARR(ntc_p)
eryth_p= DBLARR(np_p)
eryth_q= DBLARR(nq_p)
eryth_logpr= DBLARR(npr_p)
eryth_val= DBLARR(nq_p, np_p, ntc_p, nrs_p, npr_p, noz_p)


OPENR,2,'eryth3.bin',/f77
i2in=0
i4in=0L

READU,2,i2in &  eryth_noz=i2in
r4in=fltarr(i2in)
READU,2,r4in &  eryth_oz= r4in

READU,2,i2in &  eryth_npr=i2in
r4in=fltarr(i2in)
READU,2,r4in &  eryth_pr= r4in

READU,2,i2in &  eryth_nrs=i2in
r4in=fltarr(i2in)
READU,2,r4in &  eryth_rs= r4in

READU,2,i2in &  eryth_ntc=i2in
r4in=fltarr(i2in)
READU,2,r4in &  eryth_tc= r4in

READU,2,i2in &  eryth_np=i2in
r4in=fltarr(i2in)
READU,2,r4in &  eryth_p= r4in

READU,2,i2in &  eryth_nq=i2in
r4in=fltarr(i2in)
READU,2,r4in &  eryth_q= r4in

READU,2,eryth_val
CLOSE,2

END