def METHOD_NAME(fname):
    print("Reading and plotting file: {}".format(fname))
    name = fname
    f = Dataset(fname,'r')
    yr = f.variables['daysSinceStart'][:]/365.0
    yr = yr-yr[0]
    dt = f.variables['deltat'][:]/3.15e7
    print(yr.max())
    vol = f.variables['totalIceVolume'][:] / scaleVol
    if options.plotChange:
        vol = vol - vol[0]
    axVol.plot(yr, vol, label=name)
    VAF = f.variables['volumeAboveFloatation'][:] / scaleVol       
    if options.plotChange:
        VAF = VAF - VAF[0]
    axVAF.plot(yr, VAF, label=name)
    addSeaLevAx(axVAF)
    volGround = f.variables['groundedIceVolume'][:] / scaleVol
    if options.plotChange:
        volGround = volGround - volGround[0]
    axVolGround.plot(yr, volGround, label=name)
    volFloat = f.variables['floatingIceVolume'][:] / scaleVol
    if options.plotChange:
        volFloat = volFloat - volFloat[0]
    axVolFloat.plot(yr, volFloat, label=name)
    areaGrd = f.variables['groundedIceArea'][:] / 1000.0**2
    if options.plotChange:
        areaGrd = areaGrd - areaGrd[0]
    axGrdArea.plot(yr, areaGrd, label=name)
    areaFlt = f.variables['floatingIceArea'][:] / 1000.0**2
    if options.plotChange:
        areaFlt = areaFlt - areaFlt[0]
    axFltArea.plot(yr, areaFlt, label=name)
    GLflux = f.variables['groundingLineFlux'][:]
    axGLflux.plot(yr, GLflux, label=name)
    calvFlux = f.variables['totalCalvingFlux'][:]
    axCalvFlux.plot(yr, calvFlux, label=name)
    f.close()