FloPy

Demo of netCDF and shapefile export capabilities within the flopy export module.


In [1]:
from __future__ import print_function
import os
import flopy

Load our old friend...the Freyberg model


In [2]:
nam_file = "freyberg.nam"
model_ws = os.path.join("data")
ml = flopy.modflow.Modflow.load(nam_file,model_ws=model_ws,verbose=True)
print(ml.dis.sr)


Creating new model with name: freyberg
--------------------------------------------------

Parsing the namefile --> data/freyberg.nam
Setting filehandles:

--------------------------------------------------
External unit dictionary:
{19: filename:data/freyberg.rch, filetype:RCH, 2: filename:data/freyberg.list, filetype:LIST, 11: filename:data/freyberg.dis, filetype:DIS, 13: filename:data/freyberg.bas, filetype:BAS6, 14: filename:data/freyberg.oc, filetype:OC, 15: filename:data/freyberg.lpf, filetype:LPF, 18: filename:data/freyberg.riv, filetype:RIV, 51: filename:data/freyberg.hds, filetype:DATA(BINARY), 20: filename:data/freyberg.wel, filetype:WEL, 53: filename:data/freyberg.cbc, filetype:DATA(BINARY), 52: filename:data/freyberg.ddn, filetype:DATA(BINARY), 27: filename:data/freyberg.pcg, filetype:PCG}
--------------------------------------------------

ModflowBas6 free format:True

loading dis package file...
   Loading dis package with:
      1 layers, 40 rows, 20 columns, and 1 stress periods
   loading laycbd...
   loading delr...
   loading delc...
   loading top...
   loading botm...
   loading stress period data...
adding Package:  DIS
   DIS  package load...success
loading bas6 package file...
adding Package:  BAS6
   BAS6 package load...success
loading rch package file...
   loading rech stress period   1...
adding Package:  RCH
   RCH  package load...success
   LIST package load...skipped
loading oc package file...
adding Package:  OC
   OC   package load...success
loading lpf package file...
   loading IBCFCB, HDRY, NPLPF...
   loading LAYTYP...
   loading LAYAVG...
   loading CHANI...
   loading LAYVKA...
   loading LAYWET...
   loading hk layer   1...
   loading vka layer   1...
adding Package:  LPF
   LPF  package load...success
loading riv package file...
   loading <class 'flopy.modflow.mfriv.ModflowRiv'> for kper     1
adding Package:  RIV
   RIV  package load...success
   DATA(BINARY) file load...skipped
      freyberg.hds
loading wel package file...
   loading <class 'flopy.modflow.mfwel.ModflowWel'> for kper     1
adding Package:  WEL
   WEL  package load...success
   DATA(BINARY) file load...skipped
      freyberg.cbc
   DATA(BINARY) file load...skipped
      freyberg.ddn
loading pcg package file...
adding Package:  PCG
   PCG  package load...success


   The following 8 packages were successfully loaded.
      freyberg.dis
      freyberg.bas
      freyberg.rch
      freyberg.oc
      freyberg.lpf
      freyberg.riv
      freyberg.wel
      freyberg.pcg
   The following 1 packages were not loaded.
      freyberg.list



freyberg MODEL DATA VALIDATION SUMMARY:
  No errors or warnings encountered.

xul:0; yul:10000; rotation:0; proj4_str:+init=EPSG:4326; units:meters; lenuni:2

Some netCDF export capabilities:

Export a 2d array


In [3]:
# export a 2d array
ml.dis.top.export(os.path.join(model_ws, 'top.nc'))
ml.dis.top.export(os.path.join(model_ws, 'top.shp'))


2016-06-27 16:29:50.671378 WARNING: removing existing nc file: data/top.nc

2016-06-27 16:29:50.672511 WARNING: start datetime:1970-01-01T00:00:00Z

2016-06-27 16:29:50.672540 starting: initializing attributes

2016-06-27 16:29:50.672676 finished: initializing attributes, took: 0:00:00.000136

2016-06-27 16:29:50.672697 starting: initializing file

2016-06-27 16:29:50.672725 starting: initializing geometry

2016-06-27 16:29:50.674280 starting: building grid crs using proj4 string: +init=EPSG:4326

2016-06-27 16:29:50.674417 finished: building grid crs using proj4 string: +init=EPSG:4326, took: 0:00:00.000137

2016-06-27 16:29:50.674701 starting: projecting grid cell center arrays from +units=m +init=epsg:4326 to +units=m +init=epsg:4326 

2016-06-27 16:29:50.675098 starting: projecting grid cell center arrays from <pyproj.Proj object at 0x1117b45e8> to <pyproj.Proj object at 0x1117b4648>

2016-06-27 16:29:50.675297 finished: initializing geometry, took: 0:00:00.002572

2016-06-27 16:29:50.729406 starting: setting standard attributes

2016-06-27 16:29:50.731390 finished: setting standard attributes, took: 0:00:00.001984

2016-06-27 16:29:50.731412 starting: creating dimensions

2016-06-27 16:29:50.731985 finished: creating dimensions, took: 0:00:00.000573

2016-06-27 16:29:50.732006 starting: setting CRS info

2016-06-27 16:29:50.732381 finished: setting CRS info, took: 0:00:00.000375

2016-06-27 16:29:50.732413 starting: creating variable: time

2016-06-27 16:29:50.732525 finished: creating variable: time, took: 0:00:00.000112

2016-06-27 16:29:50.732824 WARNING: time_values:[ 10.]

2016-06-27 16:29:50.740023 starting: creating variable: elevation

2016-06-27 16:29:50.740257 finished: creating variable: elevation, took: 0:00:00.000234

2016-06-27 16:29:50.741414 starting: creating variable: longitude

2016-06-27 16:29:50.741615 finished: creating variable: longitude, took: 0:00:00.000201

2016-06-27 16:29:50.742423 starting: creating longitude var

2016-06-27 16:29:50.742460 starting: creating latitude var

2016-06-27 16:29:50.742485 starting: creating variable: latitude

2016-06-27 16:29:50.742577 finished: creating variable: latitude, took: 0:00:00.000092

2016-06-27 16:29:50.743116 starting: creating layer var

2016-06-27 16:29:50.743161 starting: creating variable: layer

2016-06-27 16:29:50.743448 finished: creating variable: layer, took: 0:00:00.000287

2016-06-27 16:29:50.744155 finished: creating layer var, took: 0:00:00.001039

2016-06-27 16:29:50.744193 starting: creating variable: delc

2016-06-27 16:29:50.744296 finished: creating variable: delc, took: 0:00:00.000103

2016-06-27 16:29:50.744704 starting: creating variable: delr

2016-06-27 16:29:50.744822 finished: creating variable: delr, took: 0:00:00.000118

2016-06-27 16:29:50.745368 finished: initializing file, took: 0:00:00.072671

2016-06-27 16:29:50.745576 starting: creating variable: model_top

2016-06-27 16:29:50.745682 finished: creating variable: model_top, took: 0:00:00.000106

Export a 3d array


In [4]:
#export a 3d array
ml.lpf.hk.export(os.path.join(model_ws, 'hk.nc'))
ml.lpf.hk.export(os.path.join(model_ws, 'hk.shp'))


2016-06-27 16:29:50.807230 WARNING: removing existing nc file: data/hk.nc

2016-06-27 16:29:50.808304 WARNING: start datetime:1970-01-01T00:00:00Z

2016-06-27 16:29:50.808345 starting: initializing attributes

2016-06-27 16:29:50.808550 finished: initializing attributes, took: 0:00:00.000205

2016-06-27 16:29:50.808571 starting: initializing file

2016-06-27 16:29:50.808600 starting: initializing geometry

2016-06-27 16:29:50.808632 starting: building grid crs using proj4 string: +init=EPSG:4326

2016-06-27 16:29:50.808672 finished: building grid crs using proj4 string: +init=EPSG:4326, took: 0:00:00.000040

2016-06-27 16:29:50.808789 starting: projecting grid cell center arrays from +units=m +init=epsg:4326 to +units=m +init=epsg:4326 

2016-06-27 16:29:50.808837 starting: projecting grid cell center arrays from <pyproj.Proj object at 0x114fd5ca8> to <pyproj.Proj object at 0x114fd5348>

2016-06-27 16:29:50.808882 finished: initializing geometry, took: 0:00:00.000282

2016-06-27 16:29:50.809271 starting: setting standard attributes

2016-06-27 16:29:50.809545 finished: setting standard attributes, took: 0:00:00.000274

2016-06-27 16:29:50.809560 starting: creating dimensions

2016-06-27 16:29:50.809666 finished: creating dimensions, took: 0:00:00.000106

2016-06-27 16:29:50.809679 starting: setting CRS info

2016-06-27 16:29:50.809766 finished: setting CRS info, took: 0:00:00.000087

2016-06-27 16:29:50.809786 starting: creating variable: time

2016-06-27 16:29:50.809881 finished: creating variable: time, took: 0:00:00.000095

2016-06-27 16:29:50.810143 WARNING: time_values:[ 10.]

2016-06-27 16:29:50.810859 starting: creating variable: elevation

2016-06-27 16:29:50.810993 finished: creating variable: elevation, took: 0:00:00.000134

2016-06-27 16:29:50.811522 starting: creating variable: longitude

2016-06-27 16:29:50.811634 finished: creating variable: longitude, took: 0:00:00.000112

2016-06-27 16:29:50.812519 starting: creating longitude var

2016-06-27 16:29:50.812581 starting: creating latitude var

2016-06-27 16:29:50.812630 starting: creating variable: latitude

2016-06-27 16:29:50.812790 finished: creating variable: latitude, took: 0:00:00.000160

2016-06-27 16:29:50.813461 starting: creating layer var

2016-06-27 16:29:50.813538 starting: creating variable: layer

2016-06-27 16:29:50.813889 finished: creating variable: layer, took: 0:00:00.000351

2016-06-27 16:29:50.814642 finished: creating layer var, took: 0:00:00.001181

2016-06-27 16:29:50.814686 starting: creating variable: delc

2016-06-27 16:29:50.814806 finished: creating variable: delc, took: 0:00:00.000120

2016-06-27 16:29:50.815262 starting: creating variable: delr

2016-06-27 16:29:50.815394 finished: creating variable: delr, took: 0:00:00.000132

2016-06-27 16:29:50.815804 finished: initializing file, took: 0:00:00.007233

2016-06-27 16:29:50.815962 starting: creating variable: hk

2016-06-27 16:29:50.816072 finished: creating variable: hk, took: 0:00:00.000110

Export a number of things to the same netCDF file


In [5]:
# export lots of things to the same nc file
fnc = ml.dis.botm.export(os.path.join(model_ws, 'test.nc'))
ml.lpf.hk.export(fnc)
ml.dis.top.export(fnc)

# export transient 2d
ml.rch.rech.export(fnc)


2016-06-27 16:29:50.868864 WARNING: removing existing nc file: data/test.nc

2016-06-27 16:29:50.869701 WARNING: start datetime:1970-01-01T00:00:00Z

2016-06-27 16:29:50.869741 starting: initializing attributes

2016-06-27 16:29:50.870039 finished: initializing attributes, took: 0:00:00.000298

2016-06-27 16:29:50.870089 starting: initializing file

2016-06-27 16:29:50.870128 starting: initializing geometry

2016-06-27 16:29:50.870165 starting: building grid crs using proj4 string: +init=EPSG:4326

2016-06-27 16:29:50.870205 finished: building grid crs using proj4 string: +init=EPSG:4326, took: 0:00:00.000040

2016-06-27 16:29:50.870321 starting: projecting grid cell center arrays from +units=m +init=epsg:4326 to +units=m +init=epsg:4326 

2016-06-27 16:29:50.870371 starting: projecting grid cell center arrays from <pyproj.Proj object at 0x114fd5d08> to <pyproj.Proj object at 0x114fd5348>

2016-06-27 16:29:50.870416 finished: initializing geometry, took: 0:00:00.000288

2016-06-27 16:29:50.870768 starting: setting standard attributes

2016-06-27 16:29:50.871082 finished: setting standard attributes, took: 0:00:00.000314

2016-06-27 16:29:50.871098 starting: creating dimensions

2016-06-27 16:29:50.871219 finished: creating dimensions, took: 0:00:00.000121

2016-06-27 16:29:50.871233 starting: setting CRS info

2016-06-27 16:29:50.871322 finished: setting CRS info, took: 0:00:00.000089

2016-06-27 16:29:50.871344 starting: creating variable: time

2016-06-27 16:29:50.871436 finished: creating variable: time, took: 0:00:00.000092

2016-06-27 16:29:50.871723 WARNING: time_values:[ 10.]

2016-06-27 16:29:50.872619 starting: creating variable: elevation

2016-06-27 16:29:50.872858 finished: creating variable: elevation, took: 0:00:00.000239

2016-06-27 16:29:50.873651 starting: creating variable: longitude

2016-06-27 16:29:50.873780 finished: creating variable: longitude, took: 0:00:00.000129

2016-06-27 16:29:50.874694 starting: creating longitude var

2016-06-27 16:29:50.874839 starting: creating latitude var

2016-06-27 16:29:50.874908 starting: creating variable: latitude

2016-06-27 16:29:50.875168 finished: creating variable: latitude, took: 0:00:00.000260

2016-06-27 16:29:50.876072 starting: creating layer var

2016-06-27 16:29:50.876138 starting: creating variable: layer

2016-06-27 16:29:50.876432 finished: creating variable: layer, took: 0:00:00.000294

2016-06-27 16:29:50.877481 finished: creating layer var, took: 0:00:00.001409

2016-06-27 16:29:50.877596 starting: creating variable: delc

2016-06-27 16:29:50.877838 finished: creating variable: delc, took: 0:00:00.000242

2016-06-27 16:29:50.878620 starting: creating variable: delr

2016-06-27 16:29:50.878749 finished: creating variable: delr, took: 0:00:00.000129

2016-06-27 16:29:50.879184 finished: initializing file, took: 0:00:00.009095

2016-06-27 16:29:50.879340 starting: creating variable: botm

2016-06-27 16:29:50.879444 finished: creating variable: botm, took: 0:00:00.000104

2016-06-27 16:29:50.880185 starting: creating variable: hk

2016-06-27 16:29:50.880318 finished: creating variable: hk, took: 0:00:00.000133

2016-06-27 16:29:50.881097 starting: creating variable: model_top

2016-06-27 16:29:50.881215 finished: creating variable: model_top, took: 0:00:00.000118

2016-06-27 16:29:50.881929 starting: creating variable: rech

2016-06-27 16:29:50.882038 finished: creating variable: rech, took: 0:00:00.000109

Out[5]:
<flopy.export.netcdf.NetCdf at 0x114fecb00>

Export whole packages to a netCDF file


In [6]:
# export mflist
fnc = ml.wel.export(os.path.join(model_ws, 'packages.nc'))
ml.lpf.export(fnc)
fnc.nc


2016-06-27 16:29:50.894565 WARNING: removing existing nc file: data/packages.nc

2016-06-27 16:29:50.895520 WARNING: start datetime:1970-01-01T00:00:00Z

2016-06-27 16:29:50.895545 starting: initializing attributes

2016-06-27 16:29:50.895666 finished: initializing attributes, took: 0:00:00.000121

2016-06-27 16:29:50.895687 starting: initializing file

2016-06-27 16:29:50.895717 starting: initializing geometry

2016-06-27 16:29:50.895754 starting: building grid crs using proj4 string: +init=EPSG:4326

2016-06-27 16:29:50.895802 finished: building grid crs using proj4 string: +init=EPSG:4326, took: 0:00:00.000048

2016-06-27 16:29:50.895970 starting: projecting grid cell center arrays from +units=m +init=epsg:4326 to +units=m +init=epsg:4326 

2016-06-27 16:29:50.896040 starting: projecting grid cell center arrays from <pyproj.Proj object at 0x114fd5348> to <pyproj.Proj object at 0x114fd5ee8>

2016-06-27 16:29:50.896093 finished: initializing geometry, took: 0:00:00.000376

2016-06-27 16:29:50.896623 starting: setting standard attributes

2016-06-27 16:29:50.897111 finished: setting standard attributes, took: 0:00:00.000488

2016-06-27 16:29:50.897130 starting: creating dimensions

2016-06-27 16:29:50.897260 finished: creating dimensions, took: 0:00:00.000130

2016-06-27 16:29:50.897274 starting: setting CRS info

2016-06-27 16:29:50.897410 finished: setting CRS info, took: 0:00:00.000136

2016-06-27 16:29:50.897460 starting: creating variable: time

2016-06-27 16:29:50.897679 finished: creating variable: time, took: 0:00:00.000219

2016-06-27 16:29:50.897942 WARNING: time_values:[ 10.]

2016-06-27 16:29:50.898701 starting: creating variable: elevation

2016-06-27 16:29:50.898847 finished: creating variable: elevation, took: 0:00:00.000146

2016-06-27 16:29:50.899406 starting: creating variable: longitude

2016-06-27 16:29:50.899517 finished: creating variable: longitude, took: 0:00:00.000111

2016-06-27 16:29:50.900055 starting: creating longitude var

2016-06-27 16:29:50.900090 starting: creating latitude var

2016-06-27 16:29:50.900114 starting: creating variable: latitude

2016-06-27 16:29:50.900205 finished: creating variable: latitude, took: 0:00:00.000091

2016-06-27 16:29:50.900627 starting: creating layer var

2016-06-27 16:29:50.900665 starting: creating variable: layer

2016-06-27 16:29:50.900855 finished: creating variable: layer, took: 0:00:00.000190

2016-06-27 16:29:50.901383 finished: creating layer var, took: 0:00:00.000756

2016-06-27 16:29:50.901417 starting: creating variable: delc

2016-06-27 16:29:50.901520 finished: creating variable: delc, took: 0:00:00.000103

2016-06-27 16:29:50.901923 starting: creating variable: delr

2016-06-27 16:29:50.902026 finished: creating variable: delr, took: 0:00:00.000103

2016-06-27 16:29:50.902450 finished: initializing file, took: 0:00:00.006763

2016-06-27 16:29:50.902834 starting: processing flux attribute

2016-06-27 16:29:50.902911 starting: creating variable: wel_flux

2016-06-27 16:29:50.903011 finished: creating variable: wel_flux, took: 0:00:00.000100

2016-06-27 16:29:50.903703 finished: processing flux attribute, took: 0:00:00.000869

2016-06-27 16:29:50.903748 starting: processing iface attribute

2016-06-27 16:29:50.903810 starting: creating variable: wel_iface

2016-06-27 16:29:50.903906 finished: creating variable: wel_iface, took: 0:00:00.000096

2016-06-27 16:29:50.904528 finished: processing iface attribute, took: 0:00:00.000780

2016-06-27 16:29:50.904783 starting: creating variable: hani

2016-06-27 16:29:50.904895 finished: creating variable: hani, took: 0:00:00.000112

2016-06-27 16:29:50.905717 starting: creating variable: hk

2016-06-27 16:29:50.905924 finished: creating variable: hk, took: 0:00:00.000207

2016-06-27 16:29:50.907233 starting: creating variable: ss

2016-06-27 16:29:50.907463 finished: creating variable: ss, took: 0:00:00.000230

2016-06-27 16:29:50.908792 starting: creating variable: sy

2016-06-27 16:29:50.908924 finished: creating variable: sy, took: 0:00:00.000132

2016-06-27 16:29:50.909626 starting: creating variable: vka

2016-06-27 16:29:50.909743 finished: creating variable: vka, took: 0:00:00.000117

2016-06-27 16:29:50.910460 starting: creating variable: vkcb

2016-06-27 16:29:50.910567 finished: creating variable: vkcb, took: 0:00:00.000107

2016-06-27 16:29:50.911308 starting: creating variable: wetdry

2016-06-27 16:29:50.911413 finished: creating variable: wetdry, took: 0:00:00.000105

Out[6]:
<class 'netCDF4._netCDF4.Dataset'>
root group (NETCDF4 data model, file format UNDEFINED):
    Conventions: CF-1.6
    date_created: 2016-06-27T06:29:00Z
    geospatial_vertical_positive: up
    geospatial_vertical_min: -27.3325004578
    geospatial_vertical_max: -17.3986301422
    geospatial_vertical_resolution: variable
    featureType: Grid
    origin_x: 0.0
    origin_y: 10000.0
    origin_crs: +init=EPSG:4326
    grid_rotation_from_origin: 0.0
    namefile: freyberg.nam
    solver_flux_tolerance: 0.001
    model_ws: data
    flopy_sr_xul: 0.0
    modflow_version: mf2005
    flopy_sr_yul: 10000.0
    create_hostname: IGSAAAHMLTW272M.local
    create_directory: /Users/jdhughes/Documents/Development/flopy_git/examples/Notebooks
    start_datetime: 1/1/1970
    solver_head_tolerance: 1e-05
    flopy_sr_proj4_str: +init=EPSG:4326
    create_platform: Darwin
    exe_name: mf2005.exe
    flopy_sr_rotation: 0.0
    dimensions(sizes): time(1), layer(1), y(40), x(20)
    variables(dimensions): int32 crs(), float64 time(time), float64 elevation(layer,y,x), float64 longitude(y,x), float64 latitude(y,x), float32 layer(layer), float32 delc(y), float32 delr(x), |S1 VerticalTransform(), float32 wel_flux(time,layer,y,x), float32 wel_iface(time,layer,y,x), float32 hani(layer,y,x), float32 hk(layer,y,x), float32 ss(layer,y,x), float32 sy(layer,y,x), float32 vka(layer,y,x), float32 vkcb(layer,y,x), float32 wetdry(layer,y,x)
    groups: 

Export the whole model to a netCDF and a shapefile


In [7]:
fnc = ml.export(os.path.join(model_ws, 'model.nc'))
ml.export(os.path.join(model_ws, 'model.shp'))
fnc.nc


2016-06-27 16:29:50.916410 WARNING: removing existing nc file: data/model.nc

2016-06-27 16:29:50.917103 WARNING: start datetime:1970-01-01T00:00:00Z

2016-06-27 16:29:50.917125 starting: initializing attributes

2016-06-27 16:29:50.917241 finished: initializing attributes, took: 0:00:00.000116

2016-06-27 16:29:50.917256 starting: initializing file

2016-06-27 16:29:50.917278 starting: initializing geometry

2016-06-27 16:29:50.917307 starting: building grid crs using proj4 string: +init=EPSG:4326

2016-06-27 16:29:50.917341 finished: building grid crs using proj4 string: +init=EPSG:4326, took: 0:00:00.000034

2016-06-27 16:29:50.917494 starting: projecting grid cell center arrays from +units=m +init=epsg:4326 to +units=m +init=epsg:4326 

2016-06-27 16:29:50.917562 starting: projecting grid cell center arrays from <pyproj.Proj object at 0x114fd5f48> to <pyproj.Proj object at 0x114fd5ee8>

2016-06-27 16:29:50.917613 finished: initializing geometry, took: 0:00:00.000335

2016-06-27 16:29:50.918091 starting: setting standard attributes

2016-06-27 16:29:50.918451 finished: setting standard attributes, took: 0:00:00.000360

2016-06-27 16:29:50.918475 starting: creating dimensions

2016-06-27 16:29:50.918603 finished: creating dimensions, took: 0:00:00.000128

2016-06-27 16:29:50.918616 starting: setting CRS info

2016-06-27 16:29:50.918711 finished: setting CRS info, took: 0:00:00.000095

2016-06-27 16:29:50.918740 starting: creating variable: time

2016-06-27 16:29:50.918845 finished: creating variable: time, took: 0:00:00.000105

2016-06-27 16:29:50.919175 WARNING: time_values:[ 10.]

2016-06-27 16:29:50.920052 starting: creating variable: elevation

2016-06-27 16:29:50.920186 finished: creating variable: elevation, took: 0:00:00.000134

2016-06-27 16:29:50.920731 starting: creating variable: longitude

2016-06-27 16:29:50.920838 finished: creating variable: longitude, took: 0:00:00.000107

2016-06-27 16:29:50.921393 starting: creating longitude var

2016-06-27 16:29:50.921427 starting: creating latitude var

2016-06-27 16:29:50.921451 starting: creating variable: latitude

2016-06-27 16:29:50.921542 finished: creating variable: latitude, took: 0:00:00.000091

2016-06-27 16:29:50.921966 starting: creating layer var

2016-06-27 16:29:50.922004 starting: creating variable: layer

2016-06-27 16:29:50.922200 finished: creating variable: layer, took: 0:00:00.000196

2016-06-27 16:29:50.922764 finished: creating layer var, took: 0:00:00.000798

2016-06-27 16:29:50.922803 starting: creating variable: delc

2016-06-27 16:29:50.922936 finished: creating variable: delc, took: 0:00:00.000133

2016-06-27 16:29:50.923357 starting: creating variable: delr

2016-06-27 16:29:50.923465 finished: creating variable: delr, took: 0:00:00.000108

2016-06-27 16:29:50.923920 finished: initializing file, took: 0:00:00.006664

2016-06-27 16:29:50.924126 starting: creating variable: botm

2016-06-27 16:29:50.924233 finished: creating variable: botm, took: 0:00:00.000107

2016-06-27 16:29:50.925051 starting: creating variable: thickness

2016-06-27 16:29:50.925178 finished: creating variable: thickness, took: 0:00:00.000127

2016-06-27 16:29:50.925997 starting: creating variable: model_top

2016-06-27 16:29:50.926179 finished: creating variable: model_top, took: 0:00:00.000182

2016-06-27 16:29:50.927344 starting: creating variable: ibound

2016-06-27 16:29:50.927881 finished: creating variable: ibound, took: 0:00:00.000537

2016-06-27 16:29:50.928831 starting: creating variable: strt

2016-06-27 16:29:50.928959 finished: creating variable: strt, took: 0:00:00.000128

2016-06-27 16:29:50.929736 starting: creating variable: rech

2016-06-27 16:29:50.929845 finished: creating variable: rech, took: 0:00:00.000109

2016-06-27 16:29:50.930709 starting: creating variable: hani

2016-06-27 16:29:50.930836 finished: creating variable: hani, took: 0:00:00.000127

2016-06-27 16:29:50.931812 starting: creating variable: hk

2016-06-27 16:29:50.931957 finished: creating variable: hk, took: 0:00:00.000145

2016-06-27 16:29:50.932698 starting: creating variable: ss

2016-06-27 16:29:50.932807 finished: creating variable: ss, took: 0:00:00.000109

2016-06-27 16:29:50.933654 starting: creating variable: sy

2016-06-27 16:29:50.933784 finished: creating variable: sy, took: 0:00:00.000130

2016-06-27 16:29:50.934541 starting: creating variable: vka

2016-06-27 16:29:50.934649 finished: creating variable: vka, took: 0:00:00.000108

2016-06-27 16:29:50.935407 starting: creating variable: vkcb

2016-06-27 16:29:50.935545 finished: creating variable: vkcb, took: 0:00:00.000138

2016-06-27 16:29:50.936442 starting: creating variable: wetdry

2016-06-27 16:29:50.936575 finished: creating variable: wetdry, took: 0:00:00.000133

2016-06-27 16:29:50.939622 starting: processing rbot attribute

2016-06-27 16:29:50.939734 starting: creating variable: riv_rbot

2016-06-27 16:29:50.939883 finished: creating variable: riv_rbot, took: 0:00:00.000149

2016-06-27 16:29:50.940907 finished: processing rbot attribute, took: 0:00:00.001285

2016-06-27 16:29:50.940985 starting: processing cond attribute

2016-06-27 16:29:50.941101 starting: creating variable: riv_cond

2016-06-27 16:29:50.941287 finished: creating variable: riv_cond, took: 0:00:00.000186

2016-06-27 16:29:50.942212 finished: processing cond attribute, took: 0:00:00.001227

2016-06-27 16:29:50.942268 starting: processing stage attribute

2016-06-27 16:29:50.942345 starting: creating variable: riv_stage

2016-06-27 16:29:50.942480 finished: creating variable: riv_stage, took: 0:00:00.000135

2016-06-27 16:29:50.943439 finished: processing stage attribute, took: 0:00:00.001171

2016-06-27 16:29:50.943630 starting: processing iface attribute

2016-06-27 16:29:50.943765 starting: creating variable: riv_iface

2016-06-27 16:29:50.943939 finished: creating variable: riv_iface, took: 0:00:00.000174

2016-06-27 16:29:50.944982 finished: processing iface attribute, took: 0:00:00.001352

2016-06-27 16:29:50.945472 starting: processing flux attribute

2016-06-27 16:29:50.945579 starting: creating variable: wel_flux

2016-06-27 16:29:50.945774 finished: creating variable: wel_flux, took: 0:00:00.000195

2016-06-27 16:29:50.946763 finished: processing flux attribute, took: 0:00:00.001291

2016-06-27 16:29:50.946809 starting: processing iface attribute

2016-06-27 16:29:50.946875 starting: creating variable: wel_iface

2016-06-27 16:29:50.946973 finished: creating variable: wel_iface, took: 0:00:00.000098

2016-06-27 16:29:50.947894 finished: processing iface attribute, took: 0:00:00.001085

Out[7]:
<class 'netCDF4._netCDF4.Dataset'>
root group (NETCDF4 data model, file format UNDEFINED):
    Conventions: CF-1.6
    date_created: 2016-06-27T06:29:00Z
    geospatial_vertical_positive: up
    geospatial_vertical_min: -27.3325004578
    geospatial_vertical_max: -17.3986301422
    geospatial_vertical_resolution: variable
    featureType: Grid
    origin_x: 0.0
    origin_y: 10000.0
    origin_crs: +init=EPSG:4326
    grid_rotation_from_origin: 0.0
    namefile: freyberg.nam
    solver_flux_tolerance: 0.001
    model_ws: data
    flopy_sr_xul: 0.0
    modflow_version: mf2005
    flopy_sr_yul: 10000.0
    create_hostname: IGSAAAHMLTW272M.local
    create_directory: /Users/jdhughes/Documents/Development/flopy_git/examples/Notebooks
    start_datetime: 1/1/1970
    solver_head_tolerance: 1e-05
    flopy_sr_proj4_str: +init=EPSG:4326
    create_platform: Darwin
    exe_name: mf2005.exe
    flopy_sr_rotation: 0.0
    dimensions(sizes): time(1), layer(1), y(40), x(20)
    variables(dimensions): int32 crs(), float64 time(time), float64 elevation(layer,y,x), float64 longitude(y,x), float64 latitude(y,x), float32 layer(layer), float32 delc(y), float32 delr(x), |S1 VerticalTransform(), float32 botm(layer,y,x), float32 thickness(layer,y,x), float32 model_top(y,x), int32 ibound(layer,y,x), float32 strt(layer,y,x), float32 rech(time,layer,y,x), float32 hani(layer,y,x), float32 hk(layer,y,x), float32 ss(layer,y,x), float32 sy(layer,y,x), float32 vka(layer,y,x), float32 vkcb(layer,y,x), float32 wetdry(layer,y,x), float32 riv_rbot(time,layer,y,x), float32 riv_cond(time,layer,y,x), float32 riv_stage(time,layer,y,x), float32 riv_iface(time,layer,y,x), float32 wel_flux(time,layer,y,x), float32 wel_iface(time,layer,y,x)
    groups: 

In [ ]: