ASKO PAR notes Ankur, The data in the file is in millivolts, to convert it to micromoles per m^2 per sec, use the following formula: if site='mhw' then do; v45=((v45/8.311)*199.20*1000)/(0.604*1000); end; if site='mrp' then do; v45=((v45/8.311)*158.93*1000)/(0.604*1000); end; if site='ihw' then do; v45=((v45/8.311)*199.20*1000)/(0.604*1000); end; if site='irp' then do; v45=((v45/8.311)*203.51*1000)/(0.604*1000); end; if site='pb' then do; v45=((v45/8.311)*207.47*1000)/(0.604*1000); end; To fill in gaps I use the met station that we have at the UW Ag Res Stn in Ashland. The screening commands are as follows: if agPAR<0 or agPAR>3000 then agPARbad=1; if 0 le agPAR<3000 then agPARbad=0; if v45<0 or v45>3000 then v45bad=1; if 0 le v45<3000 then v45bad=0; if (v45=. or v45bad=1) and agPARbad=1 then PAR=.; if v45 ne . and v45bad=0 then PAR=v45; if (v45=. or v45bad=1) and agPARbad=0 then PAR=agPAR; where v45 is PAR variable at each individual tower and agPAR is PAR at the Ag Stn. I have attached the Ag Stn data file. The complete list of variables in this file is: program code year day hour agts1 agts2 agts3 agta agrh agM agPAR agU agP agsnow agM0 > P.S. I am using the fluxes from the file IHW_data2003.r1w.xls file (rotation 1 with WPL) - column 8 CO2 flux - is this right? What happened to LI800_WPL storage flux - you had it for 2002 - do you have storage fluxes for 2003? Yes, you are right. 'Rotation 1" is the file to use and column 8 is CO2 flux. We have storage data for MHW and MRP and later in 2003 for IHW. The data is in a separate file. I attached the time series of storage to this e-mail. However, I have not completely screened it. The variables and units are as follows: site DOY hour co2 - co2 concentration, ppm d_800 - change in co2 concentration over 30-minute period, ppm d_800ugm - storage flux, micrograms of CO2 per m^3 of air d_800vol - storage flux, micrograms of CO2 per m^2 of ground Let me know if it makes sense. Asko