; 
; PURPOSE: Calculate outward field lines 
; WRITTEN: May 2000	Xuepu Zhao
;
pro oflfss_gaha06,ga,ha,iph,jth,oflf,Rss=rss,RBOTTOM=Rbottom,DSTEP=dstep

    if N_params() lt 1 then begin
       print,'oflfss_gaha06,ga,ha,iph,jth,oflf,rss=,rbottom=,dstep='
       return
    endif
    zgrid06,iph,jth,phd,thd,lad,cth,sth
    get_lun,unit
    openw,unit,oflf
;   iph1=iph/2-iph/4 & iph2=iph/2+iph/4
;   for i=iph1,iph2 do begin
    for i=0,iph-1 do begin
        phri=phd(i)*!DTOR
        for j=0,jth-1 do begin
            thrj=thd(j)*!DTOR
            if keyword_set(rss) then begin
               if keyword_set(rbottom) then begin
                  if keyword_set(dstep) then begin
                     oflijss_gaha06,ga,ha,phri,thrj,np,xv,yv,zv,rv,tvv,$
                                pv,rss=rss,rbottom=rbottom,dstep=dstep
                  endif else begin
                     oflijss_gaha06,ga,ha,phri,thrj,np,xv,yv,zv,rv,tvv,$
                                pv,rss=rss,rbottom=rbottom
                  endelse
               endif else begin
                  oflijss_gaha06,ga,ha,phri,thrj,np,xv,yv,zv,rv,tvv,$
                                pv,rss=rss
               endelse
            endif else begin
               oflijss_gaha06,ga,ha,phri,thrj,np,xv,yv,zv,rv,tvv,pv
            endelse
            printf,unit,np
            anp=abs(np)
            for k=0,anp-1 do printf,unit,xv(k),yv(k),zv(k),rv(k),tvv(k),pv(k)
        endfor
    endfor
    free_lun,unit

end