Hi,
When I load an existing modflow model with "Modflow.load", sf1 and sf2 in the bcf package are respectively specific yield and storage coefficient, the online manual says the opposite (https://flopy.readthedocs.io/en/latest/source/flopy.modflow.mfbcf.html).