external/esmf_time_f90/Test1.F90

   1 !
   2 ! Sub-system tests for esmf_time_f90
   3 !
   4 ! Someday, switch over to funit!
   5 !
   6
   7 MODULE my_tests
   8   USE ESMF_Mod
   9   IMPLICIT NONE
  10
  11   ! Set this to .TRUE. to make wrf_error_fatal3() print a message on failure
  12   ! instead of stopping the program.  Use for testing only (since we cannot
  13   ! catch exceptions in Fortran90!!)
  14   LOGICAL :: WRF_ERROR_FATAL_PRINT = .FALSE.
  15
  16 CONTAINS
  17
  18   ! Test printing of an ESMF_Time or ESMF_TimeInterval object.
  19   !
  20   ! Correct results are also passed in through this interface and compared
  21   ! with computed results.  PASS/FAIL messages are printed.
  22   !
  23   SUBROUTINE test_print(  t_yy,  t_mm,  t_dd,  t_h,  t_m,  t_s, t_sn, t_sd, &
  24                          ti_yy, ti_mm, ti_dd, ti_h, ti_m, ti_s, ti_sn, ti_sd, &
  25                          res_str, testname, expect_error )
  26     INTEGER, INTENT(IN), OPTIONAL :: t_YY
  27     INTEGER, INTENT(IN), OPTIONAL :: t_MM  ! month
  28     INTEGER, INTENT(IN), OPTIONAL :: t_DD  ! day of month
  29     INTEGER, INTENT(IN), OPTIONAL :: t_H
  30     INTEGER, INTENT(IN), OPTIONAL :: t_M
  31     INTEGER, INTENT(IN), OPTIONAL :: t_S
  32     INTEGER, INTENT(IN), OPTIONAL :: t_Sn
  33     INTEGER, INTENT(IN), OPTIONAL :: t_Sd
  34     INTEGER, INTENT(IN), OPTIONAL :: ti_YY
  35     INTEGER, INTENT(IN), OPTIONAL :: ti_MM  ! month
  36     INTEGER, INTENT(IN), OPTIONAL :: ti_DD  ! day of month
  37     INTEGER, INTENT(IN), OPTIONAL :: ti_H
  38     INTEGER, INTENT(IN), OPTIONAL :: ti_M
  39     INTEGER, INTENT(IN), OPTIONAL :: ti_S
  40     INTEGER, INTENT(IN), OPTIONAL :: ti_Sn
  41     INTEGER, INTENT(IN), OPTIONAL :: ti_Sd
  42     CHARACTER (LEN=*), INTENT(IN) :: res_str
  43     CHARACTER (LEN=*), INTENT(IN), OPTIONAL :: testname
  44     LOGICAL, OPTIONAL, INTENT(IN) :: expect_error
  45     ! locals
  46     INTEGER :: it_YY
  47     INTEGER :: it_MM  ! month
  48     INTEGER :: it_DD  ! day of month
  49     INTEGER :: it_H
  50     INTEGER :: it_M
  51     INTEGER :: it_S
  52     INTEGER :: it_Sn
  53     INTEGER :: it_Sd
  54     INTEGER :: iti_YY
  55     INTEGER :: iti_MM  ! month
  56     INTEGER :: iti_DD  ! day of month
  57     INTEGER :: iti_H
  58     INTEGER :: iti_M
  59     INTEGER :: iti_S
  60     INTEGER :: iti_Sn
  61     INTEGER :: iti_Sd
  62     LOGICAL :: is_t
  63     LOGICAL :: is_ti
  64     CHARACTER (LEN=512) :: itestname
  65     LOGICAL :: iexpect_error
  66     INTEGER rc
  67     TYPE(ESMF_Time)           :: t
  68     TYPE(ESMF_TimeInterval)   :: ti
  69     CHARACTER(LEN=ESMF_MAXSTR) :: str, computed_str, frac_str
  70     CHARACTER(LEN=17) :: type_str
  71     INTEGER :: res_len, computed_len, Sn, Sd
  72     LOGICAL :: test_passed
  73
  74 !  PRINT *,'DEBUG:  BEGIN test_print()'
  75     it_YY = 0
  76     it_MM = 1
  77     it_DD = 1
  78     it_H = 0
  79     it_M = 0
  80     it_S = 0
  81     it_Sn = 0
  82     it_Sd = 0
  83     iti_YY = 0
  84     iti_MM = 0
  85     iti_DD = 0
  86     iti_H = 0
  87     iti_M = 0
  88     iti_S = 0
  89     iti_Sn = 0
  90     iti_Sd = 0
  91     itestname = ''
  92     iexpect_error = .FALSE.
  93
  94     IF ( PRESENT( t_YY ) ) it_YY = t_YY
  95     IF ( PRESENT( t_MM ) ) it_MM = t_MM
  96     IF ( PRESENT( t_DD ) ) it_DD = t_DD
  97     IF ( PRESENT( t_H ) ) it_H = t_H
  98     IF ( PRESENT( t_M ) ) it_M = t_M
  99     IF ( PRESENT( t_S ) ) it_S = t_S
 100     IF ( PRESENT( t_Sn ) ) it_Sn = t_Sn
 101     IF ( PRESENT( t_Sd ) ) it_Sd = t_Sd
 102     IF ( PRESENT( ti_YY ) ) iti_YY = ti_YY
 103     IF ( PRESENT( ti_MM ) ) iti_MM = ti_MM
 104     IF ( PRESENT( ti_DD ) ) iti_DD = ti_DD
 105     IF ( PRESENT( ti_H ) ) iti_H = ti_H
 106     IF ( PRESENT( ti_M ) ) iti_M = ti_M
 107     IF ( PRESENT( ti_S ) ) iti_S = ti_S
 108     IF ( PRESENT( ti_Sn ) ) iti_Sn = ti_Sn
 109     IF ( PRESENT( ti_Sd ) ) iti_Sd = ti_Sd
 110     IF ( PRESENT( testname ) ) itestname = TRIM(testname)
 111     IF ( PRESENT( expect_error ) ) iexpect_error = expect_error
 112
 113     ! Ensure that optional arguments are consistent...
 114     is_t = ( PRESENT( t_YY ) .OR. PRESENT( t_MM ) .OR. &
 115              PRESENT( t_DD ) .OR. PRESENT( t_H ) .OR.  &
 116              PRESENT( t_M )  .OR. PRESENT( t_S ) .OR.  &
 117              PRESENT( t_Sn )  .OR. PRESENT( t_Sd ) )
 118     is_ti = ( PRESENT( ti_YY ) .OR. PRESENT( ti_MM ) .OR. &
 119               PRESENT( ti_DD ) .OR. PRESENT( ti_H ) .OR.  &
 120               PRESENT( ti_M )  .OR. PRESENT( ti_S ) .OR.  &
 121               PRESENT( ti_Sn )  .OR. PRESENT( ti_Sd ) )
 122     IF ( is_t .EQV. is_ti ) THEN
 123       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
 124         'ERROR test_print:  inconsistent args' )
 125     ENDIF
 126
 127 !PRINT *,'DEBUG:  test_print():  init objects'
 128     ! Initialize object to be tested
 129     ! modify behavior of wrf_error_fatal3 for tests expected to fail
 130     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .TRUE.
 131     Sn = 0
 132     Sd = 0
 133     IF ( is_t ) THEN
 134       type_str = 'ESMF_Time'
 135 !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeSet()'
 136 !PRINT *,'DEBUG:  test_print():  YY,MM,DD,H,M,S,Sn,Sd = ', it_YY,it_MM,it_DD,it_H,it_M,it_S,it_Sn,it_Sd
 137       CALL ESMF_TimeSet( t, YY=it_YY, MM=it_MM, DD=it_DD , &
 138                              H=it_H, M=it_M, S=it_S, Sn=it_Sn, Sd=it_Sd, rc=rc )
 139 !PRINT *,'DEBUG:  test_print():  back from ESMF_TimeSet()'
 140       CALL test_check_error( ESMF_SUCCESS, rc, &
 141                              TRIM(itestname)//'ESMF_TimeSet() ', &
 142                              __FILE__ , &
 143                              __LINE__  )
 144 !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeGet()'
 145       CALL ESMF_TimeGet( t, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
 146       CALL test_check_error( ESMF_SUCCESS, rc, &
 147                             TRIM(itestname)//'ESMF_TimeGet() ', &
 148                             __FILE__ , &
 149                             __LINE__  )
 150 !PRINT *,'DEBUG:  test_print():  back from ESMF_TimeGet(), computed_str = ',TRIM(computed_str)
 151     ELSE
 152       type_str = 'ESMF_TimeInterval'
 153 !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeIntervalSet()'
 154       CALL ESMF_TimeIntervalSet( ti, YY=iti_YY, MM=iti_MM, &
 155                                       D=iti_DD ,           &
 156                                       H=iti_H, M=iti_M,    &
 157                                       S=iti_S, Sn=iti_Sn, Sd=iti_Sd, rc=rc )
 158       CALL test_check_error( ESMF_SUCCESS, rc, &
 159                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
 160                              __FILE__ , &
 161                              __LINE__  )
 162 !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeIntervalGet()'
 163       CALL ESMF_TimeIntervalGet( ti, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
 164       CALL test_check_error( ESMF_SUCCESS, rc, &
 165                             TRIM(itestname)//'ESMF_TimeGet() ', &
 166                             __FILE__ , &
 167                             __LINE__  )
 168     ENDIF
 169     ! handle fractions
 170     IF ( Sd > 0 ) THEN
 171       IF ( Sn > 0 ) THEN
 172         WRITE(frac_str,FMT="('+',I2.2,'/',I2.2)") abs(Sn), Sd
 173       ELSE IF ( Sn < 0 ) THEN
 174         WRITE(frac_str,FMT="('-',I2.2,'/',I2.2)") abs(Sn), Sd
 175       ELSE
 176         frac_str = ''
 177       ENDIF
 178       computed_str = TRIM(computed_str)//TRIM(frac_str)
 179     ENDIF
 180     ! restore default behavior of wrf_error_fatal3
 181     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .FALSE.
 182 !PRINT *,'DEBUG:  test_print():  done init objects'
 183
 184 !PRINT *,'DEBUG:  test_print():  check result'
 185     ! check result
 186     test_passed = .FALSE.
 187     res_len = LEN_TRIM(res_str)
 188     computed_len = LEN_TRIM(computed_str)
 189     IF ( res_len == computed_len ) THEN
 190       IF ( computed_str(1:computed_len) == res_str(1:res_len) ) THEN
 191         test_passed = .TRUE.
 192       ENDIF
 193     ENDIF
 194     IF ( test_passed ) THEN
 195       WRITE(*,FMT='(A)') 'PASS:  '//TRIM(itestname)
 196     ELSE
 197       WRITE(*,'(9A)') 'FAIL:  ',TRIM(itestname),':  printing ',TRIM(type_str), &
 198         '  expected <', TRIM(res_str),'>  but computed <',TRIM(computed_str),'>'
 199     ENDIF
 200 !PRINT *,'DEBUG:  END test_print()'
 201
 202   END SUBROUTINE test_print
 203
 204
 205
 206   ! Test the following arithmetic operations on ESMF_Time and
 207   ! ESMF_TimeInterval objects:
 208   !  ESMF_Time         = ESMF_Time         + ESMF_TimeInterval
 209   !  ESMF_Time         = ESMF_TimeInterval + ESMF_Time
 210   !  ESMF_Time         = ESMF_Time         - ESMF_TimeInterval
 211   !  ESMF_TimeInterval = ESMF_Time         - ESMF_Time
 212   !  ESMF_TimeInterval = ESMF_TimeInterval + ESMF_TimeInterval
 213   !  ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
 214   !  ESMF_TimeInterval = ESMF_TimeInterval * INTEGER
 215   !  ESMF_TimeInterval = ESMF_TimeInterval / INTEGER
 216   !
 217   ! Correct results are also passed in through this interface and compared
 218   ! with computed results.  PASS/FAIL messages are printed.
 219   !
 220   ! Operations are expressed as res = op1 +|- op2
 221   !
 222   SUBROUTINE test_arithmetic( add_op, multiply_op,                                       &
 223      op1_t_yy,  op1_t_mm,  op1_t_dd,  op1_t_h,  op1_t_m,  op1_t_s,  op1_t_sn,  op1_t_sd, &
 224     op1_ti_yy, op1_ti_mm, op1_ti_dd, op1_ti_h, op1_ti_m, op1_ti_s, op1_ti_sn, op1_ti_sd, &
 225      op2_t_yy,  op2_t_mm,  op2_t_dd,  op2_t_h,  op2_t_m,  op2_t_s,  op2_t_sn,  op2_t_sd, &
 226     op2_ti_yy, op2_ti_mm, op2_ti_dd, op2_ti_h, op2_ti_m, op2_ti_s, op2_ti_sn, op2_ti_sd, &
 227     op2_int,                                                                             &
 228      res_t_yy,  res_t_mm,  res_t_dd,  res_t_h,  res_t_m,  res_t_s,  res_t_sn,  res_t_sd, &
 229     res_ti_yy, res_ti_mm, res_ti_dd, res_ti_h, res_ti_m, res_ti_s, res_ti_sn, res_ti_sd, &
 230     res_int, testname, expect_error )
 231     LOGICAL, INTENT(IN), OPTIONAL :: add_op      ! .TRUE.=add, .FALSE.=subtract
 232     LOGICAL, INTENT(IN), OPTIONAL :: multiply_op ! .TRUE.=multiply, .FALSE.=divide
 233     INTEGER, INTENT(IN), OPTIONAL :: op1_t_YY
 234     INTEGER, INTENT(IN), OPTIONAL :: op1_t_MM  ! month
 235     INTEGER, INTENT(IN), OPTIONAL :: op1_t_DD  ! day of month
 236     INTEGER, INTENT(IN), OPTIONAL :: op1_t_H
 237     INTEGER, INTENT(IN), OPTIONAL :: op1_t_M
 238     INTEGER, INTENT(IN), OPTIONAL :: op1_t_S
 239     INTEGER, INTENT(IN), OPTIONAL :: op1_t_Sn
 240     INTEGER, INTENT(IN), OPTIONAL :: op1_t_Sd
 241     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_YY
 242     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_MM  ! month
 243     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_DD  ! day of month
 244     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_H
 245     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_M
 246     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_S
 247     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_Sn
 248     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_Sd
 249     INTEGER, INTENT(IN), OPTIONAL :: op2_t_YY
 250     INTEGER, INTENT(IN), OPTIONAL :: op2_t_MM  ! month
 251     INTEGER, INTENT(IN), OPTIONAL :: op2_t_DD  ! day of month
 252     INTEGER, INTENT(IN), OPTIONAL :: op2_t_H
 253     INTEGER, INTENT(IN), OPTIONAL :: op2_t_M
 254     INTEGER, INTENT(IN), OPTIONAL :: op2_t_S
 255     INTEGER, INTENT(IN), OPTIONAL :: op2_t_Sn
 256     INTEGER, INTENT(IN), OPTIONAL :: op2_t_Sd
 257     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_YY
 258     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_MM  ! month
 259     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_DD  ! day of month
 260     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_H
 261     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_M
 262     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_S
 263     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_Sn
 264     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_Sd
 265     INTEGER, INTENT(IN), OPTIONAL :: op2_int
 266     INTEGER, INTENT(IN), OPTIONAL :: res_t_YY
 267     INTEGER, INTENT(IN), OPTIONAL :: res_t_MM  ! month
 268     INTEGER, INTENT(IN), OPTIONAL :: res_t_DD  ! day of month
 269     INTEGER, INTENT(IN), OPTIONAL :: res_t_H
 270     INTEGER, INTENT(IN), OPTIONAL :: res_t_M
 271     INTEGER, INTENT(IN), OPTIONAL :: res_t_S
 272     INTEGER, INTENT(IN), OPTIONAL :: res_t_Sn
 273     INTEGER, INTENT(IN), OPTIONAL :: res_t_Sd
 274     INTEGER, INTENT(IN), OPTIONAL :: res_ti_YY
 275     INTEGER, INTENT(IN), OPTIONAL :: res_ti_MM  ! month
 276     INTEGER, INTENT(IN), OPTIONAL :: res_ti_DD  ! day of month
 277     INTEGER, INTENT(IN), OPTIONAL :: res_ti_H
 278     INTEGER, INTENT(IN), OPTIONAL :: res_ti_M
 279     INTEGER, INTENT(IN), OPTIONAL :: res_ti_S
 280     INTEGER, INTENT(IN), OPTIONAL :: res_ti_Sn
 281     INTEGER, INTENT(IN), OPTIONAL :: res_ti_Sd
 282     INTEGER, INTENT(IN), OPTIONAL :: res_int
 283     CHARACTER (LEN=*), OPTIONAL, INTENT(IN) :: testname
 284     LOGICAL, OPTIONAL, INTENT(IN) :: expect_error
 285     ! locals
 286     LOGICAL :: iadd_op
 287     LOGICAL :: isubtract_op
 288     LOGICAL :: imultiply_op
 289     LOGICAL :: idivide_op
 290     INTEGER :: iop1_t_YY
 291     INTEGER :: iop1_t_MM  ! month
 292     INTEGER :: iop1_t_DD  ! day of month
 293     INTEGER :: iop1_t_H
 294     INTEGER :: iop1_t_M
 295     INTEGER :: iop1_t_S
 296     INTEGER :: iop1_t_Sn
 297     INTEGER :: iop1_t_Sd
 298     INTEGER :: iop1_ti_YY
 299     INTEGER :: iop1_ti_MM  ! month
 300     INTEGER :: iop1_ti_DD  ! day of month
 301     INTEGER :: iop1_ti_H
 302     INTEGER :: iop1_ti_M
 303     INTEGER :: iop1_ti_S
 304     INTEGER :: iop1_ti_Sn
 305     INTEGER :: iop1_ti_Sd
 306     INTEGER :: iop2_t_YY
 307     INTEGER :: iop2_t_MM  ! month
 308     INTEGER :: iop2_t_DD  ! day of month
 309     INTEGER :: iop2_t_H
 310     INTEGER :: iop2_t_M
 311     INTEGER :: iop2_t_S
 312     INTEGER :: iop2_t_Sn
 313     INTEGER :: iop2_t_Sd
 314     INTEGER :: iop2_ti_YY
 315     INTEGER :: iop2_ti_MM  ! month
 316     INTEGER :: iop2_ti_DD  ! day of month
 317     INTEGER :: iop2_ti_H
 318     INTEGER :: iop2_ti_M
 319     INTEGER :: iop2_ti_S
 320     INTEGER :: iop2_ti_Sn
 321     INTEGER :: iop2_ti_Sd
 322     INTEGER :: ires_t_YY
 323     INTEGER :: ires_t_MM  ! month
 324     INTEGER :: ires_t_DD  ! day of month
 325     INTEGER :: ires_t_H
 326     INTEGER :: ires_t_M
 327     INTEGER :: ires_t_S
 328     INTEGER :: ires_t_Sn
 329     INTEGER :: ires_t_Sd
 330     INTEGER :: ires_ti_YY
 331     INTEGER :: ires_ti_MM  ! month
 332     INTEGER :: ires_ti_DD  ! day of month
 333     INTEGER :: ires_ti_H
 334     INTEGER :: ires_ti_M
 335     INTEGER :: ires_ti_S
 336     INTEGER :: ires_ti_Sn
 337     INTEGER :: ires_ti_Sd
 338     LOGICAL :: op1_is_t , op2_is_t , res_is_t
 339     LOGICAL :: op1_is_ti, op2_is_ti, res_is_ti, op2_is_int
 340     LOGICAL :: res_is_int
 341     INTEGER :: num_ops, num_op1, num_op2, num_res
 342     LOGICAL :: unsupported_op, test_passed
 343     CHARACTER (LEN=512) :: itestname
 344     LOGICAL :: iexpect_error
 345     INTEGER :: rc
 346     INTEGER :: computed_int, Sn, Sd
 347     TYPE(ESMF_Time)           :: op1_t , op2_t , res_t, computed_t
 348     TYPE(ESMF_TimeInterval)   :: op1_ti, op2_ti, res_ti, computed_ti
 349     CHARACTER(LEN=ESMF_MAXSTR) :: str, op1_str, op2_str, res_str, computed_str, frac_str
 350     CHARACTER(LEN=1) :: op_str
 351     CHARACTER(LEN=17) :: op1_type_str, op2_type_str, res_type_str
 352
 353     iadd_op = .FALSE.
 354     isubtract_op = .FALSE.
 355     imultiply_op = .FALSE.
 356     idivide_op = .FALSE.
 357     iop1_t_YY = 0
 358     iop1_t_MM = 1
 359     iop1_t_DD = 1
 360     iop1_t_H = 0
 361     iop1_t_M = 0
 362     iop1_t_S = 0
 363     iop1_t_Sn = 0
 364     iop1_t_Sd = 0
 365     iop1_ti_YY = 0
 366     iop1_ti_MM = 0
 367     iop1_ti_DD = 0
 368     iop1_ti_H = 0
 369     iop1_ti_M = 0
 370     iop1_ti_S = 0
 371     iop1_ti_Sn = 0
 372     iop1_ti_Sd = 0
 373     iop2_t_YY = 0
 374     iop2_t_MM = 1
 375     iop2_t_DD = 1
 376     iop2_t_H = 0
 377     iop2_t_M = 0
 378     iop2_t_S = 0
 379     iop2_t_Sn = 0
 380     iop2_t_Sd = 0
 381     iop2_ti_YY = 0
 382     iop2_ti_MM = 0
 383     iop2_ti_DD = 0
 384     iop2_ti_H = 0
 385     iop2_ti_M = 0
 386     iop2_ti_S = 0
 387     iop2_ti_Sn = 0
 388     iop2_ti_Sd = 0
 389     ires_t_YY = 0
 390     ires_t_MM = 1
 391     ires_t_DD = 1
 392     ires_t_H = 0
 393     ires_t_M = 0
 394     ires_t_S = 0
 395     ires_t_Sn = 0
 396     ires_t_Sd = 0
 397     ires_ti_YY = 0
 398     ires_ti_MM = 0
 399     ires_ti_DD = 0
 400     ires_ti_H = 0
 401     ires_ti_M = 0
 402     ires_ti_S = 0
 403     ires_ti_Sn = 0
 404     ires_ti_Sd = 0
 405     itestname = ''
 406     iexpect_error = .FALSE.
 407
 408     IF ( PRESENT( add_op ) ) THEN
 409       iadd_op = add_op
 410       isubtract_op = ( .NOT. add_op )
 411     ENDIF
 412     IF ( PRESENT( multiply_op ) ) THEN
 413       imultiply_op = multiply_op
 414       idivide_op = ( .NOT. multiply_op )
 415     ENDIF
 416     num_ops = 0
 417     IF ( iadd_op )      num_ops = num_ops + 1
 418     IF ( isubtract_op ) num_ops = num_ops + 1
 419     IF ( imultiply_op ) num_ops = num_ops + 1
 420     IF ( idivide_op )   num_ops = num_ops + 1
 421     IF ( num_ops /= 1 ) THEN
 422       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
 423         'ERROR test_arithmetic:  inconsistent operation' )
 424     ENDIF
 425     IF ( PRESENT( op1_t_YY ) ) iop1_t_YY = op1_t_YY
 426     IF ( PRESENT( op1_t_MM ) ) iop1_t_MM = op1_t_MM
 427     IF ( PRESENT( op1_t_DD ) ) iop1_t_DD = op1_t_DD
 428     IF ( PRESENT( op1_t_H ) ) iop1_t_H = op1_t_H
 429     IF ( PRESENT( op1_t_M ) ) iop1_t_M = op1_t_M
 430     IF ( PRESENT( op1_t_S ) ) iop1_t_S = op1_t_S
 431     IF ( PRESENT( op1_t_Sn ) ) iop1_t_Sn = op1_t_Sn
 432     IF ( PRESENT( op1_t_Sd ) ) iop1_t_Sd = op1_t_Sd
 433     IF ( PRESENT( op1_ti_YY ) ) iop1_ti_YY = op1_ti_YY
 434     IF ( PRESENT( op1_ti_MM ) ) iop1_ti_MM = op1_ti_MM
 435     IF ( PRESENT( op1_ti_DD ) ) iop1_ti_DD = op1_ti_DD
 436     IF ( PRESENT( op1_ti_H ) ) iop1_ti_H = op1_ti_H
 437     IF ( PRESENT( op1_ti_M ) ) iop1_ti_M = op1_ti_M
 438     IF ( PRESENT( op1_ti_S ) ) iop1_ti_S = op1_ti_S
 439     IF ( PRESENT( op1_ti_Sn ) ) iop1_ti_Sn = op1_ti_Sn
 440     IF ( PRESENT( op1_ti_Sd ) ) iop1_ti_Sd = op1_ti_Sd
 441     IF ( PRESENT( op2_t_YY ) ) iop2_t_YY = op2_t_YY
 442     IF ( PRESENT( op2_t_MM ) ) iop2_t_MM = op2_t_MM
 443     IF ( PRESENT( op2_t_DD ) ) iop2_t_DD = op2_t_DD
 444     IF ( PRESENT( op2_t_H ) ) iop2_t_H = op2_t_H
 445     IF ( PRESENT( op2_t_M ) ) iop2_t_M = op2_t_M
 446     IF ( PRESENT( op2_t_S ) ) iop2_t_S = op2_t_S
 447     IF ( PRESENT( op2_t_Sn ) ) iop2_t_Sn = op2_t_Sn
 448     IF ( PRESENT( op2_t_Sd ) ) iop2_t_Sd = op2_t_Sd
 449     IF ( PRESENT( op2_ti_YY ) ) iop2_ti_YY = op2_ti_YY
 450     IF ( PRESENT( op2_ti_MM ) ) iop2_ti_MM = op2_ti_MM
 451     IF ( PRESENT( op2_ti_DD ) ) iop2_ti_DD = op2_ti_DD
 452     IF ( PRESENT( op2_ti_H ) ) iop2_ti_H = op2_ti_H
 453     IF ( PRESENT( op2_ti_M ) ) iop2_ti_M = op2_ti_M
 454     IF ( PRESENT( op2_ti_S ) ) iop2_ti_S = op2_ti_S
 455     IF ( PRESENT( op2_ti_Sn ) ) iop2_ti_Sn = op2_ti_Sn
 456     IF ( PRESENT( op2_ti_Sd ) ) iop2_ti_Sd = op2_ti_Sd
 457     IF ( PRESENT( res_t_YY ) ) ires_t_YY = res_t_YY
 458     IF ( PRESENT( res_t_MM ) ) ires_t_MM = res_t_MM
 459     IF ( PRESENT( res_t_DD ) ) ires_t_DD = res_t_DD
 460     IF ( PRESENT( res_t_H ) ) ires_t_H = res_t_H
 461     IF ( PRESENT( res_t_M ) ) ires_t_M = res_t_M
 462     IF ( PRESENT( res_t_S ) ) ires_t_S = res_t_S
 463     IF ( PRESENT( res_t_Sn ) ) ires_t_Sn = res_t_Sn
 464     IF ( PRESENT( res_t_Sd ) ) ires_t_Sd = res_t_Sd
 465     IF ( PRESENT( res_ti_YY ) ) ires_ti_YY = res_ti_YY
 466     IF ( PRESENT( res_ti_MM ) ) ires_ti_MM = res_ti_MM
 467     IF ( PRESENT( res_ti_DD ) ) ires_ti_DD = res_ti_DD
 468     IF ( PRESENT( res_ti_H ) ) ires_ti_H = res_ti_H
 469     IF ( PRESENT( res_ti_M ) ) ires_ti_M = res_ti_M
 470     IF ( PRESENT( res_ti_S ) ) ires_ti_S = res_ti_S
 471     IF ( PRESENT( res_ti_Sn ) ) ires_ti_Sn = res_ti_Sn
 472     IF ( PRESENT( res_ti_Sd ) ) ires_ti_Sd = res_ti_Sd
 473     IF ( PRESENT( testname ) ) itestname = TRIM(testname)
 474     IF ( PRESENT( expect_error ) ) iexpect_error = expect_error
 475
 476     ! Ensure that optional arguments are consistent...
 477     op1_is_t = ( PRESENT( op1_t_YY ) .OR. PRESENT( op1_t_MM ) .OR. &
 478                  PRESENT( op1_t_DD ) .OR. PRESENT( op1_t_H ) .OR.  &
 479                  PRESENT( op1_t_M )  .OR. PRESENT( op1_t_S ) .OR.  &
 480                  PRESENT( op1_t_Sn )  .OR. PRESENT( op1_t_Sd ) )
 481     op1_is_ti = ( PRESENT( op1_ti_YY ) .OR. PRESENT( op1_ti_MM ) .OR. &
 482                   PRESENT( op1_ti_DD ) .OR. PRESENT( op1_ti_H ) .OR.  &
 483                   PRESENT( op1_ti_M )  .OR. PRESENT( op1_ti_S ) .OR.  &
 484                   PRESENT( op1_ti_Sn )  .OR. PRESENT( op1_ti_Sd ) )
 485     op2_is_t = ( PRESENT( op2_t_YY ) .OR. PRESENT( op2_t_MM ) .OR. &
 486                  PRESENT( op2_t_DD ) .OR. PRESENT( op2_t_H ) .OR.  &
 487                  PRESENT( op2_t_M )  .OR. PRESENT( op2_t_S ) .OR.  &
 488                  PRESENT( op2_t_Sn )  .OR. PRESENT( op2_t_Sd ) )
 489     op2_is_ti = ( PRESENT( op2_ti_YY ) .OR. PRESENT( op2_ti_MM ) .OR. &
 490                   PRESENT( op2_ti_DD ) .OR. PRESENT( op2_ti_H ) .OR.  &
 491                   PRESENT( op2_ti_M )  .OR. PRESENT( op2_ti_S ) .OR.  &
 492                   PRESENT( op2_ti_Sn )  .OR. PRESENT( op2_ti_Sd ) )
 493     op2_is_int = ( PRESENT( op2_int ) )
 494     res_is_t = ( PRESENT( res_t_YY ) .OR. PRESENT( res_t_MM ) .OR. &
 495                  PRESENT( res_t_DD ) .OR. PRESENT( res_t_H ) .OR.  &
 496                  PRESENT( res_t_M )  .OR. PRESENT( res_t_S ) .OR.  &
 497                  PRESENT( res_t_Sn )  .OR. PRESENT( res_t_Sd ) )
 498     res_is_ti = ( PRESENT( res_ti_YY ) .OR. PRESENT( res_ti_MM ) .OR. &
 499                   PRESENT( res_ti_DD ) .OR. PRESENT( res_ti_H ) .OR.  &
 500                   PRESENT( res_ti_M )  .OR. PRESENT( res_ti_S ) .OR.  &
 501                   PRESENT( res_ti_Sn )  .OR. PRESENT( res_ti_Sd ) )
 502     res_is_int = ( PRESENT( res_int ) )
 503     num_op1 = 0
 504     IF ( op1_is_t   ) num_op1 = num_op1 + 1
 505     IF ( op1_is_ti  ) num_op1 = num_op1 + 1
 506     IF ( num_op1 /= 1 ) THEN
 507       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
 508         'ERROR test_arithmetic:  inconsistent args for op1' )
 509     ENDIF
 510     num_op2 = 0
 511     IF ( op2_is_t   ) num_op2 = num_op2 + 1
 512     IF ( op2_is_ti  ) num_op2 = num_op2 + 1
 513     IF ( op2_is_int ) num_op2 = num_op2 + 1
 514     IF ( num_op2 /= 1 ) THEN
 515       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
 516         'ERROR test_arithmetic:  inconsistent args for op2' )
 517     ENDIF
 518     num_res = 0
 519     IF ( res_is_t   ) num_res = num_res + 1
 520     IF ( res_is_ti  ) num_res = num_res + 1
 521     IF ( res_is_int ) num_res = num_res + 1
 522     IF ( num_res /= 1 ) THEN
 523       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
 524         'ERROR test_arithmetic:  inconsistent args for result' )
 525     ENDIF
 526
 527     ! Initialize op1
 528     IF ( op1_is_t ) THEN
 529       op1_type_str = 'ESMF_Time'
 530       CALL ESMF_TimeSet( op1_t, YY=iop1_t_YY, MM=iop1_t_MM, DD=iop1_t_DD , &
 531                                  H=iop1_t_H, M=iop1_t_M, S=iop1_t_S, Sn=iop1_t_Sn, Sd=iop1_t_Sd, rc=rc )
 532       CALL test_check_error( ESMF_SUCCESS, rc, &
 533                              TRIM(itestname)//'ESMF_TimeSet() ', &
 534                              __FILE__ , &
 535                              __LINE__  )
 536       CALL ESMF_TimeGet( op1_t, timeString=op1_str, Sn=Sn, Sd=Sd, rc=rc )
 537       CALL test_check_error( ESMF_SUCCESS, rc, &
 538                             TRIM(itestname)//'ESMF_TimeGet() ', &
 539                             __FILE__ , &
 540                             __LINE__  )
 541       ! handle fractions
 542       CALL fraction_to_string( Sn, Sd, frac_str )
 543       op1_str = TRIM(op1_str)//TRIM(frac_str)
 544     ELSE
 545       op1_type_str = 'ESMF_TimeInterval'
 546       CALL ESMF_TimeIntervalSet( op1_ti, YY=iop1_ti_YY, MM=iop1_ti_MM, &
 547                                           D=iop1_ti_DD ,               &
 548                                           H=iop1_ti_H, M=iop1_ti_M,    &
 549                                           S=iop1_ti_S, Sn=iop1_ti_Sn, Sd=iop1_ti_Sd, rc=rc )
 550       CALL test_check_error( ESMF_SUCCESS, rc, &
 551                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
 552                              __FILE__ , &
 553                              __LINE__  )
 554       CALL ESMF_TimeIntervalGet( op1_ti, timeString=op1_str, Sn=Sn, Sd=Sd, rc=rc )
 555       CALL test_check_error( ESMF_SUCCESS, rc, &
 556                             TRIM(itestname)//'ESMF_TimeGet() ', &
 557                             __FILE__ , &
 558                             __LINE__  )
 559       ! handle fractions
 560       CALL fraction_to_string( Sn, Sd, frac_str )
 561       op1_str = TRIM(op1_str)//TRIM(frac_str)
 562     ENDIF
 563     ! Initialize op2
 564     IF ( op2_is_t ) THEN
 565       op2_type_str = 'ESMF_Time'
 566       CALL ESMF_TimeSet( op2_t, YY=iop2_t_YY, MM=iop2_t_MM, DD=iop2_t_DD , &
 567                                  H=iop2_t_H, M=iop2_t_M, S=iop2_t_S, Sn=iop2_t_Sn, Sd=iop2_t_Sd, rc=rc )
 568       CALL test_check_error( ESMF_SUCCESS, rc, &
 569                              TRIM(itestname)//'ESMF_TimeSet() ', &
 570                              __FILE__ , &
 571                              __LINE__  )
 572       CALL ESMF_TimeGet( op2_t, timeString=op2_str, Sn=Sn, Sd=Sd, rc=rc )
 573       CALL test_check_error( ESMF_SUCCESS, rc, &
 574                             TRIM(itestname)//'ESMF_TimeGet() ', &
 575                             __FILE__ , &
 576                             __LINE__  )
 577       ! handle fractions
 578       CALL fraction_to_string( Sn, Sd, frac_str )
 579       op2_str = TRIM(op2_str)//TRIM(frac_str)
 580     ELSE IF ( op2_is_ti ) THEN
 581       op2_type_str = 'ESMF_TimeInterval'
 582       CALL ESMF_TimeIntervalSet( op2_ti, YY=iop2_ti_YY, MM=iop2_ti_MM, &
 583                                           D=iop2_ti_DD ,               &
 584                                           H=iop2_ti_H, M=iop2_ti_M,    &
 585                                           S=iop2_ti_S, Sn=iop2_ti_Sn, Sd=iop2_ti_Sd, rc=rc )
 586       CALL test_check_error( ESMF_SUCCESS, rc, &
 587                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
 588                              __FILE__ , &
 589                              __LINE__  )
 590       CALL ESMF_TimeIntervalGet( op2_ti, timeString=op2_str, Sn=Sn, Sd=Sd, rc=rc )
 591       CALL test_check_error( ESMF_SUCCESS, rc, &
 592                             TRIM(itestname)//'ESMF_TimeGet() ', &
 593                             __FILE__ , &
 594                             __LINE__  )
 595       ! handle fractions
 596       CALL fraction_to_string( Sn, Sd, frac_str )
 597       op2_str = TRIM(op2_str)//TRIM(frac_str)
 598     ELSE
 599       op2_type_str = 'INTEGER'
 600       IF ( op2_int > 0 ) THEN
 601         WRITE(op2_str,FMT="('+',I8.8)") ABS(op2_int)
 602       ELSE
 603         WRITE(op2_str,FMT="('-',I8.8)") ABS(op2_int)
 604       ENDIF
 605     ENDIF
 606     ! Initialize res
 607     IF ( res_is_t ) THEN  ! result is ESMF_Time
 608       res_type_str = 'ESMF_Time'
 609       CALL ESMF_TimeSet( res_t, YY=ires_t_YY, MM=ires_t_MM, DD=ires_t_DD , &
 610                                  H=ires_t_H, M=ires_t_M, S=ires_t_S, Sn=ires_t_Sn, Sd=ires_t_Sd, rc=rc )
 611       CALL test_check_error( ESMF_SUCCESS, rc, &
 612                              TRIM(itestname)//'ESMF_TimeSet() ', &
 613                              __FILE__ , &
 614                              __LINE__  )
 615       CALL ESMF_TimeGet( res_t, timeString=res_str, Sn=Sn, Sd=Sd, rc=rc )
 616       CALL test_check_error( ESMF_SUCCESS, rc, &
 617                             TRIM(itestname)//'ESMF_TimeGet() ', &
 618                             __FILE__ , &
 619                             __LINE__  )
 620       ! handle fractions
 621       CALL fraction_to_string( Sn, Sd, frac_str )
 622       res_str = TRIM(res_str)//TRIM(frac_str)
 623     ELSE IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
 624       res_type_str = 'ESMF_TimeInterval'
 625       CALL ESMF_TimeIntervalSet( res_ti, YY=ires_ti_YY, MM=ires_ti_MM, &
 626                                           D=ires_ti_DD ,               &
 627                                           H=ires_ti_H, M=ires_ti_M,    &
 628                                           S=ires_ti_S, Sn=ires_ti_Sn, Sd=ires_ti_Sd, rc=rc )
 629       CALL test_check_error( ESMF_SUCCESS, rc, &
 630                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
 631                              __FILE__ , &
 632                              __LINE__  )
 633       CALL ESMF_TimeIntervalGet( res_ti, timeString=res_str, Sn=Sn, Sd=Sd, rc=rc )
 634       CALL test_check_error( ESMF_SUCCESS, rc, &
 635                             TRIM(itestname)//'ESMF_TimeGet() ', &
 636                             __FILE__ , &
 637                             __LINE__  )
 638       ! handle fractions
 639       CALL fraction_to_string( Sn, Sd, frac_str )
 640       res_str = TRIM(res_str)//TRIM(frac_str)
 641     ELSE  ! result is INTEGER
 642       res_type_str = 'INTEGER'
 643       IF ( res_int > 0 ) THEN
 644         WRITE(res_str,FMT="('+',I8.8)") ABS(res_int)
 645       ELSE
 646         WRITE(res_str,FMT="('-',I8.8)") ABS(res_int)
 647       ENDIF
 648     ENDIF
 649
 650     ! perform requested operation
 651     unsupported_op = .FALSE.
 652     ! modify behavior of wrf_error_fatal3 for operator being tested
 653     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .TRUE.
 654     ! add
 655     IF ( iadd_op ) THEN
 656       op_str = '+'
 657       IF ( res_is_t ) THEN  ! result is ESMF_Time
 658         IF ( op1_is_t .AND. op2_is_ti ) THEN
 659           !  ESMF_Time         = ESMF_Time         + ESMF_TimeInterval
 660           computed_t = op1_t + op2_ti
 661         ELSE IF ( op1_is_ti .AND. op2_is_t ) THEN
 662           !  ESMF_Time         = ESMF_TimeInterval + ESMF_Time
 663           computed_t = op1_ti + op2_t
 664         ELSE
 665           unsupported_op = .TRUE.
 666         ENDIF
 667       ELSE  ! result is ESMF_TimeInterval
 668         IF ( op1_is_ti .AND. op2_is_ti ) THEN
 669           !  ESMF_TimeInterval = ESMF_TimeInterval + ESMF_TimeInterval
 670           computed_ti = op1_ti + op2_ti
 671         ELSE
 672           unsupported_op = .TRUE.
 673         ENDIF
 674       ENDIF
 675     ! subtract
 676     ELSE  IF ( isubtract_op ) THEN
 677       op_str = '-'
 678       IF ( res_is_t ) THEN  ! result is ESMF_Time
 679         IF ( op1_is_t .AND. op2_is_ti ) THEN
 680           !  ESMF_Time         = ESMF_Time         - ESMF_TimeInterval
 681           computed_t = op1_t - op2_ti
 682         ELSE
 683           unsupported_op = .TRUE.
 684         ENDIF
 685       ELSE  ! result is ESMF_TimeInterval
 686         IF ( op1_is_t .AND. op2_is_t ) THEN
 687           !  ESMF_TimeInterval = ESMF_Time         - ESMF_Time
 688           computed_ti = op1_t - op2_t
 689         ELSE IF ( op1_is_ti .AND. op2_is_ti ) THEN
 690           !  ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
 691           computed_ti = op1_ti - op2_ti
 692         ELSE
 693           unsupported_op = .TRUE.
 694         ENDIF
 695       ENDIF
 696     ELSE  IF ( imultiply_op ) THEN
 697       op_str = '*'
 698       IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
 699         IF ( op1_is_ti .AND. op2_is_int ) THEN
 700           !  ESMF_TimeInterval = ESMF_TimeInterval * INTEGER
 701           computed_ti = op1_ti * op2_int
 702         ELSE
 703           unsupported_op = .TRUE.
 704         ENDIF
 705       ENDIF
 706     ELSE  IF ( idivide_op ) THEN
 707       op_str = '/'
 708       IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
 709         IF ( op1_is_ti .AND. op2_is_int ) THEN
 710           !  ESMF_TimeInterval = ESMF_TimeInterval / INTEGER
 711           computed_ti = op1_ti / op2_int
 712         ELSE
 713           unsupported_op = .TRUE.
 714         ENDIF
 715       ELSE IF ( res_is_int ) THEN  ! result is INTEGER
 716         IF ( op1_is_ti .AND. op2_is_ti ) THEN
 717           !  INTEGER = ESMF_TimeInterval / ESMF_TimeInterval
 718           ! number of whole time intervals
 719           computed_int = ESMF_TimeIntervalDIVQuot( op1_ti , op2_ti )
 720         ELSE
 721           unsupported_op = .TRUE.
 722         ENDIF
 723       ENDIF
 724     ENDIF
 725     ! restore default behavior of wrf_error_fatal3
 726     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .FALSE.
 727     IF ( unsupported_op ) THEN
 728       WRITE(str,*) 'ERROR test_arithmetic ',TRIM(itestname), &
 729         ':  unsupported operation (',                           &
 730         TRIM(res_type_str),' = ',TRIM(op1_type_str),' ',TRIM(op_str),' ', &
 731         TRIM(op2_type_str),')'
 732       CALL wrf_error_fatal3( __FILE__ , __LINE__ , str )
 733     ENDIF
 734
 735     ! check result
 736     test_passed = .FALSE.
 737     IF ( res_is_t ) THEN  ! result is ESMF_Time
 738       IF ( computed_t == res_t ) THEN
 739         test_passed = .TRUE.
 740       ELSE
 741         CALL ESMF_TimeGet( computed_t, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
 742         CALL test_check_error( ESMF_SUCCESS, rc, &
 743                               TRIM(itestname)//'ESMF_TimeGet() ', &
 744                               __FILE__ , &
 745                               __LINE__  )
 746         ! handle fractions
 747         CALL fraction_to_string( Sn, Sd, frac_str )
 748         computed_str = TRIM(computed_str)//TRIM(frac_str)
 749       ENDIF
 750     ELSE IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
 751       IF ( computed_ti == res_ti ) THEN
 752         test_passed = .TRUE.
 753       ELSE
 754         CALL ESMF_TimeIntervalGet( computed_ti, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
 755         CALL test_check_error( ESMF_SUCCESS, rc, &
 756                               TRIM(itestname)//'ESMF_TimeGet() ', &
 757                               __FILE__ , &
 758                               __LINE__  )
 759         ! handle fractions
 760         CALL fraction_to_string( Sn, Sd, frac_str )
 761         computed_str = TRIM(computed_str)//TRIM(frac_str)
 762       ENDIF
 763     ELSE  ! result is INTEGER
 764       IF ( computed_int == res_int ) THEN
 765         test_passed = .TRUE.
 766       ELSE
 767         IF ( computed_int > 0 ) THEN
 768           WRITE(computed_str,FMT="('+',I8.8)") ABS(computed_int)
 769         ELSE
 770           WRITE(computed_str,FMT="('-',I8.8)") ABS(computed_int)
 771         ENDIF
 772       ENDIF
 773     ENDIF
 774     IF ( test_passed ) THEN
 775       WRITE(*,FMT='(A)') 'PASS:  '//TRIM(itestname)
 776     ELSE
 777       WRITE(*,*) 'FAIL:  ',TRIM(itestname),':  (',                        &
 778         TRIM(res_type_str),' = ',TRIM(op1_type_str),' ',TRIM(op_str),' ', &
 779         TRIM(op2_type_str),')  expected ',                                &
 780         TRIM(res_str),' = ',TRIM(op1_str),' ',TRIM(op_str),' ',           &
 781         TRIM(op2_str),'  but computed ',TRIM(computed_str)
 782     ENDIF
 783
 784   END SUBROUTINE test_arithmetic
 785
 786
 787
 788   ! simple clock creation and advance with add-subtract tests thrown in
 789   ! no self checks (yet)
 790   SUBROUTINE test_clock_advance(                                              &
 791     start_yy, start_mm, start_dd, start_h, start_m, start_s,                  &
 792      stop_yy,  stop_mm,  stop_dd,  stop_h,  stop_m,  stop_s,                  &
 793     timestep_d, timestep_h, timestep_m, timestep_s, timestep_sn, timestep_sd, &
 794     testname, increment_S, increment_Sn, increment_Sd )
 795     INTEGER, INTENT(IN), OPTIONAL :: start_YY
 796     INTEGER, INTENT(IN), OPTIONAL :: start_MM  ! month
 797     INTEGER, INTENT(IN), OPTIONAL :: start_DD  ! day of month
 798     INTEGER, INTENT(IN), OPTIONAL :: start_H
 799     INTEGER, INTENT(IN), OPTIONAL :: start_M
 800     INTEGER, INTENT(IN), OPTIONAL :: start_S
 801     INTEGER, INTENT(IN), OPTIONAL :: stop_YY
 802     INTEGER, INTENT(IN), OPTIONAL :: stop_MM  ! month
 803     INTEGER, INTENT(IN), OPTIONAL :: stop_DD  ! day of month
 804     INTEGER, INTENT(IN), OPTIONAL :: stop_H
 805     INTEGER, INTENT(IN), OPTIONAL :: stop_M
 806     INTEGER, INTENT(IN), OPTIONAL :: stop_S
 807     INTEGER, INTENT(IN), OPTIONAL :: timestep_D  ! day
 808     INTEGER, INTENT(IN), OPTIONAL :: timestep_H
 809     INTEGER, INTENT(IN), OPTIONAL :: timestep_M
 810     INTEGER, INTENT(IN), OPTIONAL :: timestep_S
 811     INTEGER, INTENT(IN), OPTIONAL :: timestep_Sn
 812     INTEGER, INTENT(IN), OPTIONAL :: timestep_Sd
 813     CHARACTER (LEN=*), OPTIONAL, INTENT(IN) :: testname
 814     INTEGER, INTENT(IN), OPTIONAL :: increment_S  ! add and subtract this
 815     INTEGER, INTENT(IN), OPTIONAL :: increment_Sn ! value each time step
 816     INTEGER, INTENT(IN), OPTIONAL :: increment_Sd
 817
 818     ! locals
 819     INTEGER :: istart_YY
 820     INTEGER :: istart_MM  ! month
 821     INTEGER :: istart_DD  ! day of month
 822     INTEGER :: istart_H
 823     INTEGER :: istart_M
 824     INTEGER :: istart_S
 825     INTEGER :: istop_YY
 826     INTEGER :: istop_MM  ! month
 827     INTEGER :: istop_DD  ! day of month
 828     INTEGER :: istop_H
 829     INTEGER :: istop_M
 830     INTEGER :: istop_S
 831     INTEGER :: itimestep_D  ! day
 832     INTEGER :: itimestep_H
 833     INTEGER :: itimestep_M
 834     INTEGER :: itimestep_S
 835     INTEGER :: itimestep_Sn
 836     INTEGER :: itimestep_Sd
 837     CHARACTER (LEN=512) :: itestname, itestfullname
 838     INTEGER :: iincrement_S
 839     INTEGER :: iincrement_Sn
 840     INTEGER :: iincrement_Sd
 841     INTEGER :: Sn, Sd
 842     INTEGER rc
 843     TYPE(ESMF_Time)           :: start_time, stop_time, current_time
 844     TYPE(ESMF_Clock), POINTER :: domain_clock
 845     TYPE(ESMF_TimeInterval)   :: timestep, increment
 846     TYPE(ESMF_Time)           :: add_time, subtract_time
 847     INTEGER :: itimestep
 848     REAL(ESMF_KIND_R8) :: dayr8
 849     CHARACTER(LEN=ESMF_MAXSTR) :: str, frac_str
 850
 851     istart_YY = 0
 852     istart_MM = 1
 853     istart_DD = 1
 854     istart_H = 0
 855     istart_M = 0
 856     istart_S = 0
 857     istop_YY = 0
 858     istop_MM = 1
 859     istop_DD = 1
 860     istop_H = 0
 861     istop_M = 0
 862     istop_S = 0
 863     itimestep_D = 0
 864     itimestep_H = 0
 865     itimestep_M = 0
 866     itimestep_S = 0
 867     itimestep_Sn = 0
 868     itimestep_Sd = 0
 869     itestname = ''
 870     iincrement_S = 0
 871     iincrement_Sn = 0
 872     iincrement_Sd = 0
 873
 874     IF ( PRESENT( start_YY ) ) istart_YY = start_YY
 875     IF ( PRESENT( start_MM ) ) istart_MM = start_MM
 876     IF ( PRESENT( start_DD ) ) istart_DD = start_DD
 877     IF ( PRESENT( start_H ) ) istart_H = start_H
 878     IF ( PRESENT( start_M ) ) istart_M = start_M
 879     IF ( PRESENT( start_S ) ) istart_S = start_S
 880     IF ( PRESENT( stop_YY ) ) istop_YY = stop_YY
 881     IF ( PRESENT( stop_MM ) ) istop_MM = stop_MM
 882     IF ( PRESENT( stop_DD ) ) istop_DD = stop_DD
 883     IF ( PRESENT( stop_H ) ) istop_H = stop_H
 884     IF ( PRESENT( stop_M ) ) istop_M = stop_M
 885     IF ( PRESENT( stop_S ) ) istop_S = stop_S
 886     IF ( PRESENT( timestep_D ) ) itimestep_D = timestep_D
 887     IF ( PRESENT( timestep_H ) ) itimestep_H = timestep_H
 888     IF ( PRESENT( timestep_M ) ) itimestep_M = timestep_M
 889     IF ( PRESENT( timestep_S ) ) itimestep_S = timestep_S
 890     IF ( PRESENT( timestep_Sn ) ) itimestep_Sn = timestep_Sn
 891     IF ( PRESENT( timestep_Sd ) ) itimestep_Sd = timestep_Sd
 892     IF ( PRESENT( testname ) ) itestname = TRIM(testname)//'_'
 893     IF ( PRESENT( increment_S ) ) iincrement_S = increment_S
 894     IF ( PRESENT( increment_Sn ) ) iincrement_Sn = increment_Sn
 895     IF ( PRESENT( increment_Sd ) ) iincrement_Sd = increment_Sd
 896
 897     ! Initialize start time, stop time, time step, clock for simple case.
 898     itestfullname = TRIM(itestname)//'SETUP'
 899     CALL ESMF_TimeSet( start_time, YY=istart_YY, MM=istart_MM, DD=istart_DD , &
 900                                    H=istart_H, M=istart_M, S=istart_S, rc=rc )
 901     CALL test_check_error( ESMF_SUCCESS, rc, &
 902                           TRIM(itestfullname)//'ESMF_TimeSet() ', &
 903                           __FILE__ , &
 904                           __LINE__  )
 905
 906     CALL ESMF_TimeGet( start_time, timeString=str, rc=rc )
 907     CALL test_check_error( ESMF_SUCCESS, rc, &
 908                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
 909                           __FILE__ , &
 910                           __LINE__  )
 911     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  start_time = <',TRIM(str),'>'
 912
 913     CALL ESMF_TimeSet( stop_time, YY=istop_YY, MM=istop_MM, DD=istop_DD , &
 914                                    H=istop_H, M=istop_M, S=istop_S, rc=rc )
 915     CALL test_check_error( ESMF_SUCCESS, rc, &
 916                           TRIM(itestfullname)//'ESMF_TimeSet() ', &
 917                           __FILE__ , &
 918                           __LINE__  )
 919
 920     CALL ESMF_TimeGet( stop_time, timeString=str, rc=rc )
 921     CALL test_check_error( ESMF_SUCCESS, rc, &
 922                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
 923                           __FILE__ , &
 924                           __LINE__  )
 925     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  stop_time = <',TRIM(str),'>'
 926
 927     CALL ESMF_TimeIntervalSet( timestep, D=itimestep_D, H=itimestep_H, &
 928                                          M=itimestep_M, S=itimestep_S, &
 929                                          Sn=itimestep_Sn, Sd=itimestep_Sd, rc=rc )
 930     CALL test_check_error( ESMF_SUCCESS, rc, &
 931                           TRIM(itestfullname)//'ESMF_TimeIntervalSet() ', &
 932                           __FILE__ , &
 933                           __LINE__  )
 934
 935     CALL ESMF_TimeIntervalGet( timestep, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
 936     CALL test_check_error( ESMF_SUCCESS, rc, &
 937                           TRIM(itestfullname)//'ESMF_TimeIntervalGet() ', &
 938                           __FILE__ , &
 939                           __LINE__  )
 940     ! handle fractions
 941     CALL fraction_to_string( Sn, Sd, frac_str )
 942     str = TRIM(str)//TRIM(frac_str)
 943     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  timestep = <',TRIM(str),'>'
 944
 945     CALL ESMF_TimeIntervalSet( increment, S=iincrement_S, &
 946                                Sn=iincrement_Sn, Sd=iincrement_Sd, rc=rc )
 947     CALL test_check_error( ESMF_SUCCESS, rc, &
 948                           TRIM(itestfullname)//'ESMF_TimeIntervalSet() ', &
 949                           __FILE__ , &
 950                           __LINE__  )
 951
 952     CALL ESMF_TimeIntervalGet( increment, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
 953     CALL test_check_error( ESMF_SUCCESS, rc, &
 954                           TRIM(itestfullname)//'ESMF_TimeIntervalGet() ', &
 955                           __FILE__ , &
 956                           __LINE__  )
 957     ! handle fractions
 958     CALL fraction_to_string( Sn, Sd, frac_str )
 959     str = TRIM(str)//TRIM(frac_str)
 960     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  increment = <',TRIM(str),'>'
 961
 962     ALLOCATE( domain_clock )
 963     domain_clock = ESMF_ClockCreate( TimeStep= timestep,  &
 964                                      StartTime=start_time, &
 965                                      StopTime= stop_time,  &
 966                                      rc=rc )
 967     CALL test_check_error( ESMF_SUCCESS, rc, &
 968                           TRIM(itestfullname)//'ESMF_ClockCreate() ', &
 969                           __FILE__ , &
 970                           __LINE__  )
 971
 972     CALL ESMF_ClockGet( domain_clock, CurrTime=current_time, &
 973                         rc=rc )
 974     CALL test_check_error( ESMF_SUCCESS, rc, &
 975                           TRIM(itestfullname)//'ESMF_ClockGet() ', &
 976                           __FILE__ , &
 977                           __LINE__  )
 978
 979     CALL ESMF_TimeGet( current_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
 980     CALL test_check_error( ESMF_SUCCESS, rc, &
 981                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
 982                           __FILE__ , &
 983                           __LINE__  )
 984     CALL fraction_to_string( Sn, Sd, frac_str )
 985     str = TRIM(str)//TRIM(frac_str)
 986     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  clock current_time = <',TRIM(str),'>'
 987
 988     CALL ESMF_TimeGet( current_time, dayOfYear_r8=dayr8, rc=rc )
 989     CALL test_check_error( ESMF_SUCCESS, rc, &
 990                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
 991                           __FILE__ , &
 992                           __LINE__  )
 993     WRITE(*,FMT='(A,A,F10.6,A)') TRIM(itestfullname),':  current_time dayOfYear_r8 = < ',dayr8,' >'
 994
 995     subtract_time = current_time - increment
 996     CALL ESMF_TimeGet( subtract_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
 997     CALL test_check_error( ESMF_SUCCESS, rc, &
 998                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
 999                           __FILE__ , &
1000                           __LINE__  )
1001     CALL fraction_to_string( Sn, Sd, frac_str )
1002     str = TRIM(str)//TRIM(frac_str)
1003     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time-increment = <',TRIM(str),'>'
1004
1005     add_time = current_time + increment
1006     CALL ESMF_TimeGet( add_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
1007     CALL test_check_error( ESMF_SUCCESS, rc, &
1008                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
1009                           __FILE__ , &
1010                           __LINE__  )
1011     CALL fraction_to_string( Sn, Sd, frac_str )
1012     str = TRIM(str)//TRIM(frac_str)
1013     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time+increment = <',TRIM(str),'>'
1014
1015     ! Advance clock.
1016     itestfullname = TRIM(itestname)//'ADVANCE'
1017     itimestep = 0
1018     DO WHILE ( .NOT. ESMF_ClockIsStopTime(domain_clock ,rc=rc) )
1019       CALL test_check_error( ESMF_SUCCESS, rc, &
1020                             TRIM(itestfullname)//'ESMF_ClockIsStopTime() ', &
1021                             __FILE__ , &
1022                             __LINE__  )
1023       itimestep = itimestep + 1
1024
1025       CALL ESMF_ClockAdvance( domain_clock, rc=rc )
1026       CALL test_check_error( ESMF_SUCCESS, rc, &
1027                             TRIM(itestfullname)//'ESMF_ClockAdvance() ', &
1028                             __FILE__ , &
1029                             __LINE__  )
1030
1031       CALL ESMF_ClockGet( domain_clock, CurrTime=current_time, &
1032                           rc=rc )
1033       CALL test_check_error( ESMF_SUCCESS, rc, &
1034                             TRIM(itestfullname)//'ESMF_ClockGet() ', &
1035                             __FILE__ , &
1036                             __LINE__  )
1037
1038       CALL ESMF_TimeGet( current_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
1039       CALL test_check_error( ESMF_SUCCESS, rc, &
1040                             TRIM(itestfullname)//'ESMF_TimeGet() ', &
1041                             __FILE__ , &
1042                             __LINE__  )
1043       CALL fraction_to_string( Sn, Sd, frac_str )
1044       str = TRIM(str)//TRIM(frac_str)
1045       WRITE(*,FMT='(A,A,I6.6,A,A,A)') TRIM(itestfullname),':  count = ', &
1046         itimestep,'  current_time = <',TRIM(str),'>'
1047
1048       subtract_time = current_time - increment
1049       CALL ESMF_TimeGet( subtract_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
1050       CALL test_check_error( ESMF_SUCCESS, rc, &
1051                             TRIM(itestfullname)//'ESMF_TimeGet() ', &
1052                             __FILE__ , &
1053                             __LINE__  )
1054       CALL fraction_to_string( Sn, Sd, frac_str )
1055       str = TRIM(str)//TRIM(frac_str)
1056       WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time-increment = <',TRIM(str),'>'
1057
1058       add_time = current_time + increment
1059       CALL ESMF_TimeGet( add_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
1060       CALL test_check_error( ESMF_SUCCESS, rc, &
1061                             TRIM(itestfullname)//'ESMF_TimeGet() ', &
1062                             __FILE__ , &
1063                             __LINE__  )
1064       CALL fraction_to_string( Sn, Sd, frac_str )
1065       str = TRIM(str)//TRIM(frac_str)
1066       WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time+increment = <',TRIM(str),'>'
1067
1068     ENDDO
1069
1070     DEALLOCATE( domain_clock )
1071
1072   END SUBROUTINE test_clock_advance
1073
1074 END MODULE my_tests
1075
1076
1077 #if defined( TIME_F90_ONLY )
1078
1079 ! TBH:  Improve the build of Test1.exe to use WRF versions of these
1080 ! TBH:  routines and remove these hacked-in duplicates!!
1081
1082 SUBROUTINE wrf_abort
1083   IMPLICIT NONE
1084 #if defined( DM_PARALLEL ) && ! defined( STUBMPI )
1085   INCLUDE 'mpif.h'
1086   INTEGER ierr
1087   CALL mpi_abort(MPI_COMM_WORLD,1,ierr)
1088 #else
1089   STOP
1090 #endif
1091 END SUBROUTINE wrf_abort
1092
1093 SUBROUTINE wrf_message( str )
1094   IMPLICIT NONE
1095   CHARACTER*(*) str
1096 #if defined( DM_PARALLEL ) && ! defined( STUBMPI)
1097   write(0,*) str
1098 #endif
1099   print*, str
1100 END SUBROUTINE wrf_message
1101
1102 ! intentionally write to stderr only
1103 SUBROUTINE wrf_message2( str )
1104   IMPLICIT NONE
1105   CHARACTER*(*) str
1106   write(0,*) str
1107 END SUBROUTINE wrf_message2
1108
1109 SUBROUTINE wrf_error_fatal3( file_str, line, str )
1110   USE my_tests
1111   IMPLICIT NONE
1112   CHARACTER*(*) file_str
1113   INTEGER , INTENT (IN) :: line  ! only print file and line if line > 0
1114   CHARACTER*(*) str
1115   CHARACTER*256 :: line_str
1116   write(line_str,'(i6)') line
1117   ! special behavior for testing since Fortran cannot catch exceptions
1118  IF ( WRF_ERROR_FATAL_PRINT ) THEN
1119   ! just print message and continue
1120   CALL wrf_message( 'ERROR IN FILE:  '//TRIM(file_str)//'  LINE:  '//TRIM(line_str) )
1121  ELSE
1122   ! normal behavior
1123 #if defined( DM_PARALLEL ) && ! defined( STUBMPI )
1124   CALL wrf_message( '-------------- FATAL CALLED ---------------' )
1125   ! only print file and line if line is positive
1126   IF ( line > 0 ) THEN
1127     CALL wrf_message( 'FATAL CALLED FROM FILE:  '//file_str//'  LINE:  '//TRIM(line_str) )
1128   ENDIF
1129   CALL wrf_message( str )
1130   CALL wrf_message( '-------------------------------------------' )
1131 #else
1132   CALL wrf_message2( '-------------- FATAL CALLED ---------------' )
1133   ! only print file and line if line is positive
1134   IF ( line > 0 ) THEN
1135     CALL wrf_message( 'FATAL CALLED FROM FILE:  '//file_str//'  LINE:  '//TRIM(line_str) )
1136   ENDIF
1137   CALL wrf_message2( str )
1138   CALL wrf_message2( '-------------------------------------------' )
1139 #endif
1140   CALL wrf_abort
1141  ENDIF
1142 END SUBROUTINE wrf_error_fatal3
1143
1144 SUBROUTINE wrf_error_fatal( str )
1145   IMPLICIT NONE
1146   CHARACTER*(*) str
1147   CALL wrf_error_fatal3 ( ' ', 0, str )
1148 END SUBROUTINE wrf_error_fatal
1149
1150 #endif
1151
1152
1153 ! Check to see if expected value == actual value
1154 ! If not, print message and exit.
1155 SUBROUTINE test_check_error( expected, actual, str, file_str, line )
1156   IMPLICIT NONE
1157   INTEGER , INTENT (IN) :: expected
1158   INTEGER , INTENT (IN) :: actual
1159   CHARACTER*(*) str
1160   CHARACTER*(*) file_str
1161   INTEGER , INTENT (IN) :: line
1162   CHARACTER (LEN=512)   :: rc_str
1163   CHARACTER (LEN=512)   :: str_with_rc
1164   IF ( expected .ne. actual ) THEN
1165     WRITE (rc_str,*) '  Routine returned error code = ',actual
1166     str_with_rc = 'FAIL:  '//TRIM(str)//TRIM(rc_str)
1167     CALL wrf_error_fatal3( file_str, line, str_with_rc )
1168   ENDIF
1169 END SUBROUTINE test_check_error
1170
1171
1172
1173 PROGRAM time_manager_test
1174   USE ESMF_Mod
1175   USE my_tests
1176   IMPLICIT NONE
1177   INTEGER :: rc
1178
1179   PRINT *,'BEGIN TEST SUITE'
1180
1181   CALL ESMF_Initialize( defaultCalendar=ESMF_CAL_GREGORIAN, rc=rc )
1182   CALL test_check_error( ESMF_SUCCESS, rc, &
1183                         'ESMF_Initialize() ', &
1184                         __FILE__ , &
1185                         __LINE__  )
1186 !  PRINT *,'DEBUG:  back from ESMF_Initialize(), rc = ',rc
1187
1188 !  CALL test_print(  t_yy,  t_mm,  t_dd,  t_h,  t_m,  t_s, &
1189 !                   ti_yy, ti_mm, ti_dd, ti_h, ti_m, ti_s, &
1190 !                   res_str, testname )
1191
1192   ! Print times
1193   ! "vanilla" tests
1194 !  PRINT *,'DEBUG:  calling 1st test_print()'
1195   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
1196     res_str='2001-12-03_01:20:10', testname='printT_1' )
1197 !  PRINT *,'DEBUG:  back from 1st test_print()'
1198   CALL test_print( t_yy=0,  t_mm=1,  t_dd=1,  t_h=0,  t_m=0,  t_s=0, &
1199     res_str='0000-01-01_00:00:00', testname='printT_2' )
1200   CALL test_print( t_yy=2003,  t_mm=12,  t_dd=30,  t_h=23,  t_m=59,  t_s=50, &
1201     res_str='2003-12-30_23:59:50', testname='printT_3' )
1202   CALL test_print( t_yy=2003,  t_mm=12,  t_dd=31,  t_h=23,  t_m=59,  t_s=50, &
1203     res_str='2003-12-31_23:59:50', testname='printT_4' )
1204   CALL test_print( t_yy=2004,  t_mm=12,  t_dd=30,  t_h=23,  t_m=59,  t_s=50, &
1205     res_str='2004-12-30_23:59:50', testname='printT_5' )
1206   CALL test_print( t_yy=2004,  t_mm=12,  t_dd=31,  t_h=23,  t_m=59,  t_s=50, &
1207     res_str='2004-12-31_23:59:50', testname='printT_6' )
1208 !$$$  NOTE that this fails -- need to fix up output string for negative year
1209 !  CALL test_print( t_yy=-2004,  t_mm=12,  t_dd=31,  t_h=23,  t_m=59,  t_s=50, &
1210 !    res_str='-2004-12-31_23:59:50', testname='printT_6' )
1211
1212   ! these test default behavior of test harness
1213   CALL test_print( t_s=0, &
1214     res_str='0000-01-01_00:00:00', testname='printT_D1' )
1215   CALL test_print( t_yy=0, &
1216     res_str='0000-01-01_00:00:00', testname='printT_D2' )
1217
1218   ! fractions
1219   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
1220     t_sn=1, t_sd=3, &
1221     res_str='2001-12-03_01:20:10+01/03', testname='printT_F1' )
1222   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
1223     t_sn=4, t_sd=3, &
1224     res_str='2001-12-03_01:20:11+01/03', testname='printT_F2' )
1225   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
1226     t_sn=12, t_sd=3, &
1227     res_str='2001-12-03_01:20:14', testname='printT_F3' )
1228   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
1229     t_sn=-1, t_sd=3, &
1230     res_str='2001-12-03_01:20:09+02/03', testname='printT_F4' )
1231
1232   ! ERROR, MM out of range
1233 !$$$here...  fix so this just prints "ERROR:  <testname>" in failure case
1234 !$$$here...  also need "expect_fail" to reverse sense of PASS/FAIL message for
1235 !$$$here...  tests that should fail
1236 !  CALL test_print( t_yy=2001,  t_mm=13,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
1237 !    res_str='2002-01-03_01:20:10', testname='printT_E1', expect_error=.TRUE. )
1238
1239   ! Print time intervals
1240   ! "vanilla" tests
1241   CALL test_print( ti_yy=0,  ti_mm=0,  ti_dd=0,  ti_h=0,  ti_m=0,  ti_s=0, &
1242     res_str='0000000000_000:000:000', testname='printTI_1' )
1243   CALL test_print( ti_yy=0,  ti_mm=0,  ti_dd=500,  ti_h=0,  ti_m=0,  ti_s=7270, &
1244     res_str='0000000500_002:001:010', testname='printTI_2' )
1245
1246   ! these test default behavior of test harness
1247   CALL test_print( ti_s=0, &
1248     res_str='0000000000_000:000:000', testname='printTI_D1' )
1249   CALL test_print( ti_yy=0, &
1250     res_str='0000000000_000:000:000', testname='printTI_D2' )
1251
1252   ! these test negative values
1253   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=-1,  ti_m=-20,  ti_s=-10, &
1254     res_str='-0000000003_001:020:010', testname='printTI_N1' )
1255
1256   ! these test mixed values
1257   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=1,  ti_m=20,  ti_s=10, &
1258     res_str='-0000000002_022:039:050', testname='printTI_M1' )
1259
1260   ! fractions
1261   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
1262     ti_sn=1, ti_sd=3, &
1263     res_str='0000000003_001:020:010+01/03', testname='printTI_F1' )
1264   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
1265     ti_sn=5, ti_sd=3, &
1266     res_str='0000000003_001:020:011+02/03', testname='printTI_F2' )
1267   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=-1,  ti_m=-20,  ti_s=-10, &
1268     ti_sn=-1, ti_sd=3, &
1269     res_str='-0000000003_001:020:010-01/03', testname='printTI_F3' )
1270   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=-1,  ti_m=-20,  ti_s=-10, &
1271     ti_sn=1, ti_sd=3, &
1272     res_str='-0000000003_001:020:009-02/03', testname='printTI_F4' )
1273
1274   ! these test non-normalized values
1275 !  CALL test_print( ti_yy=2001,  ti_mm=1,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
1276 !    res_str='02001-001-003_001:020:010', testname='printTI_NN1', expect_error=.TRUE. )
1277 !  CALL test_print( ti_yy=2001,  ti_mm=12,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
1278 !    res_str='02002-000-003_001:020:010', testname='printTI_NN2', expect_error=.TRUE. )
1279 !  CALL test_print( ti_yy=2002,  ti_mm=5,  ti_dd=500,  ti_h=0,  ti_m=0,  ti_s=7270, &
1280 !    res_str='02002-005-500_002:001:010', testname='printTI_NN3', expect_error=.TRUE. )
1281
1282   ! Addition tests
1283   ! ESMF_Time = ESMF_Time + ESMF_TimeInterval
1284   CALL test_arithmetic( add_op=.TRUE.,                                             &
1285      op1_t_yy=2001,  op1_t_mm=12,  op1_t_dd=3,  op1_t_h=1,  op1_t_m=20,  op1_t_s=10, &
1286     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=3, op2_ti_m=10, op2_ti_s=10, &
1287      res_t_yy=2001,  res_t_mm=12,  res_t_dd=3,  res_t_h=4,  res_t_m=30,  res_t_s=20, &
1288     testname='AddT_T_TI1' )
1289   CALL test_arithmetic( add_op=.TRUE.,                                             &
1290      op1_t_yy=2001,  op1_t_mm=12,  op1_t_dd=31,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1291     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1292      res_t_yy=2002,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1293     testname='AddT_T_TI2' )
1294   CALL test_arithmetic( add_op=.TRUE.,                                             &
1295      op1_t_yy=2003,  op1_t_mm=12,  op1_t_dd=31,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1296     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1297      res_t_yy=2004,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1298     testname='AddT_T_TI3' )
1299   CALL test_arithmetic( add_op=.TRUE.,                                             &
1300      op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=31,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1301     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1302      res_t_yy=2005,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1303     testname='AddT_T_TI4' )
1304   ! this case hung after the CCSM contribution
1305   CALL test_arithmetic( add_op=.TRUE.,                                             &
1306      op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1307     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1308      res_t_yy=2004,  res_t_mm=12,  res_t_dd=31,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1309     testname='AddT_T_TI5' )
1310 ! NOTE:  CCSM folks need to decide what it means to add "1 month" to Feb. 29.  And all the
1311 !        other very similar cases.  Then, write this unit test!
1312 !  CALL test_arithmetic( add_op=.TRUE.,                                             &
1313 !     op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=31,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1314 !    op2_ti_yy=   2, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1315 !     res_t_yy=2007,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1316 !    testname='AddT_T_TI6' )
1317   CALL test_arithmetic( add_op=.TRUE.,                                             &
1318      op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1319     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=365, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1320      res_t_yy=2005,  res_t_mm=12,  res_t_dd=30,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1321     testname='AddT_T_TI7' )
1322   CALL test_arithmetic( add_op=.TRUE.,                                             &
1323      op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1324     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=367, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1325      res_t_yy=2006,  res_t_mm=01,  res_t_dd=01,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1326     testname='AddT_T_TI8' )
1327   CALL test_arithmetic( add_op=.TRUE.,                                             &
1328      op1_t_yy=2003,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1329     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=365, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1330      res_t_yy=2004,  res_t_mm=12,  res_t_dd=29,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1331     testname='AddT_T_TI9' )
1332   CALL test_arithmetic( add_op=.TRUE.,                                             &
1333      op1_t_yy=2003,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1334     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=366, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1335      res_t_yy=2004,  res_t_mm=12,  res_t_dd=30,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1336     testname='AddT_T_TI10' )
1337   CALL test_arithmetic( add_op=.TRUE.,                                             &
1338      op1_t_yy=2003,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1339     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=367, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1340      res_t_yy=2004,  res_t_mm=12,  res_t_dd=31,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1341     testname='AddT_T_TI11' )
1342   CALL test_arithmetic( add_op=.TRUE.,                                             &
1343      op1_t_yy=2003,  op1_t_mm=12,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1344     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=368, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1345      res_t_yy=2005,  res_t_mm=01,  res_t_dd=01,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1346     testname='AddT_T_TI12' )
1347   CALL test_arithmetic( add_op=.TRUE.,                                             &
1348      op1_t_yy=2004,  op1_t_mm=03,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1349     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=365, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1350      res_t_yy=2005,  res_t_mm=03,  res_t_dd=30,  res_t_h=8,  res_t_m=40,  res_t_s=10, &
1351     testname='AddT_T_TI13' )
1352   CALL test_arithmetic( add_op=.TRUE.,                                             &
1353      op1_t_yy=2004,  op1_t_mm=03,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1354     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=365, op2_ti_h=22, op2_ti_m=10, op2_ti_s=10, &
1355      res_t_yy=2005,  res_t_mm=03,  res_t_dd=31,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1356     testname='AddT_T_TI14' )
1357   CALL test_arithmetic( add_op=.TRUE.,                                             &
1358      op1_t_yy=2004,  op1_t_mm=03,  op1_t_dd=30,  op1_t_h=4,  op1_t_m=30,  op1_t_s=00, &
1359     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=366, op2_ti_h=22, op2_ti_m=10, op2_ti_s=10, &
1360      res_t_yy=2005,  res_t_mm=04,  res_t_dd=01,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1361     testname='AddT_T_TI15' )
1362   ! ESMF_Time = ESMF_Time + ESMF_TimeInterval with fractions
1363   CALL test_arithmetic( add_op=.TRUE.,                                             &
1364      op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=31,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1365      op1_t_sn=01,  op1_t_sd=03, &
1366     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1367     op2_ti_sn=01, op2_ti_sd=03, &
1368      res_t_yy=2005,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1369      res_t_sn=02,  res_t_sd=03, &
1370     testname='AddT_T_TI_F1' )
1371   ! this should fail (and does)
1372 !  CALL test_arithmetic( add_op=.TRUE.,                                             &
1373 !     op1_t_yy=2004,  op1_t_mm=12,  op1_t_dd=31,  op1_t_h=22,  op1_t_m=30,  op1_t_s=00, &
1374 !     op1_t_sn=01,  op1_t_sd=03, &
1375 !    op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1376 !    op2_ti_sn=01, op2_ti_sd=03, &
1377 !     res_t_yy=2005,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1378 !     res_t_sn=01,  res_t_sd=03, &
1379 !    testname='AddT_T_TI_F2' )
1380   ! ESMF_Time = ESMF_TimeInterval + ESMF_Time
1381   CALL test_arithmetic( add_op=.TRUE.,                                             &
1382     op1_ti_yy=   0, op1_ti_mm= 0, op1_ti_dd=0, op1_ti_h=3, op1_ti_m=10, op1_ti_s=10, &
1383      op2_t_yy=2001,  op2_t_mm=12,  op2_t_dd=3,  op2_t_h=1,  op2_t_m=20,  op2_t_s=10, &
1384      res_t_yy=2001,  res_t_mm=12,  res_t_dd=3,  res_t_h=4,  res_t_m=30,  res_t_s=20, &
1385     testname='AddT_TI_T1' )
1386   CALL test_arithmetic( add_op=.TRUE.,                                             &
1387     op1_ti_yy=   0, op1_ti_mm= 0, op1_ti_dd=0, op1_ti_h=4, op1_ti_m=10, op1_ti_s=10, &
1388      op2_t_yy=2001,  op2_t_mm=12,  op2_t_dd=31,  op2_t_h=22,  op2_t_m=30,  op2_t_s=00, &
1389      res_t_yy=2002,  res_t_mm= 1,  res_t_dd=1,  res_t_h=2,  res_t_m=40,  res_t_s=10, &
1390     testname='AddT_TI_T2' )
1391   ! ESMF_TimeInterval = ESMF_TimeInterval + ESMF_TimeInterval
1392   CALL test_arithmetic( add_op=.TRUE.,                                             &
1393     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=3, op1_ti_h=1, op1_ti_m=20, op1_ti_s=10, &
1394     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=1, op2_ti_h=1, op2_ti_m=10, op2_ti_s=10, &
1395     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=4, res_ti_h=2, res_ti_m=30, res_ti_s=20, &
1396     testname='AddTI_TI_TI1' )
1397   CALL test_arithmetic( add_op=.TRUE.,                                             &
1398     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=-3, op1_ti_h=-1, op1_ti_m=-20, op1_ti_s=-10, &
1399     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=1, op2_ti_h=1, op2_ti_m=10, op2_ti_s=10, &
1400     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=-2, res_ti_h=0, res_ti_m=-10, res_ti_s=00, &
1401     testname='AddTI_TI_TI2' )
1402   CALL test_arithmetic( add_op=.TRUE.,                                             &
1403     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=-3, op1_ti_h=-1, op1_ti_m=-20, op1_ti_s=-10, &
1404     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=-1, op2_ti_h=-1, op2_ti_m=-10, op2_ti_s=-10, &
1405     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=-4, res_ti_h=-2, res_ti_m=-30, res_ti_s=-20, &
1406     testname='AddTI_TI_TI3' )
1407
1408   ! Subtraction tests
1409   ! ESMF_Time = ESMF_Time - ESMF_TimeInterval
1410   CALL test_arithmetic( add_op=.FALSE.,                                            &
1411      op1_t_yy=2001,  op1_t_mm=12,  op1_t_dd=3,  op1_t_h=1,  op1_t_m=20,  op1_t_s=10, &
1412     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=3, op2_ti_m=10, op2_ti_s=10, &
1413      res_t_yy=2001,  res_t_mm=12,  res_t_dd=2,  res_t_h=22, res_t_m=10,  res_t_s=0,  &
1414     testname='SubtractT_T_TI1' )
1415   CALL test_arithmetic( add_op=.FALSE.,                                            &
1416      op1_t_yy=2005,  op1_t_mm=1,   op1_t_dd=1,  op1_t_h=0,  op1_t_m=00,  op1_t_s=0,  &
1417     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=0, op2_ti_m=00, op2_ti_s=10, &
1418      res_t_yy=2004,  res_t_mm=12,  res_t_dd=31, res_t_h=23, res_t_m=59,  res_t_s=50, &
1419     testname='SubtractT_T_TI2' )
1420   CALL test_arithmetic( add_op=.FALSE.,                                            &
1421      op1_t_yy=2004,  op1_t_mm=1,   op1_t_dd=1,  op1_t_h=0,  op1_t_m=00,  op1_t_s=0,  &
1422     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=0, op2_ti_m=00, op2_ti_s=10, &
1423      res_t_yy=2003,  res_t_mm=12,  res_t_dd=31, res_t_h=23, res_t_m=59,  res_t_s=50, &
1424     testname='SubtractT_T_TI3' )
1425   CALL test_arithmetic( add_op=.FALSE.,                                            &
1426      op1_t_yy=2003,  op1_t_mm=1,   op1_t_dd=1,  op1_t_h=0,  op1_t_m=00,  op1_t_s=0,  &
1427     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=0, op2_ti_m=00, op2_ti_s=10, &
1428      res_t_yy=2002,  res_t_mm=12,  res_t_dd=31, res_t_h=23, res_t_m=59,  res_t_s=50, &
1429     testname='SubtractT_T_TI4' )
1430   CALL test_arithmetic( add_op=.FALSE.,                                             &
1431      op1_t_yy=2005,  op1_t_mm=04,  op1_t_dd=01,  op1_t_h=2,  op1_t_m=40,  op1_t_s=10, &
1432     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=366, op2_ti_h=22, op2_ti_m=10, op2_ti_s=10, &
1433      res_t_yy=2004,  res_t_mm=03,  res_t_dd=30,  res_t_h=4,  res_t_m=30,  res_t_s=00, &
1434     testname='SubtractT_T_TI5' )
1435   CALL test_arithmetic( add_op=.FALSE.,                                             &
1436      op1_t_yy=2006,  op1_t_mm=01,  op1_t_dd=01,  op1_t_h=8,  op1_t_m=40,  op1_t_s=10, &
1437     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=367, op2_ti_h=4, op2_ti_m=10, op2_ti_s=10, &
1438      res_t_yy=2004,  res_t_mm=12,  res_t_dd=30,  res_t_h=4,  res_t_m=30,  res_t_s=00, &
1439     testname='SubtractT_T_TI6' )
1440   ! ESMF_Time = ESMF_Time - ESMF_TimeInterval with fractions
1441   CALL test_arithmetic( add_op=.FALSE.,                                             &
1442      op1_t_yy=2005,  op1_t_mm=01,  op1_t_dd=01,  op1_t_h=00,  op1_t_m=00,  op1_t_s=00, &
1443      op1_t_sn=00,  op1_t_sd=00, &
1444     op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=0, op2_ti_m=00, op2_ti_s=01, &
1445     op2_ti_sn=01, op2_ti_sd=03, &
1446      res_t_yy=2004,  res_t_mm=12,  res_t_dd=31,  res_t_h=23,  res_t_m=59,  res_t_s=58, &
1447      res_t_sn=02,  res_t_sd=03, &
1448     testname='SubtractT_T_TI_F1' )
1449   ! ESMF_TimeInterval = ESMF_Time - ESMF_Time
1450   CALL test_arithmetic( add_op=.FALSE.,                                            &
1451      op1_t_yy=2001,  op1_t_mm=12,  op1_t_dd=3,  op1_t_h=1,  op1_t_m=20,  op1_t_s=10, &
1452      op2_t_yy=2001,  op2_t_mm=12,  op2_t_dd=1,  op2_t_h=1,  op2_t_m=10,  op2_t_s=10, &
1453     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=2, res_ti_h=0, res_ti_m=10, res_ti_s=0,  &
1454     testname='SubtractTI_T_T1' )
1455   CALL test_arithmetic( add_op=.FALSE.,                                            &
1456      op1_t_yy=2002,  op1_t_mm=1,   op1_t_dd=1,  op1_t_h=0,  op1_t_m=00,  op1_t_s=00, &
1457      op2_t_yy=2001,  op2_t_mm=12,  op2_t_dd=31, op2_t_h=23, op2_t_m=59,  op2_t_s=50, &
1458     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=0, res_ti_h=0, res_ti_m=00, res_ti_s=10, &
1459     testname='SubtractTI_T_T2' )
1460   CALL test_arithmetic( add_op=.FALSE.,                                            &
1461      op1_t_yy=2005,  op1_t_mm=1,   op1_t_dd=1,  op1_t_h=0,  op1_t_m=00,  op1_t_s=00, &
1462      op2_t_yy=2004,  op2_t_mm=12,  op2_t_dd=31, op2_t_h=23, op2_t_m=59,  op2_t_s=50, &
1463     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=0, res_ti_h=0, res_ti_m=00, res_ti_s=10, &
1464     testname='SubtractTI_T_T3' )
1465   CALL test_arithmetic( add_op=.FALSE.,                                            &
1466      op1_t_yy=2003,  op1_t_mm=03,  op1_t_dd=01, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1467      op2_t_yy=2003,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=23, op2_t_m=59,  op2_t_s=50, &
1468     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=0, res_ti_h=0, res_ti_m=00, res_ti_s=10, &
1469     testname='SubtractTI_T_T4' )
1470   CALL test_arithmetic( add_op=.FALSE.,                                            &
1471      op1_t_yy=2004,  op1_t_mm=03,  op1_t_dd=01, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1472      op2_t_yy=2004,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=23, op2_t_m=59,  op2_t_s=50, &
1473     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=1, res_ti_h=0, res_ti_m=00, res_ti_s=10, &
1474     testname='SubtractTI_T_T5' )
1475   CALL test_arithmetic( add_op=.FALSE.,                                            &
1476      op1_t_yy=2002,  op1_t_mm=02,  op1_t_dd=28, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1477      op2_t_yy=2002,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1478     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=0, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1479     testname='SubtractTI_T_T6' )
1480   CALL test_arithmetic( add_op=.FALSE.,                                            &
1481      op1_t_yy=2003,  op1_t_mm=02,  op1_t_dd=28, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1482      op2_t_yy=2002,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1483     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=365, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1484     testname='SubtractTI_T_T7' )
1485   CALL test_arithmetic( add_op=.FALSE.,                                            &
1486      op1_t_yy=2004,  op1_t_mm=02,  op1_t_dd=28, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1487      op2_t_yy=2003,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1488     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=365, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1489     testname='SubtractTI_T_T8' )
1490   CALL test_arithmetic( add_op=.FALSE.,                                            &
1491      op1_t_yy=2005,  op1_t_mm=02,  op1_t_dd=28, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1492      op2_t_yy=2004,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1493     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=366, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1494     testname='SubtractTI_T_T9' )
1495   CALL test_arithmetic( add_op=.FALSE.,                                            &
1496      op1_t_yy=2003,  op1_t_mm=03,  op1_t_dd=01, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1497      op2_t_yy=2002,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1498     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=366, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1499     testname='SubtractTI_T_T10' )
1500   CALL test_arithmetic( add_op=.FALSE.,                                            &
1501      op1_t_yy=2005,  op1_t_mm=03,  op1_t_dd=01, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1502      op2_t_yy=2004,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1503     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=367, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1504     testname='SubtractTI_T_T11' )
1505   CALL test_arithmetic( add_op=.FALSE.,                                            &
1506      op1_t_yy=2005,  op1_t_mm=03,  op1_t_dd=01, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1507      op2_t_yy=2004,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=23, op2_t_m=59,  op2_t_s=50, &
1508     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=366, res_ti_h=0, res_ti_m=00, res_ti_s=10, &
1509     testname='SubtractTI_T_T12' )
1510   CALL test_arithmetic( add_op=.FALSE.,                                            &
1511      op1_t_yy=2004,  op1_t_mm=02,  op1_t_dd=28, op1_t_h=23, op1_t_m=59,  op1_t_s=50, &
1512      op2_t_yy=2005,  op2_t_mm=03,  op2_t_dd=01, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1513     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=-366, res_ti_h=0, res_ti_m=00, res_ti_s=-10, &
1514     testname='SubtractTI_T_T13' )
1515   CALL test_arithmetic( add_op=.FALSE.,                                            &
1516      op1_t_yy=-2002,  op1_t_mm=02,  op1_t_dd=28, op1_t_h=00, op1_t_m=00,  op1_t_s=00, &
1517      op2_t_yy=-2002,  op2_t_mm=02,  op2_t_dd=28, op2_t_h=00, op2_t_m=00,  op2_t_s=00, &
1518     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=0, res_ti_h=0, res_ti_m=00, res_ti_s=00, &
1519     testname='SubtractTI_T_T14' )
1520   ! ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
1521   CALL test_arithmetic( add_op=.FALSE.,                                            &
1522     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=3, op1_ti_h=1, op1_ti_m=20, op1_ti_s=10, &
1523     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=1, op2_ti_h=1, op2_ti_m=10, op2_ti_s=10, &
1524     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=2, res_ti_h=0, res_ti_m=10, res_ti_s=0,  &
1525     testname='SubtractTI_TI_TI1' )
1526   CALL test_arithmetic( add_op=.FALSE.,                                            &
1527     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=3, op1_ti_h=1, op1_ti_m=20, op1_ti_s=10, &
1528     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=-1, op2_ti_h=-1, op2_ti_m=-10, op2_ti_s=-10, &
1529     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=4, res_ti_h=2, res_ti_m=30, res_ti_s=20,  &
1530     testname='SubtractTI_TI_TI2' )
1531   CALL test_arithmetic( add_op=.FALSE.,                                            &
1532     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=-1, op1_ti_h=-1, op1_ti_m=-10, op1_ti_s=-10, &
1533     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=-3, op2_ti_h=-1, op2_ti_m=-20, op2_ti_s=-10, &
1534     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=2, res_ti_h=0, res_ti_m=10, res_ti_s=00,  &
1535     testname='SubtractTI_TI_TI3' )
1536   ! Negative result ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
1537   CALL test_arithmetic( add_op=.FALSE.,                                            &
1538     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=1, op1_ti_h=1, op1_ti_m=10, op1_ti_s=10, &
1539     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=3, op2_ti_h=1, op2_ti_m=20, op2_ti_s=10, &
1540     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=-2, res_ti_h=0, res_ti_m=-10, res_ti_s=0,  &
1541     testname='SubtractTI_TI_TIN1' )
1542   CALL test_arithmetic( add_op=.FALSE.,                                            &
1543     op1_ti_yy=0000, op1_ti_mm=00, op1_ti_dd=-1, op1_ti_h=-1, op1_ti_m=-10, op1_ti_s=-10, &
1544     op2_ti_yy=0000, op2_ti_mm=00, op2_ti_dd=3, op2_ti_h=1, op2_ti_m=20, op2_ti_s=10, &
1545     res_ti_yy=0000, res_ti_mm=00, res_ti_dd=-4, res_ti_h=-2, res_ti_m=-30, res_ti_s=-20,  &
1546     testname='SubtractTI_TI_TIN2' )
1547
1548   ! Un-normalized ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
1549   ! this is an error
1550 !  CALL test_arithmetic( add_op=.FALSE.,                                            &
1551 !    op1_ti_yy=2001, op1_ti_mm=11, op1_ti_dd=3, op1_ti_h=1, op1_ti_m=20, op1_ti_s=10, &
1552 !    op2_ti_yy=2001, op2_ti_mm=11, op2_ti_dd=1, op2_ti_h=1, op2_ti_m=10, op2_ti_s=10, &
1553 !    res_ti_yy=0000, res_ti_mm=00, res_ti_dd=2, res_ti_h=0, res_ti_m=10, res_ti_s=0,  &
1554 !    testname='SubtractTI_TI_TIU1', expect_error=.TRUE. )
1555
1556   ! this one should FAIL, and does
1557 !  CALL test_arithmetic( add_op=.TRUE.,                                             &
1558 !     op1_t_yy=2001,  op1_t_mm=12,  op1_t_dd=3,  op1_t_h=1,  op1_t_m=20,  op1_t_s=10, &
1559 !    op2_ti_yy=   0, op2_ti_mm= 0, op2_ti_dd=0, op2_ti_h=3, op2_ti_m=10, op2_ti_s=10, &
1560 !     res_t_yy=2002,  res_t_mm=12,  res_t_dd=3,  res_t_h=4,  res_t_m=30,  res_t_s=20, &
1561 !    testname='AddTT1' )
1562
1563   ! Multiplication tests
1564   ! ESMF_TimeInterval = ESMF_TimeInterval * INTEGER
1565   CALL test_arithmetic( multiply_op=.TRUE.,                &
1566     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
1567     op2_int=2,                                             &
1568     res_ti_dd=6,  res_ti_h=24, res_ti_m=37,  res_ti_s=06,  &
1569     testname='MultiplyTI_TI_INT1' )
1570   CALL test_arithmetic( multiply_op=.TRUE.,                &
1571     op1_ti_dd=350,  op1_ti_h=23,  op1_ti_m=50,  op1_ti_s=50, &
1572     op2_int=2,                                             &
1573     res_ti_dd=701,  res_ti_h=23, res_ti_m=41,  res_ti_s=40,&
1574     testname='MultiplyTI_TI_INT2' )
1575   CALL test_arithmetic( multiply_op=.TRUE.,                &
1576     op1_ti_s=01, op1_ti_sn=03, op1_ti_sd=04,               &
1577     op2_int=8,                                             &
1578     res_ti_s=14,                                           &
1579     testname='MultiplyTI_TI_INT3' )
1580
1581   ! Division tests
1582   ! ESMF_TimeInterval = ESMF_TimeInterval / INTEGER
1583   CALL test_arithmetic( multiply_op=.FALSE.,               &
1584     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
1585     op2_int=3,                                             &
1586     res_ti_dd=1,  res_ti_h=04, res_ti_m=06,  res_ti_s=11,  &
1587     testname='DivideTI_TI_INT1' )
1588   CALL test_arithmetic( multiply_op=.FALSE.,               &
1589     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
1590     op2_int=4,                                             &
1591     res_ti_dd=0,  res_ti_h=21, res_ti_m=04,  res_ti_s=38,  &
1592     res_ti_sn=1,  res_ti_sd=4,                             &
1593     testname='DivideTI_TI_INT2' )
1594   CALL test_arithmetic( multiply_op=.FALSE.,               &
1595     op1_ti_s=01, op1_ti_sn=03, op1_ti_sd=04,               &
1596     op2_int=5,                                             &
1597     res_ti_s=0, res_ti_sn=7,  res_ti_sd=20,                &
1598     testname='DivideTI_TI_INT3' )
1599   ! INTEGER = ESMF_TimeInterval / ESMF_TimeInterval
1600   ! this operator truncates to whole integers
1601   CALL test_arithmetic( multiply_op=.FALSE.,               &
1602     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
1603     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
1604     res_int=1,                                             &
1605     testname='DivideINT_TI_TI1' )
1606   CALL test_arithmetic( multiply_op=.FALSE.,               &
1607     op1_ti_dd=6,  op1_ti_h=24,  op1_ti_m=36,  op1_ti_s=66, &
1608     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
1609     res_int=2,                                             &
1610     testname='DivideINT_TI_TI2' )
1611   CALL test_arithmetic( multiply_op=.FALSE.,               &
1612     op1_ti_dd=0,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
1613     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
1614     res_int=0,                                             &
1615     testname='DivideINT_TI_TI3' )
1616   CALL test_arithmetic( multiply_op=.FALSE.,               &
1617     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
1618     op2_ti_dd=0,  op2_ti_h=01,  op2_ti_m=00,  op2_ti_s=00, &
1619     res_int=24,                                            &
1620     testname='DivideINT_TI_TI4' )
1621   CALL test_arithmetic( multiply_op=.FALSE.,               &
1622     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
1623     op2_ti_dd=0,  op2_ti_h=00,  op2_ti_m=01,  op2_ti_s=00, &
1624     res_int=1440,                                          &
1625     testname='DivideINT_TI_TI5' )
1626   CALL test_arithmetic( multiply_op=.FALSE.,               &
1627     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
1628     op2_ti_dd=0,  op2_ti_h=00,  op2_ti_m=00,  op2_ti_s=01, &
1629     res_int=86400,                                         &
1630     testname='DivideINT_TI_TI6' )
1631   ! rounding
1632   CALL test_arithmetic( multiply_op=.FALSE.,               &
1633     op1_ti_dd=0,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=03, &
1634     op2_ti_dd=0,  op2_ti_h=00,  op2_ti_m=00,  op2_ti_s=02, &
1635     res_int=1,                                             &
1636     testname='DivideINT_TI_TIR1' )
1637   CALL test_arithmetic( multiply_op=.FALSE.,               &
1638     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=02, &
1639     op2_ti_dd=1,  op2_ti_h=00,  op2_ti_m=00,  op2_ti_s=03, &
1640     res_int=0,                                             &
1641     testname='DivideINT_TI_TIR2' )
1642   ! fractional operands
1643   CALL test_arithmetic( multiply_op=.FALSE.,               &
1644     op1_ti_m=00,  op1_ti_s=00, op1_ti_sn=03, op1_ti_sd=04, &
1645     op2_ti_m=00,  op2_ti_s=00, op2_ti_sn=03, op2_ti_sd=04, &
1646     res_int=1,                                             &
1647     testname='DivideINT_TI_TIF1' )
1648   CALL test_arithmetic( multiply_op=.FALSE.,               &
1649     op1_ti_m=00,  op1_ti_s=00, op1_ti_sn=06, op1_ti_sd=08, &
1650     op2_ti_m=00,  op2_ti_s=00, op2_ti_sn=03, op2_ti_sd=04, &
1651     res_int=1,                                             &
1652     testname='DivideINT_TI_TIF2' )
1653   CALL test_arithmetic( multiply_op=.FALSE.,               &
1654     op1_ti_m=00,  op1_ti_s=00, op1_ti_sn=03, op1_ti_sd=04, &
1655     op2_ti_m=00,  op2_ti_s=00, op2_ti_sn=04, op2_ti_sd=03, &
1656     res_int=0,                                             &
1657     testname='DivideINT_TI_TIF3' )
1658   CALL test_arithmetic( multiply_op=.FALSE.,               &
1659     op1_ti_m=00,  op1_ti_s=02, op1_ti_sn=03, op1_ti_sd=04, &
1660     op2_ti_m=00,  op2_ti_s=01, op2_ti_sn=01, op2_ti_sd=03, &
1661     res_int=2,                                             &
1662     testname='DivideINT_TI_TIF4' )
1663   ! negative operands
1664   CALL test_arithmetic( multiply_op=.FALSE.,               &
1665     op1_ti_dd=-6,  op1_ti_h=-24,  op1_ti_m=-36,  op1_ti_s=-66, &
1666     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
1667     res_int=-2,                                             &
1668     testname='DivideINT_TI_TIN1' )
1669   CALL test_arithmetic( multiply_op=.FALSE.,               &
1670     op1_ti_dd=6,  op1_ti_h=24,  op1_ti_m=36,  op1_ti_s=66, &
1671     op2_ti_dd=-3,  op2_ti_h=-12,  op2_ti_m=-18,  op2_ti_s=-33, &
1672     res_int=-2,                                             &
1673     testname='DivideINT_TI_TIN2' )
1674   CALL test_arithmetic( multiply_op=.FALSE.,               &
1675     op1_ti_dd=-6,  op1_ti_h=-24,  op1_ti_m=-36,  op1_ti_s=-66, &
1676     op2_ti_dd=-3,  op2_ti_h=-12,  op2_ti_m=-18,  op2_ti_s=-33, &
1677     res_int=2,                                             &
1678     testname='DivideINT_TI_TIN3' )
1679
1680 !$$$here...  modify these to add self-test PASS/FAIL output
1681   CALL test_clock_advance(                                                    &
1682     start_yy=2002, start_mm=12, start_dd=27, start_h=3, start_m=0, start_s=0, &
1683      stop_yy=2002,  stop_mm=12,  stop_dd=28,  stop_h=8,  stop_m=0,  stop_s=0, &
1684     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=600,                 &
1685     testname="SimpleClockAdvance" )
1686
1687   CALL test_clock_advance(                                                    &
1688     start_yy=2003, start_mm=12, start_dd=29, start_h=9, start_m=0, start_s=0, &
1689      stop_yy=2004,  stop_mm=1,   stop_dd=2,   stop_h=9,  stop_m=0,  stop_s=0, &
1690     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=3600,                &
1691     testname="StdYearClockAdvance", increment_S=10 )
1692
1693   CALL test_clock_advance(                                                    &
1694     start_yy=2004, start_mm=12, start_dd=29, start_h=9, start_m=0, start_s=0, &
1695      stop_yy=2005,  stop_mm=1,   stop_dd=2,   stop_h=9,  stop_m=0,  stop_s=0, &
1696     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=3600,                &
1697     testname="LeapYearClockAdvance", increment_S=10 )
1698
1699   ! NRCM domain 3 case:  120 seconds / 9
1700   ! 18 timesteps through end of leap year
1701   CALL test_clock_advance(                                                    &
1702     start_yy=2004, start_mm=12, start_dd=31, start_h=23, start_m=58, start_s=0,&
1703      stop_yy=2005,  stop_mm=1,   stop_dd=1,   stop_h=0,  stop_m=2,  stop_s=0, &
1704     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=13,                  &
1705     timestep_sn=1, timestep_sd=3,                                             &
1706     testname="LeapYearFractionClockAdvance",                                  &
1707     increment_S=1, increment_Sn=1, increment_Sd=3 )
1708
1709   CALL ESMF_Finalize( rc=rc )
1710   CALL test_check_error( ESMF_SUCCESS, rc, &
1711                         'ESMF_Finalize() ', &
1712                         __FILE__ , &
1713                         __LINE__  )
1714
1715   PRINT *,'END TEST SUITE'
1716
1717 END PROGRAM time_manager_test
1718