"""Test :meth:`~vacumm.data.misc.dataset.OceanDataset.get_dens` in MENOR""" # Inits ncfile = "menor.nc" # Imports from vcmq import DS, os, map2, data_sample, code_file_name # Read data ds = DS(data_sample(ncfile), 'mars', logger_level='critical') dens = ds.get_dens(squeeze=True) # Plot surface density figfile = code_file_name(ext='png') if os.path.exists(figfile): os.remove(figfile) map2(dens[-1], savefig=figfile, close=True, vmin=1025., show=False) # For unittest result = {'files': figfile}
"""Test :meth:`~vacumm.data.misc.dataset.OceanDataset.get_dens` in MENOR""" # Inits ncfile = "menor.nc" # Imports from vcmq import DS, os, map2, data_sample, code_file_name # Read data ds = DS(data_sample(ncfile), 'mars', logger_level='critical') dens = ds.get_dens(squeeze=True) # Plot surface density figfile = code_file_name(ext='png') if os.path.exists(figfile): os.remove(figfile) map2(dens[-1], savefig=figfile, close=True, vmin=1025.,show=False) # For unittest result = {'files':figfile}
# Inits ncfile = "menor.nc" lon0 = 3.31 lat0 = 43.19 lon1 = 4.78 lat1 = 42.01 # Imports from vcmq import DS, data_sample, mixed_layer_depth, N, transect, section2, curve2, code_file_name, os from vacumm.misc.plot import add_map_lines # Read data ds = DS(data_sample(ncfile), 'mars', logger_level='critical') temp = ds.get_temp(squeeze=True) dens = ds.get_dens(squeeze=True, potential=True) depth = ds.get_depth(squeeze=True) # Compute MLD mld = mixed_layer_depth(dens, depth=depth, mode='deltadens', format_axes=True) del dens # Compute transect tlons = (lon0, lon1) tlats = (lat0, lat1) tlons = N.concatenate( [N.linspace(lon0, lon1, 100.), N.linspace(lon1, lon1, 100.)]) tlats = N.concatenate( [N.linspace(lat0, lat1, 100.), N.linspace(lat1, lat0, 100.)])