pro mkparams
njobs = 30
nparams = 3
sain_params = fltarr(nparams,njobs)
sain_params = [                $
   0.8905,   0.8154,   0.4781, $
   0.0621,   0.6018,   0.2052, $
   0.0046,   0.1438,   0.4031, $
   0.2488,   0.6380,   0.6197, $
   0.2845,   0.5442,   0.8886, $
   0.1482,   0.3481,   0.9990, $
   0.4216,   0.2553,   0.7807, $
   0.2191,   0.5159,   0.0938, $
   0.3199,   0.4626,   0.3731, $
   0.9847,   0.4091,   0.1332, $
   0.1999,   0.8674,   0.3590, $
   0.0993,   0.4845,   0.1987, $
   0.4681,   0.0496,   0.7527, $
   0.3501,   0.3931,   0.9163, $
   0.7471,   0.6984,   0.2519, $
   0.5464,   0.7610,   0.6791, $
   0.5244,   0.9994,   0.7140, $
   0.6522,   0.5741,   0.6612, $
   0.6871,   0.0908,   0.0097, $
   0.9392,   0.9522,   0.8316, $
   0.4524,   0.7141,   0.4384, $
   0.1333,   0.7814,   0.5360, $
   0.5935,   0.2927,   0.5306, $
   0.7838,   0.2259,   0.1337, $
   0.8095,   0.3281,   0.3113, $
   0.9045,   0.1314,   0.8611, $
   0.7255,   0.0252,   0.5919, $
   0.6245,   0.1677,   0.0634, $
   0.3925,   0.8514,   0.9561, $
   0.8391,   0.9303,   0.2917  ]
elden_min = fltarr(njobs)
amp   = fltarr(njobs)
phase = fltarr(njobs)
amp_max = 36000.
phase_max = 12.
print,'job   amplitude       phase   elden_min'
print,'---------------------------------------'
for ijob=0,njobs-1 do begin
  i = (ijob+1)*3
  amp[ijob]   = amp_max*sain_params[i-2] 
  phase[ijob] = phase_max*sain_params[i-1] 
  elden_min[ijob] = 10.^(3.+sain_params[i-3])
  print,format="(i3,3f12.3)",ijob+1,amp[ijob],phase[ijob],elden_min[ijob]
endfor
print,format="('set amplitude = ( \',/,(6f12.3,' \'))",amp
print,format="('set phase = ( \',/,(6f12.3,' \'))",phase
print,format="('set elden_min = ( \',/,(6f12.3,' \'))",elden_min
end