pro norh_koshix,rawfile,szfile,frame $ ,imgdir=imgdir $ ,outfile=outfile,rlphdiffile=rlphdiffile $ ,version=version,freq=freq $ ,cenfnl=cenfnl,cenunit=cenunit,size=size $ ,nfrcal=nfrcal,nfrint=nfrint,clcra=clcra,clcrs=clcrs $ ,cendty=cendty,beamcor=beamcor $ ,lpfilt=lpfilt,clgain_sky=clgain_sky,clgain=clgain $ ,cltrim=cltrim,clnega=clnega $ ,jit_as_int=jit_as_int,pro_res=pro_res if (n_elements(outfile) eq 0) then outfile=curdir()+'/'+'norh_synth_input' if (n_elements(imgdir) eq 0) then imgdir=curdir()+'/' if (n_elements(version) eq 0) then version='-1' if (version eq '-1') then version='6.1' if (n_elements(freq) eq 0) then freq=17 if (n_elements(size) eq 0) then size=512 if (n_elements(cenunit) eq 0) then cenunit=1 if (n_elements(cenfnl) eq 0) then cenfnl=[0,0] if (n_elements(nfrcal) eq 0) then nfrcal=1 if (n_elements(nfrint) eq 0) then nfrint=10 if (n_elements(clcra) eq 0) then clcra=3000 if (n_elements(clcrs) eq 0) then clcrs=3000 if (n_elements(cendty) eq 0) then cendty=[999,999] if (n_elements(beamcor) eq 0) then beamcor=1 if (n_elements(cendtyl) eq 0) then cendtyl=cendty if (n_elements(cendtyr) eq 0) then cendtyr=cendty if (n_elements(lpfilt) eq 0) then lpfilt=0.1 if (n_elements(clgain_sky) eq 0) then clgain_sky=0.02 if (n_elements(clgain) eq 0) then clgain=0.2 if (n_elements(cltrim) eq 0) then cltrim=0.8 if (n_elements(clnega) eq 0) then clnega=1 if (n_elements(jit_as_int) eq 0) then jit_as_int=1 if (n_elements(pro_res) eq 0) then pro_res=1 nfrina=1 solr=1.0125 mfrx=n_elements(frame) case 1 of (version eq '5.1') : begin if (freq ne 17) then message,'Sorry, this version only supports 17GHz data' openw,unit,outfile,/get_lun printf,unit,rawfile printf,unit,szfile printf,unit,rlphdiffile printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,nfrina printf,unit,clcra printf,unit,clcrs printf,unit,solr printf,unit,clgain_sky printf,unit,clgain printf,unit,cltrim printf,unit,lpfilt printf,unit,clnega printf,unit,cendtyr printf,unit,cendtyl printf,unit,cenfnl printf,unit,jit_as_int printf,unit,pro_res printf,unit,beamcor printf,unit,imgdir close,unit free_lun,unit end (version eq '6.0') : begin if (freq ne 17) then message,'Sorry, this version only supports 17GHz data' openw,unit,outfile,/get_lun printf,unit,rawfile printf,unit,szfile printf,unit,rlphdiffile printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,nfrina printf,unit,clcra printf,unit,clcrs printf,unit,solr printf,unit,clgain_sky printf,unit,clgain printf,unit,cltrim printf,unit,lpfilt printf,unit,clnega printf,unit,cendtyr printf,unit,cendtyl printf,unit,cenunit,cenfnl printf,unit,jit_as_int printf,unit,pro_res printf,unit,beamcor printf,unit,size printf,unit,imgdir close,unit free_lun,unit end (version eq '6.1') : begin if (freq ne 17) then message,'Sorry, this version only supports 17GHz data' openw,unit,outfile,/get_lun printf,unit,imgdir printf,unit,rawfile printf,unit,szfile printf,unit,size printf,unit,cenunit,cenfnl printf,unit,0 printf,unit,mfrx printf,unit,format='(10(a,1x))',frame printf,unit,nfrcal printf,unit,nfrint printf,unit,clcra,clcrs close,unit free_lun,unit end endcase return end