diff --git a/radiation.c b/radiation.c index 1a1631f31..64e02b562 100755 --- a/radiation.c +++ b/radiation.c @@ -166,13 +166,15 @@ radiation (p, ds) freq = 0.5 * (freq_inner + freq_outer); - + + /* calculate free-free, compton and ind-compton opacities + note that we also call these with the average frequency along ds */ kappa_tot = frac_ff = kappa_ff (xplasma, freq); /* Add ff opacity */ kappa_tot += frac_comp = kappa_comp (xplasma, freq); /* 70 NSH 1108 calculate compton opacity, store it in kappa_comp and also add it to kappa_tot, the total opacity for the photon path */ kappa_tot += frac_ind_comp = kappa_ind_comp (xplasma, freq); frac_tot = frac_z = 0; /* 59a - ksl - Moved this line out of loop to avoid warning, but notes - indicate this is all disagnostic and might be removed */ + indicate this is all disagnostic and might be removed */ if (freq > phot_freq_min) diff --git a/version.h b/version.h index 261f331a3..ab3c3210d 100755 --- a/version.h +++ b/version.h @@ -1,2 +1,2 @@ -#define VERSION "ff" +#define VERSION "77a_dev" #define CHOICE 1 // Compress plasma as much as possible