@@ -152,14 +152,14 @@ def add_gwf_model(sim, gwfname, newton=False):
152152            modelname = gwfname ,
153153            newtonoptions = "NEWTON" ,
154154            save_flows = True ,
155-             model_nam_file = "{ }.nam". format ( gwfname ) ,
155+             model_nam_file = f" { gwfname }  .nam" ,
156156        )
157157    else :
158158        gwf  =  flopy .mf6 .ModflowGwf (
159159            sim ,
160160            modelname = gwfname ,
161161            save_flows = True ,
162-             model_nam_file = "{ }.nam". format ( gwfname ) ,
162+             model_nam_file = f" { gwfname }  .nam" ,
163163        )
164164
165165    # Instantiating MODFLOW 6 solver for flow model 
@@ -176,7 +176,7 @@ def add_gwf_model(sim, gwfname, newton=False):
176176        scaling_method = "NONE" ,
177177        reordering_method = "NONE" ,
178178        relaxation_factor = relax ,
179-         filename = "{ }.ims". format ( gwfname ) ,
179+         filename = f" { gwfname }  .ims" ,
180180    )
181181    sim .register_ims_package (imsgwf , [gwfname ])
182182
@@ -193,7 +193,7 @@ def add_gwf_model(sim, gwfname, newton=False):
193193        botm = bot ,
194194        idomain = 1 ,
195195        pname = "DIS" ,
196-         filename = "{ }.dis". format ( gwfname ) ,
196+         filename = f" { gwfname }  .dis" ,
197197    )
198198
199199    # Instantiating MODFLOW 6 storage package 
@@ -205,7 +205,7 @@ def add_gwf_model(sim, gwfname, newton=False):
205205        steady_state = steady ,
206206        transient = transient ,
207207        pname = "STO" ,
208-         filename = "{ }.sto". format ( gwfname ) ,
208+         filename = f" { gwfname }  .sto" ,
209209    )
210210
211211    # Instantiating MODFLOW 6 node-property flow package 
@@ -217,22 +217,22 @@ def add_gwf_model(sim, gwfname, newton=False):
217217        k33 = k33 ,
218218        save_specific_discharge = True ,
219219        pname = "NPF" ,
220-         filename = "{ }.npf". format ( gwfname ) ,
220+         filename = f" { gwfname }  .npf" ,
221221    )
222222
223223    # Instantiating MODFLOW 6 initial conditions package for flow model 
224224    flopy .mf6 .ModflowGwfic (
225225        gwf ,
226226        strt = strthd ,
227-         filename = "{ }.ic". format ( gwfname ) ,
227+         filename = f" { gwfname }  .ic" ,
228228    )
229229
230230    # Instantiating MODFLOW 6 output control package for flow model 
231231    flopy .mf6 .ModflowGwfoc (
232232        gwf ,
233233        pname = "OC" ,
234-         head_filerecord = "{ }.hds". format ( gwfname ) ,
235-         budget_filerecord = "{ }.cbc". format ( gwfname ) ,
234+         head_filerecord = f" { gwfname }  .hds" ,
235+         budget_filerecord = f" { gwfname }  .cbc" ,
236236        headprintrecord = [("COLUMNS" , 10 , "WIDTH" , 15 , "DIGITS" , 6 , "GENERAL" )],
237237        saverecord = [("HEAD" , "ALL" ), ("BUDGET" , "ALL" )],
238238        printrecord = [("HEAD" , "ALL" ), ("BUDGET" , "ALL" )],
@@ -259,7 +259,7 @@ def add_gwf_model(sim, gwfname, newton=False):
259259        connectiondata = connectiondata ,
260260        perioddata = sfr_perioddata ,
261261        pname = pname ,
262-         filename = "{ }.sfr". format ( gwfname ) ,
262+         filename = f" { gwfname }  .sfr" ,
263263    )
264264    fname  =  f"{ gwfname }  .sfr.obs" 
265265    sfr_obs  =  {
@@ -274,9 +274,7 @@ def add_gwf_model(sim, gwfname, newton=False):
274274
275275
276276def  add_gwe_model (sim , gwename ):
277-     gwe  =  flopy .mf6 .ModflowGwe (
278-         sim , modelname = gwename , model_nam_file = "{}.nam" .format (gwename )
279-     )
277+     gwe  =  flopy .mf6 .ModflowGwe (sim , modelname = gwename , model_nam_file = f"{ gwename }  .nam" )
280278    gwe .name_file .save_flows  =  True 
281279
282280    imsgwe  =  flopy .mf6 .ModflowIms (
@@ -292,7 +290,7 @@ def add_gwe_model(sim, gwename):
292290        scaling_method = "NONE" ,
293291        reordering_method = "NONE" ,
294292        relaxation_factor = relax ,
295-         filename = "{ }.ims". format ( gwename ) ,
293+         filename = f" { gwename }  .ims" ,
296294    )
297295    sim .register_ims_package (imsgwe , [gwe .name ])
298296
@@ -309,18 +307,14 @@ def add_gwe_model(sim, gwename):
309307        botm = bot ,
310308        idomain = 1 ,
311309        pname = "DIS" ,
312-         filename = "{ }.dis". format ( gwename ) ,
310+         filename = f" { gwename }  .dis" ,
313311    )
314312
315313    # Instantiating MODFLOW 6 transport initial concentrations 
316-     flopy .mf6 .ModflowGweic (
317-         gwe , strt = strt_temp , pname = "IC" , filename = "{}.ic" .format (gwename )
318-     )
314+     flopy .mf6 .ModflowGweic (gwe , strt = strt_temp , pname = "IC" , filename = f"{ gwename }  .ic" )
319315
320316    # Instantiating MODFLOW 6 transport advection package 
321-     flopy .mf6 .ModflowGweadv (
322-         gwe , scheme = scheme , pname = "ADV" , filename = "{}.adv" .format (gwename )
323-     )
317+     flopy .mf6 .ModflowGweadv (gwe , scheme = scheme , pname = "ADV" , filename = f"{ gwename }  .adv" )
324318
325319    # Instantiating MODFLOW 6 transport dispersion package 
326320    flopy .mf6 .ModflowGwecnd (
@@ -331,10 +325,11 @@ def add_gwe_model(sim, gwename):
331325        ktw = ktw  *  86400 ,
332326        kts = kts  *  86400 ,
333327        pname = "CND" ,
334-         filename = "{ }.cnd". format ( gwename ) ,
328+         filename = f" { gwename }  .cnd" ,
335329    )
336330
337-     # Instantiating MODFLOW 6 transport mass storage package (formerly "reaction" package in MT3DMS) 
331+     # Instantiating MODFLOW 6 transport mass storage package 
332+     # (formerly "reaction" package in MT3DMS) 
338333    flopy .mf6 .ModflowGweest (
339334        gwe ,
340335        save_flows = True ,
@@ -345,7 +340,7 @@ def add_gwe_model(sim, gwename):
345340        heat_capacity_solid = cps ,
346341        density_solid = rhos ,
347342        pname = "EST" ,
348-         filename = "{ }.est". format ( gwename ) ,
343+         filename = f" { gwename }  .est" ,
349344    )
350345
351346    # Instantiating MODFLOW 6 source/sink mixing package for dealing with 
@@ -355,18 +350,16 @@ def add_gwe_model(sim, gwename):
355350        gwe ,
356351        sources = sourcerecarray ,
357352        pname = "SSM" ,
358-         filename = "{ }.ssm". format ( gwename ) ,
353+         filename = f" { gwename }  .ssm" ,
359354    )
360355
361356    # Instantiate MODFLOW 6 heat transport output control package 
362357    flopy .mf6 .ModflowGweoc (
363358        gwe ,
364359        pname = "OC" ,
365-         budget_filerecord = "{}.cbc" .format (gwename ),
366-         temperature_filerecord = "{}.ucn" .format (gwename ),
367-         temperatureprintrecord = [
368-             ("COLUMNS" , 10 , "WIDTH" , 15 , "DIGITS" , 6 , "GENERAL" )
369-         ],
360+         budget_filerecord = f"{ gwename }  .cbc" ,
361+         temperature_filerecord = f"{ gwename }  .ucn" ,
362+         temperatureprintrecord = [("COLUMNS" , 10 , "WIDTH" , 15 , "DIGITS" , 6 , "GENERAL" )],
370363        saverecord = [("TEMPERATURE" , "ALL" ), ("BUDGET" , "ALL" )],
371364        printrecord = [("TEMPERATURE" , "ALL" ), ("BUDGET" , "ALL" )],
372365    )
@@ -393,7 +386,7 @@ def add_gwe_model(sim, gwename):
393386        reachperioddata = sfe_perioddata ,
394387        flow_package_name = flwpckname ,
395388        pname = "SFE" ,
396-         filename = "{ }.sfe". format ( gwename ) ,
389+         filename = f" { gwename }  .sfe" ,
397390    )
398391
399392    return  sim 
@@ -404,7 +397,7 @@ def build_models(idx, test):
404397    ws  =  test .workspace 
405398    name  =  cases [idx ]
406399
407-     print ("Building MF6 model...{}"  . format ( name ) )
400+     print (f "Building MF6 model...{ name } "  )
408401
409402    # generate names for each model 
410403    gwfname1  =  "gwf-"  +  name  +  "nwt1" 
@@ -417,9 +410,7 @@ def build_models(idx, test):
417410    )
418411
419412    # Instantiating MODFLOW 6 time discretization 
420-     flopy .mf6 .ModflowTdis (
421-         sim , nper = nper , perioddata = tdis_rc , time_units = time_units 
422-     )
413+     flopy .mf6 .ModflowTdis (sim , nper = nper , perioddata = tdis_rc , time_units = time_units )
423414
424415    # Build two flow models, one with NWT, one without 
425416    sim  =  add_gwf_model (sim , gwfname1 , newton = True )
@@ -436,7 +427,7 @@ def build_models(idx, test):
436427        exgmnamea = gwfname1 ,
437428        exgmnameb = gwename1 ,
438429        pname = "GWFGWE1" ,
439-         filename = "{ }.gwfgwe1". format ( gwename1 ) ,
430+         filename = f" { gwename1 }  .gwfgwe1" ,
440431    )
441432
442433    flopy .mf6 .ModflowGwfgwe (
@@ -445,7 +436,7 @@ def build_models(idx, test):
445436        exgmnamea = gwfname2 ,
446437        exgmnameb = gwename2 ,
447438        pname = "GWFGWE2" ,
448-         filename = "{ }.gwfgwe2". format ( gwename2 ) ,
439+         filename = f" { gwename2 }  .gwfgwe2" ,
449440    )
450441
451442    return  sim , None 
0 commit comments