SAS code for analyzing water temperature data

This code may be used to fit linear models with multivariate random effects and heterogeneous measurement-level residual variances. The code as written may be used to estimate associations between water temperature ('temp') and continuous year ('yearctr'), study reach (or field station; 'fs'), log-transformed mean July water discharge (in 1000 cms units; 'logmeanJulycms1000'), number of days from a central sampling date (for a given year, days from a standard month and day; 'jdatectr'), time of sampling (in fractions of hours from noon on the given date; 'timectrhr') and interactions thereof.