rsi0 = Clear_sky_downward_solar_flux;  %y(:,4);  %clear sky net solar
rli0 = Clear_sky_downward_long_wave_flux;  %y(:,3);  %clear sky net IR
rsi = Net_short_wave_radiation_surface;  %y(:,19);  %net solar
rli = Net_long_wave_radiation_surface;  %y(:,18);  %net IR
hsi = Sensible_heat_flux;  %y(:,22);  %sensible heat
hli = Latent_heat_flux;  %y(:,13);  %latent heat
lowcld = Low_level_cloud_cover;  %y(:,15);
midcld = Mid_level_cloud_cover;  %y(:,17);
tauni = Meridional_component_of_momentum_flux;  %y(:,16);  %meridional stress N/m^2
tauwi = Zonal_component_of_momentum_flux;  %y(:,26);  %zonal stress N/m^2
wcloud = Cloud_water*1e3;  %y(:,6)*1e3;  %cloud liq water in gm/m^2
wvap = Precipitable_water/10;  %y(:,20)/10;  %precip water in cm
P = Pressure/100;  %y(:,21)/100;  %pressure, mb
tsf = Temperature_surface-tdk;  %y(:,24)-tdk;  %surf temp, C
tdew = Dew_point_temperature-tdk;  %y(:,27)-tdk;  %dew pt temp z = 2m, C
ta = Temperature-tdk;  %y(:,29)-tdk;  %air temp z = 2m, C
ts = Soil_temperature-tdk;  %y(:,28)-tdk;  %soil temp, C
ux = u_wind;  %y(:,30);  %u_wind(units: = "m/s") z = 10m
uy = v_wind;  %y(:,31);  %v_wind(units: = "m/s") z = 10m