00001
00002
00003
00004
00005
00006
00007 #include "cddefines.h"
00008 #include "physconst.h"
00009 #include "hydro_bauman.h"
00010 #include "iso.h"
00011 #include "helike.h"
00012 #include "dense.h"
00013 #include "opacity.h"
00014 #include "atmdat.h"
00015 #include "taulines.h"
00016
00017
00018
00019
00020
00021
00022 double H_cross_section( double EgammaRyd , double EthRyd, long n, long l, long nelem )
00023 {
00024 double cs;
00025 double rel_photon_energy;
00026
00027
00028
00029
00030 rel_photon_energy = EgammaRyd / EthRyd;
00031 rel_photon_energy = MAX2( rel_photon_energy , 1. + FLT_EPSILON*2. );
00032
00033 cs = H_photo_cs(rel_photon_energy , n, l, nelem + 1 );
00034
00035 ASSERT( cs > 0. && cs < 1.E-8 );
00036
00037 return cs;
00038 }