#include "radio.hxx"
#include <simgear/scene/material/mat.hxx>
#include <Scenery/scenery.hxx>
+#include <boost/scoped_array.hpp>
#define WITH_POINT_TO_POINT 1
#include "itm.cpp"
double FGRadioTransmission::ITM_calculate_attenuation(SGGeod pos, double freq, int transmission_type) {
-
+ if((freq < 40.0) || (freq > 20000.0)) // frequency out of recommended range
+ return -1;
/** ITM default parameters
TODO: take them from tile materials (especially for sea)?
**/
int max_points = (int)floor(distance_m / point_distance);
- double delta_last = fmod(distance_m, point_distance);
+ //double delta_last = fmod(distance_m, point_distance);
deque<double> elevations;
deque<string*> materials;
while (elevations.size() <= e_size) {
probe_distance += point_distance;
SGGeod probe = SGGeod::fromGeoc(center.advanceRadM( course, probe_distance ));
- const SGMaterial *mat = 0;
+ const simgear::BVHMaterial *material = 0;
double elevation_m = 0.0;
- if (scenery->get_elevation_m( probe, elevation_m, &mat )) {
+ if (scenery->get_elevation_m( probe, elevation_m, &material )) {
+ const SGMaterial *mat;
+ mat = dynamic_cast<const SGMaterial*>(material);
if((transmission_type == 3) || (transmission_type == 4)) {
elevations.push_back(elevation_m);
if(mat) {
elevations.push_front(num_points -1);
int size = elevations.size();
- double *itm_elev;
- itm_elev = new double[size];
+ boost::scoped_array<double> itm_elev( new double[size] );
for(int i=0;i<size;i++) {
itm_elev[i]=elevations[i];
if((transmission_type == 3) || (transmission_type == 4)) {
// the sender and receiver roles are switched
- ITM::point_to_point(itm_elev, receiver_height, transmitter_height,
+ ITM::point_to_point(itm_elev.get(), receiver_height, transmitter_height,
eps_dielect, sgm_conductivity, eno, frq_mhz, radio_climate,
pol, conf, rel, dbloss, strmode, p_mode, horizons, errnum);
if( _root_node->getBoolValue( "use-clutter-attenuation", false ) )
- calculate_clutter_loss(frq_mhz, itm_elev, materials, receiver_height, transmitter_height, p_mode, horizons, clutter_loss);
+ calculate_clutter_loss(frq_mhz, itm_elev.get(), materials, receiver_height, transmitter_height, p_mode, horizons, clutter_loss);
}
else {
- ITM::point_to_point(itm_elev, transmitter_height, receiver_height,
+ ITM::point_to_point(itm_elev.get(), transmitter_height, receiver_height,
eps_dielect, sgm_conductivity, eno, frq_mhz, radio_climate,
pol, conf, rel, dbloss, strmode, p_mode, horizons, errnum);
if( _root_node->getBoolValue( "use-clutter-attenuation", false ) )
- calculate_clutter_loss(frq_mhz, itm_elev, materials, transmitter_height, receiver_height, p_mode, horizons, clutter_loss);
+ calculate_clutter_loss(frq_mhz, itm_elev.get(), materials, transmitter_height, receiver_height, p_mode, horizons, clutter_loss);
}
double pol_loss = 0.0;
//_root_node->setDoubleValue("station[0]/tx-pattern-gain", tx_pattern_gain);
//_root_node->setDoubleValue("station[0]/rx-pattern-gain", rx_pattern_gain);
- delete[] itm_elev;
for (unsigned i =0; i < materials.size(); i++) {
delete materials[i];
}