Rsoft7 part2

import sys
import os
from math import *

import numarray
import numarray.fft

import rsoft.rspytools as rspy

if len（sys.argv）<2:
  print ‘usage: filter ‘
  sys.exit（1）

fname=sys.argv[1]
if not os.path.exists（fname）:
  print “Can‘t find %s“%fname
  sys.exit（1）

prefix=fname[:-4]

#load .tmn file
tmn_data=rspy.file_to_numarray（fname）

v_t=tmn_data[:0]
v_ey=tmn_data[:1]
clight=2.9979e14

#trim off start of .tmn file containing CW ramp
npts=len（v_t）
highpow2=int（floor（log（npts*1.0）/log（2.0）））
goodpts=int（pow（2 highpow2））

#apply raised cosine window function
dt=（v_t[1]-v_t[0]）/clight
v_ey=v_ey[-goodpts:]
cos_domain=numarray.arange（goodpts typecode=‘d‘）/（goodpts-1）*pi-pi/2
filt_ey=v_ey*numarray.cos（cos_domain）
  
#make fft
v_ey_fft=numarray.fft.fft（filt_ey）
v_ey_fft_pos=numarray.abs（v_ey_fft[0:goodpts/2]）/goodpts
v_nu=numarray.arange（goodpts/2 typecode=‘d‘）/（goodpts*dt）/1e12  # freq in THz

#output windowed fft data
rspy.numarray_to_file（‘%s.fft2‘%prefix v_ey_fft_pos x=v_nu）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    9298928  2006-08-03 11:00  RSOFT7docsandsolve.pdf

     文件    3543151  2009-11-10 23:56  RSOFT7docseamprop.pdf

     文件    1272279  2006-08-02 11:00  RSOFT7docsdiffractmod.pdf

     文件     896311  2006-08-03 11:00  RSOFT7docsfemsim.pdf

     文件    2557856  2010-03-11 11:39  RSOFT7docsfullwave.pdf

     文件    1299668  2006-08-21 11:00  RSOFT7docsgratingmod.pdf

     文件     174591  2006-08-21 11:00  RSOFT7docsinstall.pdf

     文件    3921843  2006-04-05 13:48  RSOFT7docslasermod.pdf

     文件    4198886  2006-08-03 11:00  RSOFT7docsmost.pdf

     文件     540510  2002-04-03 11:00  RSOFT7docs
plot.pdf

     文件     889326  2004-10-20 11:00  RSOFT7docs
sdatabrowser.pdf

     文件    1965925  2006-08-03 11:00  RSOFT7docs
softcad.pdf

     文件        288  2006-03-30 11:00  RSOFT7etccadw32.reg

     文件        766  2003-04-15 11:00  RSOFT7etcsetup.ico

     文件       2252  2006-03-30 11:00  RSOFT7etcwinplot.reg

     文件        304  2006-03-30 10:00  RSOFT7etcxpcadw32.reg

     文件       2268  2006-03-30 10:00  RSOFT7etcxpwinplot.reg

     文件        133  2005-08-05 17:30  RSOFT7exampleBandSOLVEApplicationNotesmakefile

     文件    3927144  2008-09-17 16:19  RSOFT7exampleBandSOLVEApplicationNotespcslabsstmp_domain.vtk

     文件      15318  2005-06-15 17:51  RSOFT7exampleBandSOLVEApplicationNotespcslabsholes_solid_cladding.ind

     文件        573  2005-08-05 17:28  RSOFT7exampleBandSOLVEApplicationNotespcslabsmakefile

     文件      15868  2003-12-10 21:56  RSOFT7exampleBandSOLVEApplicationNotespcslabs
osymmetry.ind

     文件       1304  2006-04-03 14:23  RSOFT7exampleBandSOLVEApplicationNotespcslabs
osymmetry.kpts

     文件      16012  2005-06-15 17:48  RSOFT7exampleBandSOLVEApplicationNotespcslabs
ods_periodic_cladding.ind

     文件      10990  2005-05-12 13:38  RSOFT7exampleBandSOLVEApplicationNotespcswaveguidescubiccircref_clad.ind

     文件       9358  2005-05-12 13:38  RSOFT7exampleBandSOLVEApplicationNotespcswaveguidescubicclad_lightline.ind

     文件       9414  2005-08-15 11:18  RSOFT7exampleBandSOLVEApplicationNotespcswaveguidescubicclad_lightline_full.ind

     文件        343  2005-08-12 04:28  RSOFT7exampleBandSOLVEApplicationNotespcswaveguidescubiccombined.pcs

     文件       1334  2005-08-15 12:29  RSOFT7exampleBandSOLVEApplicationNotespcswaveguidescubiccombined_scan.pcs

     文件       9406  2005-05-12 13:38  RSOFT7exampleBandSOLVEApplicationNotespcswaveguidescubiccubicslab.ind

............此处省略1121个文件信息