 //From AAG Cloudwatcher formula. Need to improve futher.
 //http://www.aagware.eu/aag/cloudwatcher700/WebHelp/index.htm#page=Operational%20Aspects/23-TemperatureFactor-.htm
 //https://azug.minpet.unibas.ch/wikiobsvermes/index.php/AAG_cloud_sensor#Snow_on_the_sky_temperature_sensor
+//Td  = (K1 / 100) * (Ta - K2 / 10) + (K3 / 100) * (Exp (K4 / 1000*Ta)) ^ (K5 / 100)
+//Tsky=Ts – Td 
 double skyTemp() {
   //Constant defined above
   double Td = (K1 / 100.) * (T - K2 / 10) + (K3 / 100.) * pow((exp (K4 / 1000.* T)) , (K5 / 100.));
-  double Tsky=IR-Td;
-  return Tsky;
+  double Tsky=IR-Td; // in °C
+  return Tsky; //in °C
 }
 
 

} else {
  frezzing=0;
}

#ifdef USE_WIND_SENSOR
  Rotations = 0; // Set Rotations count to 0 ready for calculations

               break;
       case 11:     
               result=(T22ext+273)*20;
               break;               
       case 12:     
               result=Hr22int*100;
               break;       
       case 13:     
               result=(T22int+273)*20;
               break;                  
       case 14:     
               result=(DewInt+273)*20;
               break; */    
               
       default:
