Index: issm/trunk-jpl/src/c/classes/Elements/Element.h
===================================================================
--- issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 24207)
+++ issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 24208)
@@ -219,6 +219,4 @@
 		virtual void       ComputeStressTensor(void)=0;
 		virtual void       ComputeEsaStrainAndVorticity(void)=0;
-	   virtual void		 ComputeRadarAttenuation(void){_error_("not implemented yet");}; 
-		virtual void       ComputeRadarPower(void)=0;		
 		virtual void       Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters)=0;
 		virtual void       ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N,int M)=0;
Index: issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 24207)
+++ issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 24208)
@@ -784,197 +784,4 @@
 	/*Clean up and return*/
 	delete gauss;
-}
-/*}}}*/
-void       Penta::ComputeRadarAttenuation(){/*{{{*/
-   IssmDouble  xyz_list[NUMVERTICES][3];
-	IssmDouble  eps0=8.85418782e-12;
-	IssmDouble  eps_ice=3.17;
-	IssmDouble  e=2.7183;
-	IssmDouble  c=299792458;
-	IssmDouble  k=1.3806488e-23;
-	IssmDouble  eVtoJ=1.6e-19;
-	IssmDouble  Tr_M07=252.1500;
-	IssmDouble  Tr_W97=258.1500; 
-	IssmDouble  sig_ice_M07=9.2;
-	IssmDouble	sig_ice_W97=9; 
-	IssmDouble  cond_H_M07=3.2; 
-	IssmDouble  cond_H_W97=4; 
-	IssmDouble  cond_Cl_M07=0.43; 
-	IssmDouble  cond_Cl_W97=0.55; 
-	IssmDouble  cond_NH_M07=0.8; 
-	IssmDouble  cond_NH_W97=1; 
-	IssmDouble  E_M07=8.1600e-20;
-	IssmDouble  E_W97=9.2800e-20; 
-	IssmDouble  E_H_M07=3.2000e-20;
-	IssmDouble  E_H_W97=3.3600e-20;
-	IssmDouble  E_Cl_M07=3.0400e-20;
-	IssmDouble  E_Cl_W97=3.6800e-20; 
-	IssmDouble  E_NH=3.6800e-20;
-	/*IssmDouble  mol_H_hol=1.6; 
-	IssmDouble  mol_Cl_hol=0.4; 
-	IssmDouble  mol_NH_hol=0.5; 
-	IssmDouble  mol_H_lgp=0.2; 
-	IssmDouble  mol_Cl_lgp=1.8; 
-	IssmDouble  mol_NH_lgp=0.4; */
-	IssmDouble  radar_ice_period; 
-	IssmDouble  mol_H, mol_Cl, mol_NH; 
-	IssmDouble  attenuation_rate_macgregor[NUMVERTICES];
-	IssmDouble  attenuation_rate_wolff[NUMVERTICES]; 
-	IssmDouble  temperature, attenuation_rate_M07_pureice, attenuation_rate_M07_H, attenuation_rate_M07_Cl, attenuation_rate_M07_NH;
-	IssmDouble	attenuation_rate_W97_pureice, attenuation_rate_W97_H, attenuation_rate_W97_Cl, attenuation_rate_W97_NH;
-	IssmDouble m1, m2, m3, m4, m5, m6, m7, m8;
-	IssmDouble  w2, w3, w4, w5, w6, w7, w8; 
-	GaussPenta* gauss=NULL; 
-
-   /* Get node coordinates and dof list: */
-	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-
-	/*Retrieve all inputs we will be needing: */
-	Input* temp_input=inputs->GetInput(TemperatureEnum); _assert_(temp_input);
-	//Input* ice_period_input=inputs->GetInput(RadarIcePeriodEnum); _assert_(ice_period_input); 
-	
-	/* Start looping on the number of vertices: */
-	gauss=new GaussPenta();
-	for (int iv=0;iv<NUMVERTICES;iv++){
-		gauss->GaussVertex(iv);
-
-   /*Get ice temperature: */
-		temp_input->GetInputValue(&temperature,gauss);
-//		ice_period_input->GetInputValue(&radar_ice_period,gauss); 
-
-	/*Ice period condition: */
-	//	if (radar_ice_period>0.){
-			mol_H=1.6;
-			mol_Cl=0.4;
-			mol_NH=0.5;
-	//	}
-	//	else{
-	//		mol_H=0.2; 
-	//		mol_Cl=1.8;
-	//		mol_NH=0.4;
-	//	}
-
-	/*Compute M07 radar conductivity constant: */
-		m1=(10*log10(e))/(1000*eps0*sqrt(eps_ice)*c);
-		m2=E_M07/k;
-		m3=cond_H_M07*mol_H;
-		m4=E_H_M07/k;
-		m5=cond_Cl_M07*mol_Cl;
-		m6=E_Cl_M07/k;
-		m7=cond_NH_M07*mol_NH;
-		m8=E_NH/k;
-
-	/*Compute MacGregor attenuation rate: */
-		attenuation_rate_M07_pureice=m1*sig_ice_M07*exp(m2*((1/Tr_M07)-(1/temperature))); 
-		attenuation_rate_M07_H=m3*exp(m4*((1/Tr_M07)-(1/temperature)));
-		attenuation_rate_M07_Cl=m5*exp(m6*((1/Tr_M07)-(1/temperature)));
-		attenuation_rate_M07_NH=m7*exp(m8*((1/Tr_M07)-(1/temperature)));
-		attenuation_rate_macgregor[iv]=attenuation_rate_M07_pureice+attenuation_rate_M07_H+attenuation_rate_M07_Cl+attenuation_rate_M07_NH; 
-
-	/*Compute W97 radar conductivity constant: */ 
-		w2=E_W97/k;
-		w3=cond_H_W97*mol_H;
-		w4=E_H_W97/k;
-		w5=cond_Cl_W97*mol_Cl;
-		w6=E_Cl_W97/k;
-		w7=cond_NH_W97*mol_NH;
-		w8=E_NH/k;
-
-	/*Compute Wolff attenuation rate: */
-    attenuation_rate_W97_pureice=m1*sig_ice_W97*exp(w2*((1/Tr_W97)-(1/temperature)));
-	 attenuation_rate_W97_H=w3*exp(w4*((1/Tr_W97)-(1/temperature)));
-	 attenuation_rate_W97_Cl=w5*exp(w6*((1/Tr_W97)-(1/temperature)));
-	 attenuation_rate_W97_NH=w7*exp(w8*((1/Tr_W97)-(1/temperature)));
-	 attenuation_rate_wolff[iv]=attenuation_rate_W97_pureice+attenuation_rate_W97_H+attenuation_rate_W97_Cl+attenuation_rate_W97_NH;
-	}
-
-   /*Add Attenuation rate results into inputs*/
-	this->inputs->AddInput(new PentaInput(RadarAttenuationMacGregorEnum,&attenuation_rate_macgregor[0],P1Enum));
-	this->inputs->AddInput(new PentaInput(RadarAttenuationWolffEnum,&attenuation_rate_wolff[0],P1Enum));
-
-   /*Clean up and return*/
-	delete gauss;
-}
-/*}}}*/
-void       Penta::ComputeRadarPower(){/*{{{*/   
-	IssmDouble  xyz_list[NUMVERTICES][3];
-	IssmDouble  power_M07[NUMVERTICES]; 
-	IssmDouble  power_W97[NUMVERTICES]; 
-	IssmDouble	depth[NUMVERTICES];
-	IssmDouble  aircraft_elev=0.5; 
-	IssmDouble  eps_ice=3.15;
-	IssmDouble  t_tp=273.15;			/* triple point temperature [K] */
-	IssmDouble  p_tp=611.73;			/* water pressure [Pa] */
-	IssmDouble  gamma=7.4200e-07;	/* Clausius-Clapeyron constant [K/kPa] */
-	IssmDouble  attenuation_rate_macgregor, attenuation_rate_wolff, attenuation_total_M07, attenuation_total_W97; 
-	IssmDouble  thickness, surface, z, temperature, geometric_loss, reflectivity; 
-	IssmDouble  rho_ice, gravity, pressure, pressure_melting_pt, frozen_temp, basal_temp, basal_pmp;  
-	GaussPenta* gauss=NULL;
-
-   /* Get node coordinates and dof list: */
-	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-
-	/*Retrieve all inputs we will be needing: */
-   Input* atten_input_M07=inputs->GetInput(RadarAttenuationMacGregorEnum); _assert_(atten_input_M07);
-	Input* atten_input_W97=inputs->GetInput(RadarAttenuationWolffEnum); _assert_(atten_input_W97);
-	Input* surf_input=inputs->GetInput(SurfaceEnum); _assert_(surf_input);
-	Input* thick_input=inputs->GetInput(ThicknessEnum); _assert_(thick_input);
-	Input* temp_input=inputs->GetInput(TemperatureEnum); _assert_(temp_input);
-
-	/* Start looping on the number of vertices: */
-   gauss=new GaussPenta();
-	for (int iv=0;iv<NUMVERTICES;iv++){
-		gauss->GaussVertex(iv);
-	
-	/*Get all the inputs: */
-		atten_input_M07->GetInputValue(&attenuation_rate_macgregor,gauss);
-		atten_input_W97->GetInputValue(&attenuation_rate_wolff,gauss);
-		thick_input->GetInputValue(&thickness,gauss);
-		temp_input->GetInputValue(&temperature,gauss);
-		surf_input->GetInputValue(&surface,gauss);
-
-	/*Compute depth below the ice surface: */
-		z=xyz_list[iv][2];
-		depth[iv]=(surface-z)/1e3;
-	
-	/*Compute total attenuation: */ 
-		attenuation_total_M07=attenuation_rate_macgregor*depth[iv]; 
-		attenuation_total_W97=attenuation_rate_wolff*depth[iv];
-	
-	/*Compute geometric loss: */ 
-		geometric_loss=10*log10((depth[iv]+aircraft_elev)/sqrt(eps_ice)); 
-	
-	/*Compute radar power: */
-		power_M07[iv]=-geometric_loss-attenuation_total_M07;
-		power_W97[iv]=-geometric_loss-attenuation_total_W97; 
-	
-		/*Identify basal elements: */
-		if(this->IsOnBase() && iv<NUMVERTICES2D){
-
-			/*Compute pressure melting point: */ 
-			rho_ice=FindParam(MaterialsRhoIceEnum);
-			gravity=FindParam(ConstantsGEnum);
-			pressure=rho_ice*gravity*thickness; 	
-			pressure_melting_pt=t_tp-gamma*(pressure-p_tp); 
-			if((temperature-pressure_melting_pt)<=-1){
-				reflectivity=-40;
-			}
-			else if((temperature-pressure_melting_pt)>-1 && (temperature-pressure_melting_pt)<0){ 
-				reflectivity=0;
-			}
-			else{
-				reflectivity=70;
-			}
-			power_M07[iv]=power_M07[iv]+reflectivity;
-			power_W97[iv]=power_W97[iv]+reflectivity;
-		} 
-	} 
-	 
-	 /*Add power results into inputs*/
-		this->inputs->AddInput(new PentaInput(RadarPowerMacGregorEnum,&power_M07[0],P1Enum));
-		this->inputs->AddInput(new PentaInput(RadarPowerWolffEnum,&power_W97[0],P1Enum));
-
-		/*Clean up and return*/
-		delete gauss;
 }
 /*}}}*/
Index: issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 24207)
+++ issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 24208)
@@ -58,6 +58,4 @@
 		void           ComputeSigmaNN(){_error_("not implemented yet");};
 		void           ComputeStressTensor();
-		void				ComputeRadarAttenuation();
-		void           ComputeRadarPower(); 
 		void           Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters);
 		void           ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N,int M);
Index: issm/trunk-jpl/src/c/classes/Elements/Seg.h
===================================================================
--- issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 24207)
+++ issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 24208)
@@ -49,6 +49,4 @@
 		void        ComputeSigmaNN(){_error_("not implemented yet");};
 		void        ComputeStressTensor(){_error_("not implemented yet");};
-		void        ComputeRadarAttenuation(){_error_("not implemented yet");};
-		void        ComputeRadarPower(){_error_("not implemented yet");};
 		void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters){_error_("not implemented yet");};
 		void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N,int M){_error_("not implemented yet");};
Index: issm/trunk-jpl/src/c/classes/Elements/Tetra.h
===================================================================
--- issm/trunk-jpl/src/c/classes/Elements/Tetra.h	(revision 24207)
+++ issm/trunk-jpl/src/c/classes/Elements/Tetra.h	(revision 24208)
@@ -49,6 +49,4 @@
 		void        ComputeDeviatoricStressTensor(){_error_("not implemented yet");};
 		void        ComputeEsaStrainAndVorticity(){_error_("not implemented yet!");};
-		void        ComputeRadarAttenuation(){_error_("not implemented yet!");};
-		void        ComputeRadarPower(){_error_("not implemented yet");};
 		void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters);
 		void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N,int M){_error_("not implemented yet");};
Index: issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 24207)
+++ issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 24208)
@@ -63,6 +63,4 @@
 		void        ComputeStressTensor();
 		void        ComputeSurfaceNormalVelocity();
-		void        ComputeRadarAttenuation(){_error_("not implemented yet");};
-		void        ComputeRadarPower(){_error_("not implemented yet");};		
 		void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters);
 		void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int N, int M);
