Changeset 802 for codes/icosagcm/devel/Python/dynamico/time_step.py
- Timestamp:
- 01/21/19 11:34:04 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/Python/dynamico/time_step.py
r778 r802 2 2 from scipy import sparse as sparse 3 3 import math as math 4 5 import getargs 6 log_master, log_world = getargs.getLogger(__name__) 7 INFO, DEBUG, ERROR = log_master.info, log_master.debug, log_world.error 8 INFO_ALL, DEBUG_ALL = log_world.info, log_world.debug 4 9 5 10 class Struct: pass … … 87 92 def __init__(self,precision,bwd_fast_slow, dt,nstage, bsjl,bfjl): 88 93 if bwd_fast_slow is None: 89 print 'SIRK with Fortran-side time stepping'94 DEBUG('SIRK with Fortran-side time stepping') 90 95 precision=np.float64 91 96 else: 92 print 'SIRK with Python-side time stepping at precision ', precision97 DEBUG('SIRK with Python-side time stepping at precision %s' % precision) 93 98 94 99 self.bwd_fast_slow = bwd_fast_slow … … 98 103 self.csjl, self.cfjl = [x.astype(precision) for x in csjl, cfjl] 99 104 self.tauj = np.array([dt*bfjl[j,j] for j in range(self.nstage)], dtype=precision) 100 print 'Types of csjl, cfjl, tauj :', self.csjl.dtype, self.cfjl.dtype, self.tauj.dtype105 DEBUG('Types of csjl, cfjl, tauj : %s %s %s' % (self.csjl.dtype, self.cfjl.dtype, self.tauj.dtype) ) 101 106 def next(self,zj): 102 107 csjl, cfjl, tauj = self.csjl, self.cfjl, self.tauj
Note: See TracChangeset
for help on using the changeset viewer.