PageRenderTime 40ms CodeModel.GetById 20ms RepoModel.GetById 0ms app.codeStats 1ms

/wrfv2_fire/external/esmf_time_f90/Test1.F90

http://github.com/jbeezley/wrf-fire
FORTRAN Modern | 1718 lines | 1448 code | 91 blank | 179 comment | 11 complexity | ddb4b3d0d2d3b7c4ee321177a45049fb MD5 | raw file
Possible License(s): AGPL-1.0 
    

  1. !
    2. 

  2. ! Sub-system tests for esmf_time_f90
    3. 

  3. !
    4. 

  4. ! Someday, switch over to funit!  
    5. 

  5. !
    6. 

  6. 

  7. MODULE my_tests
    8. 

  8.   USE ESMF_Mod
    9. 

  9.   IMPLICIT NONE
    10. 

  10. 

  11.   ! Set this to .TRUE. to make wrf_error_fatal3() print a message on failure 
    12. 

  12.   ! instead of stopping the program.  Use for testing only (since we cannot 
    13. 

  13.   ! catch exceptions in Fortran90!!)  
    14. 

  14.   LOGICAL :: WRF_ERROR_FATAL_PRINT = .FALSE.
    15. 

  15. 

  16. CONTAINS
    17. 

  17. 

  18.   ! Test printing of an ESMF_Time or ESMF_TimeInterval object.  
    19. 

  19.   !
    20. 

  20.   ! Correct results are also passed in through this interface and compared 
    21. 

  21.   ! with computed results.  PASS/FAIL messages are printed.  
    22. 

  22.   !
    23. 

  23.   SUBROUTINE test_print(  t_yy,  t_mm,  t_dd,  t_h,  t_m,  t_s, t_sn, t_sd, &
    24. 

  24.                          ti_yy, ti_mm, ti_dd, ti_h, ti_m, ti_s, ti_sn, ti_sd, &
    25. 

  25.                          res_str, testname, expect_error )
    26. 

  26.     INTEGER, INTENT(IN), OPTIONAL :: t_YY
    27. 

  27.     INTEGER, INTENT(IN), OPTIONAL :: t_MM  ! month
    28. 

  28.     INTEGER, INTENT(IN), OPTIONAL :: t_DD  ! day of month
    29. 

  29.     INTEGER, INTENT(IN), OPTIONAL :: t_H
    30. 

  30.     INTEGER, INTENT(IN), OPTIONAL :: t_M
    31. 

  31.     INTEGER, INTENT(IN), OPTIONAL :: t_S
    32. 

  32.     INTEGER, INTENT(IN), OPTIONAL :: t_Sn
    33. 

  33.     INTEGER, INTENT(IN), OPTIONAL :: t_Sd
    34. 

  34.     INTEGER, INTENT(IN), OPTIONAL :: ti_YY
    35. 

  35.     INTEGER, INTENT(IN), OPTIONAL :: ti_MM  ! month
    36. 

  36.     INTEGER, INTENT(IN), OPTIONAL :: ti_DD  ! day of month
    37. 

  37.     INTEGER, INTENT(IN), OPTIONAL :: ti_H
    38. 

  38.     INTEGER, INTENT(IN), OPTIONAL :: ti_M
    39. 

  39.     INTEGER, INTENT(IN), OPTIONAL :: ti_S
    40. 

  40.     INTEGER, INTENT(IN), OPTIONAL :: ti_Sn
    41. 

  41.     INTEGER, INTENT(IN), OPTIONAL :: ti_Sd
    42. 

  42.     CHARACTER (LEN=*), INTENT(IN) :: res_str
    43. 

  43.     CHARACTER (LEN=*), INTENT(IN), OPTIONAL :: testname
    44. 

  44.     LOGICAL, OPTIONAL, INTENT(IN) :: expect_error
    45. 

  45.     ! locals
    46. 

  46.     INTEGER :: it_YY
    47. 

  47.     INTEGER :: it_MM  ! month
    48. 

  48.     INTEGER :: it_DD  ! day of month
    49. 

  49.     INTEGER :: it_H
    50. 

  50.     INTEGER :: it_M
    51. 

  51.     INTEGER :: it_S
    52. 

  52.     INTEGER :: it_Sn
    53. 

  53.     INTEGER :: it_Sd
    54. 

  54.     INTEGER :: iti_YY
    55. 

  55.     INTEGER :: iti_MM  ! month
    56. 

  56.     INTEGER :: iti_DD  ! day of month
    57. 

  57.     INTEGER :: iti_H
    58. 

  58.     INTEGER :: iti_M
    59. 

  59.     INTEGER :: iti_S
    60. 

  60.     INTEGER :: iti_Sn
    61. 

  61.     INTEGER :: iti_Sd
    62. 

  62.     LOGICAL :: is_t 
    63. 

  63.     LOGICAL :: is_ti
    64. 

  64.     CHARACTER (LEN=512) :: itestname
    65. 

  65.     LOGICAL :: iexpect_error
    66. 

  66.     INTEGER rc
    67. 

  67.     TYPE(ESMF_Time)           :: t
    68. 

  68.     TYPE(ESMF_TimeInterval)   :: ti
    69. 

  69.     CHARACTER(LEN=ESMF_MAXSTR) :: str, computed_str, frac_str
    70. 

  70.     CHARACTER(LEN=17) :: type_str
    71. 

  71.     INTEGER :: res_len, computed_len, Sn, Sd
    72. 

  72.     LOGICAL :: test_passed
    73. 

  73. 

  74. !  PRINT *,'DEBUG:  BEGIN test_print()'
    75. 

  75.     it_YY = 0
    76. 

  76.     it_MM = 1
    77. 

  77.     it_DD = 1
    78. 

  78.     it_H = 0
    79. 

  79.     it_M = 0
    80. 

  80.     it_S = 0
    81. 

  81.     it_Sn = 0
    82. 

  82.     it_Sd = 0
    83. 

  83.     iti_YY = 0
    84. 

  84.     iti_MM = 0
    85. 

  85.     iti_DD = 0
    86. 

  86.     iti_H = 0
    87. 

  87.     iti_M = 0
    88. 

  88.     iti_S = 0
    89. 

  89.     iti_Sn = 0
    90. 

  90.     iti_Sd = 0
    91. 

  91.     itestname = ''
    92. 

  92.     iexpect_error = .FALSE.
    93. 

  93. 

  94.     IF ( PRESENT( t_YY ) ) it_YY = t_YY
    95. 

  95.     IF ( PRESENT( t_MM ) ) it_MM = t_MM
    96. 

  96.     IF ( PRESENT( t_DD ) ) it_DD = t_DD
    97. 

  97.     IF ( PRESENT( t_H ) ) it_H = t_H
    98. 

  98.     IF ( PRESENT( t_M ) ) it_M = t_M
    99. 

  99.     IF ( PRESENT( t_S ) ) it_S = t_S
    100. 

  100.     IF ( PRESENT( t_Sn ) ) it_Sn = t_Sn
    101. 

  101.     IF ( PRESENT( t_Sd ) ) it_Sd = t_Sd
    102. 

  102.     IF ( PRESENT( ti_YY ) ) iti_YY = ti_YY
    103. 

  103.     IF ( PRESENT( ti_MM ) ) iti_MM = ti_MM
    104. 

  104.     IF ( PRESENT( ti_DD ) ) iti_DD = ti_DD
    105. 

  105.     IF ( PRESENT( ti_H ) ) iti_H = ti_H
    106. 

  106.     IF ( PRESENT( ti_M ) ) iti_M = ti_M
    107. 

  107.     IF ( PRESENT( ti_S ) ) iti_S = ti_S
    108. 

  108.     IF ( PRESENT( ti_Sn ) ) iti_Sn = ti_Sn
    109. 

  109.     IF ( PRESENT( ti_Sd ) ) iti_Sd = ti_Sd
    110. 

  110.     IF ( PRESENT( testname ) ) itestname = TRIM(testname)
    111. 

  111.     IF ( PRESENT( expect_error ) ) iexpect_error = expect_error
    112. 

  112. 

  113.     ! Ensure that optional arguments are consistent...
    114. 

  114.     is_t = ( PRESENT( t_YY ) .OR. PRESENT( t_MM ) .OR. &
    115. 

  115.              PRESENT( t_DD ) .OR. PRESENT( t_H ) .OR.  &
    116. 

  116.              PRESENT( t_M )  .OR. PRESENT( t_S ) .OR.  &
    117. 

  117.              PRESENT( t_Sn )  .OR. PRESENT( t_Sd ) )
    118. 

  118.     is_ti = ( PRESENT( ti_YY ) .OR. PRESENT( ti_MM ) .OR. &
    119. 

  119.               PRESENT( ti_DD ) .OR. PRESENT( ti_H ) .OR.  &
    120. 

  120.               PRESENT( ti_M )  .OR. PRESENT( ti_S ) .OR.  &
    121. 

  121.               PRESENT( ti_Sn )  .OR. PRESENT( ti_Sd ) )
    122. 

  122.     IF ( is_t .EQV. is_ti ) THEN
    123. 

  123.       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
    124. 

  124.         'ERROR test_print:  inconsistent args' )
    125. 

  125.     ENDIF
    126. 

  126. 

  127. !PRINT *,'DEBUG:  test_print():  init objects'
    128. 

  128.     ! Initialize object to be tested
    129. 

  129.     ! modify behavior of wrf_error_fatal3 for tests expected to fail
    130. 

  130.     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .TRUE.
    131. 

  131.     Sn = 0
    132. 

  132.     Sd = 0
    133. 

  133.     IF ( is_t ) THEN
    134. 

  134.       type_str = 'ESMF_Time'
    135. 

  135. !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeSet()'
    136. 

  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. 

  137.       CALL ESMF_TimeSet( t, YY=it_YY, MM=it_MM, DD=it_DD , &
    138. 

  138.                              H=it_H, M=it_M, S=it_S, Sn=it_Sn, Sd=it_Sd, rc=rc )
    139. 

  139. !PRINT *,'DEBUG:  test_print():  back from ESMF_TimeSet()'
    140. 

  140.       CALL test_check_error( ESMF_SUCCESS, rc, &
    141. 

  141.                              TRIM(itestname)//'ESMF_TimeSet() ', &
    142. 

  142.                              __FILE__ , &
    143. 

  143.                              __LINE__  )
    144. 

  144. !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeGet()'
    145. 

  145.       CALL ESMF_TimeGet( t, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
    146. 

  146.       CALL test_check_error( ESMF_SUCCESS, rc, &
    147. 

  147.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    148. 

  148.                             __FILE__ , &
    149. 

  149.                             __LINE__  )
    150. 

  150. !PRINT *,'DEBUG:  test_print():  back from ESMF_TimeGet(), computed_str = ',TRIM(computed_str)
    151. 

  151.     ELSE
    152. 

  152.       type_str = 'ESMF_TimeInterval'
    153. 

  153. !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeIntervalSet()'
    154. 

  154.       CALL ESMF_TimeIntervalSet( ti, YY=iti_YY, MM=iti_MM, &
    155. 

  155.                                       D=iti_DD ,           &
    156. 

  156.                                       H=iti_H, M=iti_M,    &
    157. 

  157.                                       S=iti_S, Sn=iti_Sn, Sd=iti_Sd, rc=rc )
    158. 

  158.       CALL test_check_error( ESMF_SUCCESS, rc, &
    159. 

  159.                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
    160. 

  160.                              __FILE__ , &
    161. 

  161.                              __LINE__  )
    162. 

  162. !PRINT *,'DEBUG:  test_print():  calling ESMF_TimeIntervalGet()'
    163. 

  163.       CALL ESMF_TimeIntervalGet( ti, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
    164. 

  164.       CALL test_check_error( ESMF_SUCCESS, rc, &
    165. 

  165.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    166. 

  166.                             __FILE__ , &
    167. 

  167.                             __LINE__  )
    168. 

  168.     ENDIF
    169. 

  169.     ! handle fractions
    170. 

  170.     IF ( Sd > 0 ) THEN
    171. 

  171.       IF ( Sn > 0 ) THEN
    172. 

  172.         WRITE(frac_str,FMT="('+',I2.2,'/',I2.2)") abs(Sn), Sd
    173. 

  173.       ELSE IF ( Sn < 0 ) THEN
    174. 

  174.         WRITE(frac_str,FMT="('-',I2.2,'/',I2.2)") abs(Sn), Sd
    175. 

  175.       ELSE
    176. 

  176.         frac_str = ''
    177. 

  177.       ENDIF
    178. 

  178.       computed_str = TRIM(computed_str)//TRIM(frac_str)
    179. 

  179.     ENDIF
    180. 

  180.     ! restore default behavior of wrf_error_fatal3
    181. 

  181.     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .FALSE.
    182. 

  182. !PRINT *,'DEBUG:  test_print():  done init objects'
    183. 

  183. 

  184. !PRINT *,'DEBUG:  test_print():  check result'
    185. 

  185.     ! check result
    186. 

  186.     test_passed = .FALSE.
    187. 

  187.     res_len = LEN_TRIM(res_str)
    188. 

  188.     computed_len = LEN_TRIM(computed_str)
    189. 

  189.     IF ( res_len == computed_len ) THEN
    190. 

  190.       IF ( computed_str(1:computed_len) == res_str(1:res_len) ) THEN
    191. 

  191.         test_passed = .TRUE.
    192. 

  192.       ENDIF
    193. 

  193.     ENDIF
    194. 

  194.     IF ( test_passed ) THEN
    195. 

  195.       WRITE(*,FMT='(A)') 'PASS:  '//TRIM(itestname)
    196. 

  196.     ELSE
    197. 

  197.       WRITE(*,'(9A)') 'FAIL:  ',TRIM(itestname),':  printing ',TRIM(type_str), &
    198. 

  198.         '  expected <', TRIM(res_str),'>  but computed <',TRIM(computed_str),'>'
    199. 

  199.     ENDIF
    200. 

  200. !PRINT *,'DEBUG:  END test_print()'
    201. 

  201. 

  202.   END SUBROUTINE test_print
    203. 

  203. 

  204. 

  205. 

  206.   ! Test the following arithmetic operations on ESMF_Time and 
    207. 

  207.   ! ESMF_TimeInterval objects:
    208. 

  208.   !  ESMF_Time         = ESMF_Time         + ESMF_TimeInterval
    209. 

  209.   !  ESMF_Time         = ESMF_TimeInterval + ESMF_Time
    210. 

  210.   !  ESMF_Time         = ESMF_Time         - ESMF_TimeInterval
    211. 

  211.   !  ESMF_TimeInterval = ESMF_Time         - ESMF_Time        
    212. 

  212.   !  ESMF_TimeInterval = ESMF_TimeInterval + ESMF_TimeInterval
    213. 

  213.   !  ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
    214. 

  214.   !  ESMF_TimeInterval = ESMF_TimeInterval * INTEGER
    215. 

  215.   !  ESMF_TimeInterval = ESMF_TimeInterval / INTEGER
    216. 

  216.   !
    217. 

  217.   ! Correct results are also passed in through this interface and compared 
    218. 

  218.   ! with computed results.  PASS/FAIL messages are printed.  
    219. 

  219.   !
    220. 

  220.   ! Operations are expressed as res = op1 +|- op2
    221. 

  221.   !
    222. 

  222.   SUBROUTINE test_arithmetic( add_op, multiply_op,                                       &
    223. 

  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. 

  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. 

  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. 

  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. 

  227.     op2_int,                                                                             &
    228. 

  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. 

  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. 

  230.     res_int, testname, expect_error )
    231. 

  231.     LOGICAL, INTENT(IN), OPTIONAL :: add_op      ! .TRUE.=add, .FALSE.=subtract
    232. 

  232.     LOGICAL, INTENT(IN), OPTIONAL :: multiply_op ! .TRUE.=multiply, .FALSE.=divide
    233. 

  233.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_YY
    234. 

  234.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_MM  ! month
    235. 

  235.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_DD  ! day of month
    236. 

  236.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_H
    237. 

  237.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_M
    238. 

  238.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_S
    239. 

  239.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_Sn
    240. 

  240.     INTEGER, INTENT(IN), OPTIONAL :: op1_t_Sd
    241. 

  241.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_YY
    242. 

  242.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_MM  ! month
    243. 

  243.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_DD  ! day of month
    244. 

  244.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_H
    245. 

  245.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_M
    246. 

  246.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_S
    247. 

  247.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_Sn
    248. 

  248.     INTEGER, INTENT(IN), OPTIONAL :: op1_ti_Sd
    249. 

  249.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_YY
    250. 

  250.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_MM  ! month
    251. 

  251.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_DD  ! day of month
    252. 

  252.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_H
    253. 

  253.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_M
    254. 

  254.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_S
    255. 

  255.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_Sn
    256. 

  256.     INTEGER, INTENT(IN), OPTIONAL :: op2_t_Sd
    257. 

  257.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_YY
    258. 

  258.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_MM  ! month
    259. 

  259.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_DD  ! day of month
    260. 

  260.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_H
    261. 

  261.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_M
    262. 

  262.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_S
    263. 

  263.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_Sn
    264. 

  264.     INTEGER, INTENT(IN), OPTIONAL :: op2_ti_Sd
    265. 

  265.     INTEGER, INTENT(IN), OPTIONAL :: op2_int
    266. 

  266.     INTEGER, INTENT(IN), OPTIONAL :: res_t_YY
    267. 

  267.     INTEGER, INTENT(IN), OPTIONAL :: res_t_MM  ! month
    268. 

  268.     INTEGER, INTENT(IN), OPTIONAL :: res_t_DD  ! day of month
    269. 

  269.     INTEGER, INTENT(IN), OPTIONAL :: res_t_H
    270. 

  270.     INTEGER, INTENT(IN), OPTIONAL :: res_t_M
    271. 

  271.     INTEGER, INTENT(IN), OPTIONAL :: res_t_S
    272. 

  272.     INTEGER, INTENT(IN), OPTIONAL :: res_t_Sn
    273. 

  273.     INTEGER, INTENT(IN), OPTIONAL :: res_t_Sd
    274. 

  274.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_YY
    275. 

  275.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_MM  ! month
    276. 

  276.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_DD  ! day of month
    277. 

  277.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_H
    278. 

  278.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_M
    279. 

  279.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_S
    280. 

  280.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_Sn
    281. 

  281.     INTEGER, INTENT(IN), OPTIONAL :: res_ti_Sd
    282. 

  282.     INTEGER, INTENT(IN), OPTIONAL :: res_int
    283. 

  283.     CHARACTER (LEN=*), OPTIONAL, INTENT(IN) :: testname
    284. 

  284.     LOGICAL, OPTIONAL, INTENT(IN) :: expect_error
    285. 

  285.     ! locals
    286. 

  286.     LOGICAL :: iadd_op
    287. 

  287.     LOGICAL :: isubtract_op
    288. 

  288.     LOGICAL :: imultiply_op
    289. 

  289.     LOGICAL :: idivide_op
    290. 

  290.     INTEGER :: iop1_t_YY
    291. 

  291.     INTEGER :: iop1_t_MM  ! month
    292. 

  292.     INTEGER :: iop1_t_DD  ! day of month
    293. 

  293.     INTEGER :: iop1_t_H
    294. 

  294.     INTEGER :: iop1_t_M
    295. 

  295.     INTEGER :: iop1_t_S
    296. 

  296.     INTEGER :: iop1_t_Sn
    297. 

  297.     INTEGER :: iop1_t_Sd
    298. 

  298.     INTEGER :: iop1_ti_YY
    299. 

  299.     INTEGER :: iop1_ti_MM  ! month
    300. 

  300.     INTEGER :: iop1_ti_DD  ! day of month
    301. 

  301.     INTEGER :: iop1_ti_H
    302. 

  302.     INTEGER :: iop1_ti_M
    303. 

  303.     INTEGER :: iop1_ti_S
    304. 

  304.     INTEGER :: iop1_ti_Sn
    305. 

  305.     INTEGER :: iop1_ti_Sd
    306. 

  306.     INTEGER :: iop2_t_YY
    307. 

  307.     INTEGER :: iop2_t_MM  ! month
    308. 

  308.     INTEGER :: iop2_t_DD  ! day of month
    309. 

  309.     INTEGER :: iop2_t_H
    310. 

  310.     INTEGER :: iop2_t_M
    311. 

  311.     INTEGER :: iop2_t_S
    312. 

  312.     INTEGER :: iop2_t_Sn
    313. 

  313.     INTEGER :: iop2_t_Sd
    314. 

  314.     INTEGER :: iop2_ti_YY
    315. 

  315.     INTEGER :: iop2_ti_MM  ! month
    316. 

  316.     INTEGER :: iop2_ti_DD  ! day of month
    317. 

  317.     INTEGER :: iop2_ti_H
    318. 

  318.     INTEGER :: iop2_ti_M
    319. 

  319.     INTEGER :: iop2_ti_S
    320. 

  320.     INTEGER :: iop2_ti_Sn
    321. 

  321.     INTEGER :: iop2_ti_Sd
    322. 

  322.     INTEGER :: ires_t_YY
    323. 

  323.     INTEGER :: ires_t_MM  ! month
    324. 

  324.     INTEGER :: ires_t_DD  ! day of month
    325. 

  325.     INTEGER :: ires_t_H
    326. 

  326.     INTEGER :: ires_t_M
    327. 

  327.     INTEGER :: ires_t_S
    328. 

  328.     INTEGER :: ires_t_Sn
    329. 

  329.     INTEGER :: ires_t_Sd
    330. 

  330.     INTEGER :: ires_ti_YY
    331. 

  331.     INTEGER :: ires_ti_MM  ! month
    332. 

  332.     INTEGER :: ires_ti_DD  ! day of month
    333. 

  333.     INTEGER :: ires_ti_H
    334. 

  334.     INTEGER :: ires_ti_M
    335. 

  335.     INTEGER :: ires_ti_S
    336. 

  336.     INTEGER :: ires_ti_Sn
    337. 

  337.     INTEGER :: ires_ti_Sd
    338. 

  338.     LOGICAL :: op1_is_t , op2_is_t , res_is_t
    339. 

  339.     LOGICAL :: op1_is_ti, op2_is_ti, res_is_ti, op2_is_int
    340. 

  340.     LOGICAL :: res_is_int
    341. 

  341.     INTEGER :: num_ops, num_op1, num_op2, num_res
    342. 

  342.     LOGICAL :: unsupported_op, test_passed
    343. 

  343.     CHARACTER (LEN=512) :: itestname
    344. 

  344.     LOGICAL :: iexpect_error
    345. 

  345.     INTEGER :: rc
    346. 

  346.     INTEGER :: computed_int, Sn, Sd
    347. 

  347.     TYPE(ESMF_Time)           :: op1_t , op2_t , res_t, computed_t
    348. 

  348.     TYPE(ESMF_TimeInterval)   :: op1_ti, op2_ti, res_ti, computed_ti
    349. 

  349.     CHARACTER(LEN=ESMF_MAXSTR) :: str, op1_str, op2_str, res_str, computed_str, frac_str
    350. 

  350.     CHARACTER(LEN=1) :: op_str
    351. 

  351.     CHARACTER(LEN=17) :: op1_type_str, op2_type_str, res_type_str
    352. 

  352. 

  353.     iadd_op = .FALSE.
    354. 

  354.     isubtract_op = .FALSE.
    355. 

  355.     imultiply_op = .FALSE.
    356. 

  356.     idivide_op = .FALSE.
    357. 

  357.     iop1_t_YY = 0
    358. 

  358.     iop1_t_MM = 1
    359. 

  359.     iop1_t_DD = 1
    360. 

  360.     iop1_t_H = 0
    361. 

  361.     iop1_t_M = 0
    362. 

  362.     iop1_t_S = 0
    363. 

  363.     iop1_t_Sn = 0
    364. 

  364.     iop1_t_Sd = 0
    365. 

  365.     iop1_ti_YY = 0
    366. 

  366.     iop1_ti_MM = 0
    367. 

  367.     iop1_ti_DD = 0
    368. 

  368.     iop1_ti_H = 0
    369. 

  369.     iop1_ti_M = 0
    370. 

  370.     iop1_ti_S = 0
    371. 

  371.     iop1_ti_Sn = 0
    372. 

  372.     iop1_ti_Sd = 0
    373. 

  373.     iop2_t_YY = 0
    374. 

  374.     iop2_t_MM = 1
    375. 

  375.     iop2_t_DD = 1
    376. 

  376.     iop2_t_H = 0
    377. 

  377.     iop2_t_M = 0
    378. 

  378.     iop2_t_S = 0
    379. 

  379.     iop2_t_Sn = 0
    380. 

  380.     iop2_t_Sd = 0
    381. 

  381.     iop2_ti_YY = 0
    382. 

  382.     iop2_ti_MM = 0
    383. 

  383.     iop2_ti_DD = 0
    384. 

  384.     iop2_ti_H = 0
    385. 

  385.     iop2_ti_M = 0
    386. 

  386.     iop2_ti_S = 0
    387. 

  387.     iop2_ti_Sn = 0
    388. 

  388.     iop2_ti_Sd = 0
    389. 

  389.     ires_t_YY = 0
    390. 

  390.     ires_t_MM = 1
    391. 

  391.     ires_t_DD = 1
    392. 

  392.     ires_t_H = 0
    393. 

  393.     ires_t_M = 0
    394. 

  394.     ires_t_S = 0
    395. 

  395.     ires_t_Sn = 0
    396. 

  396.     ires_t_Sd = 0
    397. 

  397.     ires_ti_YY = 0
    398. 

  398.     ires_ti_MM = 0
    399. 

  399.     ires_ti_DD = 0
    400. 

  400.     ires_ti_H = 0
    401. 

  401.     ires_ti_M = 0
    402. 

  402.     ires_ti_S = 0
    403. 

  403.     ires_ti_Sn = 0
    404. 

  404.     ires_ti_Sd = 0
    405. 

  405.     itestname = ''
    406. 

  406.     iexpect_error = .FALSE.
    407. 

  407. 

  408.     IF ( PRESENT( add_op ) ) THEN
    409. 

  409.       iadd_op = add_op
    410. 

  410.       isubtract_op = ( .NOT. add_op )
    411. 

  411.     ENDIF
    412. 

  412.     IF ( PRESENT( multiply_op ) ) THEN
    413. 

  413.       imultiply_op = multiply_op
    414. 

  414.       idivide_op = ( .NOT. multiply_op )
    415. 

  415.     ENDIF
    416. 

  416.     num_ops = 0
    417. 

  417.     IF ( iadd_op )      num_ops = num_ops + 1
    418. 

  418.     IF ( isubtract_op ) num_ops = num_ops + 1
    419. 

  419.     IF ( imultiply_op ) num_ops = num_ops + 1
    420. 

  420.     IF ( idivide_op )   num_ops = num_ops + 1
    421. 

  421.     IF ( num_ops /= 1 ) THEN
    422. 

  422.       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
    423. 

  423.         'ERROR test_arithmetic:  inconsistent operation' )
    424. 

  424.     ENDIF
    425. 

  425.     IF ( PRESENT( op1_t_YY ) ) iop1_t_YY = op1_t_YY
    426. 

  426.     IF ( PRESENT( op1_t_MM ) ) iop1_t_MM = op1_t_MM
    427. 

  427.     IF ( PRESENT( op1_t_DD ) ) iop1_t_DD = op1_t_DD
    428. 

  428.     IF ( PRESENT( op1_t_H ) ) iop1_t_H = op1_t_H
    429. 

  429.     IF ( PRESENT( op1_t_M ) ) iop1_t_M = op1_t_M
    430. 

  430.     IF ( PRESENT( op1_t_S ) ) iop1_t_S = op1_t_S
    431. 

  431.     IF ( PRESENT( op1_t_Sn ) ) iop1_t_Sn = op1_t_Sn
    432. 

  432.     IF ( PRESENT( op1_t_Sd ) ) iop1_t_Sd = op1_t_Sd
    433. 

  433.     IF ( PRESENT( op1_ti_YY ) ) iop1_ti_YY = op1_ti_YY
    434. 

  434.     IF ( PRESENT( op1_ti_MM ) ) iop1_ti_MM = op1_ti_MM
    435. 

  435.     IF ( PRESENT( op1_ti_DD ) ) iop1_ti_DD = op1_ti_DD
    436. 

  436.     IF ( PRESENT( op1_ti_H ) ) iop1_ti_H = op1_ti_H
    437. 

  437.     IF ( PRESENT( op1_ti_M ) ) iop1_ti_M = op1_ti_M
    438. 

  438.     IF ( PRESENT( op1_ti_S ) ) iop1_ti_S = op1_ti_S
    439. 

  439.     IF ( PRESENT( op1_ti_Sn ) ) iop1_ti_Sn = op1_ti_Sn
    440. 

  440.     IF ( PRESENT( op1_ti_Sd ) ) iop1_ti_Sd = op1_ti_Sd
    441. 

  441.     IF ( PRESENT( op2_t_YY ) ) iop2_t_YY = op2_t_YY
    442. 

  442.     IF ( PRESENT( op2_t_MM ) ) iop2_t_MM = op2_t_MM
    443. 

  443.     IF ( PRESENT( op2_t_DD ) ) iop2_t_DD = op2_t_DD
    444. 

  444.     IF ( PRESENT( op2_t_H ) ) iop2_t_H = op2_t_H
    445. 

  445.     IF ( PRESENT( op2_t_M ) ) iop2_t_M = op2_t_M
    446. 

  446.     IF ( PRESENT( op2_t_S ) ) iop2_t_S = op2_t_S
    447. 

  447.     IF ( PRESENT( op2_t_Sn ) ) iop2_t_Sn = op2_t_Sn
    448. 

  448.     IF ( PRESENT( op2_t_Sd ) ) iop2_t_Sd = op2_t_Sd
    449. 

  449.     IF ( PRESENT( op2_ti_YY ) ) iop2_ti_YY = op2_ti_YY
    450. 

  450.     IF ( PRESENT( op2_ti_MM ) ) iop2_ti_MM = op2_ti_MM
    451. 

  451.     IF ( PRESENT( op2_ti_DD ) ) iop2_ti_DD = op2_ti_DD
    452. 

  452.     IF ( PRESENT( op2_ti_H ) ) iop2_ti_H = op2_ti_H
    453. 

  453.     IF ( PRESENT( op2_ti_M ) ) iop2_ti_M = op2_ti_M
    454. 

  454.     IF ( PRESENT( op2_ti_S ) ) iop2_ti_S = op2_ti_S
    455. 

  455.     IF ( PRESENT( op2_ti_Sn ) ) iop2_ti_Sn = op2_ti_Sn
    456. 

  456.     IF ( PRESENT( op2_ti_Sd ) ) iop2_ti_Sd = op2_ti_Sd
    457. 

  457.     IF ( PRESENT( res_t_YY ) ) ires_t_YY = res_t_YY
    458. 

  458.     IF ( PRESENT( res_t_MM ) ) ires_t_MM = res_t_MM
    459. 

  459.     IF ( PRESENT( res_t_DD ) ) ires_t_DD = res_t_DD
    460. 

  460.     IF ( PRESENT( res_t_H ) ) ires_t_H = res_t_H
    461. 

  461.     IF ( PRESENT( res_t_M ) ) ires_t_M = res_t_M
    462. 

  462.     IF ( PRESENT( res_t_S ) ) ires_t_S = res_t_S
    463. 

  463.     IF ( PRESENT( res_t_Sn ) ) ires_t_Sn = res_t_Sn
    464. 

  464.     IF ( PRESENT( res_t_Sd ) ) ires_t_Sd = res_t_Sd
    465. 

  465.     IF ( PRESENT( res_ti_YY ) ) ires_ti_YY = res_ti_YY
    466. 

  466.     IF ( PRESENT( res_ti_MM ) ) ires_ti_MM = res_ti_MM
    467. 

  467.     IF ( PRESENT( res_ti_DD ) ) ires_ti_DD = res_ti_DD
    468. 

  468.     IF ( PRESENT( res_ti_H ) ) ires_ti_H = res_ti_H
    469. 

  469.     IF ( PRESENT( res_ti_M ) ) ires_ti_M = res_ti_M
    470. 

  470.     IF ( PRESENT( res_ti_S ) ) ires_ti_S = res_ti_S
    471. 

  471.     IF ( PRESENT( res_ti_Sn ) ) ires_ti_Sn = res_ti_Sn
    472. 

  472.     IF ( PRESENT( res_ti_Sd ) ) ires_ti_Sd = res_ti_Sd
    473. 

  473.     IF ( PRESENT( testname ) ) itestname = TRIM(testname)
    474. 

  474.     IF ( PRESENT( expect_error ) ) iexpect_error = expect_error
    475. 

  475. 

  476.     ! Ensure that optional arguments are consistent...
    477. 

  477.     op1_is_t = ( PRESENT( op1_t_YY ) .OR. PRESENT( op1_t_MM ) .OR. &
    478. 

  478.                  PRESENT( op1_t_DD ) .OR. PRESENT( op1_t_H ) .OR.  &
    479. 

  479.                  PRESENT( op1_t_M )  .OR. PRESENT( op1_t_S ) .OR.  &
    480. 

  480.                  PRESENT( op1_t_Sn )  .OR. PRESENT( op1_t_Sd ) )
    481. 

  481.     op1_is_ti = ( PRESENT( op1_ti_YY ) .OR. PRESENT( op1_ti_MM ) .OR. &
    482. 

  482.                   PRESENT( op1_ti_DD ) .OR. PRESENT( op1_ti_H ) .OR.  &
    483. 

  483.                   PRESENT( op1_ti_M )  .OR. PRESENT( op1_ti_S ) .OR.  &
    484. 

  484.                   PRESENT( op1_ti_Sn )  .OR. PRESENT( op1_ti_Sd ) )
    485. 

  485.     op2_is_t = ( PRESENT( op2_t_YY ) .OR. PRESENT( op2_t_MM ) .OR. &
    486. 

  486.                  PRESENT( op2_t_DD ) .OR. PRESENT( op2_t_H ) .OR.  &
    487. 

  487.                  PRESENT( op2_t_M )  .OR. PRESENT( op2_t_S ) .OR.  &
    488. 

  488.                  PRESENT( op2_t_Sn )  .OR. PRESENT( op2_t_Sd ) )
    489. 

  489.     op2_is_ti = ( PRESENT( op2_ti_YY ) .OR. PRESENT( op2_ti_MM ) .OR. &
    490. 

  490.                   PRESENT( op2_ti_DD ) .OR. PRESENT( op2_ti_H ) .OR.  &
    491. 

  491.                   PRESENT( op2_ti_M )  .OR. PRESENT( op2_ti_S ) .OR.  &
    492. 

  492.                   PRESENT( op2_ti_Sn )  .OR. PRESENT( op2_ti_Sd ) )
    493. 

  493.     op2_is_int = ( PRESENT( op2_int ) )
    494. 

  494.     res_is_t = ( PRESENT( res_t_YY ) .OR. PRESENT( res_t_MM ) .OR. &
    495. 

  495.                  PRESENT( res_t_DD ) .OR. PRESENT( res_t_H ) .OR.  &
    496. 

  496.                  PRESENT( res_t_M )  .OR. PRESENT( res_t_S ) .OR.  &
    497. 

  497.                  PRESENT( res_t_Sn )  .OR. PRESENT( res_t_Sd ) )
    498. 

  498.     res_is_ti = ( PRESENT( res_ti_YY ) .OR. PRESENT( res_ti_MM ) .OR. &
    499. 

  499.                   PRESENT( res_ti_DD ) .OR. PRESENT( res_ti_H ) .OR.  &
    500. 

  500.                   PRESENT( res_ti_M )  .OR. PRESENT( res_ti_S ) .OR.  &
    501. 

  501.                   PRESENT( res_ti_Sn )  .OR. PRESENT( res_ti_Sd ) )
    502. 

  502.     res_is_int = ( PRESENT( res_int ) )
    503. 

  503.     num_op1 = 0
    504. 

  504.     IF ( op1_is_t   ) num_op1 = num_op1 + 1
    505. 

  505.     IF ( op1_is_ti  ) num_op1 = num_op1 + 1
    506. 

  506.     IF ( num_op1 /= 1 ) THEN
    507. 

  507.       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
    508. 

  508.         'ERROR test_arithmetic:  inconsistent args for op1' )
    509. 

  509.     ENDIF
    510. 

  510.     num_op2 = 0
    511. 

  511.     IF ( op2_is_t   ) num_op2 = num_op2 + 1
    512. 

  512.     IF ( op2_is_ti  ) num_op2 = num_op2 + 1
    513. 

  513.     IF ( op2_is_int ) num_op2 = num_op2 + 1
    514. 

  514.     IF ( num_op2 /= 1 ) THEN
    515. 

  515.       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
    516. 

  516.         'ERROR test_arithmetic:  inconsistent args for op2' )
    517. 

  517.     ENDIF
    518. 

  518.     num_res = 0
    519. 

  519.     IF ( res_is_t   ) num_res = num_res + 1
    520. 

  520.     IF ( res_is_ti  ) num_res = num_res + 1
    521. 

  521.     IF ( res_is_int ) num_res = num_res + 1
    522. 

  522.     IF ( num_res /= 1 ) THEN
    523. 

  523.       CALL wrf_error_fatal3( __FILE__ , __LINE__ , &
    524. 

  524.         'ERROR test_arithmetic:  inconsistent args for result' )
    525. 

  525.     ENDIF
    526. 

  526. 

  527.     ! Initialize op1
    528. 

  528.     IF ( op1_is_t ) THEN
    529. 

  529.       op1_type_str = 'ESMF_Time'
    530. 

  530.       CALL ESMF_TimeSet( op1_t, YY=iop1_t_YY, MM=iop1_t_MM, DD=iop1_t_DD , &
    531. 

  531.                                  H=iop1_t_H, M=iop1_t_M, S=iop1_t_S, Sn=iop1_t_Sn, Sd=iop1_t_Sd, rc=rc )
    532. 

  532.       CALL test_check_error( ESMF_SUCCESS, rc, &
    533. 

  533.                              TRIM(itestname)//'ESMF_TimeSet() ', &
    534. 

  534.                              __FILE__ , &
    535. 

  535.                              __LINE__  )
    536. 

  536.       CALL ESMF_TimeGet( op1_t, timeString=op1_str, Sn=Sn, Sd=Sd, rc=rc )
    537. 

  537.       CALL test_check_error( ESMF_SUCCESS, rc, &
    538. 

  538.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    539. 

  539.                             __FILE__ , &
    540. 

  540.                             __LINE__  )
    541. 

  541.       ! handle fractions
    542. 

  542.       CALL fraction_to_string( Sn, Sd, frac_str )
    543. 

  543.       op1_str = TRIM(op1_str)//TRIM(frac_str)
    544. 

  544.     ELSE
    545. 

  545.       op1_type_str = 'ESMF_TimeInterval'
    546. 

  546.       CALL ESMF_TimeIntervalSet( op1_ti, YY=iop1_ti_YY, MM=iop1_ti_MM, &
    547. 

  547.                                           D=iop1_ti_DD ,               &
    548. 

  548.                                           H=iop1_ti_H, M=iop1_ti_M,    &
    549. 

  549.                                           S=iop1_ti_S, Sn=iop1_ti_Sn, Sd=iop1_ti_Sd, rc=rc )
    550. 

  550.       CALL test_check_error( ESMF_SUCCESS, rc, &
    551. 

  551.                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
    552. 

  552.                              __FILE__ , &
    553. 

  553.                              __LINE__  )
    554. 

  554.       CALL ESMF_TimeIntervalGet( op1_ti, timeString=op1_str, Sn=Sn, Sd=Sd, rc=rc )
    555. 

  555.       CALL test_check_error( ESMF_SUCCESS, rc, &
    556. 

  556.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    557. 

  557.                             __FILE__ , &
    558. 

  558.                             __LINE__  )
    559. 

  559.       ! handle fractions
    560. 

  560.       CALL fraction_to_string( Sn, Sd, frac_str )
    561. 

  561.       op1_str = TRIM(op1_str)//TRIM(frac_str)
    562. 

  562.     ENDIF
    563. 

  563.     ! Initialize op2
    564. 

  564.     IF ( op2_is_t ) THEN
    565. 

  565.       op2_type_str = 'ESMF_Time'
    566. 

  566.       CALL ESMF_TimeSet( op2_t, YY=iop2_t_YY, MM=iop2_t_MM, DD=iop2_t_DD , &
    567. 

  567.                                  H=iop2_t_H, M=iop2_t_M, S=iop2_t_S, Sn=iop2_t_Sn, Sd=iop2_t_Sd, rc=rc )
    568. 

  568.       CALL test_check_error( ESMF_SUCCESS, rc, &
    569. 

  569.                              TRIM(itestname)//'ESMF_TimeSet() ', &
    570. 

  570.                              __FILE__ , &
    571. 

  571.                              __LINE__  )
    572. 

  572.       CALL ESMF_TimeGet( op2_t, timeString=op2_str, Sn=Sn, Sd=Sd, rc=rc )
    573. 

  573.       CALL test_check_error( ESMF_SUCCESS, rc, &
    574. 

  574.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    575. 

  575.                             __FILE__ , &
    576. 

  576.                             __LINE__  )
    577. 

  577.       ! handle fractions
    578. 

  578.       CALL fraction_to_string( Sn, Sd, frac_str )
    579. 

  579.       op2_str = TRIM(op2_str)//TRIM(frac_str)
    580. 

  580.     ELSE IF ( op2_is_ti ) THEN
    581. 

  581.       op2_type_str = 'ESMF_TimeInterval'
    582. 

  582.       CALL ESMF_TimeIntervalSet( op2_ti, YY=iop2_ti_YY, MM=iop2_ti_MM, &
    583. 

  583.                                           D=iop2_ti_DD ,               &
    584. 

  584.                                           H=iop2_ti_H, M=iop2_ti_M,    &
    585. 

  585.                                           S=iop2_ti_S, Sn=iop2_ti_Sn, Sd=iop2_ti_Sd, rc=rc )
    586. 

  586.       CALL test_check_error( ESMF_SUCCESS, rc, &
    587. 

  587.                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
    588. 

  588.                              __FILE__ , &
    589. 

  589.                              __LINE__  )
    590. 

  590.       CALL ESMF_TimeIntervalGet( op2_ti, timeString=op2_str, Sn=Sn, Sd=Sd, rc=rc )
    591. 

  591.       CALL test_check_error( ESMF_SUCCESS, rc, &
    592. 

  592.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    593. 

  593.                             __FILE__ , &
    594. 

  594.                             __LINE__  )
    595. 

  595.       ! handle fractions
    596. 

  596.       CALL fraction_to_string( Sn, Sd, frac_str )
    597. 

  597.       op2_str = TRIM(op2_str)//TRIM(frac_str)
    598. 

  598.     ELSE
    599. 

  599.       op2_type_str = 'INTEGER'
    600. 

  600.       IF ( op2_int > 0 ) THEN
    601. 

  601.         WRITE(op2_str,FMT="('+',I8.8)") ABS(op2_int)
    602. 

  602.       ELSE
    603. 

  603.         WRITE(op2_str,FMT="('-',I8.8)") ABS(op2_int)
    604. 

  604.       ENDIF
    605. 

  605.     ENDIF
    606. 

  606.     ! Initialize res
    607. 

  607.     IF ( res_is_t ) THEN  ! result is ESMF_Time
    608. 

  608.       res_type_str = 'ESMF_Time'
    609. 

  609.       CALL ESMF_TimeSet( res_t, YY=ires_t_YY, MM=ires_t_MM, DD=ires_t_DD , &
    610. 

  610.                                  H=ires_t_H, M=ires_t_M, S=ires_t_S, Sn=ires_t_Sn, Sd=ires_t_Sd, rc=rc )
    611. 

  611.       CALL test_check_error( ESMF_SUCCESS, rc, &
    612. 

  612.                              TRIM(itestname)//'ESMF_TimeSet() ', &
    613. 

  613.                              __FILE__ , &
    614. 

  614.                              __LINE__  )
    615. 

  615.       CALL ESMF_TimeGet( res_t, timeString=res_str, Sn=Sn, Sd=Sd, rc=rc )
    616. 

  616.       CALL test_check_error( ESMF_SUCCESS, rc, &
    617. 

  617.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    618. 

  618.                             __FILE__ , &
    619. 

  619.                             __LINE__  )
    620. 

  620.       ! handle fractions
    621. 

  621.       CALL fraction_to_string( Sn, Sd, frac_str )
    622. 

  622.       res_str = TRIM(res_str)//TRIM(frac_str)
    623. 

  623.     ELSE IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
    624. 

  624.       res_type_str = 'ESMF_TimeInterval'
    625. 

  625.       CALL ESMF_TimeIntervalSet( res_ti, YY=ires_ti_YY, MM=ires_ti_MM, &
    626. 

  626.                                           D=ires_ti_DD ,               &
    627. 

  627.                                           H=ires_ti_H, M=ires_ti_M,    &
    628. 

  628.                                           S=ires_ti_S, Sn=ires_ti_Sn, Sd=ires_ti_Sd, rc=rc )
    629. 

  629.       CALL test_check_error( ESMF_SUCCESS, rc, &
    630. 

  630.                              TRIM(itestname)//'ESMF_TimeIntervalSet() ', &
    631. 

  631.                              __FILE__ , &
    632. 

  632.                              __LINE__  )
    633. 

  633.       CALL ESMF_TimeIntervalGet( res_ti, timeString=res_str, Sn=Sn, Sd=Sd, rc=rc )
    634. 

  634.       CALL test_check_error( ESMF_SUCCESS, rc, &
    635. 

  635.                             TRIM(itestname)//'ESMF_TimeGet() ', &
    636. 

  636.                             __FILE__ , &
    637. 

  637.                             __LINE__  )
    638. 

  638.       ! handle fractions
    639. 

  639.       CALL fraction_to_string( Sn, Sd, frac_str )
    640. 

  640.       res_str = TRIM(res_str)//TRIM(frac_str)
    641. 

  641.     ELSE  ! result is INTEGER
    642. 

  642.       res_type_str = 'INTEGER'
    643. 

  643.       IF ( res_int > 0 ) THEN
    644. 

  644.         WRITE(res_str,FMT="('+',I8.8)") ABS(res_int)
    645. 

  645.       ELSE
    646. 

  646.         WRITE(res_str,FMT="('-',I8.8)") ABS(res_int)
    647. 

  647.       ENDIF
    648. 

  648.     ENDIF
    649. 

  649. 

  650.     ! perform requested operation
    651. 

  651.     unsupported_op = .FALSE.
    652. 

  652.     ! modify behavior of wrf_error_fatal3 for operator being tested
    653. 

  653.     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .TRUE.
    654. 

  654.     ! add
    655. 

  655.     IF ( iadd_op ) THEN
    656. 

  656.       op_str = '+'
    657. 

  657.       IF ( res_is_t ) THEN  ! result is ESMF_Time
    658. 

  658.         IF ( op1_is_t .AND. op2_is_ti ) THEN
    659. 

  659.           !  ESMF_Time         = ESMF_Time         + ESMF_TimeInterval
    660. 

  660.           computed_t = op1_t + op2_ti
    661. 

  661.         ELSE IF ( op1_is_ti .AND. op2_is_t ) THEN
    662. 

  662.           !  ESMF_Time         = ESMF_TimeInterval + ESMF_Time
    663. 

  663.           computed_t = op1_ti + op2_t
    664. 

  664.         ELSE
    665. 

  665.           unsupported_op = .TRUE.
    666. 

  666.         ENDIF
    667. 

  667.       ELSE  ! result is ESMF_TimeInterval
    668. 

  668.         IF ( op1_is_ti .AND. op2_is_ti ) THEN
    669. 

  669.           !  ESMF_TimeInterval = ESMF_TimeInterval + ESMF_TimeInterval
    670. 

  670.           computed_ti = op1_ti + op2_ti
    671. 

  671.         ELSE
    672. 

  672.           unsupported_op = .TRUE.
    673. 

  673.         ENDIF
    674. 

  674.       ENDIF
    675. 

  675.     ! subtract
    676. 

  676.     ELSE  IF ( isubtract_op ) THEN
    677. 

  677.       op_str = '-'
    678. 

  678.       IF ( res_is_t ) THEN  ! result is ESMF_Time
    679. 

  679.         IF ( op1_is_t .AND. op2_is_ti ) THEN
    680. 

  680.           !  ESMF_Time         = ESMF_Time         - ESMF_TimeInterval
    681. 

  681.           computed_t = op1_t - op2_ti
    682. 

  682.         ELSE
    683. 

  683.           unsupported_op = .TRUE.
    684. 

  684.         ENDIF
    685. 

  685.       ELSE  ! result is ESMF_TimeInterval
    686. 

  686.         IF ( op1_is_t .AND. op2_is_t ) THEN
    687. 

  687.           !  ESMF_TimeInterval = ESMF_Time         - ESMF_Time        
    688. 

  688.           computed_ti = op1_t - op2_t
    689. 

  689.         ELSE IF ( op1_is_ti .AND. op2_is_ti ) THEN
    690. 

  690.           !  ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
    691. 

  691.           computed_ti = op1_ti - op2_ti
    692. 

  692.         ELSE
    693. 

  693.           unsupported_op = .TRUE.
    694. 

  694.         ENDIF
    695. 

  695.       ENDIF
    696. 

  696.     ELSE  IF ( imultiply_op ) THEN
    697. 

  697.       op_str = '*'
    698. 

  698.       IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
    699. 

  699.         IF ( op1_is_ti .AND. op2_is_int ) THEN
    700. 

  700.           !  ESMF_TimeInterval = ESMF_TimeInterval * INTEGER
    701. 

  701.           computed_ti = op1_ti * op2_int
    702. 

  702.         ELSE
    703. 

  703.           unsupported_op = .TRUE.
    704. 

  704.         ENDIF
    705. 

  705.       ENDIF
    706. 

  706.     ELSE  IF ( idivide_op ) THEN
    707. 

  707.       op_str = '/'
    708. 

  708.       IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
    709. 

  709.         IF ( op1_is_ti .AND. op2_is_int ) THEN
    710. 

  710.           !  ESMF_TimeInterval = ESMF_TimeInterval / INTEGER
    711. 

  711.           computed_ti = op1_ti / op2_int
    712. 

  712.         ELSE
    713. 

  713.           unsupported_op = .TRUE.
    714. 

  714.         ENDIF
    715. 

  715.       ELSE IF ( res_is_int ) THEN  ! result is INTEGER
    716. 

  716.         IF ( op1_is_ti .AND. op2_is_ti ) THEN
    717. 

  717.           !  INTEGER = ESMF_TimeInterval / ESMF_TimeInterval
    718. 

  718.           ! number of whole time intervals
    719. 

  719.           computed_int = ESMF_TimeIntervalDIVQuot( op1_ti , op2_ti )
    720. 

  720.         ELSE
    721. 

  721.           unsupported_op = .TRUE.
    722. 

  722.         ENDIF
    723. 

  723.       ENDIF
    724. 

  724.     ENDIF
    725. 

  725.     ! restore default behavior of wrf_error_fatal3
    726. 

  726.     IF ( iexpect_error ) WRF_ERROR_FATAL_PRINT = .FALSE.
    727. 

  727.     IF ( unsupported_op ) THEN
    728. 

  728.       WRITE(str,*) 'ERROR test_arithmetic ',TRIM(itestname), &
    729. 

  729.         ':  unsupported operation (',                           &
    730. 

  730.         TRIM(res_type_str),' = ',TRIM(op1_type_str),' ',TRIM(op_str),' ', &
    731. 

  731.         TRIM(op2_type_str),')'
    732. 

  732.       CALL wrf_error_fatal3( __FILE__ , __LINE__ , str )
    733. 

  733.     ENDIF
    734. 

  734. 

  735.     ! check result
    736. 

  736.     test_passed = .FALSE.
    737. 

  737.     IF ( res_is_t ) THEN  ! result is ESMF_Time
    738. 

  738.       IF ( computed_t == res_t ) THEN
    739. 

  739.         test_passed = .TRUE.
    740. 

  740.       ELSE
    741. 

  741.         CALL ESMF_TimeGet( computed_t, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
    742. 

  742.         CALL test_check_error( ESMF_SUCCESS, rc, &
    743. 

  743.                               TRIM(itestname)//'ESMF_TimeGet() ', &
    744. 

  744.                               __FILE__ , &
    745. 

  745.                               __LINE__  )
    746. 

  746.         ! handle fractions
    747. 

  747.         CALL fraction_to_string( Sn, Sd, frac_str )
    748. 

  748.         computed_str = TRIM(computed_str)//TRIM(frac_str)
    749. 

  749.       ENDIF
    750. 

  750.     ELSE IF ( res_is_ti ) THEN  ! result is ESMF_TimeInterval
    751. 

  751.       IF ( computed_ti == res_ti ) THEN
    752. 

  752.         test_passed = .TRUE.
    753. 

  753.       ELSE
    754. 

  754.         CALL ESMF_TimeIntervalGet( computed_ti, timeString=computed_str, Sn=Sn, Sd=Sd, rc=rc )
    755. 

  755.         CALL test_check_error( ESMF_SUCCESS, rc, &
    756. 

  756.                               TRIM(itestname)//'ESMF_TimeGet() ', &
    757. 

  757.                               __FILE__ , &
    758. 

  758.                               __LINE__  )
    759. 

  759.         ! handle fractions
    760. 

  760.         CALL fraction_to_string( Sn, Sd, frac_str )
    761. 

  761.         computed_str = TRIM(computed_str)//TRIM(frac_str)
    762. 

  762.       ENDIF
    763. 

  763.     ELSE  ! result is INTEGER
    764. 

  764.       IF ( computed_int == res_int ) THEN
    765. 

  765.         test_passed = .TRUE.
    766. 

  766.       ELSE
    767. 

  767.         IF ( computed_int > 0 ) THEN
    768. 

  768.           WRITE(computed_str,FMT="('+',I8.8)") ABS(computed_int)
    769. 

  769.         ELSE
    770. 

  770.           WRITE(computed_str,FMT="('-',I8.8)") ABS(computed_int)
    771. 

  771.         ENDIF
    772. 

  772.       ENDIF
    773. 

  773.     ENDIF
    774. 

  774.     IF ( test_passed ) THEN
    775. 

  775.       WRITE(*,FMT='(A)') 'PASS:  '//TRIM(itestname)
    776. 

  776.     ELSE
    777. 

  777.       WRITE(*,*) 'FAIL:  ',TRIM(itestname),':  (',                        &
    778. 

  778.         TRIM(res_type_str),' = ',TRIM(op1_type_str),' ',TRIM(op_str),' ', &
    779. 

  779.         TRIM(op2_type_str),')  expected ',                                &
    780. 

  780.         TRIM(res_str),' = ',TRIM(op1_str),' ',TRIM(op_str),' ',           &
    781. 

  781.         TRIM(op2_str),'  but computed ',TRIM(computed_str)
    782. 

  782.     ENDIF
    783. 

  783. 

  784.   END SUBROUTINE test_arithmetic
    785. 

  785. 

  786. 

  787. 

  788.   ! simple clock creation and advance with add-subtract tests thrown in
    789. 

  789.   ! no self checks (yet)
    790. 

  790.   SUBROUTINE test_clock_advance(                                              &
    791. 

  791.     start_yy, start_mm, start_dd, start_h, start_m, start_s,                  &
    792. 

  792.      stop_yy,  stop_mm,  stop_dd,  stop_h,  stop_m,  stop_s,                  &
    793. 

  793.     timestep_d, timestep_h, timestep_m, timestep_s, timestep_sn, timestep_sd, &
    794. 

  794.     testname, increment_S, increment_Sn, increment_Sd )
    795. 

  795.     INTEGER, INTENT(IN), OPTIONAL :: start_YY
    796. 

  796.     INTEGER, INTENT(IN), OPTIONAL :: start_MM  ! month
    797. 

  797.     INTEGER, INTENT(IN), OPTIONAL :: start_DD  ! day of month
    798. 

  798.     INTEGER, INTENT(IN), OPTIONAL :: start_H
    799. 

  799.     INTEGER, INTENT(IN), OPTIONAL :: start_M
    800. 

  800.     INTEGER, INTENT(IN), OPTIONAL :: start_S
    801. 

  801.     INTEGER, INTENT(IN), OPTIONAL :: stop_YY
    802. 

  802.     INTEGER, INTENT(IN), OPTIONAL :: stop_MM  ! month
    803. 

  803.     INTEGER, INTENT(IN), OPTIONAL :: stop_DD  ! day of month
    804. 

  804.     INTEGER, INTENT(IN), OPTIONAL :: stop_H
    805. 

  805.     INTEGER, INTENT(IN), OPTIONAL :: stop_M
    806. 

  806.     INTEGER, INTENT(IN), OPTIONAL :: stop_S
    807. 

  807.     INTEGER, INTENT(IN), OPTIONAL :: timestep_D  ! day
    808. 

  808.     INTEGER, INTENT(IN), OPTIONAL :: timestep_H
    809. 

  809.     INTEGER, INTENT(IN), OPTIONAL :: timestep_M
    810. 

  810.     INTEGER, INTENT(IN), OPTIONAL :: timestep_S
    811. 

  811.     INTEGER, INTENT(IN), OPTIONAL :: timestep_Sn
    812. 

  812.     INTEGER, INTENT(IN), OPTIONAL :: timestep_Sd
    813. 

  813.     CHARACTER (LEN=*), OPTIONAL, INTENT(IN) :: testname
    814. 

  814.     INTEGER, INTENT(IN), OPTIONAL :: increment_S  ! add and subtract this 
    815. 

  815.     INTEGER, INTENT(IN), OPTIONAL :: increment_Sn ! value each time step
    816. 

  816.     INTEGER, INTENT(IN), OPTIONAL :: increment_Sd
    817. 

  817. 

  818.     ! locals
    819. 

  819.     INTEGER :: istart_YY
    820. 

  820.     INTEGER :: istart_MM  ! month
    821. 

  821.     INTEGER :: istart_DD  ! day of month
    822. 

  822.     INTEGER :: istart_H
    823. 

  823.     INTEGER :: istart_M
    824. 

  824.     INTEGER :: istart_S
    825. 

  825.     INTEGER :: istop_YY
    826. 

  826.     INTEGER :: istop_MM  ! month
    827. 

  827.     INTEGER :: istop_DD  ! day of month
    828. 

  828.     INTEGER :: istop_H
    829. 

  829.     INTEGER :: istop_M
    830. 

  830.     INTEGER :: istop_S
    831. 

  831.     INTEGER :: itimestep_D  ! day
    832. 

  832.     INTEGER :: itimestep_H
    833. 

  833.     INTEGER :: itimestep_M
    834. 

  834.     INTEGER :: itimestep_S
    835. 

  835.     INTEGER :: itimestep_Sn
    836. 

  836.     INTEGER :: itimestep_Sd
    837. 

  837.     CHARACTER (LEN=512) :: itestname, itestfullname
    838. 

  838.     INTEGER :: iincrement_S
    839. 

  839.     INTEGER :: iincrement_Sn
    840. 

  840.     INTEGER :: iincrement_Sd
    841. 

  841.     INTEGER :: Sn, Sd
    842. 

  842.     INTEGER rc
    843. 

  843.     TYPE(ESMF_Time)           :: start_time, stop_time, current_time
    844. 

  844.     TYPE(ESMF_Clock), POINTER :: domain_clock
    845. 

  845.     TYPE(ESMF_TimeInterval)   :: timestep, increment
    846. 

  846.     TYPE(ESMF_Time)           :: add_time, subtract_time
    847. 

  847.     INTEGER :: itimestep
    848. 

  848.     REAL(ESMF_KIND_R8) :: dayr8
    849. 

  849.     CHARACTER(LEN=ESMF_MAXSTR) :: str, frac_str
    850. 

  850. 

  851.     istart_YY = 0
    852. 

  852.     istart_MM = 1
    853. 

  853.     istart_DD = 1
    854. 

  854.     istart_H = 0
    855. 

  855.     istart_M = 0
    856. 

  856.     istart_S = 0
    857. 

  857.     istop_YY = 0
    858. 

  858.     istop_MM = 1
    859. 

  859.     istop_DD = 1
    860. 

  860.     istop_H = 0
    861. 

  861.     istop_M = 0
    862. 

  862.     istop_S = 0
    863. 

  863.     itimestep_D = 0
    864. 

  864.     itimestep_H = 0
    865. 

  865.     itimestep_M = 0
    866. 

  866.     itimestep_S = 0
    867. 

  867.     itimestep_Sn = 0
    868. 

  868.     itimestep_Sd = 0
    869. 

  869.     itestname = ''
    870. 

  870.     iincrement_S = 0
    871. 

  871.     iincrement_Sn = 0
    872. 

  872.     iincrement_Sd = 0
    873. 

  873. 

  874.     IF ( PRESENT( start_YY ) ) istart_YY = start_YY
    875. 

  875.     IF ( PRESENT( start_MM ) ) istart_MM = start_MM
    876. 

  876.     IF ( PRESENT( start_DD ) ) istart_DD = start_DD
    877. 

  877.     IF ( PRESENT( start_H ) ) istart_H = start_H
    878. 

  878.     IF ( PRESENT( start_M ) ) istart_M = start_M
    879. 

  879.     IF ( PRESENT( start_S ) ) istart_S = start_S
    880. 

  880.     IF ( PRESENT( stop_YY ) ) istop_YY = stop_YY
    881. 

  881.     IF ( PRESENT( stop_MM ) ) istop_MM = stop_MM
    882. 

  882.     IF ( PRESENT( stop_DD ) ) istop_DD = stop_DD
    883. 

  883.     IF ( PRESENT( stop_H ) ) istop_H = stop_H
    884. 

  884.     IF ( PRESENT( stop_M ) ) istop_M = stop_M
    885. 

  885.     IF ( PRESENT( stop_S ) ) istop_S = stop_S
    886. 

  886.     IF ( PRESENT( timestep_D ) ) itimestep_D = timestep_D
    887. 

  887.     IF ( PRESENT( timestep_H ) ) itimestep_H = timestep_H
    888. 

  888.     IF ( PRESENT( timestep_M ) ) itimestep_M = timestep_M
    889. 

  889.     IF ( PRESENT( timestep_S ) ) itimestep_S = timestep_S
    890. 

  890.     IF ( PRESENT( timestep_Sn ) ) itimestep_Sn = timestep_Sn
    891. 

  891.     IF ( PRESENT( timestep_Sd ) ) itimestep_Sd = timestep_Sd
    892. 

  892.     IF ( PRESENT( testname ) ) itestname = TRIM(testname)//'_'
    893. 

  893.     IF ( PRESENT( increment_S ) ) iincrement_S = increment_S
    894. 

  894.     IF ( PRESENT( increment_Sn ) ) iincrement_Sn = increment_Sn
    895. 

  895.     IF ( PRESENT( increment_Sd ) ) iincrement_Sd = increment_Sd
    896. 

  896. 

  897.     ! Initialize start time, stop time, time step, clock for simple case. 
    898. 

  898.     itestfullname = TRIM(itestname)//'SETUP'
    899. 

  899.     CALL ESMF_TimeSet( start_time, YY=istart_YY, MM=istart_MM, DD=istart_DD , &
    900. 

  900.                                    H=istart_H, M=istart_M, S=istart_S, rc=rc )
    901. 

  901.     CALL test_check_error( ESMF_SUCCESS, rc, &
    902. 

  902.                           TRIM(itestfullname)//'ESMF_TimeSet() ', &
    903. 

  903.                           __FILE__ , &
    904. 

  904.                           __LINE__  )
    905. 

  905. 

  906.     CALL ESMF_TimeGet( start_time, timeString=str, rc=rc )
    907. 

  907.     CALL test_check_error( ESMF_SUCCESS, rc, &
    908. 

  908.                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
    909. 

  909.                           __FILE__ , &
    910. 

  910.                           __LINE__  )
    911. 

  911.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  start_time = <',TRIM(str),'>'
    912. 

  912. 

  913.     CALL ESMF_TimeSet( stop_time, YY=istop_YY, MM=istop_MM, DD=istop_DD , &
    914. 

  914.                                    H=istop_H, M=istop_M, S=istop_S, rc=rc )
    915. 

  915.     CALL test_check_error( ESMF_SUCCESS, rc, &
    916. 

  916.                           TRIM(itestfullname)//'ESMF_TimeSet() ', &
    917. 

  917.                           __FILE__ , &
    918. 

  918.                           __LINE__  )
    919. 

  919. 

  920.     CALL ESMF_TimeGet( stop_time, timeString=str, rc=rc )
    921. 

  921.     CALL test_check_error( ESMF_SUCCESS, rc, &
    922. 

  922.                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
    923. 

  923.                           __FILE__ , &
    924. 

  924.                           __LINE__  )
    925. 

  925.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  stop_time = <',TRIM(str),'>'
    926. 

  926. 

  927.     CALL ESMF_TimeIntervalSet( timestep, D=itimestep_D, H=itimestep_H, &
    928. 

  928.                                          M=itimestep_M, S=itimestep_S, &
    929. 

  929.                                          Sn=itimestep_Sn, Sd=itimestep_Sd, rc=rc )
    930. 

  930.     CALL test_check_error( ESMF_SUCCESS, rc, &
    931. 

  931.                           TRIM(itestfullname)//'ESMF_TimeIntervalSet() ', &
    932. 

  932.                           __FILE__ , &
    933. 

  933.                           __LINE__  )
    934. 

  934. 

  935.     CALL ESMF_TimeIntervalGet( timestep, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    936. 

  936.     CALL test_check_error( ESMF_SUCCESS, rc, &
    937. 

  937.                           TRIM(itestfullname)//'ESMF_TimeIntervalGet() ', &
    938. 

  938.                           __FILE__ , &
    939. 

  939.                           __LINE__  )
    940. 

  940.     ! handle fractions
    941. 

  941.     CALL fraction_to_string( Sn, Sd, frac_str )
    942. 

  942.     str = TRIM(str)//TRIM(frac_str)
    943. 

  943.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  timestep = <',TRIM(str),'>'
    944. 

  944. 

  945.     CALL ESMF_TimeIntervalSet( increment, S=iincrement_S, &
    946. 

  946.                                Sn=iincrement_Sn, Sd=iincrement_Sd, rc=rc )
    947. 

  947.     CALL test_check_error( ESMF_SUCCESS, rc, &
    948. 

  948.                           TRIM(itestfullname)//'ESMF_TimeIntervalSet() ', &
    949. 

  949.                           __FILE__ , &
    950. 

  950.                           __LINE__  )
    951. 

  951. 

  952.     CALL ESMF_TimeIntervalGet( increment, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    953. 

  953.     CALL test_check_error( ESMF_SUCCESS, rc, &
    954. 

  954.                           TRIM(itestfullname)//'ESMF_TimeIntervalGet() ', &
    955. 

  955.                           __FILE__ , &
    956. 

  956.                           __LINE__  )
    957. 

  957.     ! handle fractions
    958. 

  958.     CALL fraction_to_string( Sn, Sd, frac_str )
    959. 

  959.     str = TRIM(str)//TRIM(frac_str)
    960. 

  960.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  increment = <',TRIM(str),'>'
    961. 

  961. 

  962.     ALLOCATE( domain_clock )
    963. 

  963.     domain_clock = ESMF_ClockCreate( TimeStep= timestep,  &
    964. 

  964.                                      StartTime=start_time, &
    965. 

  965.                                      StopTime= stop_time,  &
    966. 

  966.                                      rc=rc )
    967. 

  967.     CALL test_check_error( ESMF_SUCCESS, rc, &
    968. 

  968.                           TRIM(itestfullname)//'ESMF_ClockCreate() ', &
    969. 

  969.                           __FILE__ , &
    970. 

  970.                           __LINE__  )
    971. 

  971. 

  972.     CALL ESMF_ClockGet( domain_clock, CurrTime=current_time, &
    973. 

  973.                         rc=rc )
    974. 

  974.     CALL test_check_error( ESMF_SUCCESS, rc, &
    975. 

  975.                           TRIM(itestfullname)//'ESMF_ClockGet() ', &
    976. 

  976.                           __FILE__ , &
    977. 

  977.                           __LINE__  )
    978. 

  978. 

  979.     CALL ESMF_TimeGet( current_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    980. 

  980.     CALL test_check_error( ESMF_SUCCESS, rc, &
    981. 

  981.                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
    982. 

  982.                           __FILE__ , &
    983. 

  983.                           __LINE__  )
    984. 

  984.     CALL fraction_to_string( Sn, Sd, frac_str )
    985. 

  985.     str = TRIM(str)//TRIM(frac_str)
    986. 

  986.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  clock current_time = <',TRIM(str),'>'
    987. 

  987. 

  988.     CALL ESMF_TimeGet( current_time, dayOfYear_r8=dayr8, rc=rc )
    989. 

  989.     CALL test_check_error( ESMF_SUCCESS, rc, &
    990. 

  990.                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
    991. 

  991.                           __FILE__ , &
    992. 

  992.                           __LINE__  )
    993. 

  993.     WRITE(*,FMT='(A,A,F10.6,A)') TRIM(itestfullname),':  current_time dayOfYear_r8 = < ',dayr8,' >'
    994. 

  994. 

  995.     subtract_time = current_time - increment
    996. 

  996.     CALL ESMF_TimeGet( subtract_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    997. 

  997.     CALL test_check_error( ESMF_SUCCESS, rc, &
    998. 

  998.                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
    999. 

  999.                           __FILE__ , &
    1000. 

  1000.                           __LINE__  )
    1001. 

  1001.     CALL fraction_to_string( Sn, Sd, frac_str )
    1002. 

  1002.     str = TRIM(str)//TRIM(frac_str)
    1003. 

  1003.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time-increment = <',TRIM(str),'>'
    1004. 

  1004. 

  1005.     add_time = current_time + increment
    1006. 

  1006.     CALL ESMF_TimeGet( add_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    1007. 

  1007.     CALL test_check_error( ESMF_SUCCESS, rc, &
    1008. 

  1008.                           TRIM(itestfullname)//'ESMF_TimeGet() ', &
    1009. 

  1009.                           __FILE__ , &
    1010. 

  1010.                           __LINE__  )
    1011. 

  1011.     CALL fraction_to_string( Sn, Sd, frac_str )
    1012. 

  1012.     str = TRIM(str)//TRIM(frac_str)
    1013. 

  1013.     WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time+increment = <',TRIM(str),'>'
    1014. 

  1014. 

  1015.     ! Advance clock.  
    1016. 

  1016.     itestfullname = TRIM(itestname)//'ADVANCE'
    1017. 

  1017.     itimestep = 0
    1018. 

  1018.     DO WHILE ( .NOT. ESMF_ClockIsStopTime(domain_clock ,rc=rc) )
    1019. 

  1019.       CALL test_check_error( ESMF_SUCCESS, rc, &
    1020. 

  1020.                             TRIM(itestfullname)//'ESMF_ClockIsStopTime() ', &
    1021. 

  1021.                             __FILE__ , &
    1022. 

  1022.                             __LINE__  )
    1023. 

  1023.       itimestep = itimestep + 1
    1024. 

  1024. 

  1025.       CALL ESMF_ClockAdvance( domain_clock, rc=rc )
    1026. 

  1026.       CALL test_check_error( ESMF_SUCCESS, rc, &
    1027. 

  1027.                             TRIM(itestfullname)//'ESMF_ClockAdvance() ', &
    1028. 

  1028.                             __FILE__ , &
    1029. 

  1029.                             __LINE__  )
    1030. 

  1030. 

  1031.       CALL ESMF_ClockGet( domain_clock, CurrTime=current_time, &
    1032. 

  1032.                           rc=rc )
    1033. 

  1033.       CALL test_check_error( ESMF_SUCCESS, rc, &
    1034. 

  1034.                             TRIM(itestfullname)//'ESMF_ClockGet() ', &
    1035. 

  1035.                             __FILE__ , &
    1036. 

  1036.                             __LINE__  )
    1037. 

  1037. 

  1038.       CALL ESMF_TimeGet( current_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    1039. 

  1039.       CALL test_check_error( ESMF_SUCCESS, rc, &
    1040. 

  1040.                             TRIM(itestfullname)//'ESMF_TimeGet() ', &
    1041. 

  1041.                             __FILE__ , &
    1042. 

  1042.                             __LINE__  )
    1043. 

  1043.       CALL fraction_to_string( Sn, Sd, frac_str )
    1044. 

  1044.       str = TRIM(str)//TRIM(frac_str)
    1045. 

  1045.       WRITE(*,FMT='(A,A,I6.6,A,A,A)') TRIM(itestfullname),':  count = ', &
    1046. 

  1046.         itimestep,'  current_time = <',TRIM(str),'>'
    1047. 

  1047. 

  1048.       subtract_time = current_time - increment
    1049. 

  1049.       CALL ESMF_TimeGet( subtract_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    1050. 

  1050.       CALL test_check_error( ESMF_SUCCESS, rc, &
    1051. 

  1051.                             TRIM(itestfullname)//'ESMF_TimeGet() ', &
    1052. 

  1052.                             __FILE__ , &
    1053. 

  1053.                             __LINE__  )
    1054. 

  1054.       CALL fraction_to_string( Sn, Sd, frac_str )
    1055. 

  1055.       str = TRIM(str)//TRIM(frac_str)
    1056. 

  1056.       WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time-increment = <',TRIM(str),'>'
    1057. 

  1057. 

  1058.       add_time = current_time + increment
    1059. 

  1059.       CALL ESMF_TimeGet( add_time, timeString=str, Sn=Sn, Sd=Sd, rc=rc )
    1060. 

  1060.       CALL test_check_error( ESMF_SUCCESS, rc, &
    1061. 

  1061.                             TRIM(itestfullname)//'ESMF_TimeGet() ', &
    1062. 

  1062.                             __FILE__ , &
    1063. 

  1063.                             __LINE__  )
    1064. 

  1064.       CALL fraction_to_string( Sn, Sd, frac_str )
    1065. 

  1065.       str = TRIM(str)//TRIM(frac_str)
    1066. 

  1066.       WRITE(*,FMT='(A,A,A,A)') TRIM(itestfullname),':  current_time+increment = <',TRIM(str),'>'
    1067. 

  1067. 

  1068.     ENDDO
    1069. 

  1069. 

  1070.     DEALLOCATE( domain_clock )
    1071. 

  1071.   
    1072. 

  1072.   END SUBROUTINE test_clock_advance
    1073. 

  1073. 

  1074. END MODULE my_tests
    1075. 

  1075. 

  1076. 

  1077. #if defined( TIME_F90_ONLY ) 
    1078. 

  1078. 

  1079. ! TBH:  Improve the build of Test1.exe to use WRF versions of these 
    1080. 

  1080. ! TBH:  routines and remove these hacked-in duplicates!!  
    1081. 

  1081. 

  1082. SUBROUTINE wrf_abort
    1083. 

  1083.   IMPLICIT NONE
    1084. 

  1084. #if defined( DM_PARALLEL ) && ! defined( STUBMPI )
    1085. 

  1085.   INCLUDE 'mpif.h'
    1086. 

  1086.   INTEGER ierr
    1087. 

  1087.   CALL mpi_abort(MPI_COMM_WORLD,1,ierr)
    1088. 

  1088. #else
    1089. 

  1089.   STOP
    1090. 

  1090. #endif
    1091. 

  1091. END SUBROUTINE wrf_abort
    1092. 

  1092. 

  1093. SUBROUTINE wrf_message( str )
    1094. 

  1094.   IMPLICIT NONE
    1095. 

  1095.   CHARACTER*(*) str
    1096. 

  1096. #if defined( DM_PARALLEL ) && ! defined( STUBMPI)
    1097. 

  1097.   write(0,*) str
    1098. 

  1098. #endif
    1099. 

  1099.   print*, str
    1100. 

  1100. END SUBROUTINE wrf_message
    1101. 

  1101. 

  1102. ! intentionally write to stderr only
    1103. 

  1103. SUBROUTINE wrf_message2( str )
    1104. 

  1104.   IMPLICIT NONE
    1105. 

  1105.   CHARACTER*(*) str
    1106. 

  1106.   write(0,*) str
    1107. 

  1107. END SUBROUTINE wrf_message2
    1108. 

  1108. 

  1109. SUBROUTINE wrf_error_fatal3( file_str, line, str )
    1110. 

  1110.   USE my_tests
    1111. 

  1111.   IMPLICIT NONE
    1112. 

  1112.   CHARACTER*(*) file_str
    1113. 

  1113.   INTEGER , INTENT (IN) :: line  ! only print file and line if line > 0
    1114. 

  1114.   CHARACTER*(*) str
    1115. 

  1115.   CHARACTER*256 :: line_str
    1116. 

  1116.   write(line_str,'(i6)') line
    1117. 

  1117.   ! special behavior for testing since Fortran cannot catch exceptions
    1118. 

  1118.  IF ( WRF_ERROR_FATAL_PRINT ) THEN
    1119. 

  1119.   ! just print message and continue
    1120. 

  1120.   CALL wrf_message( 'ERROR IN FILE:  '//TRIM(file_str)//'  LINE:  '//TRIM(line_str) )
    1121. 

  1121.  ELSE
    1122. 

  1122.   ! normal behavior
    1123. 

  1123. #if defined( DM_PARALLEL ) && ! defined( STUBMPI )
    1124. 

  1124.   CALL wrf_message( '-------------- FATAL CALLED ---------------' )
    1125. 

  1125.   ! only print file and line if line is positive
    1126. 

  1126.   IF ( line > 0 ) THEN
    1127. 

  1127.     CALL wrf_message( 'FATAL CALLED FROM FILE:  '//file_str//'  LINE:  '//TRIM(line_str) )
    1128. 

  1128.   ENDIF
    1129. 

  1129.   CALL wrf_message( str )
    1130. 

  1130.   CALL wrf_message( '-------------------------------------------' )
    1131. 

  1131. #else
    1132. 

  1132.   CALL wrf_message2( '-------------- FATAL CALLED ---------------' )
    1133. 

  1133.   ! only print file and line if line is positive
    1134. 

  1134.   IF ( line > 0 ) THEN
    1135. 

  1135.     CALL wrf_message( 'FATAL CALLED FROM FILE:  '//file_str//'  LINE:  '//TRIM(line_str) )
    1136. 

  1136.   ENDIF
    1137. 

  1137.   CALL wrf_message2( str )
    1138. 

  1138.   CALL wrf_message2( '-------------------------------------------' )
    1139. 

  1139. #endif
    1140. 

  1140.   CALL wrf_abort
    1141. 

  1141.  ENDIF
    1142. 

  1142. END SUBROUTINE wrf_error_fatal3
    1143. 

  1143. 

  1144. SUBROUTINE wrf_error_fatal( str )
    1145. 

  1145.   IMPLICIT NONE
    1146. 

  1146.   CHARACTER*(*) str
    1147. 

  1147.   CALL wrf_error_fatal3 ( ' ', 0, str )
    1148. 

  1148. END SUBROUTINE wrf_error_fatal
    1149. 

  1149. 

  1150. #endif
    1151. 

  1151. 

  1152. 

  1153. ! Check to see if expected value == actual value
    1154. 

  1154. ! If not, print message and exit.
    1155. 

  1155. SUBROUTINE test_check_error( expected, actual, str, file_str, line )
    1156. 

  1156.   IMPLICIT NONE
    1157. 

  1157.   INTEGER , INTENT (IN) :: expected
    1158. 

  1158.   INTEGER , INTENT (IN) :: actual
    1159. 

  1159.   CHARACTER*(*) str
    1160. 

  1160.   CHARACTER*(*) file_str
    1161. 

  1161.   INTEGER , INTENT (IN) :: line
    1162. 

  1162.   CHARACTER (LEN=512)   :: rc_str
    1163. 

  1163.   CHARACTER (LEN=512)   :: str_with_rc
    1164. 

  1164.   IF ( expected .ne. actual ) THEN
    1165. 

  1165.     WRITE (rc_str,*) '  Routine returned error code = ',actual
    1166. 

  1166.     str_with_rc = 'FAIL:  '//TRIM(str)//TRIM(rc_str)
    1167. 

  1167.     CALL wrf_error_fatal3( file_str, line, str_with_rc )
    1168. 

  1168.   ENDIF
    1169. 

  1169. END SUBROUTINE test_check_error
    1170. 

  1170. 

  1171. 

  1172. 

  1173. PROGRAM time_manager_test
    1174. 

  1174.   USE ESMF_Mod
    1175. 

  1175.   USE my_tests
    1176. 

  1176.   IMPLICIT NONE
    1177. 

  1177.   INTEGER :: rc
    1178. 

  1178. 

  1179.   PRINT *,'BEGIN TEST SUITE'
    1180. 

  1180. 

  1181.   CALL ESMF_Initialize( defaultCalendar=ESMF_CAL_GREGORIAN, rc=rc )
    1182. 

  1182.   CALL test_check_error( ESMF_SUCCESS, rc, &
    1183. 

  1183.                         'ESMF_Initialize() ', &
    1184. 

  1184.                         __FILE__ , &
    1185. 

  1185.                         __LINE__  )
    1186. 

  1186. !  PRINT *,'DEBUG:  back from ESMF_Initialize(), rc = ',rc
    1187. 

  1187. 

  1188. !  CALL test_print(  t_yy,  t_mm,  t_dd,  t_h,  t_m,  t_s, &
    1189. 

  1189. !                   ti_yy, ti_mm, ti_dd, ti_h, ti_m, ti_s, &
    1190. 

  1190. !                   res_str, testname )
    1191. 

  1191. 

  1192.   ! Print times
    1193. 

  1193.   ! "vanilla" tests
    1194. 

  1194. !  PRINT *,'DEBUG:  calling 1st test_print()'
    1195. 

  1195.   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
    1196. 

  1196.     res_str='2001-12-03_01:20:10', testname='printT_1' )
    1197. 

  1197. !  PRINT *,'DEBUG:  back from 1st test_print()'
    1198. 

  1198.   CALL test_print( t_yy=0,  t_mm=1,  t_dd=1,  t_h=0,  t_m=0,  t_s=0, &
    1199. 

  1199.     res_str='0000-01-01_00:00:00', testname='printT_2' )
    1200. 

  1200.   CALL test_print( t_yy=2003,  t_mm=12,  t_dd=30,  t_h=23,  t_m=59,  t_s=50, &
    1201. 

  1201.     res_str='2003-12-30_23:59:50', testname='printT_3' )
    1202. 

  1202.   CALL test_print( t_yy=2003,  t_mm=12,  t_dd=31,  t_h=23,  t_m=59,  t_s=50, &
    1203. 

  1203.     res_str='2003-12-31_23:59:50', testname='printT_4' )
    1204. 

  1204.   CALL test_print( t_yy=2004,  t_mm=12,  t_dd=30,  t_h=23,  t_m=59,  t_s=50, &
    1205. 

  1205.     res_str='2004-12-30_23:59:50', testname='printT_5' )
    1206. 

  1206.   CALL test_print( t_yy=2004,  t_mm=12,  t_dd=31,  t_h=23,  t_m=59,  t_s=50, &
    1207. 

  1207.     res_str='2004-12-31_23:59:50', testname='printT_6' )
    1208. 

  1208. !$$$  NOTE that this fails -- need to fix up output string for negative year
    1209. 

  1209. !  CALL test_print( t_yy=-2004,  t_mm=12,  t_dd=31,  t_h=23,  t_m=59,  t_s=50, &
    1210. 

  1210. !    res_str='-2004-12-31_23:59:50', testname='printT_6' )
    1211. 

  1211. 

  1212.   ! these test default behavior of test harness
    1213. 

  1213.   CALL test_print( t_s=0, &
    1214. 

  1214.     res_str='0000-01-01_00:00:00', testname='printT_D1' )
    1215. 

  1215.   CALL test_print( t_yy=0, &
    1216. 

  1216.     res_str='0000-01-01_00:00:00', testname='printT_D2' )
    1217. 

  1217. 

  1218.   ! fractions
    1219. 

  1219.   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
    1220. 

  1220.     t_sn=1, t_sd=3, &
    1221. 

  1221.     res_str='2001-12-03_01:20:10+01/03', testname='printT_F1' )
    1222. 

  1222.   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
    1223. 

  1223.     t_sn=4, t_sd=3, &
    1224. 

  1224.     res_str='2001-12-03_01:20:11+01/03', testname='printT_F2' )
    1225. 

  1225.   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
    1226. 

  1226.     t_sn=12, t_sd=3, &
    1227. 

  1227.     res_str='2001-12-03_01:20:14', testname='printT_F3' )
    1228. 

  1228.   CALL test_print( t_yy=2001,  t_mm=12,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
    1229. 

  1229.     t_sn=-1, t_sd=3, &
    1230. 

  1230.     res_str='2001-12-03_01:20:09+02/03', testname='printT_F4' )
    1231. 

  1231. 

  1232.   ! ERROR, MM out of range
    1233. 

  1233. !$$$here...  fix so this just prints "ERROR:  <testname>" in failure case
    1234. 

  1234. !$$$here...  also need "expect_fail" to reverse sense of PASS/FAIL message for 
    1235. 

  1235. !$$$here...  tests that should fail
    1236. 

  1236. !  CALL test_print( t_yy=2001,  t_mm=13,  t_dd=3,  t_h=1,  t_m=20,  t_s=10, &
    1237. 

  1237. !    res_str='2002-01-03_01:20:10', testname='printT_E1', expect_error=.TRUE. )
    1238. 

  1238. 

  1239.   ! Print time intervals
    1240. 

  1240.   ! "vanilla" tests
    1241. 

  1241.   CALL test_print( ti_yy=0,  ti_mm=0,  ti_dd=0,  ti_h=0,  ti_m=0,  ti_s=0, &
    1242. 

  1242.     res_str='0000000000_000:000:000', testname='printTI_1' )
    1243. 

  1243.   CALL test_print( ti_yy=0,  ti_mm=0,  ti_dd=500,  ti_h=0,  ti_m=0,  ti_s=7270, &
    1244. 

  1244.     res_str='0000000500_002:001:010', testname='printTI_2' )
    1245. 

  1245. 

  1246.   ! these test default behavior of test harness
    1247. 

  1247.   CALL test_print( ti_s=0, &
    1248. 

  1248.     res_str='0000000000_000:000:000', testname='printTI_D1' )
    1249. 

  1249.   CALL test_print( ti_yy=0, &
    1250. 

  1250.     res_str='0000000000_000:000:000', testname='printTI_D2' )
    1251. 

  1251. 

  1252.   ! these test negative values
    1253. 

  1253.   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=-1,  ti_m=-20,  ti_s=-10, &
    1254. 

  1254.     res_str='-0000000003_001:020:010', testname='printTI_N1' )
    1255. 

  1255. 

  1256.   ! these test mixed values
    1257. 

  1257.   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=1,  ti_m=20,  ti_s=10, &
    1258. 

  1258.     res_str='-0000000002_022:039:050', testname='printTI_M1' )
    1259. 

  1259. 

  1260.   ! fractions
    1261. 

  1261.   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
    1262. 

  1262.     ti_sn=1, ti_sd=3, &
    1263. 

  1263.     res_str='0000000003_001:020:010+01/03', testname='printTI_F1' )
    1264. 

  1264.   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
    1265. 

  1265.     ti_sn=5, ti_sd=3, &
    1266. 

  1266.     res_str='0000000003_001:020:011+02/03', testname='printTI_F2' )
    1267. 

  1267.   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=-1,  ti_m=-20,  ti_s=-10, &
    1268. 

  1268.     ti_sn=-1, ti_sd=3, &
    1269. 

  1269.     res_str='-0000000003_001:020:010-01/03', testname='printTI_F3' )
    1270. 

  1270.   CALL test_print( ti_yy=0000,  ti_mm=0,  ti_dd=-3,  ti_h=-1,  ti_m=-20,  ti_s=-10, &
    1271. 

  1271.     ti_sn=1, ti_sd=3, &
    1272. 

  1272.     res_str='-0000000003_001:020:009-02/03', testname='printTI_F4' )
    1273. 

  1273. 

  1274.   ! these test non-normalized values
    1275. 

  1275. !  CALL test_print( ti_yy=2001,  ti_mm=1,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
    1276. 

  1276. !    res_str='02001-001-003_001:020:010', testname='printTI_NN1', expect_error=.TRUE. )
    1277. 

  1277. !  CALL test_print( ti_yy=2001,  ti_mm=12,  ti_dd=3,  ti_h=1,  ti_m=20,  ti_s=10, &
    1278. 

  1278. !    res_str='02002-000-003_001:020:010', testname='printTI_NN2', expect_error=.TRUE. )
    1279. 

  1279. !  CALL test_print( ti_yy=2002,  ti_mm=5,  ti_dd=500,  ti_h=0,  ti_m=0,  ti_s=7270, &
    1280. 

  1280. !    res_str='02002-005-500_002:001:010', testname='printTI_NN3', expect_error=.TRUE. )
    1281. 

  1281. 

  1282.   ! Addition tests
    1283. 

  1283.   ! ESMF_Time = ESMF_Time + ESMF_TimeInterval
    1284. 

  1284.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1285. 

  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. 

  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. 

  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. 

  1288.     testname='AddT_T_TI1' )
    1289. 

  1289.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1290. 

  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. 

  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. 

  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. 

  1293.     testname='AddT_T_TI2' )
    1294. 

  1294.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1295. 

  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. 

  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. 

  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. 

  1298.     testname='AddT_T_TI3' )
    1299. 

  1299.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1300. 

  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. 

  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. 

  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. 

  1303.     testname='AddT_T_TI4' )
    1304. 

  1304.   ! this case hung after the CCSM contribution
    1305. 

  1305.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1306. 

  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. 

  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. 

  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. 

  1309.     testname='AddT_T_TI5' )
    1310. 

  1310. ! NOTE:  CCSM folks need to decide what it means to add "1 month" to Feb. 29.  And all the 
    1311. 

  1311. !        other very similar cases.  Then, write this unit test!  
    1312. 

  1312. !  CALL test_arithmetic( add_op=.TRUE.,                                             &
    1313. 

  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. 

  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. 

  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. 

  1316. !    testname='AddT_T_TI6' )
    1317. 

  1317.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1318. 

  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. 

  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. 

  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. 

  1321.     testname='AddT_T_TI7' )
    1322. 

  1322.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1323. 

  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. 

  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. 

  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. 

  1326.     testname='AddT_T_TI8' )
    1327. 

  1327.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1328. 

  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. 

  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. 

  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. 

  1331.     testname='AddT_T_TI9' )
    1332. 

  1332.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1333. 

  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. 

  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. 

  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. 

  1336.     testname='AddT_T_TI10' )
    1337. 

  1337.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1338. 

  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. 

  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. 

  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. 

  1341.     testname='AddT_T_TI11' )
    1342. 

  1342.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1343. 

  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. 

  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. 

  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. 

  1346.     testname='AddT_T_TI12' )
    1347. 

  1347.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1348. 

  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. 

  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. 

  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. 

  1351.     testname='AddT_T_TI13' )
    1352. 

  1352.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1353. 

  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. 

  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. 

  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. 

  1356.     testname='AddT_T_TI14' )
    1357. 

  1357.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1358. 

  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. 

  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. 

  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. 

  1361.     testname='AddT_T_TI15' )
    1362. 

  1362.   ! ESMF_Time = ESMF_Time + ESMF_TimeInterval with fractions
    1363. 

  1363.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1364. 

  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. 

  1365.      op1_t_sn=01,  op1_t_sd=03, &
    1366. 

  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. 

  1367.     op2_ti_sn=01, op2_ti_sd=03, &
    1368. 

  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. 

  1369.      res_t_sn=02,  res_t_sd=03, &
    1370. 

  1370.     testname='AddT_T_TI_F1' )
    1371. 

  1371.   ! this should fail (and does)
    1372. 

  1372. !  CALL test_arithmetic( add_op=.TRUE.,                                             &
    1373. 

  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. 

  1374. !     op1_t_sn=01,  op1_t_sd=03, &
    1375. 

  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. 

  1376. !    op2_ti_sn=01, op2_ti_sd=03, &
    1377. 

  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. 

  1378. !     res_t_sn=01,  res_t_sd=03, &
    1379. 

  1379. !    testname='AddT_T_TI_F2' )
    1380. 

  1380.   ! ESMF_Time = ESMF_TimeInterval + ESMF_Time
    1381. 

  1381.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1382. 

  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. 

  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. 

  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. 

  1385.     testname='AddT_TI_T1' )
    1386. 

  1386.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1387. 

  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. 

  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. 

  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. 

  1390.     testname='AddT_TI_T2' )
    1391. 

  1391.   ! ESMF_TimeInterval = ESMF_TimeInterval + ESMF_TimeInterval
    1392. 

  1392.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1393. 

  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. 

  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. 

  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. 

  1396.     testname='AddTI_TI_TI1' )
    1397. 

  1397.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1398. 

  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. 

  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. 

  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. 

  1401.     testname='AddTI_TI_TI2' )
    1402. 

  1402.   CALL test_arithmetic( add_op=.TRUE.,                                             &
    1403. 

  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. 

  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. 

  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. 

  1406.     testname='AddTI_TI_TI3' )
    1407. 

  1407. 

  1408.   ! Subtraction tests
    1409. 

  1409.   ! ESMF_Time = ESMF_Time - ESMF_TimeInterval
    1410. 

  1410.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1411. 

  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. 

  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. 

  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. 

  1414.     testname='SubtractT_T_TI1' )
    1415. 

  1415.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1416. 

  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. 

  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. 

  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. 

  1419.     testname='SubtractT_T_TI2' )
    1420. 

  1420.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1421. 

  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. 

  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. 

  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. 

  1424.     testname='SubtractT_T_TI3' )
    1425. 

  1425.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1426. 

  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. 

  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. 

  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. 

  1429.     testname='SubtractT_T_TI4' )
    1430. 

  1430.   CALL test_arithmetic( add_op=.FALSE.,                                             &
    1431. 

  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. 

  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. 

  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. 

  1434.     testname='SubtractT_T_TI5' )
    1435. 

  1435.   CALL test_arithmetic( add_op=.FALSE.,                                             &
    1436. 

  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. 

  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. 

  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. 

  1439.     testname='SubtractT_T_TI6' )
    1440. 

  1440.   ! ESMF_Time = ESMF_Time - ESMF_TimeInterval with fractions
    1441. 

  1441.   CALL test_arithmetic( add_op=.FALSE.,                                             &
    1442. 

  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. 

  1443.      op1_t_sn=00,  op1_t_sd=00, &
    1444. 

  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. 

  1445.     op2_ti_sn=01, op2_ti_sd=03, &
    1446. 

  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. 

  1447.      res_t_sn=02,  res_t_sd=03, &
    1448. 

  1448.     testname='SubtractT_T_TI_F1' )
    1449. 

  1449.   ! ESMF_TimeInterval = ESMF_Time - ESMF_Time
    1450. 

  1450.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1451. 

  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. 

  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. 

  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. 

  1454.     testname='SubtractTI_T_T1' )
    1455. 

  1455.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1456. 

  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. 

  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. 

  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. 

  1459.     testname='SubtractTI_T_T2' )
    1460. 

  1460.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1461. 

  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. 

  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. 

  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. 

  1464.     testname='SubtractTI_T_T3' )
    1465. 

  1465.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1466. 

  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. 

  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. 

  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. 

  1469.     testname='SubtractTI_T_T4' )
    1470. 

  1470.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1471. 

  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. 

  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. 

  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. 

  1474.     testname='SubtractTI_T_T5' )
    1475. 

  1475.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1476. 

  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. 

  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. 

  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. 

  1479.     testname='SubtractTI_T_T6' )
    1480. 

  1480.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1481. 

  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. 

  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. 

  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. 

  1484.     testname='SubtractTI_T_T7' )
    1485. 

  1485.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1486. 

  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. 

  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. 

  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. 

  1489.     testname='SubtractTI_T_T8' )
    1490. 

  1490.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1491. 

  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. 

  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. 

  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. 

  1494.     testname='SubtractTI_T_T9' )
    1495. 

  1495.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1496. 

  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. 

  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. 

  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. 

  1499.     testname='SubtractTI_T_T10' )
    1500. 

  1500.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1501. 

  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. 

  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. 

  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. 

  1504.     testname='SubtractTI_T_T11' )
    1505. 

  1505.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1506. 

  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. 

  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. 

  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. 

  1509.     testname='SubtractTI_T_T12' )
    1510. 

  1510.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1511. 

  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. 

  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. 

  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. 

  1514.     testname='SubtractTI_T_T13' )
    1515. 

  1515.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1516. 

  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. 

  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. 

  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. 

  1519.     testname='SubtractTI_T_T14' )
    1520. 

  1520.   ! ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
    1521. 

  1521.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1522. 

  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. 

  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. 

  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. 

  1525.     testname='SubtractTI_TI_TI1' )
    1526. 

  1526.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1527. 

  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. 

  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. 

  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. 

  1530.     testname='SubtractTI_TI_TI2' )
    1531. 

  1531.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1532. 

  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. 

  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. 

  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. 

  1535.     testname='SubtractTI_TI_TI3' )
    1536. 

  1536.   ! Negative result ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
    1537. 

  1537.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1538. 

  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. 

  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. 

  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. 

  1541.     testname='SubtractTI_TI_TIN1' )
    1542. 

  1542.   CALL test_arithmetic( add_op=.FALSE.,                                            &
    1543. 

  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. 

  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. 

  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. 

  1546.     testname='SubtractTI_TI_TIN2' )
    1547. 

  1547. 

  1548.   ! Un-normalized ESMF_TimeInterval = ESMF_TimeInterval - ESMF_TimeInterval
    1549. 

  1549.   ! this is an error
    1550. 

  1550. !  CALL test_arithmetic( add_op=.FALSE.,                                            &
    1551. 

  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. 

  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. 

  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. 

  1554. !    testname='SubtractTI_TI_TIU1', expect_error=.TRUE. )
    1555. 

  1555. 

  1556.   ! this one should FAIL, and does
    1557. 

  1557. !  CALL test_arithmetic( add_op=.TRUE.,                                             &
    1558. 

  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. 

  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. 

  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. 

  1561. !    testname='AddTT1' )
    1562. 

  1562. 

  1563.   ! Multiplication tests
    1564. 

  1564.   ! ESMF_TimeInterval = ESMF_TimeInterval * INTEGER
    1565. 

  1565.   CALL test_arithmetic( multiply_op=.TRUE.,                &
    1566. 

  1566.     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
    1567. 

  1567.     op2_int=2,                                             &
    1568. 

  1568.     res_ti_dd=6,  res_ti_h=24, res_ti_m=37,  res_ti_s=06,  &
    1569. 

  1569.     testname='MultiplyTI_TI_INT1' )
    1570. 

  1570.   CALL test_arithmetic( multiply_op=.TRUE.,                &
    1571. 

  1571.     op1_ti_dd=350,  op1_ti_h=23,  op1_ti_m=50,  op1_ti_s=50, &
    1572. 

  1572.     op2_int=2,                                             &
    1573. 

  1573.     res_ti_dd=701,  res_ti_h=23, res_ti_m=41,  res_ti_s=40,&
    1574. 

  1574.     testname='MultiplyTI_TI_INT2' )
    1575. 

  1575.   CALL test_arithmetic( multiply_op=.TRUE.,                &
    1576. 

  1576.     op1_ti_s=01, op1_ti_sn=03, op1_ti_sd=04,               &
    1577. 

  1577.     op2_int=8,                                             &
    1578. 

  1578.     res_ti_s=14,                                           &
    1579. 

  1579.     testname='MultiplyTI_TI_INT3' )
    1580. 

  1580. 

  1581.   ! Division tests
    1582. 

  1582.   ! ESMF_TimeInterval = ESMF_TimeInterval / INTEGER
    1583. 

  1583.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1584. 

  1584.     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
    1585. 

  1585.     op2_int=3,                                             &
    1586. 

  1586.     res_ti_dd=1,  res_ti_h=04, res_ti_m=06,  res_ti_s=11,  &
    1587. 

  1587.     testname='DivideTI_TI_INT1' )
    1588. 

  1588.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1589. 

  1589.     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
    1590. 

  1590.     op2_int=4,                                             &
    1591. 

  1591.     res_ti_dd=0,  res_ti_h=21, res_ti_m=04,  res_ti_s=38,  &
    1592. 

  1592.     res_ti_sn=1,  res_ti_sd=4,                             &
    1593. 

  1593.     testname='DivideTI_TI_INT2' )
    1594. 

  1594.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1595. 

  1595.     op1_ti_s=01, op1_ti_sn=03, op1_ti_sd=04,               &
    1596. 

  1596.     op2_int=5,                                             &
    1597. 

  1597.     res_ti_s=0, res_ti_sn=7,  res_ti_sd=20,                &
    1598. 

  1598.     testname='DivideTI_TI_INT3' )
    1599. 

  1599.   ! INTEGER = ESMF_TimeInterval / ESMF_TimeInterval
    1600. 

  1600.   ! this operator truncates to whole integers
    1601. 

  1601.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1602. 

  1602.     op1_ti_dd=3,  op1_ti_h=12,  op1_ti_m=18,  op1_ti_s=33, &
    1603. 

  1603.     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
    1604. 

  1604.     res_int=1,                                             &
    1605. 

  1605.     testname='DivideINT_TI_TI1' )
    1606. 

  1606.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1607. 

  1607.     op1_ti_dd=6,  op1_ti_h=24,  op1_ti_m=36,  op1_ti_s=66, &
    1608. 

  1608.     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
    1609. 

  1609.     res_int=2,                                             &
    1610. 

  1610.     testname='DivideINT_TI_TI2' )
    1611. 

  1611.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1612. 

  1612.     op1_ti_dd=0,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
    1613. 

  1613.     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
    1614. 

  1614.     res_int=0,                                             &
    1615. 

  1615.     testname='DivideINT_TI_TI3' )
    1616. 

  1616.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1617. 

  1617.     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
    1618. 

  1618.     op2_ti_dd=0,  op2_ti_h=01,  op2_ti_m=00,  op2_ti_s=00, &
    1619. 

  1619.     res_int=24,                                            &
    1620. 

  1620.     testname='DivideINT_TI_TI4' )
    1621. 

  1621.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1622. 

  1622.     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
    1623. 

  1623.     op2_ti_dd=0,  op2_ti_h=00,  op2_ti_m=01,  op2_ti_s=00, &
    1624. 

  1624.     res_int=1440,                                          &
    1625. 

  1625.     testname='DivideINT_TI_TI5' )
    1626. 

  1626.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1627. 

  1627.     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=00, &
    1628. 

  1628.     op2_ti_dd=0,  op2_ti_h=00,  op2_ti_m=00,  op2_ti_s=01, &
    1629. 

  1629.     res_int=86400,                                         &
    1630. 

  1630.     testname='DivideINT_TI_TI6' )
    1631. 

  1631.   ! rounding
    1632. 

  1632.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1633. 

  1633.     op1_ti_dd=0,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=03, &
    1634. 

  1634.     op2_ti_dd=0,  op2_ti_h=00,  op2_ti_m=00,  op2_ti_s=02, &
    1635. 

  1635.     res_int=1,                                             &
    1636. 

  1636.     testname='DivideINT_TI_TIR1' )
    1637. 

  1637.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1638. 

  1638.     op1_ti_dd=1,  op1_ti_h=00,  op1_ti_m=00,  op1_ti_s=02, &
    1639. 

  1639.     op2_ti_dd=1,  op2_ti_h=00,  op2_ti_m=00,  op2_ti_s=03, &
    1640. 

  1640.     res_int=0,                                             &
    1641. 

  1641.     testname='DivideINT_TI_TIR2' )
    1642. 

  1642.   ! fractional operands
    1643. 

  1643.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1644. 

  1644.     op1_ti_m=00,  op1_ti_s=00, op1_ti_sn=03, op1_ti_sd=04, &
    1645. 

  1645.     op2_ti_m=00,  op2_ti_s=00, op2_ti_sn=03, op2_ti_sd=04, &
    1646. 

  1646.     res_int=1,                                             &
    1647. 

  1647.     testname='DivideINT_TI_TIF1' )
    1648. 

  1648.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1649. 

  1649.     op1_ti_m=00,  op1_ti_s=00, op1_ti_sn=06, op1_ti_sd=08, &
    1650. 

  1650.     op2_ti_m=00,  op2_ti_s=00, op2_ti_sn=03, op2_ti_sd=04, &
    1651. 

  1651.     res_int=1,                                             &
    1652. 

  1652.     testname='DivideINT_TI_TIF2' )
    1653. 

  1653.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1654. 

  1654.     op1_ti_m=00,  op1_ti_s=00, op1_ti_sn=03, op1_ti_sd=04, &
    1655. 

  1655.     op2_ti_m=00,  op2_ti_s=00, op2_ti_sn=04, op2_ti_sd=03, &
    1656. 

  1656.     res_int=0,                                             &
    1657. 

  1657.     testname='DivideINT_TI_TIF3' )
    1658. 

  1658.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1659. 

  1659.     op1_ti_m=00,  op1_ti_s=02, op1_ti_sn=03, op1_ti_sd=04, &
    1660. 

  1660.     op2_ti_m=00,  op2_ti_s=01, op2_ti_sn=01, op2_ti_sd=03, &
    1661. 

  1661.     res_int=2,                                             &
    1662. 

  1662.     testname='DivideINT_TI_TIF4' )
    1663. 

  1663.   ! negative operands
    1664. 

  1664.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1665. 

  1665.     op1_ti_dd=-6,  op1_ti_h=-24,  op1_ti_m=-36,  op1_ti_s=-66, &
    1666. 

  1666.     op2_ti_dd=3,  op2_ti_h=12,  op2_ti_m=18,  op2_ti_s=33, &
    1667. 

  1667.     res_int=-2,                                             &
    1668. 

  1668.     testname='DivideINT_TI_TIN1' )
    1669. 

  1669.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1670. 

  1670.     op1_ti_dd=6,  op1_ti_h=24,  op1_ti_m=36,  op1_ti_s=66, &
    1671. 

  1671.     op2_ti_dd=-3,  op2_ti_h=-12,  op2_ti_m=-18,  op2_ti_s=-33, &
    1672. 

  1672.     res_int=-2,                                             &
    1673. 

  1673.     testname='DivideINT_TI_TIN2' )
    1674. 

  1674.   CALL test_arithmetic( multiply_op=.FALSE.,               &
    1675. 

  1675.     op1_ti_dd=-6,  op1_ti_h=-24,  op1_ti_m=-36,  op1_ti_s=-66, &
    1676. 

  1676.     op2_ti_dd=-3,  op2_ti_h=-12,  op2_ti_m=-18,  op2_ti_s=-33, &
    1677. 

  1677.     res_int=2,                                             &
    1678. 

  1678.     testname='DivideINT_TI_TIN3' )
    1679. 

  1679. 

  1680. !$$$here...  modify these to add self-test PASS/FAIL output
    1681. 

  1681.   CALL test_clock_advance(                                                    &
    1682. 

  1682.     start_yy=2002, start_mm=12, start_dd=27, start_h=3, start_m=0, start_s=0, &
    1683. 

  1683.      stop_yy=2002,  stop_mm=12,  stop_dd=28,  stop_h=8,  stop_m=0,  stop_s=0, &
    1684. 

  1684.     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=600,                 &
    1685. 

  1685.     testname="SimpleClockAdvance" )
    1686. 

  1686. 

  1687.   CALL test_clock_advance(                                                    &
    1688. 

  1688.     start_yy=2003, start_mm=12, start_dd=29, start_h=9, start_m=0, start_s=0, &
    1689. 

  1689.      stop_yy=2004,  stop_mm=1,   stop_dd=2,   stop_h=9,  stop_m=0,  stop_s=0, &
    1690. 

  1690.     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=3600,                &
    1691. 

  1691.     testname="StdYearClockAdvance", increment_S=10 )
    1692. 

  1692. 

  1693.   CALL test_clock_advance(                                                    &
    1694. 

  1694.     start_yy=2004, start_mm=12, start_dd=29, start_h=9, start_m=0, start_s=0, &
    1695. 

  1695.      stop_yy=2005,  stop_mm=1,   stop_dd=2,   stop_h=9,  stop_m=0,  stop_s=0, &
    1696. 

  1696.     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=3600,                &
    1697. 

  1697.     testname="LeapYearClockAdvance", increment_S=10 )
    1698. 

  1698. 

  1699.   ! NRCM domain 3 case:  120 seconds / 9 
    1700. 

  1700.   ! 18 timesteps through end of leap year
    1701. 

  1701.   CALL test_clock_advance(                                                    &
    1702. 

  1702.     start_yy=2004, start_mm=12, start_dd=31, start_h=23, start_m=58, start_s=0,&
    1703. 

  1703.      stop_yy=2005,  stop_mm=1,   stop_dd=1,   stop_h=0,  stop_m=2,  stop_s=0, &
    1704. 

  1704.     timestep_d=0, timestep_h=0, timestep_m=0, timestep_s=13,                  &
    1705. 

  1705.     timestep_sn=1, timestep_sd=3,                                             &
    1706. 

  1706.     testname="LeapYearFractionClockAdvance",                                  &
    1707. 

  1707.     increment_S=1, increment_Sn=1, increment_Sd=3 )
    1708. 

  1708. 

  1709.   CALL ESMF_Finalize( rc=rc )
    1710. 

  1710.   CALL test_check_error( ESMF_SUCCESS, rc, &
    1711. 

  1711.                         'ESMF_Finalize() ', &
    1712. 

  1712.                         __FILE__ , &
    1713. 

  1713.                         __LINE__  )
    1714. 

  1714. 

  1715.   PRINT *,'END TEST SUITE'
    1716. 

  1716. 

  1717. END PROGRAM time_manager_test
    1718. 

  1718. 

  1719.