flopy 3.2.1__zip → 3.2.2__zip

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfdis.py

@@ -11,7 +11,7 @@ MODFLOW Guide
 import sys
 import numpy as np
 from flopy.mbase import Package
-from flopy.utils import util_2d,util_3d
+from flopy.utils import util_2d,util_3d,reference
 
 class ModflowDis(Package):
     """
@@ -92,7 +92,7 @@ class ModflowDis(Package):
     def __init__(self, model, nlay=1, nrow=2, ncol=2, nper=1, delr=1.0,
                  delc=1.0, laycbd=0, top=1, botm=0, perlen=1, nstp=1,
                  tsmult=1, steady=True, itmuni=4, lenuni=2, extension='dis',
-                 unitnumber=11):
+                 unitnumber=11, xul=None, yul=None, rotation=0.0):
 
         # Call ancestor's init to set self.parent, extension, name and unit
         # number
@@ -102,7 +102,7 @@ class ModflowDis(Package):
         self.ncol = ncol
         self.nlay = nlay
         self.nper = nper
-        
+
         # Set values of all parameters
         self.heading = '# Discretization file for MODFLOW, generated by Flopy.'
         self.laycbd = util_2d(model, (self.nlay,), np.int, laycbd,
@@ -124,12 +124,15 @@ class ModflowDis(Package):
                               name='tsmult')
         self.steady = util_2d(model, (self.nper,), np.bool,
                               steady,name='steady')
-        self.itmuni = itmuni
-        self.lenuni = lenuni
+        self.itmuni = int(itmuni)
+        self.lenuni = int(lenuni)
         self.parent.add_package(self)
         self.itmuni_dict = {0: "undefined", 1: "seconds", 2: "minutes",
                             3: "hours", 4: "days", 5: "years"}
 
+        self.sr = reference.SpatialReference(self.delr.array, self.delc.array, self.lenuni, xul=xul,
+                                             yul=yul, rotation=rotation)
+
     def checklayerthickness(self):
         """
         Check layer thickness.
@@ -516,100 +519,14 @@ class ModflowDis(Package):
             tsmult.append(a3)
             steady.append(a4)
 
-        #--create dis object instance
+        # create dis object instance
         dis = ModflowDis(model, nlay, nrow, ncol, nper, delr, delc, laycbd, 
                          top, botm, perlen, nstp, tsmult, steady, itmuni,
                          lenuni)
-        #--return dis object instance
+        # return dis object instance
         return dis
 
 
-    def plot(self):
-        try:
-            import pylab as plt
-        except Exception as e:
-            print("error importing pylab: " + str(e))
-            return
-
-        #get the bas for ibound masking
-        bas = self.parent.bas6
-        if bas is not None:
-            ibnd = bas.getibound()
-        else:
-            ibnd = np.ones((self.nlay, self.nrow, self.ncol))
-
-        cmap = plt.cm.winter
-        cmap.set_bad('w', 1.0)
-        fs = 5
-
-        #the width and height of each subplot
-        delt = 2.0
-        shape = (2, self.nlay+1)
-        fig = plt.figure(figsize=(delt+(self.nlay*delt), delt * 2.0))
-        #fig = plt.figure()
-
-        #plot the time stepping info in the upper left corner
-        ax_time = plt.subplot2grid(shape, (0, 0))
-        ax_time.set_title("time stepping information", fontsize=fs)
-        idx = np.arange((self.nper))
-        width = 0.5
-        bars = ax_time.bar(idx, self.perlen, width=width, edgecolor="none",
-                    facecolor='b')
-        for s, b in zip(self.steady, bars):
-            if s:
-                b.set_color('b')
-            else:
-                b.set_color('c')
-
-        ax_time.set_xticks(idx+(width/2.0))
-        ax_time.set_xticklabels(idx, fontsize=fs, rotation=90)
-        ax_time.set_xlabel("stress period", fontsize=fs)
-        ax_time.set_ylabel("time (" + self.itmuni_dict[self.itmuni] + ")",
-                           fontsize=fs)
-        ax_time.set_yticklabels(ax_time.get_yticks(), fontsize=fs)
-
-        #plot model top in the lower left
-        ax_top = plt.subplot2grid(shape, (1, 0),aspect='equal')
-        top = self.top.array
-        top = np.ma.masked_where(ibnd[0] == 0, top)
-        ax_top.set_xlabel("column", fontsize=fs)
-        ax_top.set_ylabel("row", fontsize=fs)
-        ax_top.imshow(top, cmap=cmap, alpha=0.7, interpolation="none")
-        ax_top.set_title("model top - max,min: {0:G},{1:G}"
-                         .format(top.max(), top.min()), fontsize=fs)
-        ax_top.set_xticklabels(ax_top.get_xticks(), fontsize=fs)
-        ax_top.set_yticklabels(ax_top.get_yticks(), fontsize=fs)
-
-
-        botm = self.botm.array
-        for k in range(self.nlay):
-            ax_botm = plt.subplot2grid(shape, (0, k+1), aspect="equal")
-            ax_thk = plt.subplot2grid(shape, (1, k+1), aspect="equal")
-            ax_thk.set_xlabel("column", fontsize=fs)
-            ax_thk.set_xlabel("row", fontsize=fs)
-            #botm for this layer
-            b = botm[k]
-            b = np.ma.masked_where(ibnd[k] == 0, b)
-            ax_botm.imshow(b,cmap=cmap, alpha=0.7,
-                           interpolation="none")
-            ax_botm.set_title("botm of layer {0:d} - max,min : {1:G},{2:G}"
-                              .format(k+1, b.max(), b.min()), fontsize=fs)
-            #thickness of this layer
-            if k == 0:
-                t = top - botm[k]
-            else:
-                t = botm[k-1] - botm[k]
-            t = np.ma.masked_where(ibnd[k] == 0, t)
-            ax_thk.imshow(t,cmap=cmap, alpha=0.7,
-                           interpolation="none")
-            ax_thk.set_title("thickness of layer {0:d} - max,min : {1:G},{2:G}"
-                             .format(k+1, t.max(), t.min()), fontsize=fs)
-            ax_thk.set_yticklabels([])
-            ax_botm.set_yticklabels([])
-            ax_botm.set_xticklabels([])
-            ax_thk.set_xlabel("column", fontsize=fs)
-            ax_thk.set_xticklabels(ax_thk.get_xticks(), fontsize=fs)
-        plt.tight_layout()
-        plt.show()
+
 
 

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfdrn.py

@@ -108,7 +108,7 @@ class ModflowDrn(Package):
             self.dtype = dtype
         else:
             self.dtype = self.get_default_dtype(structured=self.parent.structured)
-        self.stress_period_data = mflist(model, self.dtype, stress_period_data)
+        self.stress_period_data = mflist(self, stress_period_data)
         self.parent.add_package(self)
 
     def __repr__(self):

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfevt.py

@@ -190,7 +190,7 @@ class ModflowEvt(Package):
         except:
             pass
 
-        #--dataset 3 and 4 - parameters data
+        # dataset 3 and 4 - parameters data
         pak_parms = None
         if npar > 0:
             pak_parms = mfparbc.loadarray(f, npar, model.verbose)
@@ -267,9 +267,9 @@ class ModflowEvt(Package):
                     current_ievt = t
                 ievt[iper] = current_ievt
 
-        #--create evt object
+        # create evt object
         evt = ModflowEvt(model, nevtop=nevtop, ipakcb=ipakcb, 
                          surf=surf, evtr=evtr, exdp=exdp, ievt=ievt)
 
-        #--return evt object
+        # return evt object
         return evt

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfghb.py

@@ -112,7 +112,7 @@ class ModflowGhb(Package):
             self.dtype = dtype
         else:
             self.dtype = self.get_default_dtype(structured=self.parent.structured)
-        self.stress_period_data = mflist(model, self.dtype, stress_period_data)
+        self.stress_period_data = mflist(self, stress_period_data)
 
     def __repr__(self):
         return 'GHB package class'
@@ -202,6 +202,6 @@ class ModflowGhb(Package):
         """
 
         if model.verbose:
-            sys.stdout.write('loading gmg package file...\n')
+            sys.stdout.write('loading ghb package file...\n')
 
         return Package.load(model, ModflowGhb, f, nper)

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfgmg.py

@@ -4,9 +4,10 @@ from flopy.mbase import Package
 class ModflowGmg(Package):
     '''GMG package
     Only programmed to work with the default values; may need work for other options'''
-    def __init__(self, model, mxiter=50, iiter=30, iadamp=0, \
-                 hclose=1e-5, rclose=1e-5, relax=1.0, ioutgmg=0, \
-                 ism=0, isc=0, damp=1.0,dup=0.75,dlow=0.01,\
+    def __init__(self, model, mxiter=50, iiter=30, iadamp=0, 
+                 hclose=1e-5, rclose=1e-5, relax=1.0, 
+                 ioutgmg=0, iunitmhc=0, 
+                 ism=0, isc=0, damp=1.0, dup=0.75, dlow=0.01,
                  chglimit=1.0,extension='gmg', unitnumber=27):
         Package.__init__(self, model, extension, 'GMG', unitnumber) # Call ancestor's init to set self.parent, extension, name and unit number
         self.heading = '# GMG for MODFLOW, generated by Flopy.'
@@ -24,7 +25,7 @@ class ModflowGmg(Package):
         self.chglimit = chglimit
         self.damp = damp
         self.ioutgmg = ioutgmg
-        self.iunitmhc = 0
+        self.iunitmhc = iunitmhc
         self.parent.add_package(self)
 
     def __repr__( self ):
@@ -34,17 +35,20 @@ class ModflowGmg(Package):
         # Open file for writing
         f_gmg = open(self.fn_path, 'w')
         f_gmg.write('%s\n' % self.heading)        
-        #--ds0
-        f_gmg.write('{0:9.3e} {1:9d} {2:9.3e} {3:9d}\n'\
-             .format(self.rclose,self.iiter,self.hclose,self.mxiter))        
-        #--ds1
-        f_gmg.write('{0:9.3e} {1:9d} {2:9d} {3:9d}\n'\
-             .format(self.damp,self.iadamp,self.ioutgmg,self.iunitmhc))
-        #--ds2
-        f_gmg.write('{0:9d} {1:9d} {2:9.3e} {3:9.3e} {4:9.3e}\n'\
-             .format(self.ism,self.isc,self.dup,self.dlow,self.chglimit))
-        #--ds3
-        f_gmg.write('{0:10.3e}\n'.format(self.relax))
+        # dataset 0
+        f_gmg.write('{} {} {} {}\n'\
+             .format(self.rclose, self.iiter, self.hclose, self.mxiter))        
+        # dataset 1
+        f_gmg.write('{} {} {} {}\n'\
+             .format(self.damp, self.iadamp, self.ioutgmg, self.iunitmhc))
+        # dataset 2
+        f_gmg.write('{} {} '.format(self.ism, self.isc))
+        if self.iadamp == 2:
+            f_gmg.write('{} {} {}'\
+                 .format(self.dup, self.dlow, self.chglimit))
+        f_gmg.write('\n')
+        # dataset 3
+        f_gmg.write('{}\n'.format(self.relax))
         f_gmg.close()
 
 
@@ -87,11 +91,48 @@ class ModflowGmg(Package):
             filename = f
             f = open(filename, 'r')
         #dataset 0 -- header
+        while True:
+            line = f.readline()
+            if line[0] != '#':
+                break
+        #dataset 0
+        t = line.strip().split()
+        rclose = float(t[0])
+        iiter = int(t[1])
+        hclose = float(t[2])
+        mxiter = int(t[3])
+        #dataset 1
+        line = f.readline()
+        t = line.strip().split()
+        damp = float(t[0])
+        iadamp = int(t[1])
+        ioutgmg = int(t[2])
+        try:
+            iunitmhc = int(t[3])            
+        except:
+            iunitmhc = 0
+        #dataset 2
+        line = f.readline()
+        t = line.strip().split()
+        ism = int(t[0])
+        isc = int(t[1])
+        dup, dlow, chglimit = 0.75, 0.01, 1.0
+        if iadamp == 2:
+            dup = float(t[2])
+            dlow = float(t[3])
+            chglimit = float(t[4])
+        #dataset 3
+        line = f.readline()
+        t = line.strip().split()
+        relax = float(t[0])
 
-        print('   Warning: load method not completed. default gmg object created.')
-
-        #--close the open file
+        # close the open file
         f.close()
 
-        gmg = ModflowGmg(model)
+        # create the gmg object
+        gmg = ModflowGmg(model, mxiter=mxiter, iiter=iiter, iadamp=iadamp,
+                         hclose=hclose, rclose=rclose, relax=relax,
+                         ioutgmg=ioutgmg, iunitmhc=iunitmhc,
+                         ism=ism, isc=isc, damp=damp,
+                         dup=dup, dlow=dlow, chglimit=chglimit)
         return gmg

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfhfb.py

@@ -239,7 +239,7 @@ class ModflowHfb(Package):
         nphfb = int(t[0])
         mxfb = int(t[1])
         nhfbnp = int(t[2])
-        #--check for no-print suppressor
+        # check for no-print suppressor
         options = []
         aux_names = []
         if len(t) > 2:
@@ -254,11 +254,11 @@ class ModflowHfb(Package):
                     aux_names.append(t[it + 1].lower())
                     it += 1
                 it += 1
-        #--data set 2 and 3
+        # data set 2 and 3
         if nphfb > 0:
             dt = ModflowHfb.get_empty(1).dtype
             pak_parms = mfparbc.load(f, nphfb, dt, model.verbose)
-        #--data set 4
+        # data set 4
         bnd_output = None
         if nhfbnp > 0:
             specified = ModflowHfb.get_empty(nhfbnp)
@@ -269,7 +269,7 @@ class ModflowHfb(Package):
                 t = line.strip().split()
                 specified[ibnd] = tuple(t[:len(specified.dtype.names)])
 
-            #--convert indices to zero-based
+            # convert indices to zero-based
             specified['k'] -= 1
             specified['irow1'] -= 1
             specified['icol1'] -= 1
@@ -295,7 +295,7 @@ class ModflowHfb(Package):
 
                 par_current = ModflowHfb.get_empty(par_dict['nlst'])
 
-                #--
+                #
                 if model.mfpar.pval is None:
                     parval = np.float(par_dict['parval'])
                 else:
@@ -304,11 +304,11 @@ class ModflowHfb(Package):
                     except:
                         parval = np.float(par_dict['parval'])
 
-                #--fill current parameter data (par_current)
+                # fill current parameter data (par_current)
                 for ibnd, t in enumerate(data_dict):
                     par_current[ibnd] = tuple(t[:len(par_current.dtype.names)])
 
-                #--convert indices to zero-based
+                # convert indices to zero-based
                 par_current['k'] -= 1
                 par_current['irow1'] -= 1
                 par_current['icol1'] -= 1

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mflpf.py

@@ -347,13 +347,13 @@ class ModflowLpf(Package):
             t = line.strip().split()
             wetfct, iwetit, ihdwet = float(t[0]), int(t[1]), int(t[2])
 
-        #--parameters data
+        # parameters data
         par_types = []
         if nplpf > 0:
             par_types, parm_dict = mfpar.load(f, nplpf, model.verbose)
             #print parm_dict
 
-        #--non-parameter data
+        # non-parameter data
         transient = not model.get_package('DIS').steady.all()
         hk = [0] * nlay
         hani = [0] * nlay
@@ -385,17 +385,20 @@ class ModflowLpf(Package):
             if model.verbose:
                 print('   loading vka layer {0:3d}...'.format(k+1))
             if 'vka' not in par_types and 'vani' not in par_types:
-                t = util_2d.load(f, model, (nrow,ncol), np.float32, 'vka',
+                t = util_2d.load(f, model, (nrow, ncol), np.float32, 'vka',
                                  ext_unit_dict)
             else:
                 line = f.readline()
-                t = mfpar.parameter_fill(model, (nrow, ncol), 'vka', parm_dict, findlayer=k)
+                key = 'vka'
+                if 'vani' in par_types:
+                    key = 'vani'
+                t = mfpar.parameter_fill(model, (nrow, ncol), key, parm_dict, findlayer=k)
             vka[k] = t
             if transient:
                 if model.verbose:
                     print('   loading ss layer {0:3d}...'.format(k+1))
                 if 'ss' not in par_types:
-                    t = util_2d.load(f, model, (nrow,ncol), np.float32, 'ss',
+                    t = util_2d.load(f, model, (nrow, ncol), np.float32, 'ss',
                                      ext_unit_dict)
                 else:
                     line = f.readline()
@@ -405,7 +408,7 @@ class ModflowLpf(Package):
                     if model.verbose:
                         print('   loading sy layer {0:3d}...'.format(k+1))
                     if 'sy' not in par_types:
-                        t = util_2d.load(f, model, (nrow,ncol), np.float32, 'sy',
+                        t = util_2d.load(f, model, (nrow, ncol), np.float32, 'sy',
                                          ext_unit_dict)
                     else:
                         line = f.readline()
@@ -416,7 +419,7 @@ class ModflowLpf(Package):
                 if model.verbose:
                     print('   loading vkcb layer {0:3d}...'.format(k+1))
                 if 'vkcb' not in par_types:
-                    t = util_2d.load(f, model, (nrow,ncol), np.float32, 'vkcb',
+                    t = util_2d.load(f, model, (nrow, ncol), np.float32, 'vkcb',
                                      ext_unit_dict)
                 else:
                     line = f.readline()
@@ -425,13 +428,13 @@ class ModflowLpf(Package):
             if (laywet[k] != 0 and laytyp[k] != 0):
                 if model.verbose:
                     print('   loading wetdry layer {0:3d}...'.format(k+1))
-                t = util_2d.load(f, model, (nrow,ncol), np.float32, 'wetdry',
+                t = util_2d.load(f, model, (nrow, ncol), np.float32, 'wetdry',
                                  ext_unit_dict)
                 wetdry[k] = t
 
-        #--create instance of lpf class
+        # create instance of lpf class
         lpf = ModflowLpf(model, ilpfcb=ilpfcb, laytyp=laytyp, layavg=layavg, chani=chani,
-                         laywet=laywet,hdry=hdry, iwdflg=iwetdry,
+                         layvka=layvka, laywet=laywet, hdry=hdry, iwdflg=iwetdry,
                          wetfct=wetfct, iwetit=iwetit, ihdwet=ihdwet,
                          hk=hk, hani=hani, vka=vka, ss=ss, sy=sy, vkcb=vkcb,
                          wetdry=wetdry, storagecoefficient=storagecoefficient,

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfmlt.py

@@ -154,11 +154,11 @@ class ModflowMlt(Package):
                 if 'function' in t[1].lower() or 'expression' in t[1].lower():
                     readArray = False
                     kwrd = t[1].lower()
-            #--load data
+            # load data
             if readArray:
                 t = util_2d.load(f, model, (nrow, ncol), np.float32, mltnam,
                                  ext_unit_dict)
-                #--add unit number to list of external files in ext_unit_dict to remove.
+                # add unit number to list of external files in ext_unit_dict to remove.
                 if t.locat is not None:
                     model.add_pop_key_list(t.locat)
             else:
@@ -167,7 +167,7 @@ class ModflowMlt(Package):
                 t = ModflowMlt.mult_function(mult_dict, line)
             mult_dict[mltnam] = t
 
-        #--create mlt dictionary
+        # create mlt dictionary
         mlt = ModflowMlt(model, mult_dict=mult_dict)
 
         return mlt

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfnwt.py

@@ -325,7 +325,7 @@ class ModflowNwt(Package):
 
         print('  ...load method not completed. default nwt file created.')
 
-        #--close the open file
+        # close the open file
         f.close()
 
         nwt = ModflowNwt(model)

{flopy-3.2.1 → flopy-3.2.2}/flopy/modflow/mfoc.py

@@ -203,7 +203,7 @@ class ModflowOc(Package):
         f_oc = open(self.fn_path, 'w')
         f_oc.write('{}\n'.format(self.heading))
 
-        #--write options
+        # write options
 
         f_oc.write('HEAD PRINT FORMAT {0:3.0f}\n'\
                    .format(self.ihedfm))            
@@ -229,7 +229,10 @@ class ModflowOc(Package):
                        .format(self.unit_number[4]))
 
         if self.compact:
-            f_oc.write('COMPACT BUDGET FILES')
+            f_oc.write('COMPACT BUDGET FILES\n')
+        
+        # add a line separator between header and stress
+        #  period data
         f_oc.write('\n')
 
 
@@ -263,7 +266,7 @@ class ModflowOc(Package):
                     f_oc.write('\n')
                     ddnref = ''
 
-        #--close oc file
+        # close oc file
         f_oc.close()
 
     @staticmethod
@@ -344,7 +347,7 @@ class ModflowOc(Package):
             filename = f
             f = open(filename, 'r')
 
-        #--read header
+        # read header
         ipos = f.tell()
         while True:
             line = f.readline()
@@ -361,10 +364,10 @@ class ModflowOc(Package):
                 except:
                     f.seek(ipos)
                     pass
-                #--exit so the remaining data can be read 
+                # exit so the remaining data can be read
                 #  from the file based on numericformat
                 break
-            #--set pointer to current position in the OC file    
+            # set pointer to current position in the OC file
             ipos = f.tell()
                  
         
@@ -377,16 +380,16 @@ class ModflowOc(Package):
                     lnlst = line.strip().split()
                     incode, ihddfl = int(lnlst[0]), int(lnlst[1])
                     ibudfl, icbcfl = int(lnlst[2]), int(lnlst[3])
-                    #--new print and save flags are needed if incode is not 
+                    # new print and save flags are needed if incode is not
                     #  less than 0.
                     if incode >= 0:
                         lines = []
-                    #--use print options from the last time step
+                    # use print options from the last time step
                     else:
                         if len(lines) > 0:
                             stress_period_data[(iperoc, itsoc)] = list(lines)
                         continue
-                    #--set print and save budget flags
+                    # set print and save budget flags
                     if ibudfl != 0:
                         lines.append('PRINT BUDGET')
                     if icbcfl != 0:
@@ -494,13 +497,13 @@ class ModflowOc(Package):
                 elif 'PERIOD' in lnlst[0].upper():
                     if len(lines) > 0:
                         if iperoc > 0:
-                            #--create period step tuple
+                            # create period step tuple
                             kperkstp = (iperoc-1, itsoc-1)
-                            #--save data
+                            # save data
                             stress_period_data[kperkstp] = lines
-                        #--reset lines
+                        # reset lines
                         lines = []
-                    #--turn off oc if required
+                    # turn off oc if required
                     if iperoc > 0:
                         if itsoc==nstp[iperoc-1]:
                             iperoc1 = iperoc + 1
@@ -510,10 +513,13 @@ class ModflowOc(Package):
                             itsoc1 = itsoc + 1
                     else:
                         iperoc1, itsoc1 = iperoc, itsoc 
-                    #--update iperoc and itsoc
+                    # update iperoc and itsoc
                     iperoc = int(lnlst[1])
                     itsoc = int(lnlst[3])
-                    #--add a empty list if necessary
+                    # do not used data that exceeds nper
+                    if iperoc > nper:
+                        break
+                    # add a empty list if necessary
                     iempty = False
                     if iperoc != iperoc1:
                         iempty = True
@@ -535,11 +541,11 @@ class ModflowOc(Package):
     
             #store the last record in word
             if len(lines) > 0:
-                #--create period step tuple
+                # create period step tuple
                 kperkstp = (iperoc-1, itsoc-1)
-                #--save data
+                # save data
                 stress_period_data[kperkstp] = lines
-                #--add a empty list if necessary
+                # add a empty list if necessary
                 iempty = False
                 if iperoc != iperoc1:
                     iempty = True
@@ -550,7 +556,7 @@ class ModflowOc(Package):
                     kperkstp = (iperoc1-1, itsoc1-1)
                     stress_period_data[kperkstp] = []
                     
-        #--reset unit numbers
+        # reset unit numbers
         unitnumber=[14, 51, 52, 53]
         if ihedun > 0:
             model.add_pop_key_list(ihedun)
@@ -562,7 +568,7 @@ class ModflowOc(Package):
                 cboufm = True
                 
 
-        #--create instance of oc class
+        # create instance of oc class
         oc = ModflowOc(model, ihedfm=ihedfm, iddnfm=iddnfm,
                  chedfm=chedfm, cddnfm=cddnfm, cboufm=cboufm,
                  compact=compact,