Converting a MODFLOW-2005 style SFR Package to MODFLOW 6

Sometimes you may want to convert a MODFLOW-2005 or MODFLOW-NWT model to MODFLOW 6, but the USGS mf5to6 converter fails on the SFR Package. This notebook demonstrates how to convert a MODFLOW-2005 (or MODFLOW-NWT) SFR Package input file to MODFLOW 6 using SFRmaker.

[1]:

import flopy
import sfrmaker

# Load the model using Flopy
model = flopy.modflow.Modflow.load('tf_with_sfr.nam', model_ws='../tylerforks/tylerforks')

#Check that the SFR package was loaded
model.get_package_list()

[1]:

['DIS', 'BAS6', 'UPW', 'RCH', 'OC', 'NWT', 'CHD', 'SFR']

[2]:

# Instantiate an `Mf6SFR` object (with the Flopy MODFLOW-2005 style SFR instance)
mf5to6 = sfrmaker.Mf6SFR(model.sfr)

# Write the MODFLOW 6 SFR Package file
mf5to6.write_file('tf_mf6.sfr')

converting reach and segment data to package data...
converting segment data to period data...
wrote tf_mf6.sfr