source: trunk/Calendar/daysinmonth.pro @ 9

;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
;+
; NAME:daysinmonth
;
; PURPOSE:give the nomber of days in a specefic month
;
; CATEGORY: calandar
;
; CALLING SEQUENCE: ndays = daysinmonth([month[, year]])
;
; INPUTS (optionnal): month and year. Year is used only if the
;      common variable key_caltype = 'greg'. In that case, month
;      and year must have the same number of elements. if not
;      provided, we take month and year from "time" common variable.
;
; KEYWORD PARAMETERS:
;
; OUTPUTS:
;
; COMMON BLOCKS: cm_4cal
;        
; SIDE EFFECTS: return -1 in case of error
;
; RESTRICTIONS:
;
; EXAMPLE:
;
; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr)
;                       2/7/98
; update/review/english/new commons: June 2005 Sebastien Masson.
;-
;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
function daysinmonth, month, year
;------------------------------------------------------------
; include commons
@cm_4cal
IF NOT keyword_set(key_forgetold) THEN BEGIN
@updatenew
ENDIF
;------------------------------------------------------------
  IF n_elements(key_caltype) EQ 0 THEN key_caltype = 'greg'
  CASE key_caltype OF
    '360d': if n_elements(month) GT 1 THEN $
        return, replicate(30, n_elements(month)) ELSE return, 30
    'noleap':BEGIN
      days_in_mth = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31]
      IF n_elements(month) EQ 0 THEN caldat, time, month
      return, days_in_mth[month-1]
    END
    'greg':BEGIN
      CASE n_params() OF
        0:caldat, time, month, day, year
        2:IF n_elements(month) NE n_elements(year) THEN $
          return, report('month and year must have the same number of elements')
        ELSE:return, report('daysinmonth accept 0 or 2 input parameters')
      ENDCASE
      days_in_mth = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31]
      return, days_in_mth[month-1] + leapyr(year)*(month EQ 2)
    END 
    ELSE:return, report('only 3 types of calendar are accepted: greg, 360d and noleap')
  ENDCASE

END