Added a basic model for simulating the Greenland Attenuation Lengths.…

AraSim.cc

@@ -127,7 +127,6 @@ int main(int argc, char **argv) {   // read setup.txt file
     // IceModel *icemodel=new IceModel(ICE_MODEL + NOFZ*10,CONSTANTICETHICKNESS * 1000 + CONSTANTCRUST * 100 + FIXEDELEVATION * 10 + 0,MOOREBAY);// creates Antarctica ice model
     IceModel *icemodel=new IceModel(settings1->ICE_MODEL + settings1->NOFZ*10,settings1->CONSTANTICETHICKNESS * 1000 + settings1->CONSTANTCRUST * 100 + settings1->FIXEDELEVATION * 10 + 0,settings1->MOOREBAY);// creates Antarctica ice model
     // IceModel inherits from EarthModel  
-
     cout<<endl;
     cout<<"Surface at (log:0, lat:0) : "<<icemodel->Surface(0., 0.)<<endl;
     cout<<"SurfaceAboveGeoid at (log:0, lat:0) : "<<icemodel->SurfaceAboveGeoid(0., 0.)<<endl;

IceModel.cc

@@ -171,6 +171,31 @@ void IceModel::setUpIceModel(int model) {
       i++;
     }
   westlanddown.close();
+
+    //Add Greenland attenuation lengths
+  ifstream greenlandUpper("data/greenland_attenlength_bogorodsky_upper.txt");
+  if(greenlandUpper.fail())
+    {cerr << "Failed to open greenland_attenlength_bogorodsky_upper.txt";
+      exit(1);
+    }
+  i=0;
+  while(greenlandUpper>>d_greenlandUpper[i]>>l_greenlandUpper[i])
+    {
+      i++;
+    }
+  greenlandUpper.close();
+
+  ifstream greenlandLower("data/greenland_attenlength_bogorodsky_lower.txt");
+  if(greenlandLower.fail())
+    {cerr << "Failed to open greenland_attenlength_bogorodsky_lower.txt";
+      exit(1);
+    }
+  i=0;
+  while(greenlandLower>>d_greenlandLower[i]>>l_greenlandLower[i])
+    {
+      i++;
+    }
+  greenlandLower.close();
 
 
   // new ARA ice attenuation measurement values (at 300 MHz)
@@ -1378,6 +1403,7 @@ double IceModel::EffectiveAttenuationLength(const Position &pos,const int &which
 	    cerr << " wrong attenuation length " <<endl;
       	}
 
+
       if(mooreBayFlag)//if use Moore's Bay measured data for the west land
 	attenuation_length*=1.717557; //about 450 m (field attenuation length) for one whole way when assuming -3dB for the power loss at the bottom
     }
@@ -1399,6 +1425,7 @@ double IceModel::EffectiveAttenuationLength(const Position &pos,const int &which
 	 cerr << " wrong attenuation length " <<endl;
      } //else
 
+
   return attenuation_length;
 } //EffectiveAttenuationLengthUp
 
@@ -1441,6 +1468,13 @@ double IceModel::EffectiveAttenuationLength(Settings *settings1, const Position
 	  else
 	    cerr << " wrong attenuation length " <<endl;
       	}
+      //Adding Greenland Ice attenuation lengths
+      if(settings1->GREENLAND_ATTEN == 1) {
+        attenuation_length=l_greenlandUpper[depth_index]; 
+      }
+      else if(settings1->GREENLAND_ATTEN == 2) {
+        attenuation_length=l_greenlandLower[depth_index]; 
+      }
 
       //if(mooreBayFlag)//if use Moore's Bay measured data for the west land
       if(settings1->MOOREBAY)//if use Moore's Bay measured data for the west land
@@ -1463,6 +1497,13 @@ double IceModel::EffectiveAttenuationLength(Settings *settings1, const Position
        else
 	 cerr << " wrong attenuation length " <<endl;
      } //else
+      //Adding Greenland Ice attenuation lengths
+      if(settings1->GREENLAND_ATTEN == 1) {
+        attenuation_length=l_greenlandUpper[depth_index]; 
+      }
+      else if(settings1->GREENLAND_ATTEN == 2) {
+        attenuation_length=l_greenlandLower[depth_index]; 
+      }
 
   return attenuation_length;
 } //EffectiveAttenuationLengthUp

IceModel.h

@@ -94,7 +94,11 @@ void ENtoLonLat(int e_coord,
   const static int N_westlanddown=420;
   double d_westlanddown[N_westlanddown],l_westlanddown[N_westlanddown];
 
-
+  // Attenuation lengths of Greenland Ice using the Bogorodsky model from Figure 7 of doi.org/10.1017/jog.2022.40
+  const static int N_greenlandUpper=2810;
+  double d_greenlandUpper[N_greenlandUpper], l_greenlandUpper[N_greenlandUpper];
+  const static int N_greenlandLower=2810;
+  double d_greenlandLower[N_greenlandLower], l_greenlandLower[N_greenlandLower];
 
 public:
 

Settings.cc

@@ -293,6 +293,8 @@ outputdir="outputs"; // directory where outputs go
     APPLY_NOISE_FIGURE=0; // default: 0 - don't use new noise figure information
 
     CUSTOM_ELECTRONICS=0; //default: 0 -- don't use custom electronics, load regular "ARA_Electronics_TotalGain_TwoFilter.tst"
+
+    GREENLAND_ATTEN=0; //default:0 -- Uses South Pole attenuation length.  1 & 2 use Greenland attenuation lengths via RNO-G.
 
 
     /*

Settings.h

@@ -300,6 +300,10 @@ class Settings
 
 	int CUSTOM_ELECTRONICS; //0 (default): use the regular "ARA_Electronics_TotalGain_TwoFilter.txt" file
 							//1 : load a custom electronics file, stored as "custom_electronics.txt" in the `data` directory
+
+    int GREENLAND_ATTEN; // 0 (default): Use regular Ara attenuation lengths
+                         // 1:  Use upper limit of bogorodsky curve fron RNO-G
+                         // 2:  Use lower limit of bogorodsky curve fron RNO-G
 
 
 //arrays for saving read in event features in EVENT_GENERATION_MODE=1