restore,'fixshots.idl'
;xnpl=xnstre & ynpl=ynstre
fmt='(10(f6.2,1x))'
fmt1='(3(f7.3,1x))'
num=n_elements(rmx)
openw,lu,'fixshots'+'.dat',/get_lun 
printf,lu,string
for ii=0,num-1 do begin &$
;xsw=xnpl*rmx(ii)*1000.0 & ysw=ynpl*rmx(ii)*1000.0 &$
;elc=elon(ii) & nlc=nlat(ii) &$
print,ii,rmx(ii) &$
;XYTOLL,xsw,ysw,eelon,nnlat,elc,nlc &$
;eelon=reform(eelon,151,39) &$
;nnlat=reform(nnlat,151,39) &$
printf,lu,format=fmt1,rmx(ii),nlat(ii),elon(ii) &$
print,rmx(ii),nlat(ii),elon(ii) &$
for j=0,35 do printf,lu,format=fmt,usnap(*,j,ii) &$
for j=0,35 do printf,lu,format=fmt,vsnap(*,j,ii) &$
endfor
close,lu & free_lun,lu 
;
exit
