In [1]:
import iris
import iris.plot as iplt
import matplotlib.pyplot as plt
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
return f(*args, **kwds)
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
return f(*args, **kwds)
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
return f(*args, **kwds)
In [2]:
wfo_file = '/g/data1b/oi10/replicas/CMIP6/CMIP/IPSL/IPSL-CM6A-LR/piControl/r1i1p1f1/Omon/wfo/gn/v20181123/wfo_Omon_IPSL-CM6A-LR_piControl_r1i1p1f1_gn_285001-304912.nc'
wfo_cube = iris.load_cube(wfo_file, 'water_flux_into_sea_water')
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/site-packages/iris/fileformats/cf.py:798: UserWarning: Missing CF-netCDF measure variable 'areacello', referenced by netCDF variable 'wfo'
warnings.warn(message % (variable_name, nc_var_name))
In [3]:
wfo_cube
Out[3]:
Water Flux Into Sea Water (kg m-2 s-1)
time
--
--
Shape
2400
332
362
Dimension coordinates
time
x
-
-
Auxiliary coordinates
latitude
-
x
x
longitude
-
x
x
Attributes
CMIP6_CV_version
cv=6.2.3.5-2-g63b123e
Conventions
CF-1.7 CMIP-6.2
EXPID
piControl
activity_id
CMIP
branch_method
standard
branch_time_in_child
0.0
branch_time_in_parent
36524.0
contact
ipsl-cmip6@listes.ipsl.fr
creation_date
2018-09-11T13:58:03Z
data_specs_version
01.00.21
description
computed as the water flux into the ocean divided by the area of the ocean...
dr2xml_md5sum
92ddb3d0d8ce79f498d792fc8e559dcf
dr2xml_version
1.13
experiment
pre-industrial control
experiment_id
piControl
external_variables
areacello
forcing_index
1
frequency
mon
further_info_url
https://furtherinfo.es-doc.org/CMIP6.IPSL.IPSL-CM6A-LR.piControl.none....
grid
native ocean tri-polar grid with 105 k ocean cells
grid_label
gn
history
none
initialization_index
1
institution
Institut Pierre Simon Laplace, Paris 75252, France
institution_id
IPSL
interval_operation
2700 s
interval_write
1 month
license
CMIP6 model data produced by IPSL is licensed under a Creative Commons...
mip_era
CMIP6
model_version
6.1.6
name
/ccc/work/cont003/gencmip6/p86maf/IGCM_OUT/IPSLCM6/PROD/piControl/CM61...
nominal_resolution
100 km
online_operation
average
parent_activity_id
CMIP
parent_experiment_id
piControl-spinup
parent_mip_era
CMIP6
parent_source_id
IPSL-CM6A-LR
parent_time_units
days since 1750-01-01 00:00:00
parent_variant_label
r1i1p1f1
physics_index
1
product
model-output
realization_index
1
realm
ocean
source
IPSL-CM6A-LR (2017): atmos: LMDZ (NPv6, N96; 144 x 143 longitude/latitude;...
source_id
IPSL-CM6A-LR
source_type
AOGCM BGC
sub_experiment
none
sub_experiment_id
none
table_id
Omon
title
IPSL-CM6A-LR model output prepared for CMIP6 / CMIP piControl
tracking_id
hdl:21.14100/c3aa234a-ee21-4ee7-a476-4d80952a6d43
variable_id
wfo
variant_label
r1i1p1f1
Cell methods
mean where sea
area
mean
time
In [4]:
plot_cube = wfo_cube.collapsed('time', iris.analysis.MEAN)
/g/data/r87/dbi599/miniconda3/envs/ocean/lib/python3.6/importlib/_bootstrap.py:219: RuntimeWarning: numpy.dtype size changed, may indicate binary incompatibility. Expected 96, got 88
return f(*args, **kwds)
In [5]:
plot_cube
Out[5]:
Water Flux Into Sea Water (kg m-2 s-1)
--
--
Shape
332
362
Auxiliary coordinates
latitude
x
x
longitude
x
x
Scalar coordinates
time
2950-01-01 00:00:00, bound=(2850-01-01 00:00:00, 3050-01-01 00:00:00)
Attributes
CMIP6_CV_version
cv=6.2.3.5-2-g63b123e
Conventions
CF-1.7 CMIP-6.2
EXPID
piControl
activity_id
CMIP
branch_method
standard
branch_time_in_child
0.0
branch_time_in_parent
36524.0
contact
ipsl-cmip6@listes.ipsl.fr
creation_date
2018-09-11T13:58:03Z
data_specs_version
01.00.21
description
computed as the water flux into the ocean divided by the area of the ocean...
dr2xml_md5sum
92ddb3d0d8ce79f498d792fc8e559dcf
dr2xml_version
1.13
experiment
pre-industrial control
experiment_id
piControl
external_variables
areacello
forcing_index
1
frequency
mon
further_info_url
https://furtherinfo.es-doc.org/CMIP6.IPSL.IPSL-CM6A-LR.piControl.none....
grid
native ocean tri-polar grid with 105 k ocean cells
grid_label
gn
history
none
initialization_index
1
institution
Institut Pierre Simon Laplace, Paris 75252, France
institution_id
IPSL
interval_operation
2700 s
interval_write
1 month
license
CMIP6 model data produced by IPSL is licensed under a Creative Commons...
mip_era
CMIP6
model_version
6.1.6
name
/ccc/work/cont003/gencmip6/p86maf/IGCM_OUT/IPSLCM6/PROD/piControl/CM61...
nominal_resolution
100 km
online_operation
average
parent_activity_id
CMIP
parent_experiment_id
piControl-spinup
parent_mip_era
CMIP6
parent_source_id
IPSL-CM6A-LR
parent_time_units
days since 1750-01-01 00:00:00
parent_variant_label
r1i1p1f1
physics_index
1
product
model-output
realization_index
1
realm
ocean
source
IPSL-CM6A-LR (2017): atmos: LMDZ (NPv6, N96; 144 x 143 longitude/latitude;...
source_id
IPSL-CM6A-LR
source_type
AOGCM BGC
sub_experiment
none
sub_experiment_id
none
table_id
Omon
title
IPSL-CM6A-LR model output prepared for CMIP6 / CMIP piControl
tracking_id
hdl:21.14100/c3aa234a-ee21-4ee7-a476-4d80952a6d43
variable_id
wfo
variant_label
r1i1p1f1
Cell methods
mean where sea
area
mean
time
mean
time
In [6]:
vmin = -0.0003
vmax = 0.0003
#iplt.pcolormesh(plot_cube, cmap='RdBu_r', vmin=vmin, vmax=vmax)
plt.pcolormesh(plot_cube.data, cmap='RdBu_r', vmin=vmin, vmax=vmax)
cb = plt.colorbar()
cb.set_label(str(plot_cube.units))
plt.show()
In [ ]:
Content source: DamienIrving/ocean-analysis
Similar notebooks: