the "raw" data file has data exactly as I got them from MRS..the only
thing I added to these data was a time stamp in the 1st several
columns of the data file (note: the time stamp is in MST and
represents the center time of the averaging time period)...the
"mrs_c1_2004.dat" data file includes some minor data cleanup and two
extra parameters (the last 2 columns)...which are: precip_cumsum
(cumulative precipitation, starting at zero on January 1), and
precip_hourly (hourly precipitation total).

7/22/04:

I put the MRS C1 hourly data on toyon as:


  1562144 Jul 22 11:43 /projects/niwot_ridge/data_c1/mrs_c1_2002.mat
  1561976 Jul 22 11:43 /projects/niwot_ridge/data_c1/mrs_c1_2003.mat
  1562368 Jul 22 11:43 /projects/niwot_ridge/data_c1/mrs_c1_2004.mat


these will give you some basic met parameters to go along with the
sodar data.  The parameters available are:


s_Pbaro_c1            s_jd_c1               s_wd_avg_c1           
s_RH_c1               s_jd_mst_c1           s_wd_std_c1           
s_Rnet_sw_c1          s_nadp_c1             s_ws_max_c1           
s_Tair_c1             s_precip_max_c1       s_ws_max_time_c1      
s_Tpanel_c1           s_precip_max_time_c1  s_ws_scalar_c1        
s_Tsoil_c1            s_qsoil_c1            s_ws_vector_c1        
s_battery_c1          s_snow_max_c1         s_time_ymdhms_mst_c1
s_fileID_c1           s_snow_max_time_c1


where,


             s_jd_c1: Day of year (whole number)
         s_jd_mst_c1: Decimal Day of year (center of averaging time period, 
                      in Mountain Standard Time)
s_time_ymdhms_mst_c1: a Nx6 matrix...[Year Month Day Hour Minute Second]...example:


 s_time_ymdhms_mst_c1(1:10,:)


        2004           1           1           0          30           0
        2004           1           1           1          30           0
        2004           1           1           2          30           0
        2004           1           1           3          30           0
        2004           1           1           4          30           0
        2004           1           1           5          30           0
        2004           1           1           6          30           0
        2004           1           1           7          30           0
        2004           1           1           8          30           0
        2004           1           1           9          30           0



       s_Tpanel_c1:  average panel temperature (degC)
         s_Tair_c1:  average air temperature (i think at 5m or 10m?)
           s_RH_c1:  average relative humidity (percent)
        s_Pbaro_c1:  average barometric pressure (mbar)
    s_ws_scalar_c1:  scalar wind speed average (m/s)
    s_ws_vector_c1:  vector speed average (m/s)
       s_wd_avg_c1:  average wind direction (direction wind 
                     is coming FROM, deg from N, 0-360)
       s_wd_std_c1:  standard deviation of direction (deg)
       s_ws_max_c1:  maximum wind speed (m/s)
  s_ws_max_time_c1:  time of maximum Wind Speed
      s_Rnet_sw_c1:  total solar radiation (MegaJoules/m2)
        s_Tsoil_c1:  average soil temperature (C)
        s_qsoil_c1:  average soil volumetric moisture
      s_battery_c1:  average battery voltage
         s_nadp_c1:  average NADP (not sure what this is)
   s_precip_max_c1:  maximum precipitation (i think this is cumulative 
                     max precipitation in cm....might need to check on this)
s_precip_max_time_c1:  time of maximum precipitation
     s_snow_max_c1: maximum snow depth
s_snow_max_time_c1: time of maximum snow depth



these are hourly averaged data with the avg period at the top of each
hour...eg, avgs are between 1-2, 2-3, 3-4, etc.


if you want to plot the air temperature in matlab you could do:


  load /projects/niwot_ridge/data_c1/mrs_c1_2004.mat
  plot(s_jd_mst_c1,s_Tair_c1)  


What i put here are data for the entire year...if you are only
interested in data between 05/17 - 07/08 (which is between days
138-190) then do:


   r=find(s_jd_mst_c1 >= 138 & s_jd_mst_c1 <= 191);
   plot(s_jd_mst_c1(r),s_Tair_c1(r))
   set(gca,'xlim',[min(s_jd_mst_c1(r)) max(s_jd_mst_c1(r))])
   grid on;


if you have any other questions let me know...