#include "fortranfunctions.h"

#include <iostream>
#include<cstdlib>

extern "C" void calc_limb_refraction (int* isl, int* iel, float* pout,  double* appang,
	double *trueang, float* za) 
    
{

	static int ic1= -1;
	static int ic2 = 1;
	for(int ir = *isl; ir <= *iel; ++ir) {
		trueang[ir-1] = appang[ir-1] + ::refraction_( &ir, &ic1, &ic2, pout );
	} 
}