<s> from __future__ import print_function <EOL> import os <EOL> import sys <EOL> import math <EOL> import numpy as np <EOL> import flopy . modflow as mf <EOL> import flopy . mt3d as mt3 <EOL> import flopy . seawat as swt <EOL> import flopy . utils as fu <EOL> import matplotlib . pyplot as plt <EOL> updates = { '<STR_LIT>' : [ '<STR_LIT>' , '<STR_LIT>' ] , <EOL> '<STR_LIT>' : '<STR_LIT>' , <EOL> '<STR_LIT>' : <NUM_LIT:0> , <EOL> '<STR_LIT>' : <NUM_LIT> , <EOL> '<STR_LIT>' : <NUM_LIT:7> , <EOL> '<STR_LIT>' : <NUM_LIT:8> , <EOL> '<STR_LIT>' : <NUM_LIT:7> , <EOL> '<STR_LIT>' : <NUM_LIT:7> } <EOL> plt . rcParams . update ( updates ) <EOL> cleanFiles = False <EOL> skipRuns = False <EOL> fext = '<STR_LIT>' <EOL> narg = len ( sys . argv ) <EOL> iarg = <NUM_LIT:0> <EOL> if narg > <NUM_LIT:1> : <EOL> while iarg < narg - <NUM_LIT:1> : <EOL> iarg += <NUM_LIT:1> <EOL> basearg = sys . argv [ iarg ] . lower ( ) <EOL> if basearg == '<STR_LIT>' : <EOL> cleanFiles = True <EOL> elif basearg == '<STR_LIT>' : <EOL> skipRuns = True <EOL> elif basearg == '<STR_LIT>' : <EOL> fext = '<STR_LIT>' <EOL> dirs = [ os . path . join ( '<STR_LIT>' ) , os . path . join ( '<STR_LIT>' ) ] <EOL> if cleanFiles : <EOL> print ( '<STR_LIT>' ) <EOL> print ( '<STR_LIT>' ) <EOL> file_dict = { } <EOL> file_dict [ '<STR_LIT:.>' ] = os . listdir ( '<STR_LIT:.>' ) <EOL> file_dict [ dirs [ <NUM_LIT:0> ] ] = os . listdir ( dirs [ <NUM_LIT:0> ] ) <EOL> file_dict [ dirs [ <NUM_LIT:1> ] ] = os . listdir ( dirs [ <NUM_LIT:1> ] ) <EOL> for key , files in list ( file_dict . items ( ) ) : <EOL> for f in files : <EOL> if os . path . isdir ( f ) : <EOL> continue <EOL> if '<STR_LIT>' != os . path . splitext ( f ) [ <NUM_LIT:1> ] . lower ( ) : <EOL> print ( '<STR_LIT>' . format ( os . path . basename ( f ) ) ) <EOL> os . remove ( os . path . join ( key , f ) ) <EOL> for d in dirs : <EOL> if os . path . exists ( d ) : <EOL> os . rmdir ( d ) <EOL> sys . exit ( <NUM_LIT:1> ) <EOL> for d in dirs : <EOL> if not os . path . exists ( d ) : <EOL> os . mkdir ( d ) <EOL> modelname = '<STR_LIT>' <EOL> mf_name = '<STR_LIT>' <EOL> nper = <NUM_LIT:1> <EOL> perlen = <NUM_LIT> <EOL> nstp = <NUM_LIT:1000> <EOL> nlay , nrow , ncol = <NUM_LIT:1> , <NUM_LIT:1> , <NUM_LIT> <EOL> delr = <NUM_LIT> <EOL> nsurf = <NUM_LIT:2> <EOL> x = np . arange ( <NUM_LIT:0.5> * delr , ncol * delr , delr ) <EOL> xedge = np . linspace ( <NUM_LIT:0> , float ( ncol ) * delr , len ( x ) + <NUM_LIT:1> ) <EOL> ibound = np . ones ( ( nrow , ncol ) , np . int ) <EOL> ibound [ <NUM_LIT:0> , <NUM_LIT:0> ] = - <NUM_LIT:1> <EOL> z0 = np . zeros ( ( nlay , nrow , ncol ) , np . float ) <EOL> z1 = np . zeros ( ( nlay , nrow , ncol ) , np . float ) <EOL> z0 [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT:30> : <NUM_LIT> ] = np . arange ( - <NUM_LIT> , - <NUM_LIT> , - <NUM_LIT:5> ) <EOL> z0 [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT> : ] = - <NUM_LIT> <EOL> z1 [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT> : <NUM_LIT:30> ] = np . arange ( - <NUM_LIT> , - <NUM_LIT> , - <NUM_LIT:5> ) <EOL> z1 [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT:30> : ] = - <NUM_LIT> <EOL> z = [ ] <EOL> z . append ( z0 ) <EOL> z . append ( z1 ) <EOL> ssz = <NUM_LIT> <EOL> isource = np . ones ( ( nrow , ncol ) , '<STR_LIT:int>' ) <EOL> isource [ <NUM_LIT:0> , <NUM_LIT:0> ] = <NUM_LIT:2> <EOL> modelname = '<STR_LIT>' <EOL> print ( '<STR_LIT>' , modelname ) <EOL> ml = mf . Modflow ( modelname , version = '<STR_LIT>' , exe_name = mf_name , model_ws = dirs [ <NUM_LIT:0> ] ) <EOL> discret = mf . ModflowDis ( ml , nlay = <NUM_LIT:1> , ncol = ncol , nrow = nrow , delr = delr , delc = <NUM_LIT:1> , top = <NUM_LIT:0> , botm = [ - <NUM_LIT> ] , <EOL> nper = nper , perlen = perlen , nstp = nstp ) <EOL> bas = mf . ModflowBas ( ml , ibound = ibound , strt = <NUM_LIT> ) <EOL> bcf = mf . ModflowBcf ( ml , laycon = <NUM_LIT:0> , tran = <NUM_LIT:2> * <NUM_LIT> ) <EOL> swi = mf . ModflowSwi2 ( ml , nsrf = nsurf , istrat = <NUM_LIT:1> , toeslope = <NUM_LIT> , tipslope = <NUM_LIT> , nu = [ <NUM_LIT:0> , <NUM_LIT> , <NUM_LIT> ] , <EOL> zeta = z , ssz = ssz , isource = isource , nsolver = <NUM_LIT:1> ) <EOL> oc = mf . ModflowOc88 ( ml , save_head_every = <NUM_LIT:1000> ) <EOL> pcg = mf . ModflowPcg ( ml ) <EOL> ml . write_input ( ) <EOL> if not skipRuns : <EOL> m = ml . run_model ( silent = False ) <EOL> zetafile = os . path . join ( dirs [ <NUM_LIT:0> ] , '<STR_LIT>' . format ( modelname ) ) <EOL> zobj = fu . CellBudgetFile ( zetafile ) <EOL> zkstpkper = zobj . get_kstpkper ( ) <EOL> zeta = zobj . get_data ( kstpkper = zkstpkper [ - <NUM_LIT:1> ] , text = '<STR_LIT>' ) [ <NUM_LIT:0> ] <EOL> zeta2 = zobj . get_data ( kstpkper = zkstpkper [ - <NUM_LIT:1> ] , text = '<STR_LIT>' ) [ <NUM_LIT:0> ] <EOL> modelname = '<STR_LIT>' <EOL> print ( '<STR_LIT>' , modelname ) <EOL> ml = mf . Modflow ( modelname , version = '<STR_LIT>' , exe_name = mf_name , model_ws = dirs [ <NUM_LIT:0> ] ) <EOL> discret = mf . ModflowDis ( ml , nlay = <NUM_LIT:1> , ncol = ncol , nrow = nrow , delr = delr , delc = <NUM_LIT:1> , top = <NUM_LIT:0> , botm = [ - <NUM_LIT> ] , <EOL> nper = nper , perlen = perlen , nstp = nstp ) <EOL> bas = mf . ModflowBas ( ml , ibound = ibound , strt = <NUM_LIT> ) <EOL> bcf = mf . ModflowBcf ( ml , laycon = <NUM_LIT:0> , tran = <NUM_LIT:2> * <NUM_LIT> ) <EOL> swi = mf . ModflowSwi2 ( ml , nsrf = nsurf , istrat = <NUM_LIT:0> , toeslope = <NUM_LIT> , tipslope = <NUM_LIT> , nu = [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT> , <NUM_LIT> ] , <EOL> zeta = z , ssz = ssz , isource = isource , nsolver = <NUM_LIT:1> ) <EOL> oc = mf . ModflowOc88 ( ml , save_head_every = <NUM_LIT:1000> ) <EOL> pcg = mf . ModflowPcg ( ml ) <EOL> ml . write_input ( ) <EOL> if not skipRuns : <EOL> m = ml . run_model ( silent = False ) <EOL> zetafile = os . path . join ( dirs [ <NUM_LIT:0> ] , '<STR_LIT>' . format ( modelname ) ) <EOL> zobj = fu . CellBudgetFile ( zetafile ) <EOL> zkstpkper = zobj . get_kstpkper ( ) <EOL> zetavd = zobj . get_data ( kstpkper = zkstpkper [ - <NUM_LIT:1> ] , text = '<STR_LIT>' ) [ <NUM_LIT:0> ] <EOL> zetavd2 = zobj . get_data ( kstpkper = zkstpkper [ - <NUM_LIT:1> ] , text = '<STR_LIT>' ) [ <NUM_LIT:0> ] <EOL> swtexe_name = '<STR_LIT>' <EOL> modelname = '<STR_LIT>' <EOL> print ( '<STR_LIT>' , modelname ) <EOL> swt_xmax = <NUM_LIT> <EOL> swt_zmax = <NUM_LIT> <EOL> swt_delr = <NUM_LIT:1.0> <EOL> swt_delc = <NUM_LIT:1.0> <EOL> swt_delz = <NUM_LIT:0.5> <EOL> swt_ncol = int ( swt_xmax / swt_delr ) <EOL> swt_nrow = <NUM_LIT:1> <EOL> swt_nlay = int ( swt_zmax / swt_delz ) <EOL> print ( swt_nlay , swt_nrow , swt_ncol ) <EOL> swt_ibound = np . ones ( ( swt_nlay , swt_nrow , swt_ncol ) , np . int ) <EOL> swt_ibound [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT:0> ] = - <NUM_LIT:1> <EOL> swt_x = np . arange ( <NUM_LIT:0.5> * swt_delr , swt_ncol * swt_delr , swt_delr ) <EOL> swt_xedge = np . linspace ( <NUM_LIT:0> , float ( ncol ) * delr , len ( swt_x ) + <NUM_LIT:1> ) <EOL> swt_top = <NUM_LIT:0.> <EOL> z0 = swt_top <EOL> swt_botm = np . zeros ( ( swt_nlay ) , np . float ) <EOL> swt_z = np . zeros ( ( swt_nlay ) , np . float ) <EOL> zcell = - swt_delz / <NUM_LIT> <EOL> for ilay in range ( <NUM_LIT:0> , swt_nlay ) : <EOL> z0 -= swt_delz <EOL> swt_botm [ ilay ] = z0 <EOL> swt_z [ ilay ] = zcell <EOL> zcell -= swt_delz <EOL> swt_X , swt_Z = np . meshgrid ( swt_x , swt_z ) <EOL> icbund = np . ones ( ( swt_nlay , swt_nrow , swt_ncol ) , np . int ) <EOL> sconc = np . ones ( ( swt_nlay , swt_nrow , swt_ncol ) , np . float ) <EOL> sconcp = np . zeros ( ( swt_nlay , swt_ncol ) , np . float ) <EOL> xsb = <NUM_LIT> <EOL> xbf = <NUM_LIT> <EOL> for ilay in range ( <NUM_LIT:0> , swt_nlay ) : <EOL> for icol in range ( <NUM_LIT:0> , swt_ncol ) : <EOL> if swt_x [ icol ] > xsb : <EOL> sconc [ ilay , <NUM_LIT:0> , icol ] = <NUM_LIT:0.5> <EOL> if swt_x [ icol ] > xbf : <EOL> sconc [ ilay , <NUM_LIT:0> , icol ] = <NUM_LIT:0.0> <EOL> for icol in range ( <NUM_LIT:0> , swt_ncol ) : <EOL> sconcp [ ilay , icol ] = sconc [ ilay , <NUM_LIT:0> , icol ] <EOL> xsb += swt_delz <EOL> xbf += swt_delz <EOL> itype = mt3 . Mt3dSsm . itype_dict ( ) <EOL> ssm_data = { <NUM_LIT:0> : [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT> , itype [ '<STR_LIT>' ] ] } <EOL> timprs = ( np . arange ( <NUM_LIT:5> ) + <NUM_LIT:1> ) * <NUM_LIT> <EOL> nprs = len ( timprs ) <EOL> ml = [ ] <EOL> ml = mf . Modflow ( modelname , version = '<STR_LIT>' , exe_name = swtexe_name , model_ws = dirs [ <NUM_LIT:1> ] ) <EOL> discret = mf . ModflowDis ( ml , nrow = swt_nrow , ncol = swt_ncol , nlay = swt_nlay , <EOL> delr = swt_delr , delc = swt_delc , laycbd = <NUM_LIT:0> , top = swt_top , botm = swt_botm , <EOL> nper = nper , perlen = perlen , nstp = <NUM_LIT:1> , steady = False ) <EOL> bas = mf . ModflowBas ( ml , ibound = swt_ibound , strt = <NUM_LIT> ) <EOL> lpf = mf . ModflowLpf ( ml , hk = <NUM_LIT> , vka = <NUM_LIT> , ss = <NUM_LIT:0.0> , sy = <NUM_LIT:0.0> , laytyp = <NUM_LIT:0> , layavg = <NUM_LIT:0> ) <EOL> oc = mf . ModflowOc88 ( ml , save_head_every = <NUM_LIT:1> , item2 = [ [ <NUM_LIT:0> , <NUM_LIT:1> , <NUM_LIT:0> , <NUM_LIT:0> ] ] ) <EOL> pcg = mf . ModflowPcg ( ml ) <EOL> ml . write_input ( ) <EOL> mt = mt3 . Mt3dms ( modelname , '<STR_LIT>' , modflowmodel = ml , <EOL> model_ws = dirs [ <NUM_LIT:1> ] ) <EOL> adv = mt3 . Mt3dAdv ( mt , mixelm = - <NUM_LIT:1> , <EOL> percel = <NUM_LIT> , <EOL> nadvfd = <NUM_LIT:0> , <EOL> nplane = <NUM_LIT:4> , <EOL> mxpart = <NUM_LIT> , <EOL> itrack = <NUM_LIT:2> , <EOL> dceps = <NUM_LIT> , <EOL> npl = <NUM_LIT:16> , <EOL> nph = <NUM_LIT:16> , <EOL> npmin = <NUM_LIT:8> , <EOL> npmax = <NUM_LIT> ) <EOL> btn = mt3 . Mt3dBtn ( mt , icbund = <NUM_LIT:1> , prsity = ssz , sconc = sconc , ifmtcn = - <NUM_LIT:1> , <EOL> chkmas = False , nprobs = <NUM_LIT:10> , nprmas = <NUM_LIT:10> , dt0 = <NUM_LIT:0.0> , ttsmult = <NUM_LIT> , ttsmax = <NUM_LIT> , <EOL> ncomp = <NUM_LIT:1> , nprs = nprs , timprs = timprs , mxstrn = <NUM_LIT> ) <EOL> dsp = mt3 . Mt3dDsp ( mt , al = <NUM_LIT:0.> , trpt = <NUM_LIT:1.> , trpv = <NUM_LIT:1.> , dmcoef = <NUM_LIT:0.> ) <EOL> gcg = mt3 . Mt3dGcg ( mt , mxiter = <NUM_LIT:1> , iter1 = <NUM_LIT:50> , isolve = <NUM_LIT:3> , cclose = <NUM_LIT> , iprgcg = <NUM_LIT:5> ) <EOL> ssm = mt3 . Mt3dSsm ( mt , stress_period_data = ssm_data ) <EOL> mt . write_input ( ) <EOL> mswtf = swt . Seawat ( modelname , '<STR_LIT>' , modflowmodel = ml , mt3dmsmodel = mt , <EOL> exe_name = swtexe_name , model_ws = dirs [ <NUM_LIT:1> ] ) <EOL> vdf = swt . SeawatVdf ( mswtf , nswtcpl = <NUM_LIT:1> , iwtable = <NUM_LIT:0> , densemin = <NUM_LIT:0> , densemax = <NUM_LIT:0> , denseref = <NUM_LIT> , denseslp = <NUM_LIT> , firstdt = <NUM_LIT> ) <EOL> mswtf . write_input ( ) <EOL> if not skipRuns : <EOL> m = mswtf . run_model ( silent = False ) <EOL> ucnfile = os . path . join ( dirs [ <NUM_LIT:1> ] , '<STR_LIT>' ) <EOL> uobj = fu . UcnFile ( ucnfile ) <EOL> times = uobj . get_times ( ) <EOL> print ( times ) <EOL> ukstpkper = uobj . get_kstpkper ( ) <EOL> print ( ukstpkper ) <EOL> c = uobj . get_data ( totim = times [ - <NUM_LIT:1> ] ) <EOL> conc = np . zeros ( ( swt_nlay , swt_ncol ) , np . float ) <EOL> for icol in range ( <NUM_LIT:0> , swt_ncol ) : <EOL> for ilay in range ( <NUM_LIT:0> , swt_nlay ) : <EOL> conc [ ilay , icol ] = c [ ilay , <NUM_LIT:0> , icol ] <EOL> fwid = <NUM_LIT> <EOL> fhgt = <NUM_LIT> <EOL> flft = <NUM_LIT> <EOL> frgt = <NUM_LIT> <EOL> fbot = <NUM_LIT> <EOL> ftop = <NUM_LIT> <EOL> print ( '<STR_LIT>' ) <EOL> xsf = plt . figure ( figsize = ( fwid , fhgt ) , facecolor = '<STR_LIT:w>' ) <EOL> xsf . subplots_adjust ( wspace = <NUM_LIT> , hspace = <NUM_LIT> , left = flft , right = frgt , bottom = fbot , top = ftop ) <EOL> ax = xsf . add_subplot ( <NUM_LIT:3> , <NUM_LIT:1> , <NUM_LIT:1> ) <EOL> ax . text ( - <NUM_LIT> , <NUM_LIT> , '<STR_LIT:A>' , transform = ax . transAxes , va = '<STR_LIT>' , ha = '<STR_LIT>' , size = '<STR_LIT>' ) <EOL> ax . plot ( [ <NUM_LIT> , <NUM_LIT> ] , [ <NUM_LIT:0> , - <NUM_LIT> ] , '<STR_LIT:k>' ) <EOL> ax . plot ( [ <NUM_LIT> , <NUM_LIT> ] , [ <NUM_LIT:0> , - <NUM_LIT> ] , '<STR_LIT:k>' ) <EOL> ax . set_xlim ( <NUM_LIT:0> , <NUM_LIT> ) <EOL> ax . set_ylim ( - <NUM_LIT> , <NUM_LIT:0> ) <EOL> ax . set_yticks ( np . arange ( - <NUM_LIT> , <NUM_LIT:1> , <NUM_LIT:10> ) ) <EOL> ax . text ( <NUM_LIT:50> , - <NUM_LIT:20> , '<STR_LIT>' , va = '<STR_LIT>' , ha = '<STR_LIT>' ) <EOL> ax . text ( <NUM_LIT> , - <NUM_LIT:20> , '<STR_LIT>' , va = '<STR_LIT>' , ha = '<STR_LIT>' ) <EOL> ax . text ( <NUM_LIT> , - <NUM_LIT:20> , '<STR_LIT>' , va = '<STR_LIT>' , ha = '<STR_LIT>' ) <EOL> ax . set_ylabel ( '<STR_LIT>' ) <EOL> ax = xsf . add_subplot ( <NUM_LIT:3> , <NUM_LIT:1> , <NUM_LIT:2> ) <EOL> ax . text ( - <NUM_LIT> , <NUM_LIT> , '<STR_LIT:B>' , transform = ax . transAxes , va = '<STR_LIT>' , ha = '<STR_LIT>' , size = '<STR_LIT>' ) <EOL> zp = zeta [ <NUM_LIT:0> , <NUM_LIT:0> , : ] <EOL> p = ( zp < <NUM_LIT:0.0> ) & ( zp > - <NUM_LIT> ) <EOL> ax . plot ( x [ p ] , zp [ p ] , '<STR_LIT:b>' , linewidth = <NUM_LIT> , drawstyle = '<STR_LIT>' ) <EOL> zp = zeta2 [ <NUM_LIT:0> , <NUM_LIT:0> , : ] <EOL> p = ( zp < <NUM_LIT:0.0> ) & ( zp > - <NUM_LIT> ) <EOL> ax . plot ( x [ p ] , zp [ p ] , '<STR_LIT:b>' , linewidth = <NUM_LIT> , drawstyle = '<STR_LIT>' ) <EOL> cc = ax . contour ( swt_X , swt_Z , conc , levels = [ <NUM_LIT> , <NUM_LIT> ] , colors = '<STR_LIT:k>' , linestyles = '<STR_LIT>' , linewidths = <NUM_LIT> , zorder = <NUM_LIT> ) <EOL> ax . plot ( [ - <NUM_LIT> , - <NUM_LIT:100> ] , [ - <NUM_LIT> , - <NUM_LIT:100> ] , '<STR_LIT:b>' , linewidth = <NUM_LIT> , label = '<STR_LIT>' ) <EOL> ax . plot ( [ - <NUM_LIT> , - <NUM_LIT:100> ] , [ - <NUM_LIT> , - <NUM_LIT:100> ] , '<STR_LIT:k>' , linewidth = <NUM_LIT> , label = '<STR_LIT>' ) <EOL> leg = ax . legend ( loc = '<STR_LIT>' , numpoints = <NUM_LIT:1> ) <EOL> leg . _drawFrame = False <EOL> ax . set_xlim ( <NUM_LIT:0> , <NUM_LIT> ) <EOL> ax . set_ylim ( - <NUM_LIT> , <NUM_LIT:0> ) <EOL> ax . set_yticks ( np . arange ( - <NUM_LIT> , <NUM_LIT:1> , <NUM_LIT:10> ) ) <EOL> ax . set_ylabel ( '<STR_LIT>' ) <EOL> ax = xsf . add_subplot ( <NUM_LIT:3> , <NUM_LIT:1> , <NUM_LIT:3> ) <EOL> ax . text ( - <NUM_LIT> , <NUM_LIT> , '<STR_LIT:C>' , transform = ax . transAxes , va = '<STR_LIT>' , ha = '<STR_LIT>' , size = '<STR_LIT>' ) <EOL> dr = zeta [ <NUM_LIT:0> , <NUM_LIT:0> , : ] <EOL> ax . plot ( x , dr , '<STR_LIT:b>' , linewidth = <NUM_LIT> , drawstyle = '<STR_LIT>' ) <EOL> dr = zeta2 [ <NUM_LIT:0> , <NUM_LIT:0> , : ] <EOL> ax . plot ( x , dr , '<STR_LIT:b>' , linewidth = <NUM_LIT> , drawstyle = '<STR_LIT>' ) <EOL> dr = zetavd [ <NUM_LIT:0> , <NUM_LIT:0> , : ] <EOL> ax . plot ( x , dr , '<STR_LIT:r>' , linewidth = <NUM_LIT> , drawstyle = '<STR_LIT>' ) <EOL> dr = zetavd2 [ <NUM_LIT:0> , <NUM_LIT:0> , : ] <EOL> ax . plot ( x , dr , '<STR_LIT:r>' , linewidth = <NUM_LIT> , drawstyle = '<STR_LIT>' ) <EOL> ax . plot ( [ - <NUM_LIT> , - <NUM_LIT:100> ] , [ - <NUM_LIT> , - <NUM_LIT:100> ] , '<STR_LIT:b>' , linewidth = <NUM_LIT> , label = '<STR_LIT>' ) <EOL> ax . plot ( [ - <NUM_LIT> , - <NUM_LIT:100> ] , [ - <NUM_LIT> , - <NUM_LIT:100> ] , '<STR_LIT:r>' , linewidth = <NUM_LIT> , label = '<STR_LIT>' ) <EOL> leg = ax . legend ( loc = '<STR_LIT>' , numpoints = <NUM_LIT:1> ) <EOL> leg . _drawFrame = False <EOL> ax . set_xlim ( <NUM_LIT:0> , <NUM_LIT> ) <EOL> ax . set_ylim ( - <NUM_LIT> , <NUM_LIT:0> ) <EOL> ax . set_yticks ( np . arange ( - <NUM_LIT> , <NUM_LIT:1> , <NUM_LIT:10> ) ) <EOL> ax . set_xlabel ( '<STR_LIT>' ) <EOL> ax . set_ylabel ( '<STR_LIT>' ) <EOL> outfig = '<STR_LIT>' . format ( fext ) <EOL> xsf . savefig ( outfig , dpi = <NUM_LIT> ) <EOL> print ( '<STR_LIT>' , outfig ) </s>
