# -Function to return the julian day of the year
def julian(self):
    y = self.curdate.year
    start = self.datetime.datetime(y, 1, 1).toordinal()
    current = self.curdate.toordinal()
    day = current - start + 1
    return day, 1


# -Function to calculate the number of timesteps for the model run
def timesteps(self):
    nrTimeSteps = (self.enddate - self.startdate).days + 1
    print(
        "Running SPHY for "
        + str(self.startdate.day)
        + "-"
        + str(self.startdate.month)
        + "-"
        + str(self.startdate.year)
        + " through "
        + str(self.enddate.day)
        + "-"
        + str(self.enddate.month)
        + "-"
        + str(self.enddate.year)
    )
    print("with " + str(nrTimeSteps) + " daily timesteps")
    return nrTimeSteps