;-------------------------------------------------------------
;datalist_to_ingal_driver.pro
;-------------------------------------------------------------
pro datalist_to_ingal_driver_f,datalistide

  if(n_params() le 0) then begin
     print,'--------------------------------------'
     print,'datalist_to_ingal_driver_f,datalistide'
     print,'--------------------------------------'
     print,''
     print,'makes simple ingal-files for galaxies listed in datalistide'
     print,'example:'
     print,"datalist_to_ingal_driver_f,'s4g_pipeline2_test081209'"
     print,'heikki.salo@oulu.fi   081209'
     print,'datalist version=2  (sky,dsky) added 090610'
     return
  endif


;make galfit input parameter files
;for "10galaxies"
;creates b, bd, bdbarnmodel

;------------------------------------------------------
;read datalist-file containing pipeline 1&2 parameters
; eq. infile='pipe3_example.datalist'

;runfile = write idl-driver procedure for making the
;          GALFIT interations with the INGAL's created here


  infile=datalistide+'.datalist'
  print,infile

  runfile=datalistide+'_driver.pro'
  print,runfile

  close,1
  openr,1,infile

  close,2
  openw,2,runfile
  printf,2,';Automatically generated driver-procedure'
  printf,2,';'+!stime

  printf,2,';new=1 -> make decomposition iterations'
  printf,2,';new=0 -> display previously done decomposition iterations'
  printf,2,';         (if they exist; otherwise make iteration)'
  printf,2,';ps=0 --> screen'
  printf,2,';ps=2 --> ps and png files'

  printf,2,'new=1'
  printf,2,'ps=2'
  printf,2,'  '


  line=''
  readf,1,line                  ;version line

  version=0
  if(line eq 'version=1') then version=1
  if(line eq 'version=2') then version=2

  while not eof(1) do begin   
     readf,1,line
     start=strmid(line,0,1)
     
     if(strmid(line,0,1) ne '#') then begin ;not a comment 
        print,'------------------------------------------------------------'
        print,line
        print,'------------------------------------------------------------'
        res=strsplit(line,/extr)


           ide=res(0)
           data_fits=res(1)
           xc=float(res(2))-1.
           yc=float(res(3))-1.
           pas=res(4)
           ells=res(5)
           smas=res(6)

        if(version eq 0) then begin
           smult='1.3'
           skys='0.'
           dskys='0.'
           mask_fits=res(7)
           sigma_fits=res(8)
           psf=res(9)
        endif

        if(version eq 1) then begin
           smult=res(7)
           skys='0.'
           dskys='0.'
           mask_fits=res(8)
           sigma_fits=res(9)
           psf=res(10)
        endif

        if(version eq 2) then begin
           smult=res(7)
           skys=res(8)
           dskys=res(9)
           mask_fits=res(10)
           sigma_fits=res(11)
           psf=res(12)
        endif


;cut around denter
        cut=float(smas)*float(smult)

        
        rexpo_guess=float(smas)/4.


;make simple ingal files
;-------------------------------------------------------------
;bulge

        ide_use=ide+'_b'      
        datalist_to_ingal,ide_use,data_fits, sigma_fits, mask_fits, psf,$
                          xc_idl=xc,yc_idl=yc,cut=cut

        printf,2,'galfit_run_new,"'+ide_use,'",ps=ps,new=new'

;-------------------------------------------------------------
;bulge+disk

        if(strpos(ells,'zz') eq -1) then begin
           ide_use=ide+'_bd'
           disk=[1.-float(ells),float(pas),rexpo_guess]
           datalist_to_ingal,ide_use,data_fits, sigma_fits, mask_fits, psf,$
                             xc_idl=xc,yc_idl=yc,disk=disk,cut=cut

           printf,2,'galfit_run_new,"'+ide_use,'",ps=ps,new=new'

           ide_use=ide+'_bdf'
           disk=[1.-float(ells),float(pas),rexpo_guess]

           datalist_to_ingal,ide_use,data_fits, sigma_fits, mask_fits, psf,$
                             xc_idl=xc,yc_idl=yc,disk=disk,cut=cut,/free

           printf,2,';galfit_run_new,"'+ide_use,'",ps=ps,new=new'
        
        endif

;-------------------------------------------------------------
;bulge+zdisk

        if(strpos(ells,'zz') ne -1) then begin
           ide_use=ide+'_bz'
           zdisk=float(pas)
           datalist_to_ingal,ide_use,data_fits, sigma_fits, mask_fits, psf,$
                             xc_idl=xc,yc_idl=yc,zdisk=zdisk,cut=cut

           printf,2,'galfit_run_new,"'+ide_use,'",ps=ps,new=new'
        endif


;MODEL TEMPLATES FOR CUT/PASTE/EDIT

;-------------------------------------------------------------
;bulge+disk+bar+nuc  model
        if(strpos(ells,'zz') eq -1) then begin
           ide_use=ide+'_bdbarnmodel'
           disk=[1.-float(ells),float(pas),rexpo_guess]
           datalist_to_ingal,ide_use,data_fits, sigma_fits, mask_fits, psf,$
                             xc_idl=xc,yc_idl=yc,disk=disk,/bar,/nuc,cut=cut
           printf,2,';galfit_run_new,"'+ide_use,'",ps=ps,new=new'
           printf,2,';'
        endif

;-------------------------------------------------------------
;bulge+zdisk+bar+nuc model
        if(strpos(ells,'zz') ne -1) then begin
           ide_use=ide+'_bzbarnmodel'
           zdisk=float(pas)
           datalist_to_ingal,ide_use,data_fits, sigma_fits, mask_fits, psf,$
                             xc_idl=xc,yc_idl=yc,zdisk=zdisk,/bar,/nuc,cut=cut

           printf,2,';galfit_run_new,"'+ide_use,'",ps=ps,new=new'
           printf,2,';'
        endif



     endif
skip:
  endwhile

  close,1
  printf,2,'end'
  close,2

end