Feature/atrayano/history start stop

CHANGELOG.md

@@ -9,6 +9,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ### Added
 
+- Added optional start_date and start_time to control the output window for each History collection. No output will be written before then. If not specified, these default to the beginning of the experiment.
+
 - Added loggers when writing or reading weight files
 - Added new option to AGCM.rc `overwrite_checkpoint` to allow checkpoint files to be overwritten. By default still will not overwrite checkpoints
 

gridcomps/History/MAPL_HistoryCollection.F90

@@ -47,12 +47,15 @@ module MAPL_HistoryCollectionMod
      integer                            :: acc_offset
      integer                            :: ref_date
      integer                            :: ref_time
+     integer                            :: start_date
+     integer                            :: start_time
      integer                            :: end_date
      integer                            :: end_time
      integer                            :: duration
      type(ESMF_Alarm)                   :: his_alarm ! when to write file
      type(ESMF_Alarm)                   :: seg_alarm ! segment alarm controls when to write to new file
      type(ESMF_Alarm)                   :: mon_alarm
+     type(ESMF_Alarm)                   :: start_alarm
      type(ESMF_Alarm)                   :: end_alarm
      integer,pointer                    :: expSTATE (:)
      integer                            :: unit
@@ -70,6 +73,7 @@ module MAPL_HistoryCollectionMod
      integer                            :: verbose
      integer                            :: xyoffset
      logical                            :: disabled
+     logical                            :: skipWriting
      logical                            :: subVm
      logical                            :: backwards ! Adds support for clock running in reverse direction
      logical                            :: useNewFormat

gridcomps/History/MAPL_HistoryGridComp.F90

@@ -811,6 +811,17 @@ subroutine Initialize ( gc, import, dumexport, clock, rc )
                                       label=trim(string) // 'ref_time:',_RC )
        _ASSERT(is_valid_time(list(n)%ref_time),'Invalid ref_time')
 
+       call ESMF_ConfigGetAttribute ( cfg, list(n)%start_date, default=-999, &
+                                      label=trim(string) // 'start_date:',_RC )
+       if (list(n)%start_date /= -999) then
+          _ASSERT(is_valid_date(list(n)%start_date),'Invalid start_date')
+       end if
+       call ESMF_ConfigGetAttribute ( cfg, list(n)%start_time, default=-999, &
+                                      label=trim(string) // 'start_time:',_RC )
+       if (list(n)%start_time /= -999) then
+          _ASSERT(is_valid_time(list(n)%start_time),'Invalid start_time')
+       end if
+
        call ESMF_ConfigGetAttribute ( cfg, list(n)%end_date, default=-999, &
                                       label=trim(string) // 'end_date:',_RC )
        if (list(n)%end_date /= -999) then
@@ -1349,6 +1360,38 @@ subroutine Initialize ( gc, import, dumexport, clock, rc )
           intState%stampOffset(n) = Frequency ! we go to the beginning of the month
        end if
 
+! End Alarm based on start_date and start_time
+! ----------------------------------------
+       if( list(n)%start_date.ne.-999 .and. list(n)%start_time.ne.-999 ) then
+          REF_TIME(1) =     list(n)%start_date/10000
+          REF_TIME(2) = mod(list(n)%start_date,10000)/100
+          REF_TIME(3) = mod(list(n)%start_date,100)
+          REF_TIME(4) =     list(n)%start_time/10000
+          REF_TIME(5) = mod(list(n)%start_time,10000)/100
+          REF_TIME(6) = mod(list(n)%start_time,100)
+
+          call ESMF_TimeSet( RingTime, YY = REF_TIME(1), &
+                                       MM = REF_TIME(2), &
+                                       DD = REF_TIME(3), &
+                                       H  = REF_TIME(4), &
+                                       M  = REF_TIME(5), &
+                                       S  = REF_TIME(6), calendar=cal, rc=rc )
+       else
+          RingTime = CurrTime
+       end if
+       list(n)%start_alarm = ESMF_AlarmCreate( clock=clock, RingTime=RingTime, sticky=.false., _RC )
+
+       list(n)%skipWriting = .true.
+       if (RingTime == CurrTime) then
+          call  ESMF_AlarmRingerOn(list(n)%start_alarm, _RC )
+          list(n)%skipWriting = .false.
+       else
+          if (RingTime < CurrTime .NEQV. list(n)%backwards) then
+             list(n)%skipWriting = .false.
+          endif
+       end if
+
+
 ! End Alarm based on end_date and end_time
 ! ----------------------------------------
        if( list(n)%end_date.ne.-999 .and. list(n)%end_time.ne.-999 ) then
@@ -3427,6 +3470,14 @@ subroutine Run ( gc, import, export, clock, rc )
 
   ! decide if we are writing based on alarms
 
+   do n=1,nlist
+      if (list(n)%skipWriting) then
+         if (ESMF_AlarmIsRinging(list(n)%start_alarm)) then
+            list(n)%skipWriting = .false.
+         endif
+      endif
+   end do
+
    do n=1,nlist
       if (list(n)%disabled .or. ESMF_AlarmIsRinging(list(n)%end_alarm) ) then
          list(n)%disabled = .true.
@@ -3454,6 +3505,8 @@ subroutine Run ( gc, import, export, clock, rc )
          end if
       end if
 
+      if (list(n)%skipWriting) writing(n) = .false.
+
        if (writing(n) .and. .not.IntState%average(n)) then
           ! R8 to R4 copy (if needed!)
           do m=1,list(n)%field_set%nfields