;
; PURPOSE: read ga and ha from  ghfile
; NOTE: g and h in GAUSS both WSO and others
; WRITTEN:	May 2006	Xuepu Zhao
;
pro gaha_ghf06,ghf,nmax,ga,ha

    if N_params() lt 1 then begin
       print,'pro gaha_ghf06,ghf,nmax,ga,ha'
       print,'INPUT: ghf,nmax'
       return
    endif
;
;   if keyword_set(n0) then n0=0 else n0=1
    nn=nmax+1
    ga=fltarr(nn,nn)
    ha=ga
    str=' '
    get_lun,unit
    openr,unit,ghf
    readf,unit,str
    readf,unit,nmax0
  ; for n=0,Nmax do begin              
    readf,unit,ni,mi,g0,h0
    for n=1,Nmax do begin              ; 26JAN2006
        for m=0,n do begin
            readf,unit,ni,mi,g0,h0
            ga(n,m)=g0
            ha(n,m)=h0
        endfor
    endfor
    free_lun,unit

end