Weather Logistics UK

Read months.pro

From WeatherMedia

Jump to: navigation, search

Download our Climate and Energy Calculator - See How Heating Energy Demand for Household or Business is Affected By UK Temperatures

Please View Our HDD (Heating_degree_days) information page for more details on how to estimate your gas / oil bills for the winter ... Our download consists of a 5MB .xls file

Requires
Seasonal_forecast_code, thermal wind, NAOJet_Clouds_Monthly_Series.pro, read_months.pro

see: Upper Atmospheric Models (WikiPedia) & Designer Weather - How Seasonal Weather Forecasts are Made for further details

pro read_months, YEAR, dayno, month, DOM, YEARx=YEARx

create a 6 year climatology of jet stream latitides using the following date codes

days_leap = [indgen(31) + 1.0, indgen(29) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0] days = [indgen(31) + 1.0, indgen(28) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0, indgen(30) + 1.0, indgen(31) + 1.0]

months_leap = [replicate(1.0, 31), replicate(2.0, 29), replicate(3.0, 31), replicate(4.0, 30), replicate(5.0, 31), replicate(6.0, 30), replicate(7.0, 31), replicate(8.0, 31), replicate(9.0, 31), replicate(10.0, 30), replicate(11.0, 30), replicate(12.0, 31)] months = [replicate(1.0, 31), replicate(2.0, 28), replicate(3.0, 31), replicate(4.0, 30), replicate(5.0, 31), replicate(6.0, 30), replicate(7.0, 31), replicate(8.0, 31), replicate(9.0, 31), replicate(10.0, 30), replicate(11.0, 30), replicate(12.0, 31)]

if keyword_set(YEARx) eq 1.0 then begin

   if YEARx eq '2004' then begin
       YEAR=[replicate(2004, 366)]
       dayno=[indgen(366) + 1.0]
       month = [months_leap]
       DOM = [days_leap]
   endif
   if YEARx eq '2005' then begin
       YEAR=[replicate(2005, 365)]
       dayno=[indgen(365) + 1.0]
       month = [months]
       DOM = [days]
   endif
   if YEARx eq '2006' then begin
       YEAR=[replicate(2006, 365)]
       dayno=[indgen(365) + 1.0]
       month = [months]
       DOM = [days]
   endif
   if YEARx eq '2007' then begin
       YEAR=[replicate(2007, 365)]
       month = [months]
       dayno=[indgen(365) + 1.0]
       DOM = [days]
   endif
   if YEARx eq '2008' then begin
       YEAR=[replicate(2008, 366)]
       dayno=[indgen(366) + 1.0]
       month = [months_leap]
       DOM = [days_leap]
   endif
   if YEARx eq '2009' then begin
       YEAR=[replicate(2009, 365)]
       dayno=[indgen(365) + 1.0]
       month = [months]
       DOM = [days]
   endif

endif else begin

   YEAR=[replicate(2004, 366), replicate(2005, 365), replicate(2006, 365), replicate(2007, 365), replicate(2008, 366), replicate(2009, 365)]
   dayno=[indgen(366) + 1.0, indgen(365) + 1.0, indgen(365) + 1.0, indgen(365) + 1.0, indgen(366) + 1.0, indgen(365)]
   month=[months_leap, months, months, months, months_leap, months]
   DOM = [days_leap, days, days, days, days_leap, days]

endelse

end

Personal tools