file: \data\2015\Hawaii-2015-01\doc\aux_descrip.txt date: 07-Apr-2015; by: MF/MLML (rev: 09Apr2015) what: Description of auxiliary data for MOBY-refresh Resonon Blue #1 Resonon sent us a "Moby Blue wiring chart" for the 2x 5-pin bulkhead connectors on the instrument. Connector 1, pin 1,2 = Slit Thermistor, pin 3,4 = Near-camera Thermistor, pin 5 = Unused. Connector 2, pin 1 = Humidity Sensor -ve (0 V), pin 2 = Humidity Sensor Signal Out(*), pin 3 = Humidity Sensor +ve (+3.3 V), pin 4,5 = Under-cap Thermistor. * Note: No load has been bridged to ground as suggested in the Honeywell Typical Application Circuit. The 5-pin connectors/cables are from SeaCon, Hummer series, HUMG-5-BCR (bulkhead) and HUM-5-CCP (cable). Thermistors are from US Sensors, type PT103J2: 10 kOhm @ 25 degC, +/- 0.2 degC (0 to 70 degC) "Under-cap Thermistor" is next to the humidity sensor. Humidity Sensor is from Honeywell, PN: HIH-5030, +/- 3%RH, -40 to 85 degC, 0 to 100 %RH, non-condensing environment, input voltage 2.7 to 5.5 Vdc, current draw 200 uAmp Typical Application Circuit: 65 kOhm minimum load between Vout and -Ve(0V) - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - For now, I am acquiring auxiliary data by reading an Agilent 34970A DVM. Typically, each channel is sampled 10x times per loop, on a continuous loop. Chan 1 = Slit Thermistor (Ohm) [US Sensors PT103J2] Chan 2 = Near-camera Thermistor (Ohm) [US Sensors PT103J2] Chan 3 = Under-cap Thermistor (Ohm) [US Sensors PT103J2] Chan 4 = Humidity Sensor supply Voltage (Vdc) Chan 5 = Humidity Sensor output Voltage (Vdc) Chan 6 = Ambient Thermistor (Ohm) [Radio Shack 271-110A] Humidity Sensor supply Voltage via Agilent E3645A (nominal 3.3 Vdc) MATLAB ex: >> [s1,a1,b1,s2,a2,b2]=readslm('auxdata.txt',3,2); >> x=s1; clear s1 a1 b1 s2 a2 b2, whos x % [720x14] >> plot(x(:,2),'x-'),ylabel('year'),zoom on,pause,close >> plot(x(:,3),'x-'),ylabel('jday'),zoom on,pause,close >> plot(x(:,4),'x-'),ylabel('hour'),zoom on,pause,close >> o1=[]; o2=[]; o3=[]; o4=[]; vi=[]; vo=[]; % Ohm therm#1,2,3,4; Volts RH in & out >> for i=1:size(x,1) >> if x(i,1)==1, o1=[o1; x(i,4), x(i,5:end)]; % [dec hr, Ohm data...] >> elseif x(i,1)==2, o2=[o2; x(i,4), x(i,5:end)]; % [dec hr, Ohm data...] >> elseif x(i,1)==3, o3=[o3; x(i,4), x(i,5:end)]; % [dec hr, Ohm data...] >> elseif x(i,1)==4, vi=[vi; x(i,4), x(i,5:end)]; % [dec hr, Volt data...] >> elseif x(i,1)==5, vo=[vo; x(i,4), x(i,5:end)]; % [dec hr, Volt data...] >> elseif x(i,1)==6, o4=[o4; x(i,4), x(i,5:end)]; % [dec hr, Ohm data...] >> end >> end, clear x i ans - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Specifications for Honeywell HIH-5030, Low Voltage Humidity Sensors http://sensing.honeywell.com/index.php?ci_id=49692 V_out = (V_supply)(0.00636(sensor RH) + 0.1515), typical at 25 degC True RH = (Sensor RH)/(1.0546 - 0.00216T), T in degC or, for MATLAB, sensor_RH = (V_out/V_supply - 0.1515)/0.00636 True_RH = Sensor_RH/(1.0546 - 0.00216*sensor_T) sensor_T = degC - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - U.S. Sensor 10k Ohm NTC thermistor, PT103J2, Resistance vs Temperature http://www.ussensor.com/standard-precision-interchangeable-thermistors-pt103j2-table PT103J2 REV NONE (RESISTANCE VS. TEMPERATURE TABLE) Column #1: Temp (°C), Column #2: Resistance (Ω) 10,000 Ω @ +25°C Steinhart-Hart equation: 1/T = C1 + C2*(lnR) + C3*(lnR)^3 Temperature (T) is in Kelvin, and resistance (R) is in ohms (0 degC == 273.15 K) C1=1.12485e-3; C2=2.34793e-4; C3=0.85453e-7; 3rd order regress coef's (...for Labview) X=natural_logarithm( Ohm ) T=A + B*X + C*X^2 + D*X^3 ( degC ) A=508.26; B=-99.7397; C=7.0545; D=-0.20863; - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Radio Shack 10k-Ohm Thermistor, PN: 271-110A 10k ohms +/- 1% at 25 degC Steinhart-Hart equation: 1/T = C1 + C2*(lnR) + C3*(lnR)^3 Temperature (T) is in Kelvin, and resistance (R) is in ohms (0 degC == 273.15 K) C1=0.901823e-3; C2=2.48641e-4; C3=2.06308e-7; 3rd order regress coef's (...for Labview) X=natural_logarithm( Ohm ) T=A + B*X + C*X^2 + D*X^3 ( degC ) A=551.92; B=-99.654; C=5.9139; D=-0.14178; - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - MATLAB ex: >> % Convert thermistor (Ohm) to temperature (degC) >> t1=mythermcnv(o1,'\wrk\thermist\thrm_USS.dat'); t1(:,1)=o1(:,1); >> t2=mythermcnv(o2,'\wrk\thermist\thrm_USS.dat'); t2(:,1)=o2(:,1); >> t3=mythermcnv(o3,'\wrk\thermist\thrm_USS.dat'); t3(:,1)=o3(:,1); >> t4=mythermcnv(o4,'\wrk\thermist\thrm_RS.dat'); t4(:,1)=o4(:,1); whos >> >> fm; % plot degC >> x=t1; errorbar(x(:,1),nanmean(x(:,2:end)')',nanstd(x(:,2:end)')' ,'b'); clear x, hold on >> x=t2; errorbar(x(:,1),nanmean(x(:,2:end)')',nanstd(x(:,2:end)')' ,'g'); clear x, >> x=t3; errorbar(x(:,1),nanmean(x(:,2:end)')',nanstd(x(:,2:end)')' ,'r'); clear x, >> x=t4; errorbar(x(:,1),nanmean(x(:,2:end)')',nanstd(x(:,2:end)')' ,'k'); clear x >> ylabel('degC');xlabel('dec hr'); legend('T1 = Slit','T2 = Near-camera','T3 = Under-cap','T4 = Ambient',4) >> zoom on,pause,close