In [18]:
import AGBTracksUtils
import ResolvedStellarPops as rsp
import fileIO
import os
In [15]:
input_file = '../../input_agb_caf09_s_mar13.dat'
infile = rsp.fileIO.input_file(input_file, default_dict=fileIO.agb_input_defaults())
AGBTracksUtils.AGB_file_setup(infile)
In [16]:
cd ~/research/TP-AGBcalib/code/TPAGB-calib/
In [20]:
metal_dir = infile.metal_dirs[0]
In [22]:
agb_tracks = rsp.fileIO.get_files(os.path.join(infile.working_dir,
metal_dir),
infile.track_identifier)
agb_tracks.sort()
In [28]:
track = fileIO.get_numeric_data(agb_tracks[0])
In [30]:
track.Qs
rows = [q for q in track.Qs]
In [1]:
examine1TP_file = '/home/rosenfield/Dropbox/TP-AGBcal/debug/examine1TP_CAF09_S_OCT13.dat'
In [3]:
lines = open(examine1TP_file, 'r').readlines()
In [25]:
import pprint
import numpy as np
In [94]:
data = np.array([l for l in lines if l.startswith('Z')])
nind = [i for i,l in enumerate(data) if not 'COLIBRI' in l]
mass = np.array([m.strip().split()[3] for m in data], dtype=float)
inds, = np.nonzero((mass>0.55) & (mass < 5.))
missing = list(set(nind) & set(inds))
for i in missing:
if len(data[i].split()) < len(data[i-1].split()):
print data[i].strip()