Changeset 693 for IOIPSL/trunk/src

Timestamp:

07/29/09 17:49:31 (15 years ago)

Author:

bellier

Message:

Added an integer function "ioget_year_len(year)" returning the
true length of a year in days according to the used calendar.

File:

• IOIPSL/trunk/src/calendar.f90 (modified) (2 diffs)

IOIPSL/trunk/src/calendar.f90

@@ -44 +44 @@
   PRIVATE
   PUBLIC :: ymds2ju,ju2ymds,tlen2itau,isittime,ioconf_calendar, &
- &          ioget_calendar,ioget_mon_len,itau2date,ioget_timestamp, &
- &          ioconf_startdate,itau2ymds,time_diff,time_add, &
- &          lock_calendar
+ &          ioget_calendar,ioget_mon_len,ioget_year_len,itau2date, &
+ &          ioget_timestamp,ioconf_startdate,itau2ymds, &
+ &          time_diff,time_add,lock_calendar
 !-
   INTERFACE ioget_calendar
@@ -909 +909 @@
 !===
 !-
+INTEGER FUNCTION ioget_year_len (year)
+!!--------------------------------------------------------------------
+!! The "ioget_year_len" function returns
+!! the number of days in "year", in the current calendar.
+!!
+!! INTEGER FUNCTION ioget_year_len (year)
+!!
+!! INPUT
+!!
+!! (I) year  : year
+!!
+!! OUTPUT
+!!
+!! (I) ioget_year_len : number of days in the year
+!!--------------------------------------------------------------------
+  IMPLICIT NONE
+!-
+  INTEGER,INTENT(IN) :: year
+!-
+  INTEGER :: yl
+!---------------------------------------------------------------------
+  SELECT CASE(TRIM(calendar_used))
+  CASE('gregorian')
+    yl = 365
+    IF (    ((MOD(year,4) == 0).AND.(MOD(year,100) /= 0)) &
+        .OR.(MOD(year,400) == 0) ) THEN
+      yl = yl+1
+    ENDIF
+  CASE('julian')
+    yl = 365
+    IF (MOD(year,4) == 0) THEN
+      yl = yl+1
+    ENDIF
+  CASE DEFAULT
+    yl = NINT(one_year)
+  END SELECT
+  ioget_year_len = yl
+!--------------------------
+END FUNCTION ioget_year_len
+!-
+!===
+!-
 SUBROUTINE ioget_timestamp (string)
 !---------------------------------------------------------------------