src/stellar_spectrum.c

   1 /***************************************************************************
   2  * This File is a part of the CADS/UV Stellar Spectrum model software      *
   3  *   Copyright (C) 2012 by CADS/UV Software Team,                          *
   4  *                         Indian Institute of Astrophysics                *
   5  *                         Bangalore 560034                                *
   6  *                         cads_AT_iiap.res.in                             *
   7  *                                                                         *
   8  *   This program is free software; you can redistribute it and/or modify  *
   9  *   it under the terms of the GNU General Public License as published by  *
  10  *   the Free Software Foundation; either version 2 of the License, or     *
  11  *   (at your option) any later version.                                   *
  12  *                                                                         *
  13  *   This program is distributed in the hope that it will be useful,       *
  14  *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
  15  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
  16  *   GNU General Public License for more details.                          *
  17  *                                                                         *
  18  *   You should have received a copy of the GNU General Public License     *
  19  *   along with this program; if not, write to the                         *
  20  *   Free Software Foundation, Inc.,                                       *
  21  *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
  22  ***************************************************************************/
  23
  24 /*--------------------------------------------------------------------------
  25 File          : mag2flux.c
  26 Author        : Sujatha N V
  27 Date          : 12th August, 2004
  28 Version       : 1.0
  29
  30 Program Description:
  31 This program reads sp_type, magnitude & E(B-V) from an input
  32 parameter file & compares the spectral type in the  TABLE in file :
  33 MKSpTypewwww.in. Selects the corresponding KURUCZ flux file
  34 coresponding to Temp. T and gravity g of the sp. type. Then
  35 calculate scale value for Kurucz flux and calculate the flux at
  36 Earth for different wavelengths read from 'wave.dat'
  37 output file contains: Wavelength, Flux at Earth F(lambda)
  38
  39 ---------------------------------------------------------------------------*/
  40
  41 /*--------------------------------------------------------------------------
  42 Rivision            : Author, dd/mm/yyyy
  43 Revision history    :     -----
  44
  45 Reks, 12th July, 2005.
  46 1. Introduced signal handling (signal_handling.c and fun_prototypes.h)
  47 2. Modified all printf for error messages to fprintf->STDERR
  48    Those two modifications are a must to be able to run this on web...
  49
  50 Fayaz, Sept., 2005.
  51 1. Added UVS_msgs.c, to print error messages (supersedes #2 in previous
  52    revision) and other fancy things like a prologue() before starting.
  53
  54 Reks, 08th Dec. 2005
  55 1. Can take care of comments in input parameter file.
  56
  57 Reks, 07th July, 2007
  58 1. configure-make
  59 2. renamed project to mag2counts
  60 3. renamed file to mag2flux
  61 4. Better error messages and overall readability
  62 5. moved all functions into kurucz_functions.c
  63
  64 JM: 10 July, 2007
  65 1. Corrected several memory leaks in text handling
  66
  67 JM: 15 July, 2007
  68 1. Cleaned up code in main program; moved parameter reading etc to functions.
  69    (Still have to go through other functions.)
  70
  71 Reks: 21 July, 2007
  72 1. Removed SPDATA_FROOT from paramfile.
  73 2. SPDATA_DIR now points to the top level directory that contain all data
  74    files, subdirectories KURUCZ_FLUX and MkSpDATA
  75 3. shifted wave.dat and crossec0.dat into SPDATA_DIR (default location)
  76    experienced users may over-ride this from parameter file.
  77 4. added command line options to create default paramfile and help message
  78 5. Still have to fix the memory leak problems with sptype inside GET_REF1!
  79
  80 JM: 21 July, 2007
  81 1. fixed all memory leaks!
  82
  83 Reks: 01 July 2012:
  84 1. Rebranded to CADS/UVS, reset version to 1.0 in new cads git repository
  85 2. Renamed to stellar_spectrum
  86 3. read/write ops shifted to read_write.c
  87 ----------------------------------------------------------------------------*/
  88
  89 #include "stellar_spectrum.h"
  90 #include <sys/types.h>
  91
  92 /* Main program starts here*/
  93 /****************************************************************************/
  94 int main (int argc, char *argv[])
  95 {
  96
  97     FILE            *out_ptr = NULL, *wave_ptr = NULL;
  98     struct STARS     stars;
  99     struct S_TABLE   mkstTG;
 100     int              gr = 0, tst;
 101     char             OutFile[MAX_TEXT];
 102     char             ParamFile[80];
 103     int c;           /* For getopt */
 104
 105     /* Before we begin... set the signals and print a banner */
 106     set_signals();
 107     PRINT_BANNER();
 108
 109     strcpy(ParamFile, paramfile);/*paramfile is defined in header file*/
 110
 111     /* Command line options*/
 112     while ( (c = my_getopt (argc, argv, "gh") ) != EOF)
 113         switch (c)
 114         {
 115         case 'g':
 116             gen_stspec_params(ParamFile);
 117             return(EXIT_SUCCESS);
 118             break;
 119         case 'h':
 120             usage();
 121             break;
 122         default:
 123             usage();
 124             break;
 125         }
 126
 127     /* Open and read the parameter file */
 128     READ_PARAMS(paramfile, &stars, OutFile);
 129     /* try opening the wave file */
 130     if ((wave_ptr = fopen(stars.WaveFile,"r")) == NULL)
 131     {
 132         fprintf(stderr, "ERROR  : Unable to open file %s\n", stars.WaveFile);
 133         exit(EXIT_FAILURE);
 134     }
 135     /*Output file*/
 136     if ((out_ptr = fopen(OutFile, "w")) == NULL)
 137     {
 138         fprintf(stderr, "ERROR: Unable to write to file %s\n", OutFile);
 139         exit(EXIT_FAILURE);
 140     }
 141
 142     WRITE_INFO(out_ptr, stars);/*Writes program information to file*/
 143
 144     /* for all wavelengths listed in this file */
 145     while ((feof(wave_ptr))==0)
 146     {
 147
 148         tst = GET_WAVE_NAME(wave_ptr, stars.dir_root, &mkstTG);
 149         if (tst == EXIT_SUCCESS)
 150         {
 151
 152             /* Read the wavelength file*/
 153             printf("\rWorking: |");
 154             READ_TABLE(&mkstTG);
 155
 156             /* get the crossection */
 157             printf("\rWorking: /");
 158             RCross_Sec(&stars, mkstTG.w_filename);
 159
 160             /*Reads ISM Absorption CrosSection*/
 161             printf("\rWorking: -");
 162             gr  = GET_REF1(&stars, &mkstTG);
 163
 164             printf("\rWorking: \\");
 165             if (gr == 0)
 166             {
 167                 ALPHA(&stars);
 168                 if (stars.alpha != 0.0) WRITE_LINE(out_ptr, &stars);
 169             }
 170
 171         }/*End tst*/
 172     } /* end of while loop */
 173     printf("\rInfo   : Done\n");
 174     fclose(wave_ptr);
 175     fclose(out_ptr);
 176
 177     return(EXIT_SUCCESS);
 178 }