!-----  Ex4_4: Calendar   -----
      INTEGER:: n, s
      PRINT*, 'Input a year:'
      READ*, n
      s = 365*n + n/4 - n/100 + n/400
      PRINT*, n, "年12月31日は，１年１月１日から数えて", s, "日目です。"
      END