(ATCScript "DailyMonthlyFlowStats.spt"
  (Arg "DataFileObject")
  (Arg "Locations")
  (Arg "CurScenario")

  (LoadLibrary "ATCScriptGrid" "ATCScriptData" "ATCScriptMisc")

  (Unset Report) ('get rid off possible old grid)

  (Set CRLF (+ (Chr 13) (Chr 10)))
  (Set Space10 "          ")
  (Set Space20 "                    ")
  (Set Space38 "                                      ")
  (Set ReportString "")

  (If (= CurScenario "CurScenario")
    (Warn "CurScenario must be set before calling DailyMonthlyFlowStats script." "Flow Stats Report")
    (Else
      (Set LocCounter 0)
      (+= ReportString
        (+ "Daily and Monthly Flow Statistics for '" CurScenario "' scenario."  CRLF
          "   Run Made " (Format (FileAttribute (+ InputDir CurScenario ".hbn") "date") "yyyy-mm-dd hh:mm") CRLF 
          "   " gHspfRunInfo CRLF
          "   Simulation Period: " HspfnVals " " HspfTimeUnits 
                             " from " (Format HspfStartInterval "mm/dd/yy")
                             " to "   (Format HspfEndInterval   "mm/dd/yy") CRLF
          "   (Units:Inches)" CRLF CRLF
        )
      )
      (For CurLocation in Locations
        (Increment LocCounter)
        ('Get daily and monthly simulated flow)
        (Set SimDataSet        (DataFile DataFileObject Data Location CurLocation 
                                                             Scenario CurScenario
                                                             Constituent "Flow"))
        (Set SimDataSet        (DataTSer SimDataSet SubSetByDate SDate EDate))
        (Set DailySimDataSet   (DataTSer SimDataSet Aggregate 4 1 "Aver"))
        (Set MonthlySimDataSet (DataTSer SimDataSet Aggregate 5 1 "Aver"))
        ('Get daily and monthly observed flow)
        (Set ObsDataSet        (DataFile DataFileObject Data Location CurLocation 
                                                             Scenario "OBSERVED"
                                                             Constituent "Flow"))
        (Set ObsDataSet        (DataTSer ObsDataSet SubSetByDate SDate EDate))
        (Set DailyObsDataSet   (DataTSer ObsDataSet Aggregate 4 1 "Aver"))
        (Set MonthlyObsDataSet (DataTSer ObsDataSet Aggregate 5 1 "Aver"))

        (+= ReportString (+ CurLocation ":" (Format "Daily" Space10) CRLF))
        (+= ReportString (+ (Format "" Space38)
                            (Format "Simulated" Space10)
                            (Format "Observed" Space10)
                            CRLF
                         )
        )
        (+= ReportString (FlowStats DailySimDataSet DailyObsDataSet))
        (+= ReportString (+ CurLocation ":" (Format "Monthly" Space10) CRLF))
        (+= ReportString (+ (Format "" Space38)
                            (Format "Simulated" Space10)
                            (Format "Observed" Space10)
                            CRLF
                         )
        )
        (+= ReportString (FlowStats MonthlySimDataSet MonthlyObsDataSet))
      )
    )
  )
  (Return ReportString)

  (Unset ReportString)
  (Unset SimDataSet)
  (Unset ObsDataSet)
  (Unset MonthlySimDataSet)
  (Unset DailySimDataSet)
  (Unset MonthlyObsDataSet)
  (Unset DailyObsDataSet)
  (Unset LocCounter)

)