;
; PURPOSE: produce NL, FPs, and a set of loops from crns to crne
; WRITTEN: 2002Jun26 Xuepu Zhao
;
pro pnlfpfl_ss,crns,crne,obs,OUTDIR=outdir,OPFL=opfl,r1,r2,EPS=eps

    if N_params( ) lt 1 then begin
       print,'pnlfpfl_ss,crns,crne,obs,outfir=,/opfl,r1,r2,/eps'
       return
    endif
  ; INPUT files
  ; r1=1.4 & r2=1.5
  ; obs='kpk'
    CASE 1 OF
     obs eq 'wso': subdir='/scr221/scr84/xuepu/WSOBLSC/'
     obs eq 'kpk': subdir='/scr221/scr84/xuepu/KPKBRSC/'
     obs eq 'mdi': subdir='/scr221/scr84/xuepu/MDIBRSC/'
     else: print,'No input for obs'
    ENDCASE
    tcrn=crne-crns+1
  ; stcrn=strmid(strtrim(tcrn,2),0,3)
    crnt=indgen(tcrn)+crns 
    scrnt=strtrim(crnt,2)
    brf=subdir+'br25f/br_ss'+scrnt+'_2.50.fits'
    iflf=subdir+'iflf/ifl_ss'+scrnt+'_2.50.d'
    oflf=subdir+'oflf/'+'ofl_ss'+scrnt+'_9.d'
  ; if obs eq 'kpk' $
  ;    then oflf=subdir+'oflf/'+obs+'ofl'+scrnt+'.d' $
  ;    else oflf=subdir+'oflf/'+'ofl_ss'+scrnt+'_9.d'

    if keyword_set(outdir) then $
       outp=outdir else $
       outp='/scr318/zhao/TEST_2002/'+obs
    redu=tcrn MOD 24
    for ind=0,tcrn-1,24 do begin
  ; OUTPUT file
        if ind lt tcrn-redu then ind2=ind+23 else ind2=tcrn-1
        if keyword_set(eps) then begin
           if keyword_set(opfl) then ofile=outp+'nlfpfl_'+scrnt(ind)+'_'+scrnt(ind2)+'.eps' $
                                else ofile=outp+'nlfp_'+scrnt(ind)+'_'+scrnt(ind2)+'.eps'
           psout,241,ofile,/eps
        endif else begin
           if keyword_set(opfl) then ofile=outp+'nlfpfl_'+scrnt(ind)+'_'+scrnt(ind2)+'.ps' $
                                else ofile=outp+'nlfp_'+scrnt(ind)+'_'+scrnt(ind2)+'.ps
           psout,241,ofile
        endelse
        !p.multi=[0,3,8,0,1]
        for ii=ind,ind2 do begin
            inse_crn,crnt(ii),sdts=sdts,sdte=sdte
            sdts=strmid(sdts,0,14)
            sdte=strmid(sdte,0,14)
            ttl='CR'+scrnt(ii)+'  '+sdte+' <= '+sdts
            fnm=findfile(brf(ii),count=count)
            if count eq 1 then begin
               pnlfpc_brfiflf,brf(ii),iflf(ii),0.2,1.01,/sinl,ttl=ttl,/fits
               if keyword_set(opfl) then opfl_oflf,oflf(ii),r1,r2,/sinl
            endif else begin
               bra=fltarr(360,180)
               contour,bra,xrange=[0,360],yrange=[-1,1],xstyle=5,ystyle=5,title=ttl
               sc_annot,1
            endelse
        endfor
        !p.multi=0
        psdone
    endfor

end