cmd.h

Go to the documentation of this file.

#ifndef __CMD_H__
#define __CMD_H__

#ifdef __cplusplus
extern "C" {
#endif 

/*  gridded CMD parameters */
#define XSTEPS 500
#define YSTEPS 500

/*  x and y ranges: in function set_labels_ranges */

/*  end gridded CMD parameters */
/*  plot ranges for non-gridded CMD: in the routine plot_cmd */

/*  factor determining the number ratio of stars you want in your (non-gridded) CMD */
/*  e.g. if 1e6 stars are needed to reach the galaxy stellar mass and NUMFAC  */
/*  is 100 then about 1e4 stars will be used for the CMD. */
#define NUMFAC 1

/*  determines the maximum number of stars you want in the CMD. If you use this */
/*  it makes the most sense to take NUMFAC == 1. */
/*  this is also an array index so it can't be bigger than the largest possible */
/*  integer */
#define MAX_NUM_STARS 5e5

#ifdef CMD_GIRARDI

/*  number of isochrone files */
#define N_CMD_FILES 7
/* #define N_CMD_FILES 6 */

/*  max number of columns to be read from the isochrone files */
#define CMD_FILES_COLUMNS 16

/*  max number of rows to be read from the isochrone files */
#define CMD_FILES_ROWS 10000

/*  number of metallicity bins */
#define CMD_FILE_ZMET_BINS N_CMD_FILES

/*  number of age bins  */
#define CMD_FILE_AGE_BINS 64
/* #define CMD_FILE_AGE_BINS 50 */

/*  number of mass bins */
#define CMD_FILE_MASS_BINS 240

/*  the columns where the different variables are found in the files. */
/*  REMEMBER: start at 0 */
#define POS_MASS 1
#define POS_LUMINOSITY 3
#define POS_LOGTEMP 4
#define POS_AGE 0
#define POS_MB  8
#define POS_MV  9
#define POS_MR  10
#define POS_MI  11

extern double cmd_file_metallicity[];                     
#endif /*  CMD_GIRARDI */


#ifdef GRIDDEDCMD

extern struct cmd_data
{
  /*  will contain the values at the centers of the box. */
  double logtemp[XSTEPS]; /*  temperature */
  double loglumin[YSTEPS];  /*  luminosity */
  
  /*  must be float for pgplot */
  float templum[XSTEPS][YSTEPS];

  int ntempsteps;
  int nluminsteps;

  double diag_half;   /*  half the diagonal of an element */
  double logluminmin;
  double logluminmax;
  double logtempmin;
  double logtempmax;
} cmd;

/*  rough guess for the number of bins an isochrone part will */
/*  pass out mass to */
#define NUM_SEC_VALS ((int) (XSTEPS*YSTEPS/100) + 1)
  //#define NUM_SEC_VALS 200

extern int cmd_numbers[NUM_SEC_VALS][2];

/* temporarily store the values of the isochrone that go onto the grid */
extern float cmd_value[NUM_SEC_VALS];
extern int cmd_value_size;

#endif /*  GRIDDEDCMD */

#ifdef DO_CMD_TEST
/*  array that will contain the mass of each sampled star on the CMD */
extern float * cmd_masses;
#endif


/*  struct that will contain the data read from the isochrone files */
/*  cmd_file[age][metallicity] */
extern struct cmd_file_data
{
  double metallicity;                          
  double logage;                               
  double mass[CMD_FILE_MASS_BINS];
  float loglumin[CMD_FILE_MASS_BINS];
  float logtemperature[CMD_FILE_MASS_BINS];
} **cmd_file;


extern struct cmd_res_data
{
  int izmet;                                  
  int iage;                                   
  double metallicity;                         
  double logage;                              
  double massmin;
  double massmax;

  int imassmax;                               
  double norm;

  double mass[CMD_FILE_MASS_BINS];
  float logtemperature[CMD_FILE_MASS_BINS];
  float loglumin[CMD_FILE_MASS_BINS];

#ifdef GRIDDEDCMD
  double xval[CMD_FILE_MASS_BINS];            
  double yval[CMD_FILE_MASS_BINS];            
#endif

}
#ifdef GRIDDEDCMD
  cmd_res;
#else
cmd_res, *cmd_row;

extern float cmd_final[2][(int) MAX_NUM_STARS];
#endif


extern struct isochrone_data
{
  int num_isochrones;
  
  struct cmd_res_data * chrone;

}iso;


#ifndef GRIDDEDCMD
extern double * xvals;
#endif

#ifdef __cplusplus
}
#endif 

#endif /*  __CMD_H__ */

---