Prior to the 2015 NCK days icebreaker, OpenEarth.nl organizes a hands-on sprint session to explain how to access and visualize some great open access core-description and grain size datasets of the Dutch coast:
i) Sediment Atlas WaddenZee (RWS): 7500(!) distribution curves in Marsdiep-EmsDollard
ii) TNO 13.000(!) core descriptions for the Dutch continental shelf via NODC/SeaDataNet (and much, much more data available via DINOloket.nl)
iii) TNO interpolated dz10, dz50, dz90 maps with 200m resolution for Dutch Continental Shelf
iv) Along the Dutch beaches, grain sizes for dune safety assessment (RWS)
During this sprint session you’ll learn how to work with these free data in Matlab and/or python. Bring your own laptop to work with the data, and take home your data products and visualizations.
Time: March 18th, 13:30-17:00 (max 10 participants, sponsored by NCK)
Location: http://www.ijgenweisschoorl.nl (where icebreaker starts at 18:30)
Registration: Gerben J. de Boer (Van Oord) Gerben.deBoer@vanoord.com ism TNO, Deltares, TUD
Flyer: NCK days 2015 sprintsession north sea grain sizes.pdf
Course information:
Access information for the data can be found at Dataset documentation.
You can download netCDF files (to prevent wifi overload) with BITSadmin. An example is for the TNO gain size maps is:
set from=http://opendap.deltares.nl/thredds/fileServer/opendap/tno/ncp/ set into=d:\opendap\tno\ bitsadmin /transfer Job001 /download /priority normal %from%slib_juli2007.nc %into%slib_juli2007.nc bitsadmin /transfer Job001 /download /priority normal %from%grind_fbr2007.nc %into%grind_fbr2007.nc bitsadmin /transfer Job001 /download /priority normal %from%dz10_juli2007.nc %into%dz10_juli2007.nc bitsadmin /transfer Job001 /download /priority normal %from%dz50_juli2007.nc %into%dz50_juli2007.nc bitsadmin /transfer Job001 /download /priority normal %from%dz90_juli2007.nc %into%dz90_juli2007.nc
Yoou can access subsets these files with the Matlab netCDF package called ncread and with the python package netcdf4-python, see our Matlab example and
python example. A Matlab simple example how to make a publication quality figure is given here
%% load data (local cache), D contains data, M contains all metadata [D,M] = nc2struct('slib_juli2007.nc') %% load data (via opendap) L = nc2struct('http://opendap.deltares.nl/thredds/dodsC/opendap/deltares/landboundaries/northsea.nc') %% plot pcolorcorcen(D.longitude,D.latitude,D.mud) colorbarwithvtext([M.mud.long_name,'[',M.mud.units,']']) tickmap('ll') grid on axislat plot(L.lon,L.lat,'k') set(findfont,'fontsize',8) title({M.nc_global.institution,'NCP map mud [version jul 2007]'}) print2screensize('slib_juli2007')