def writeModifCFile(self,filename,base_station):
"""
write modified cgx type c file (without header)
need to check which is corrected from tides
re-write the sd as former SD/sqrt(duration)
last column is the first gravity value of the 'o' file
hence the first of the dictionnary, corrected for the tides
add an extra column with the keepdata value (0 or 1) for further
loading
should add a checkbox in the tree object to select the base station
IMPORTANT: grav is corrected from tides (as in the classical raw data),
which is different than for classical c files
and SD is SD, not SD/sqrt(dur) as for classical c files
OR an other option could be to really write 'c' files and modify the
reading function (read_c_file)
"""
offset=self.station_dic[base_station].grav[0]
print "filename: %s"%filename
file=open(filename,'w')
for station_id,sta in sorted(self.station_dic.iteritems(), key=lambda x: x[1].t[1]):
#for station_id,sta in self.station_dic.iteritems():
if sta.keepitem==1:
for i in range(len(sta.t)):
file.write("%6d %4.3f %1.3f %2d %2d %2.1f %2.1f %2.2f %1.3f %3d %3.3f %02d%02d%02d %02d%02d%02d 0 0.000 %3.4f %1d\n"%(sta.station[i],\
sta.grav[i],sta.sd[i],
sta.dur[i], sta.rej[i], sta.tiltx[i], sta.tilty[i],
sta.temp[i],sta.etc[i],sta.t[i].timetuple()[7],
sta.t[i].hour*60+sta.t[i].minute+sta.t[i].second/60,
sta.t[i].day,sta.t[i].month,sta.t[i].year-2000,
sta.t[i].hour,sta.t[i].minute,sta.t[i].second,
sta.grav[i]-offset, sta.keepdata[i]))
file.close()
评论列表
文章目录
