Ejemplo n.º 1
0
    def makeMCMaxStars(self, models):
        '''
        Set parameters for star_list taken from the MCMax database.
        
        Based on the model id of MCMax.
        
        @param models: model_ids for the MCMax db
        @type models: list(string)
        @return: The model instances 
        @rtype: list(Star())
        
        '''

        star_grid = Star.makeStars(models=models,code='MCMax',id_type='MCMax',\
                                   path=self.path)
        for star, model in zip(star_grid, models):
            filepath = os.path.join(cc.path.mout,'models',\
                                    star['LAST_MCMAX_MODEL'])
            denstemp = os.path.join(filepath, 'denstemp.dat')
            logfile = os.path.join(filepath, 'log.dat')
            grid_shape = DataIO.getMCMaxOutput(filename=denstemp,incr=1,\
                                               keyword='NGRAINS',single=0)[0]
            star.update({'NTHETA':int(grid_shape[1]),\
                         'NRAD':int(grid_shape[0]),\
                         'T_STAR':float(DataIO.getMCMaxOutput(filename=logfile,\
                                                incr=0,\
                                                keyword='STELLAR TEMPERATURE',\
                                                single=0)[0][2]),\
                         'R_STAR':float(DataIO.getMCMaxOutput(filename=logfile,\
                                                incr=0,\
                                                keyword='STELLAR RADIUS',\
                                                single=0)[0][2])})
        return star_grid
Ejemplo n.º 2
0
 def makeMCMaxStars(self,models):
     
     '''
     Set parameters for star_list taken from the MCMax database.
     
     Based on the model id of MCMax.
     
     @param models: model_ids for the MCMax db
     @type models: list(string)
     @return: The model instances 
     @rtype: list(Star())
     
     '''
     
     star_grid = Star.makeStars(models=models,code='MCMax',id_type='MCMax',\
                                path=self.path)
     for star,model in zip(star_grid,models):    
         filepath = os.path.join(cc.path.mout,'models',\
                                 star['LAST_MCMAX_MODEL'])
         denstemp = os.path.join(filepath,'denstemp.dat')
         logfile = os.path.join(filepath,'log.dat')
         grid_shape = DataIO.getMCMaxOutput(filename=denstemp,incr=1,\
                                            keyword='NGRAINS',single=0)[0]
         star.update({'NTHETA':int(grid_shape[1]),\
                      'NRAD':int(grid_shape[0]),\
                      'T_STAR':float(DataIO.getMCMaxOutput(filename=logfile,\
                                             incr=0,\
                                             keyword='STELLAR TEMPERATURE',\
                                             single=0)[0][2]),\
                      'R_STAR':float(DataIO.getMCMaxOutput(filename=logfile,\
                                             incr=0,\
                                             keyword='STELLAR RADIUS',\
                                             single=0)[0][2])})            
     return star_grid  
Ejemplo n.º 3
0
    def setModels(self,star_grid=[],models=[]):
        
        '''
        Load the models and remember them.
        
        @keyword star_grid: The parameter sets, if not given: model ids needed
        
                            (default: [])
        @type star_grid: list[Star()]
        @keyword models: the PACS ids, only relevant if star_grid == [] and if 
                         instrument == PACS. In all other cases a star_grid is
                         required.
        
                         (default: [])
        @type models: list[string]
        @keyword extra_keywords: any instrument specific keywords that you need
        
                                 (default: dict())
        @type extra_keywords: dict
        
        '''        

        #-- The unresolved-data case
        if self.instrument:
            instr = self.instrument.instrument.upper()
            print '***********************************'
            print '* Checking Sphinx models for comparison with %s.'%instr
            if instr == 'SPIRE': models = []
            if not star_grid and not models:
                raise IOError('Statistics.setModels requires either ' + \
                              'star_grid or models to be defined.')
            #-- star_grid can be reconstructed for PACS through its database
            elif not star_grid:
                star_grid = Star.makeStars(models=models,id_type='PACS',\
                                           path=path_code,code=self.code)
                self.instrument.addStarPars(star_grid)
            if set([s['MOLECULE'] and 1 or 0 for s in star_grid]) == set([0]): 
                return
            self.instrument.prepareSphinx(star_grid)
            self.star_grid = [star 
                              for star in star_grid 
                              if star['LAST_%s_MODEL'%instr] <> None]
        #-- The FREQ_RESO case
        else:
            if not star_grid:
                raise IOError('Statistics.setModels requires a ' + \
                              'star_grid to be defined for freq-resolved data.')
            if set([s['MOLECULE'] and 1 or 0 for s in star_grid]) == set([0]): 
                return
            self.star_grid = star_grid
Ejemplo n.º 4
0
    def makeStars(self, models):
        '''
        Make a Star list based on the Chemistry model ids.
        
        @param models: model_ids for the Chemistry db
        @type models: list(string)
        
        @return: the parameter sets
        @rtype: list[Star()]
        
        '''

        star_grid = Star.makeStars(models=models,\
                                   id_type='Chemistry',\
                                   code='Chemistry',path=self.path)
        [star.addCoolingPars() for star in star_grid]
        return star_grid
Ejemplo n.º 5
0
 def makeStars(self,models):
     
     '''
     Make a Star list based on the Chemistry model ids.
     
     @param models: model_ids for the Chemistry db
     @type models: list(string)
     
     @return: the parameter sets
     @rtype: list[Star()]
     
     '''
     
     star_grid = Star.makeStars(models=models,\
                                id_type='Chemistry',\
                                code='Chemistry',path=self.path)
     [star.addCoolingPars() for star in star_grid]
     return star_grid