; NAME :	saveft1
; PURPOSE : convert data of solar flare telescope to the standard
;		data format of MAGPACK2
;
; CATEGORY :
; 	solar magnetic field computation package (MAGPACK2)
; CALLING SEQUENCE :
;		.run saveft1
; INPUTS :
;		arrays h,img,q,u,v should exist in the memory
; OPTIONAL INPUT PARAMETERS : 
;	none
; KEYWORD PARAMETERS :
; OUTPUTS :
;
; COMMON BLOCKS : 	none
; SIDE EFFECTS :
; RESTRICTIONS :	X-window environment is assumed
; PROCEDURE :
;
; MODIFICATION HISTORY :
;	T.Sakurai, 1992
;
;------------------------------------------------------------------------
;
convft1, h,img,q,u,v, mparam, rdata

pol2b, mparam, rdata
curlb1, mparam, rdata, /keep_bp

if strmid(fnam0,13,2) eq '.p' then app='' else app='.'

mfile = '/solar/flare/fieldlines/FB'+strmid(fnam0, 2,12) +app+'d'

; 	mfile = '/solar/flare/fieldlines/FB'+strmid(fnam0, 2,12) +app+'D'

lfile = '/solar/flare/fieldlines/FB'+strmid(fnam0, 2,12) +app+'L'

writem, mparam, rdata, filename=mfile
print,'output file=',mfile
print,'create data file for field line computation? (y/n)'
yn=''
read,yn
if yn eq 'n' or yn eq 'N' then goto, nofile

seqno=0
countup:	seqno=seqno+1
	datafile='/home/wind/solar/data/fname'+string(seqno,format='(i2.2)')+'.dat'
	exist=findfile(datafile,count=count)
	if count gt 0 then goto, countup

print,'data file =',datafile
openw, unit, datafile, /get_lun
printf, unit, ''
printf, unit, mfile
printf, unit, lfile
close, unit
free_lun, unit
nofile:
dispm, mparam, rdata
end