Skip to content

Instantly share code, notes, and snippets.

@is
Created June 29, 2016 02:12
Show Gist options
  • Select an option

  • Save is/ab6ae0f173576e5c632fdc7256bb975b to your computer and use it in GitHub Desktop.

Select an option

Save is/ab6ae0f173576e5c632fdc7256bb975b to your computer and use it in GitHub Desktop.
check O3 emission.
import numpy as np
import netCDF4 as nc
fn = 'Daily.Combine.Surf.Column.AOD.together.d03.2016062'
nc = nc.Dataset(fn)
print(nc.variables)
o3 = nc.variables['O3']
o3arr = o3[:]
print(np.min(o3arr), np.max(o3arr), np.mean(o3arr))
print(len(np.where(o3arr < 1e-10)[0]))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment