def METHOD_NAME(args):
    '''
    Read in a netCDF file generated by Nalu and post-process.
    '''
    # Load the data
    data = ABLStatsFileClass(stats_file='abl_statistics.nc')
    '''
    Plot time-history for different heights
    '''
    # The heights to plot [m]
    heights = [70, 100, 150, 300]
    # The time vector
    time = data.time
    plt.clf()
    # Loop through all the heights
    for z in heights:
        # Velocity as function of time
        u = data.u_h(z)
        v = data.v_h(z)
        # Velocity magnitude
        u_mag = np.sqrt(u**2 + v**2)
        plt.plot(time, u_mag, label=r'$U_{\rm mag},$ $z=$'+str(z)+' [m]')
    plt.legend(loc='best')
    plt.xlabel('Time [s]')
    plt.ylabel('Velocity [m/s]')
    plt.savefig('Velocity.pdf')
    '''
    Plot time-averaged profiles for all heights
    '''
    # Velocity
    z = data.heights
    u = data.time_average(field='velocity', index=0, times=[0, 1500])
    v = data.time_average(field='velocity', index=1, times=[0, 1500])
    u_mag = np.sqrt(u**2 + v**2)
    plt.clf()
    plt.plot(u, z, '-o', label='U')
    plt.plot(v, z, '-o', label='V')
    plt.plot(u_mag, z, '-o', label='U mag')
    plt.legend(loc='best')
    plt.xlabel('Velocity [m/s]')
    plt.ylabel('z [m]')
    plt.savefig('Velocity_average.pdf')
    # Resolved Stress
    uu = data.time_average(field='resolved_stress', index=0, times=[0, 1500])
    vv = data.time_average(field='resolved_stress', index=2, times=[0, 1500])
    plt.clf()
    plt.plot(uu, z, '-o', label='U')
    plt.plot(vv, z, '-o', label='V')
    plt.legend(loc='best')
    plt.xlabel("Resolved Stress [m$^2$/s$^2$]")
    plt.ylabel('z [m]')
    plt.savefig('rs_average.pdf')