Use Case 11

Problem Definition:

As part of a project on Arctic climate change, an undergraduate student wishes to look at different ECVs on a polar stereographic projection.

Required Toolbox Features:

Access to and ingestion of CCI ECV data (e.g. sea ice, ice sheets, sea level, SST, clouds aerosol)
Access to and ingestion of ECV data from external server
Remapping to fit data onto user-chosen projection
Spatial and temporal subsetting
Gap-filling (user-chosen strategy)
Generation of scalable maps



In [ ]:

    
# To reduce bandwith, spatial and temporal subsetting is done when downloading the dataset
from cate import ops
from cate.util.monitor import ConsoleMonitor

monitor = ConsoleMonitor()



In [ ]:

    
sst_polar = ops.open_dataset(ds_id='esacci.SST.day.L4.SSTdepth.multi-sensor.multi-platform.OSTIA.1-1.r1',
                             region=(-180, 60, 180, 90),
                             time_range='2007-01-01, 2007-01-15',
                             var_names='analysed_sst, sea_ice_fraction',
                             force_local=True,
                             local_ds_id='SST_polar_2007',
                             normalize=False,
                             monitor=monitor)



In [ ]:

    
# Above creates the dataset, but errors out when trying to open it
sst_polar = ops.open_dataset('local.SST_polar_2007')



In [ ]:

    
%matplotlib inline
ops.plot_map(sst_polar, var='sea_ice_fraction', region=(-180, 60, 180, 90), projection='NorthPolarStereo',
             contour_plot=True, file='sst.svg')



In [ ]:

    
cld_polar = ops.open_dataset(ds_id='esacci.CLOUD.mon.L3C.CLD_PRODUCTS.AVHRR.multi-platform.AVHRR-AM.2-0.r1',
                             region=(-180, 60, 180, 90),
                             time_range='2007-01-01, 2007-01-31',
                             var_names='cfc, cee',
                             force_local=True,
                             local_ds_id='CLOUDS_polar_2007',
                             normalize=False,
                             monitor=monitor)



In [ ]:

    
clouds = ops.open_dataset('local.CLOUDS_polar_2007')



In [ ]:

    
ops.plot_map(clouds, var='cfc', region=(-180, 60, 180, 90), projection='NorthPolarStereo', file='polar_clouds.svg')



In [ ]:

    
aerosol_polar = ops.open_dataset(ds_id='esacci.AEROSOL.mon.L3.AAI.multi-sensor.multi-platform.ms_uvai.1-5-7.r1',
                                 region=(-180, 60, 180, 90),
                                 time_range='2007-01-01, 2007-03-31',
                                 force_local=True,
                                 local_ds_id='AEROSOL_polar_2007',
                                 monitor=monitor)



In [ ]:

    
aerosol = ops.open_dataset('local.AEROSOL_polar_2007')



In [ ]:

    
ops.plot_map(aerosol, region=(-180, 60, 180, 90), var='absorbing_aerosol_index', projection='NorthPolarStereo',
             file='aerosol_polar.svg')



In [ ]: