getHbgrid {stellaR}

HB models from ZAHB to thermal pulses

Description

Import stellar track data for a grid of HB models starting from ZAHB and ending at the onset of thermal pulses.

Usage

 getHbgrid(z, y, ml, afe,
             baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
     

Arguments

z

the initial metallicity of the model to get. Allowed values are 0.0001, 0.0002, 0.0003, 0.0004, 0.0005, 0.0006, 0.0007, 0.0008, 0.0009, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01.

y

the initial helium abundance of the model to get. Allowed values are 0.25, 0.27, 0.33, 0.38, 0.42.

ml

the mixing-length of the model to get. Allowed values are 1.7, 1.8, 1.9.

afe

the alpha-enhanchment of the model to get. It can be one of afe = 0 for [alpha/Fe] = 0.0 or afe = 1 for [alpha/Fe] = 0.3.

baseURL

the URL of the base directory of the database from where get the models.

Details

As a sanity check, the function performs a consistency test on the parameters z, y, ml, afe by a call to testComposition.

By default the function requires an Internet access for data gathering. Tracks and isochrones are downloaded from the on-line repository CDS (http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26). It is also possible to gather data from a local download of the catalog, specifying as baseURL the root directory of the download. In this case the structure of the catalog and the names of files and directory must be preserved.

Value

getHbgrid returns an object of class hbset, i.e. a list with one element of class hb for each HB model.

If baseURL points to an invalid path or CDS ftp site is unavailable the function returns NA and issues a warning.

References

M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.

See Also

getTrk, getHb, getIso, getTrkSet.

Examples

  ### slow!
  ## Not run: 
  hbgrid <- getHbgrid(0.002, 0.25, 1.7, 0)

  ### get data from local directory /data
  hbgrid <- getHbgrid(0.002, 0.25, 1.7, 0, baseURL="/data/")
## End(Not run)