%matplotlib inline
import pandas as pd
import socket
host = socket.getfqdn()
from core import load, zoom, calc, save,plots,monitor
#reload funcs after updating ./core/*.py
import importlib
importlib.reload(load)
importlib.reload(zoom)
importlib.reload(calc)
importlib.reload(save)
importlib.reload(plots)
importlib.reload(monitor)
<module 'core.monitor' from '/ccc/work/cont003/gen7420/odakatin/monitor-sedna/notebook/core/monitor.py'>
# 'month': = 'JOBID' almost month but not really,
# If you submit the job with job scheduler, above
#below are list of enviroment variable one can pass
#%env local='2"
# local : if True run dask local cluster, if not true, put number of workers
# setted in the 'local'
# if no 'local ' given, local will be setted automatically to 'True'
#%env ychunk='2'
#%env tchunk='2'
# controls chunk. 'False' sets no modification from original netcdf file's chunk.
# ychunk=10 will group the original netcdf file to 10 by 10
# tchunk=1 will chunk the time coordinate one by one
#%env control=Fluxnet
# name of control file to be used for computation/plots/save/
#%env file_exp=
# 'file_exp': Which 'experiment' name is it?
#. this corresopnds to intake catalog name without path and .yaml
#%env year=
# for Validation, this correspoinds to path/year/month 's year
# for monitoring, this corresponids to 'date' having * means do all files in the monitoring directory
# setting it as *0[0-9] &*1[0-9]& *[2-3][0-9], the job can be separated in three lots.
#%env month=
# for monitoring this corresponds to file path path-XIOS.{month}/
#
#%env save= proceed saving? True or False , Default is setted as True
#%env plot= proceed plotting? True or False , Default is setted as True
#%env calc= proceed computation? or just load computed result? True or False , Default is setted as True
#%env calc=True
%%time
# 'savefig': Do we save output in html? or not. keep it true.
savefig=True
client,cluster,control,catalog_url,month,year,daskreport,outputpath = load.set_control(host)
!mkdir -p $outputpath
!mkdir -p $daskreport
client
local True using host= irene5317.c-irene.mg1.tgcc.ccc.cea.fr starting dask cluster on local= True workers 16 10000000000 False not local in tgcc rome local cluster starting This code is running on irene5317.c-irene.mg1.tgcc.ccc.cea.fr using SEDNA_ALPHA_MONITOR file experiment, read from ../lib/SEDNA_ALPHA_MONITOR.yaml on year= * on month= 23 outputpath= ../results/SEDNA_ALPHA_MONITOR/23/ daskreport= ../results/dask/2526596irene5317.c-irene.mg1.tgcc.ccc.cea.fr_SEDNA_ALPHA_MONITOR_23Ice_quantities/ CPU times: user 324 ms, sys: 250 ms, total: 575 ms Wall time: 9.6 s
Client
|
Cluster
|
df=load.controlfile(control)
#Take out 'later' tagged computations
#df=df[~df['Value'].str.contains('later')]
df
Value | Inputs | Equation | Zone | Plot | Colourmap | MinMax | Unit | Oldname | Unnamed: 10 | |
---|---|---|---|---|---|---|---|---|---|---|
Ice_quantities | param.e1te2t,icemod.sivelo,icemod.sivolu,icemo... | calc.Ice_quant(data) | BBFG | Ice_intquant | None | (0,20) | cm s^(-1) | I-2 |
Each computation consists of
%%time
import os
calcswitch=os.environ.get('calc', 'True')
print('calcswitch=',calcswitch)
data = load.datas(catalog_url,df.Inputs,month,year,daskreport) if calcswitch=='True' else 0
data
calcswitch= True ../lib/SEDNA_ALPHA_MONITOR.yaml using param_xios reading ../lib/SEDNA_ALPHA_MONITOR.yaml using param_xios reading <bound method DataSourceBase.describe of sources: param_xios: args: combine: by_coords concat_dim: y urlpath: /ccc/work/cont003/gen7420/odakatin/CONFIGS/SEDNA/SEDNA-I/SEDNA_Domain_cfg_Tgt_20210423_tsh10m_L1/param_f32/x_*.nc xarray_kwargs: compat: override coords: minimal data_vars: minimal parallel: true description: SEDNA NEMO parameters from MPI output nav_lon lat fails driver: intake_xarray.netcdf.NetCDFSource metadata: catalog_dir: /ccc/work/cont003/gen7420/odakatin/monitor-sedna/notebook/../lib/ > {'name': 'param_xios', 'container': 'xarray', 'plugin': ['netcdf'], 'driver': ['netcdf'], 'description': 'SEDNA NEMO parameters from MPI output nav_lon lat fails', 'direct_access': 'forbid', 'user_parameters': [{'name': 'path', 'description': 'file coordinate', 'type': 'str', 'default': '/ccc/work/cont003/gen7420/odakatin/CONFIGS/SEDNA/MESH/SEDNA_mesh_mask_Tgt_20210423_tsh10m_L1/param'}], 'metadata': {}, 'args': {'urlpath': '/ccc/work/cont003/gen7420/odakatin/CONFIGS/SEDNA/SEDNA-I/SEDNA_Domain_cfg_Tgt_20210423_tsh10m_L1/param_f32/x_*.nc', 'combine': 'by_coords', 'concat_dim': 'y'}} 0 read icemod ['sivelo', 'siconc', 'sivolu'] using load_data_xios reading icemod using load_data_xios reading <bound method DataSourceBase.describe of sources: data_xios: args: combine: by_coords concat_dim: time_counter,x,y urlpath: /ccc/scratch/cont003/gen7420/talandel/ONGOING-RUNS/SEDNA-ALPHA-XIOS.23/SEDNA-ALPHA_1d_icemod_*_0[0-5][0-9][0-9].nc xarray_kwargs: compat: override coords: minimal data_vars: minimal drop_variables: !!set botpres: null deptht_bounds: null depthu_bounds: null iicestru: null iicestrv: null intstrx: null intstry: null mldkz5: null rhop_sig0: null siages: null sidive: null sisali: null sishea: null sistre: null sitemp: null snthic: null snvolu: null sometauy: null sozotaux: null time_centered_bounds: null time_counter_bounds: null utau_atmoce: null utau_iceoce: null uwspd10: null vtau_atmoce: null vtau_iceoce: null vwspd10: null parallel: true preprocess: !!python/name:core.load.prep '' description: SEDNA NEMO outputs from different xios server driver: intake_xarray.netcdf.NetCDFSource metadata: catalog_dir: /ccc/work/cont003/gen7420/odakatin/monitor-sedna/notebook/../lib/ > took 432.994909286499 seconds 0 merging icemod ['sivelo', 'siconc', 'sivolu'] param nav_lon will be included in data param e1te2t will be included in data param nav_lat will be included in data param mask2d will be included in data sum_num (13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12) start rechunking with (130, 122, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 120, 48) end of y_rechunk CPU times: user 3min 23s, sys: 30.8 s, total: 3min 54s Wall time: 7min 27s
<xarray.Dataset> Dimensions: (t: 31, x: 6560, y: 6540) Coordinates: nav_lat (y, x) float32 dask.array<chunksize=(130, 6560), meta=np.ndarray> nav_lon (y, x) float32 dask.array<chunksize=(130, 6560), meta=np.ndarray> * t (t) object 2004-07-01 12:00:00 ... 2004-07-31 12:00:00 * y (y) int64 1 2 3 4 5 6 7 8 ... 6534 6535 6536 6537 6538 6539 6540 * x (x) int64 1 2 3 4 5 6 7 8 ... 6554 6555 6556 6557 6558 6559 6560 e1te2t (y, x) float64 dask.array<chunksize=(130, 6560), meta=np.ndarray> mask2d (y, x) bool dask.array<chunksize=(130, 6560), meta=np.ndarray> Data variables: sivelo (t, y, x) float32 dask.array<chunksize=(1, 130, 6560), meta=np.ndarray> siconc (t, y, x) float32 dask.array<chunksize=(1, 130, 6560), meta=np.ndarray> sivolu (t, y, x) float32 dask.array<chunksize=(1, 130, 6560), meta=np.ndarray> Attributes: name: /ccc/scratch/cont003/gen7420/talandel/ONGOING-RU... description: ice variables title: ice variables Conventions: CF-1.6 timeStamp: 2021-Aug-04 15:32:27 GMT uuid: be0d49a1-d456-4e73-a33e-e553bcdcaa0a ibegin: 0 ni: 6560 jbegin: 0 nj: 13 DOMAIN_number_total: 544 DOMAIN_number: 0 DOMAIN_dimensions_ids: [2 3] DOMAIN_size_global: [6560 6540] DOMAIN_size_local: [6560 13] DOMAIN_position_first: [1 1] DOMAIN_position_last: [6560 13] DOMAIN_halo_size_start: [0 0] DOMAIN_halo_size_end: [0 0] DOMAIN_type: box start_date: 20030101 output_frequency: 1d CONFIG: SEDNA CASE: ALPHA
|
|
array([cftime.DatetimeNoLeap(2004, 7, 1, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 2, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 3, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 4, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 5, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 6, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 7, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 8, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 9, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 10, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 11, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 12, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 13, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 14, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 15, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 16, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 17, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 18, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 19, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 20, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 21, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 22, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 23, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 24, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 25, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 26, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 27, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 28, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 29, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 30, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 31, 12, 0, 0, 0)], dtype=object)
array([ 1, 2, 3, ..., 6538, 6539, 6540])
array([ 1, 2, 3, ..., 6558, 6559, 6560])
|
|
|
|
|
%%time
monitor.auto(df,data,savefig,daskreport,outputpath,file_exp='SEDNA'
)
#calc= True #save= True #plot= False Zone='BBFG' Value='Ice_quantities' cmap='None' clabel='cm s^(-1)' clim= (0, 20) outputpath='../results/SEDNA_ALPHA_MONITOR/23/' nc_outputpath='../nc_results/SEDNA_ALPHA_MONITOR/23/' filename='SEDNA_Ice_intquant_BBFG_Ice_quantities' #2 Zooming Data dataa= zoom.BBFG(data)
<xarray.Dataset> Dimensions: (t: 31, x: 6560, y: 5264) Coordinates: * x (x) int64 1 2 3 4 5 6 7 8 ... 6554 6555 6556 6557 6558 6559 6560 * y (y) int64 1277 1278 1279 1280 1281 ... 6536 6537 6538 6539 6540 nav_lat (y, x) float32 dask.array<chunksize=(56, 6560), meta=np.ndarray> nav_lon (y, x) float32 dask.array<chunksize=(56, 6560), meta=np.ndarray> * t (t) object 2004-07-01 12:00:00 ... 2004-07-31 12:00:00 e1te2t (y, x) float64 dask.array<chunksize=(56, 6560), meta=np.ndarray> mask2d (y, x) bool dask.array<chunksize=(56, 6560), meta=np.ndarray> Data variables: sivelo (t, y, x) float32 dask.array<chunksize=(1, 56, 6560), meta=np.ndarray> siconc (t, y, x) float32 dask.array<chunksize=(1, 56, 6560), meta=np.ndarray> sivolu (t, y, x) float32 dask.array<chunksize=(1, 56, 6560), meta=np.ndarray> Attributes: name: /ccc/scratch/cont003/gen7420/talandel/ONGOING-RU... description: ice variables title: ice variables Conventions: CF-1.6 timeStamp: 2021-Aug-04 15:32:27 GMT uuid: be0d49a1-d456-4e73-a33e-e553bcdcaa0a ibegin: 0 ni: 6560 jbegin: 0 nj: 13 DOMAIN_number_total: 544 DOMAIN_number: 0 DOMAIN_dimensions_ids: [2 3] DOMAIN_size_global: [6560 6540] DOMAIN_size_local: [6560 13] DOMAIN_position_first: [1 1] DOMAIN_position_last: [6560 13] DOMAIN_halo_size_start: [0 0] DOMAIN_halo_size_end: [0 0] DOMAIN_type: box start_date: 20030101 output_frequency: 1d CONFIG: SEDNA CASE: ALPHA
array([ 1, 2, 3, ..., 6558, 6559, 6560])
array([1277, 1278, 1279, ..., 6538, 6539, 6540])
|
|
array([cftime.DatetimeNoLeap(2004, 7, 1, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 2, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 3, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 4, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 5, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 6, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 7, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 8, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 9, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 10, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 11, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 12, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 13, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 14, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 15, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 16, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 17, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 18, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 19, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 20, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 21, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 22, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 23, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 24, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 25, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 26, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 27, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 28, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 29, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 30, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 31, 12, 0, 0, 0)], dtype=object)
|
|
|
|
|
#3 Start computing dtaa= calc.Ice_quant(data)
<xarray.Dataset> Dimensions: (t: 31) Coordinates: * t (t) object 2004-07-01 12:00:00 ... 2004-07-31 12:00:00 Data variables: Ice volume (t) float64 dask.array<chunksize=(1,), meta=np.ndarray> Ice area (t) float64 dask.array<chunksize=(1,), meta=np.ndarray> Ice extent (t) float64 dask.array<chunksize=(1,), meta=np.ndarray> Ice drift (t) float64 dask.array<chunksize=(1,), meta=np.ndarray>
array([cftime.DatetimeNoLeap(2004, 7, 1, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 2, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 3, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 4, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 5, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 6, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 7, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 8, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 9, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 10, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 11, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 12, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 13, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 14, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 15, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 16, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 17, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 18, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 19, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 20, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 21, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 22, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 23, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 24, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 25, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 26, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 27, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 28, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 29, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 30, 12, 0, 0, 0), cftime.DatetimeNoLeap(2004, 7, 31, 12, 0, 0, 0)], dtype=object)
|
|
|
|
#4 Saving SEDNA_Ice_intquant_BBFG_Ice_quantities dtaa=save.datas(data,plot=Plot,path=nc_outputpath,filename=filename) start saving data saving data in a csv file ../nc_results/SEDNA_ALPHA_MONITOR/23/SEDNA_Ice_intquant_BBFG_Ice_quantities2004-07-01_2004-07-31.nc save computed data at ../nc_results/SEDNA_ALPHA_MONITOR/23/SEDNA_Ice_intquant_BBFG_Ice_quantities2004-07-01_2004-07-31.nc completed CPU times: user 1h 38min 31s, sys: 1min 35s, total: 1h 40min 7s Wall time: 1h 39min 23s