MOi hydrodynamics (psy4v3r1 daily)

• Location of data on lorenz: /storage/shared/oceanparcels/input_data/MOi/psy4v3r1
• location of coordinates on lorenz: /storage/shared/oceanparcels/input_data/MOi/domain_ORCA0083-N006
• Simulation: NEMO Global 1/12° Analysis & Forecast psy4v3r1-daily (GLO12)
• Region: Global
• Period: 2007/01/01 - 2023/03/15 (spin-up from 2006/10/11)
• Frequency: daily
• Variables: See full list here
• Vertical levels: 50 depth levels
• Grid: ORCA0083 Native C-grid
• Source: https://tds.mercator-ocean.fr/thredds/dodsC
• Documentation: https://www.mercator-ocean.fr/wp-content/uploads/products_pdf/MOI_GLO12.pdf
• Used by: Mikael Kaandorp, Daan Reijnders, Claudio Pierard

Parcels code to load 2D fields

data_path = '/storage/shared/oceanparcels/input_data/MOi/'
ufiles = sorted(glob(data_path+'psy4v3r1/psy4v3r1-daily_U_*.nc'))
vfiles = [f.replace('_U_', '_V_') for f in ufiles]
wfiles = [f.replace('_U_', '_W_') for f in ufiles]
mesh_mask = data_path + 'domain_ORCA0083-N006/PSY4V3R1_mesh_hgr.nc'

filenames = {'U': {'lon': mesh_mask, 'lat': mesh_mask, 'depth': wfiles[0], 'data': ufiles},
             'V': {'lon': mesh_mask, 'lat': mesh_mask, 'depth': wfiles[0], 'data': vfiles}}
variables = {'U': 'vozocrtx', 'V': 'vomecrty'}
dimensions = {'U': {'lon': 'glamf', 'lat': 'gphif', 'depth': 'depthw', 'time': 'time_counter'},
              'V': {'lon': 'glamf', 'lat': 'gphif', 'depth': 'depthw', 'time': 'time_counter'}}
indices = {'depth': range(0, 2)}  # Note at least two depth layers need to be loaded
fieldset = FieldSet.from_nemo(filenames, variables, dimensions, indices=indices)

Parcels code to load 3D fields

data_path = '/storage/shared/oceanparcels/input_data/MOi/'
ufiles = sorted(glob(data_path+'psy4v3r1/psy4v3r1-daily_U_*.nc'))
vfiles = [f.replace('_U_', '_V_') for f in ufiles]
wfiles = [f.replace('_U_', '_W_') for f in ufiles]
mesh_mask = data_path + 'domain_ORCA0083-N006/PSY4V3R1_mesh_hgr.nc'

filenames = {'U': {'lon': mesh_mask, 'lat': mesh_mask, 'depth': wfiles[0], 'data': ufiles},
             'V': {'lon': mesh_mask, 'lat': mesh_mask, 'depth': wfiles[0], 'data': vfiles},
             'W': {'lon': mesh_mask, 'lat': mesh_mask, 'depth': wfiles[0], 'data': wfiles}}
variables = {'U': 'vozocrtx', 'V': 'vomecrty', 'W': 'vovecrtz'}
dimensions = {'U': {'lon': 'glamf', 'lat': 'gphif', 'depth': 'depthw', 'time': 'time_counter'},
              'V': {'lon': 'glamf', 'lat': 'gphif', 'depth': 'depthw', 'time': 'time_counter'},
              'W': {'lon': 'glamf', 'lat': 'gphif', 'depth': 'depthw', 'time': 'time_counter'}}

fieldset = FieldSet.from_nemo(filenames, variables, dimensions)