#NEWLINE# def __init__( self , model , isothm = 0 , ireact = 0 , igetsc = 1 , rhob = 1800.0 , prsity2 = 0.1 , srconc = 0.0 , sp1 = 0.0 , sp2 = 0.0 , rc1 = 0.0 , rc2 = 0.0 , extension = 'rct' , unitnumber = None ) : #NEWLINE# #INDENT# if ( unitnumber is None ) : #NEWLINE# #INDENT# unitnumber = self . unitnumber #NEWLINE# #UNINDENT# Package . __init__ ( self , model , extension , 'RCT' , unitnumber ) #NEWLINE# nrow = model . nrow #NEWLINE# ncol = model . ncol #NEWLINE# nlay = model . nlay #NEWLINE# self . heading1 = '# RCT for MT3DMS, generated by Flopy.' #NEWLINE# self . isothm = isothm #NEWLINE# self . ireact = ireact #NEWLINE# self . irctop = 2 #NEWLINE# self . igetsc = igetsc #NEWLINE# self . rhob = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , rhob , name = 'rhob' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . prsity2 = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , prsity2 , name = 'prsity2' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . srconc = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , srconc , name = 'srconc' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . sp1 = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , sp1 , name = 'sp1' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . sp2 = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , sp2 , name = 'sp2' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . rc1 = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , rc1 , name = 'rc1' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . rc2 = Util3d ( model , ( nlay , nrow , ncol ) , np . float32 , rc2 , name = 'rc2' , locat = self . unit_number [ 0 ] ) #NEWLINE# self . parent . add_package ( self ) #NEWLINE# return