Radiative losses H, Ca II, Mg II database
=========================================

Calculated by the MALI code by Petr Heinzel.

The results are in the directory 'Radiative_Losses'.
This directory contains 6 subdirectories VAL_A - VAL_F according to the initial models.

Each subdirectory contains 2989 files 'RADLOS?????.DAT' (00000-02988), 
corresponding to the scaled models. The ordering of the RADLOS files
is the same as for the model files (see 'model_files.txt').

Each file has 103 rows (103 depths z) and 6 columns:
   -----------------------------
   id, z, col_mass, LH, LCa, LMg  [cgs units]
   -----------------------------
LH are the radiative losses at the depth z for hydrogen, LCa for Ca II, and LMg for Mg II.
Total radiative losses = LH + LCa + LMg

LH includes Lyman and Balmer lines and hydrogen continua,
LCa includes the Ca II H & K lines and the Ca II infrared triplet,
LMg includes the Mg II h & k lines.

The units are [erg/sec/cm^3], divide them by 1000 to convert to [W/m^2].

 The conversion of z to geometrical height h (km, 0 at tau_500=1) is
     for i=0,102 do h[i] = (z[102]-7.5e6-z[i])/1.e5,
 because tau_500=1 is near 7.5e6 cm (75 km) in all models.

----------------------------------------------------------------
READRADLOS_ID.PRO
PRO readradlos_id,init_model,id,radtab,flag
; 
; Reading radiative losses from the database
; 
; INPUT:
;  init_model (string) in the form 'C' defines the folder name
;  id (integer) - ID number of the model in the folder 'Radiative_Losses/VAL_'+init_model
; OUTPUTS:
;  radtab - radiative-losses array DOUBLE(6,103)
;     [h,m,LH,LCa,LMg,Ltot] - height [km], col.mass [g/cm^2], losses for H, CaII, MgII,
;                             and total losses [erg/s/cm^3]
;  flag = 1 when the output is ok, flag = 0 when the output contains NaN
; Makes a plot of LH (red), LCa (violet), LMg (green), Ltot (black) vs. height

READRADLOS_TM.PRO
PRO readradlos_tm,init_model,tp,mp,radtab,flag
; 
; Reading radiative losses from the database
; 
; INPUT:
;  init_model (string) in the form 'C' defines the folder name
;  tp - temperature parameter in the range -0.03 to 0.03
;  mp - column-mass parameter in the range 0.2 to 5
; OUTPUTS:
;  radtab - radiative-losses array DOUBLE(6,103)
;     [h,m,LH,LCa,LMg,Ltot] - height [km], col.mass [g/cm^2], losses for H, CaII, MgII,
;                             and total losses [erg/s/cm^3]
;  flag = 1 when the output is ok, flag = 0 when the output contains NaN
; Makes a plot of LH (red), LCa (violet), LMg (green), Ltot (black) vs. height