;
; read gw and hw from  ghfile
;
pro gaha_ghfcsss,ghf,nmax,gw,hw

    if N_params() lt 1 then begin
       print,'pro gaha_ghfcsss,ghf,nmax,gw,hw'
       print,'INPUT: ghf,nmax'
       return
    endif

    nn=nmax+1
    gw=fltarr(nn,nn)
    hw=gw

    str=' '
    get_lun,unit
    openr,unit,ghf
    readf,unit,str
    readf,unit,n90
    for n=0,Nmax do begin
        for m=0,n do begin
            readf,unit,ni,mi,g0,h0
            gw(n,m)=g0
            hw(n,m)=h0
        endfor
    endfor
    free_lun,unit

end