climate-commits mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From huiky...@apache.org
Subject [1/3] climate git commit: CLIMATE-861 - An example to plot aerosol optical depth climatology from NASA MISR instrument
Date Thu, 25 Aug 2016 19:32:00 GMT
Repository: climate
Updated Branches:
  refs/heads/master 1c85945fa -> 9deb819f7


CLIMATE-861 - An example to plot aerosol optical depth climatology from NASA MISR instrument

- example/draw_climatology_map_MISR_AOD.py has been added.


Project: http://git-wip-us.apache.org/repos/asf/climate/repo
Commit: http://git-wip-us.apache.org/repos/asf/climate/commit/992ab591
Tree: http://git-wip-us.apache.org/repos/asf/climate/tree/992ab591
Diff: http://git-wip-us.apache.org/repos/asf/climate/diff/992ab591

Branch: refs/heads/master
Commit: 992ab591b8be01061340697750bf66d508adc301
Parents: 818784e
Author: huikyole <huikyole@argo.jpl.nasa.gov>
Authored: Wed Aug 24 16:28:54 2016 -0700
Committer: huikyole <huikyole@argo.jpl.nasa.gov>
Committed: Wed Aug 24 16:28:54 2016 -0700

----------------------------------------------------------------------
 examples/draw_climatology_map_MISR_AOD.py | 34 ++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)
----------------------------------------------------------------------


http://git-wip-us.apache.org/repos/asf/climate/blob/992ab591/examples/draw_climatology_map_MISR_AOD.py
----------------------------------------------------------------------
diff --git a/examples/draw_climatology_map_MISR_AOD.py b/examples/draw_climatology_map_MISR_AOD.py
new file mode 100644
index 0000000..1af834d
--- /dev/null
+++ b/examples/draw_climatology_map_MISR_AOD.py
@@ -0,0 +1,34 @@
+import ocw.dataset as ds
+import ocw.data_source.local as local
+import ocw.dataset_processor as dsp
+import ocw.plotter as plotter
+
+import numpy as np
+import numpy.ma as ma
+
+
+''' data source: https://dx.doi.org/10.6084/m9.figshare.3753321.v1'''
+dataset = local.load_file('/home/huikyole/climate/examples/AOD_monthly_2000-MAR_2016-FEB_from_MISR_L3_JOINT.nc',

+                          'nonabsorbing_ave')
+''' Subset the data for East Asia'''
+Bounds = ds.Bounds(lat_min=20, lat_max=57.7, lon_min=90, lon_max=150)
+dataset = dsp.subset(dataset, Bounds)
+
+'''The original dataset includes nonabsorbing AOD values between March 2000 and February
2015. 
+dsp.temporal_subset will extract data in September-October-November.'''
+dataset_SON = dsp.temporal_subset(dataset, month_start=9, month_end=11, average_each_year=True)
+
+ny, nx = dataset_SON.values.shape[1:]
+
+# multi-year mean aod
+clim_aod = ma.zeros([3, ny, nx])
+
+clim_aod[0,:] = ma.mean(dataset_SON.values, axis=0) # 16-year mean
+clim_aod[1,:] = ma.mean(dataset_SON.values[-5:,:], axis=0) # the last 5-year mean
+clim_aod[2,:] = dataset_SON.values[-1,:] # the last year's value
+
+# plot clim_aod (3 subplots)
+plotter.draw_contour_map(clim_aod, dataset_SON.lats, dataset_SON.lons, 
+                         fname='nonabsorbing_AOD_clim_East_Asia_Sep-Nov',
+                         gridshape=[1,3],subtitles=['2000-2015: 16 years','2011-2015: 5 years',
'2015: 1 year'], 
+                         clevs=np.arange(21)*0.02)


Mime
View raw message