ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/JSOC/proj/globalhs/scripts/idl/mkpar1_asym.idl
Revision: 1.1
Committed: Fri Apr 23 09:16:06 2010 UTC (13 years, 5 months ago) by tplarson
Branch: MAIN
CVS Tags: Ver_6-0, Ver_6-1, Ver_6-2, Ver_6-3, Ver_6-4, Ver_5-14, Ver_5-13, Ver_5-12, Ver_5-11, Ver_5-10, Ver_5-9, Ver_8-8, Ver_8-2, Ver_8-3, Ver_8-0, Ver_8-1, Ver_8-6, Ver_8-7, Ver_8-4, Ver_8-5, Ver_7-1, Ver_7-0, Ver_8-10, Ver_8-11, Ver_8-12, Ver_9-0
Log Message:
idl functions, procedures, scripts, and input files needed by doretidl, dogfidl, and dopkgn

File Contents

# Content
1
2 coef=[-3.768, 12.808, -17.242, 11.684, -4.180, 0.750, -0.0529]
3 q=readx('par0',nlin=nl)
4 f0=q[2,*]
5 x=f0[*]/1000
6 diff=poly(x,coef)
7 w1=where(f0 lt 1000)
8 diff[w1]=0
9 w2=where(f0 gt 4000)
10 diff[w2]=0.5
11 f=1000*x+diff
12
13 c=[1.55,-0.07,-0.01]
14 x=f/1000
15 a=mypoly(x)
16 asym=a#c
17 asym[w1]=1.55
18 asym[w2]=0.0
19
20 openw,1,'par1'
21 for i=0,nl-1 do printf,1,q[0:1,i],f[i],q[3:11,i],asym[i],format='(i2,i4,f11.4,99f10.6)'
22 close,1