Exemplo n.º 1
0
 def modeltime(self):
     # build model time
     data_frame = {'perlen': self.dis.perlen.array,
                   'nstp': self.dis.nstp.array,
                   'tsmult': self.dis.tsmult.array}
     self._model_time = ModelTime(data_frame,
                                  self.dis.itmuni_dict[self.dis.itmuni],
                                  self.dis.start_datetime, self.dis.steady)
     return self._model_time
Exemplo n.º 2
0
Arquivo: mt.py Projeto: tdeyster/flopy
 def modeltime(self):
     # build model time
     data_frame = {
         "perlen": self.mf.dis.perlen.array,
         "nstp": self.mf.dis.nstp.array,
         "tsmult": self.mf.dis.tsmult.array,
     }
     self._model_time = ModelTime(
         data_frame,
         self.mf.dis.itmuni_dict[self.mf.dis.itmuni],
         self.dis.start_datetime,
         self.dis.steady.array,
     )
     return self._model_time
Exemplo n.º 3
0
 def modeltime(self):
     # build model time
     tdis = self.simulation.get_package('tdis')
     itmuni = tdis.time_units.get_data()
     start_date_time = tdis.start_date_time.get_data()
     if itmuni is None:
         itmuni = 0
     if start_date_time is None:
         start_date_time = '01-01-1970'
     period_data = tdis.perioddata.get_data()
     data_frame = {'perlen': period_data['perlen'],
                   'nstp': period_data['nstp'],
                   'tsmult': period_data['tsmult']}
     self._model_time = ModelTime(data_frame, itmuni, start_date_time)
     return self._model_time
Exemplo n.º 4
0
    def modeltime(self):
        tdis = self.simulation.get_package('tdis')
        period_data = tdis.perioddata.get_data()

        # build steady state data
        sto = self.get_package('sto')
        if sto is None:
            steady = np.full((len(period_data['perlen'])), True, dtype=bool)
        else:
            steady = np.full((len(period_data['perlen'])), False, dtype=bool)
            ss_periods = sto.steady_state.get_active_key_dict()
            tr_periods = sto.transient.get_active_key_dict()
            if ss_periods:
                last_ss_value = False
                # loop through steady state array
                for index, value in enumerate(steady):
                    # resolve if current index is steady state or transient
                    if index in ss_periods:
                        last_ss_value = True
                    elif index in tr_periods:
                        last_ss_value = False
                    if last_ss_value == True:
                        steady[index] = True

        # build model time
        itmuni = tdis.time_units.get_data()
        start_date_time = tdis.start_date_time.get_data()
        if itmuni is None:
            itmuni = 0
        if start_date_time is None:
            start_date_time = '01-01-1970'
        data_frame = {
            'perlen': period_data['perlen'],
            'nstp': period_data['nstp'],
            'tsmult': period_data['tsmult']
        }
        self._model_time = ModelTime(data_frame, itmuni, start_date_time,
                                     steady)
        return self._model_time