diff --git a/Forecast_tools/download_mercator_frcst_python.m b/Forecast_tools/download_mercator_frcst_python.m
index 0402779f8c60e4cb3eb7e10052f18a4f5c00d61d..324c5a576ce95be48021810e29026ffcec5c2770 100644
--- a/Forecast_tools/download_mercator_frcst_python.m
+++ b/Forecast_tools/download_mercator_frcst_python.m
@@ -1,10 +1,9 @@
-function mercator_name=download_mercator_frcst_python(pathMotu,user,password,mercator_type, ...
- motu_url, service_id,product_id, ...
+function mercator_name=download_mercator_frcst_python(pathCMC,user,password,mercator_type, ...
+ product_id, ...
lh,lf, ...
lonmin,lonmax,latmin,latmax,zmax, ...
FRCST_dir,FRCST_prefix,raw_mercator_name,Yorig)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%
% Extract a subgrid from mercator to get a CROCO forcing
% Store that into monthly files.
% Take care of the Greenwitch Meridian.
@@ -42,6 +41,9 @@ function mercator_name=download_mercator_frcst_python(pathMotu,user,password,mer
download_raw_data=1;
convert_raw2crocotools=1; % convert -> crocotools format data
%
+disp(['Making output data directory ',FRCST_dir]) % create directory
+eval(['!mkdir ',FRCST_dir])
+%
% Set variable names according to mercator type data
%
if mercator_type==1,
@@ -90,54 +92,42 @@ disp([' '])
if download_raw_data
%
% Get data (now in 4 different files)
- %
- % temporarily removing Matlab lib path to avoid conflict with Python,
- % (mandatory with python 2.7.X)
- %
- pathdyld=getenv('DYLD_LIBRARY_PATH');
- setenv('DYLD_LIBRARY_PATH','');
- pathld=getenv('LD_LIBRARY_PATH');
- setenv('LD_LIBRARY_PATH','');
-
- get_file_python_mercator(pathMotu,mercator_type, ... % SSH
- motu_url,service_id,product_id{1}, ...
- vars(1), ...
- [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
- {datestr(tiempo_inicial,'yyyy-mm-dd') ...
- datestr(tiempo_final, 'yyyy-mm-dd')}, ...
- {user password}, ...
- [raw_mercator_name(1:end-3),'_z.nc']);
-
- get_file_python_mercator(pathMotu,mercator_type, ... % U,V
- motu_url,service_id,product_id{2}, ...
- vars(2:3), ...
- [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
- {datestr(tiempo_inicial,'yyyy-mm-dd') ...
- datestr(tiempo_final, 'yyyy-mm-dd')}, ...
- {user password}, ...
- [raw_mercator_name(1:end-3),'_u.nc']);
-
- get_file_python_mercator(pathMotu,mercator_type, ... % T
- motu_url,service_id,product_id{3}, ...
- vars(4), ...
- [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
- {datestr(tiempo_inicial,'yyyy-mm-dd') ...
- datestr(tiempo_final, 'yyyy-mm-dd')}, ...
- {user password}, ...
- [raw_mercator_name(1:end-3),'_t.nc']);
-
- get_file_python_mercator(pathMotu,mercator_type, ... % S
- motu_url,service_id,product_id{4}, ...
- vars(5), ...
- [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
- {datestr(tiempo_inicial,'yyyy-mm-dd') ...
- datestr(tiempo_final, 'yyyy-mm-dd')}, ...
- {user password}, ...
- [raw_mercator_name(1:end-3),'_s.nc']);
-
- setenv('DYLD_LIBRARY_PATH',pathdyld); % set back Matlab libs path
- setenv('LD_LIBRARY_PATH',pathld); % set back Matlab libs path
-
+ get_file_python_mercator(pathCMC,mercator_type, ... %SSH
+ product_id{1}, ...
+ vars(1), ...
+ [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
+ {datestr(tiempo_inicial,'yyyy-mm-dd') ...
+ datestr(tiempo_final, 'yyyy-mm-dd')}, ...
+ {user password}, ...
+ [raw_mercator_name(1:end-3),'_z.nc']);
+
+ get_file_python_mercator(pathCMC,mercator_type, ... %U, V
+ product_id{2}, ...
+ vars(2:3), ...
+ [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
+ {datestr(tiempo_inicial,'yyyy-mm-dd') ...
+ datestr(tiempo_final, 'yyyy-mm-dd')}, ...
+ {user password}, ...
+ [raw_mercator_name(1:end-3),'_u.nc']);
+
+ get_file_python_mercator(pathCMC,mercator_type, ... % T
+ product_id{3}, ...
+ vars(4), ...
+ [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
+ {datestr(tiempo_inicial,'yyyy-mm-dd') ...
+ datestr(tiempo_final, 'yyyy-mm-dd')}, ...
+ {user password}, ...
+ [raw_mercator_name(1:end-3),'_t.nc']);
+
+ get_file_python_mercator(pathCMC,mercator_type, ... % S
+ product_id{4}, ...
+ vars(5), ...
+ [lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
+ {datestr(tiempo_inicial,'yyyy-mm-dd') ...
+ datestr(tiempo_final, 'yyyy-mm-dd')}, ...
+ {user password}, ...
+ [raw_mercator_name(1:end-3),'_s.nc']);
+
end % download_raw_data
if convert_raw2crocotools
@@ -145,9 +135,6 @@ if convert_raw2crocotools
% Convert data format and write in a more CROCOTOOLS
% compatible input file
%
- disp(['Making output data directory ',FRCST_dir]) % create directory
- eval(['!mkdir ',FRCST_dir])
- %
mercator_name=[FRCST_dir,FRCST_prefix,num2str(rundate),'.cdf'];
if exist(mercator_name)
disp('Mercator file already exist => overwrite it')
diff --git a/Forecast_tools/make_OGCM_frcst.m b/Forecast_tools/make_OGCM_frcst.m
index 3c11f8e40125b08cb6a518dc362f948d6164d432..a9bdac3e78f9e54a2a61898c6b4d51bca7ff7cdb 100644
--- a/Forecast_tools/make_OGCM_frcst.m
+++ b/Forecast_tools/make_OGCM_frcst.m
@@ -33,14 +33,14 @@
% Updated 14-Oct-2020 by P. Marchesiello
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%clear all
-%close all
+clear all
+close all
tic
%%%%%%%%%%%%%%%%%%%%% USERS DEFINED VARIABLES %%%%%%%%%%%%%%%%%%%%%%%%
%
% Common parameters
%
-crocotools_param
+crocotools_param;
%
% Specific to forecast
%
@@ -95,7 +95,6 @@ pm=nc{'pm'}(:);
pn=nc{'pn'}(:);
rmask=nc{'mask_rho'}(:);
close(nc)
-
%---------------------------------------------------------------
% Extract data from the Internet
%---------------------------------------------------------------
@@ -121,9 +120,15 @@ if Download_data
% Use Motu python
%
disp('Download data...')
+ % eval(['OGCM_name=download_', ...
+ % OGCM,'_frcst_python(pathMotu,user,password,mercator_type,', ...
+ % 'motu_url_fcst,service_id_fcst,product_id_fcst,', ...
+ % 'hdays,fdays,', ...
+ % 'lonmin,lonmax,latmin,latmax,hmax,', ...
+ % 'FRCST_dir,FRCST_prefix,raw_mercator_name,Yorig);'])
eval(['OGCM_name=download_', ...
- OGCM,'_frcst_python(pathMotu,user,password,mercator_type,', ...
- 'motu_url_fcst,service_id_fcst,product_id_fcst,', ...
+ OGCM,'_frcst_python(pathCMC,user,password,mercator_type,', ...
+ 'product_id_fcst,', ...
'hdays,fdays,', ...
'lonmin,lonmax,latmin,latmax,hmax,', ...
'FRCST_dir,FRCST_prefix,raw_mercator_name,Yorig);'])
@@ -217,7 +222,7 @@ for tndx=1:length(time)
disp([' Time step : ',num2str(tndx),' of ',num2str(length(time)),' :'])
interp_OGCM_frcst(OGCM_name,Roa,interp_method,...
lonU,latU,lonV,latV,lonT,latT,Z,trange(tndx),...
- nc_clm,nc_bry,lon,lat,angle,h,pm,pn,rmask, ...
+ nc_clm,nc_bry,lon,lat,angle,h,pm,pn,rmask, ...
tndx,vtransform,obc)
end
diff --git a/Oforc_OGCM/Copernicus_Marine_Toolbox_installation.md b/Oforc_OGCM/Copernicus_Marine_Toolbox_installation.md
new file mode 100644
index 0000000000000000000000000000000000000000..484c5f3316631c2117a3932ab538993abc80beed
--- /dev/null
+++ b/Oforc_OGCM/Copernicus_Marine_Toolbox_installation.md
@@ -0,0 +1,15 @@
+# Copernicus Marine Toolbox installation
+
+The Copernicus Marine Toolbox is used to download and extract data from the Copernicus Marine Data Store (Global Ocean Physics Reanalysis, Global Ocean Physics Analysis and Forecast, ...) to create oceanic open boundary and initial conditions. [ https://data.marine.copernicus.eu/products ]
+
+All the instruction can be found here https://help.marine.copernicus.eu/en/collections/4060068-copernicus-marine-toolbox. It require to install a dedicated conda env, by default named cmt_1.0.
+Note that you can use micromamba instead of mamba to install the python environement.
+
+To use it , you need to activate the cmt_1.0 : `mamba activate cmt_1.0`
+
+You will have access to the copernicusmarine executable with various sub-command, very useful to get, extract and download data from the Copernicus Marine Data Store.
+
+To be used in the Matlab croco_tools, in the crocotools_param.m, you need to :
+
+- define the path to the copernicusmarine executable (`pathCMS`) in crocotools_param.m at *section 7, Option for make_OGCM or make_OGCM_mercator*
+- define your copernicus marine `login` and `password` in crocotools_param.m at *section 7, Options for for make_OGCM or make_OGCM_mercator*
diff --git a/Oforc_OGCM/README b/Oforc_OGCM/README
deleted file mode 100644
index e6c1fb9025401f548d2fd508d1c2042509bafab5..0000000000000000000000000000000000000000
--- a/Oforc_OGCM/README
+++ /dev/null
@@ -1,24 +0,0 @@
-==
-make_OGCM.m is a way to get initial and boundary conditions
-for CROCO simulations from different global oceanic models.
-for the moment it works with SODA and ECCO.
-it uses DODS to extract the subgrids.
-
-Pierrick 2006/1/19.
-
-==
-Update:
-Now it use the ECCO2 data, every 3 days
-Gildas 2014/2/14
-
-==
-24-11-2020 : Update by Gildas Cambon
-Now you can process oceanic frocing from the global daily glorys12 reanalysis
-at 1/12 degree resolution.
-It is distributed by copernicus/mecator => https://resources.marine.copernicus.eu/
-
-It use the python-motu client to download a space and time
-extraction of the reanalysis.
-Have a look at the routines :
-- download_mercator_frcst_python.m
-- get_file_python_mercator.m
diff --git a/Oforc_OGCM/README.md b/Oforc_OGCM/README.md
new file mode 100644
index 0000000000000000000000000000000000000000..31f4a7caa65eb28e987352656969cbe0d0a60ef5
--- /dev/null
+++ b/Oforc_OGCM/README.md
@@ -0,0 +1,20 @@
+## 23-02-2024 : Update to use the Copernicus Marine Toolbox
+
+ - see Copernicus_Marine_Toolbox_installation.md
+
+## 24-11-2020 : Update by Gildas Cambon
+Now you can process oceanic forcing from the global daily glorys12 reanalysis
+at 1/12 degree resolution.
+It is distributed by copernicus/mecator => https://resources.marine.copernicus.eu/
+
+It use the python-motu client to download a space and time
+extraction of the reanalysis.
+Have a look at the routines :
+- download_mercator_frcst_python.m
+- get_file_python_mercator.m
+
+## 14-02-2014 : Gildas Cambon
+Now it use the ECCO2 data, every 3 days
+
+## 01-19-2006 : Pierrick Penven
+make_OGCM.m is a way to get initial and boundary conditions. For CROCO simulations from different global oceanic models. For the moment it works with SODA and ECCO. it uses DODS to extract the subgrids
diff --git a/Oforc_OGCM/download_mercator_python.m b/Oforc_OGCM/download_mercator_python.m
index eb98f398503619748aaeda2822c4cda7e63cea6c..6eada2b9d0847b726726189c04a7cd0bf5f0c4ce 100644
--- a/Oforc_OGCM/download_mercator_python.m
+++ b/Oforc_OGCM/download_mercator_python.m
@@ -1,5 +1,5 @@
-function download_mercator_python(pathMotu,user,password,mercator_type, ...
- motu_url,service_id,product_id, ...
+function download_mercator_python(pathCMC,user,password,mercator_type, ...
+ product_id, ...
Y,M, ...
lonmin,lonmax,latmin,latmax,zmax, ...
OGCM_dir,OGCM_prefix,thedatemonth,Yorig)
@@ -37,7 +37,6 @@ function download_mercator_python(pathMotu,user,password,mercator_type, ...
% Updated 20-Aug-2008 by Matthieu Caillaud & P. Marchesiello
% Update 23-Oct-2020 by Gildas Cambon
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% pathMotu is a deprecated parameter !
download_raw_data=1;
convert_raw2crocotools=1; % convert -> crocotools format data
%
@@ -84,28 +83,16 @@ raw_mercator_name=[OGCM_dir,'raw_motu_',OGCM_prefix,thedatemonth,'.nc'];
if download_raw_data
%
- % Get data
- % (temporarily removing Matlab lib path to avoid conflict with Python, mandatory with python 2.7.X)
- % for example problem like :
- % Execution failed: /usr/lib/python2.7/lib-dynload/pyexpat.x86_64-linux-gnu.so:
- % undefined symbol: XML_SetHashSalt" )
-
- pathdyld=getenv('DYLD_LIBRARY_PATH');
- setenv('DYLD_LIBRARY_PATH','');
- pathld=getenv('LD_LIBRARY_PATH');
- setenv('LD_LIBRARY_PATH','');
-
- get_file_python_mercator(pathMotu,mercator_type, ...
- motu_url,service_id,product_id, ...
+ % Get data
+ get_file_python_mercator(pathCMC,mercator_type, ...
+ product_id, ...
vars, ...
[lonmin-1 lonmax+1 latmin-1 latmax+1 0 zmax], ...
{datestr(tiempo_inicial,'yyyy-mm-dd') ...
datestr(tiempo_final, 'yyyy-mm-dd')}, ...
{user password}, ...
- raw_mercator_name);
-
- setenv('DYLD_LIBRARY_PATH',pathdyld); % set back Matlab libs path
- setenv('LD_LIBRARY_PATH',pathld); % set back Matlab libs path
+ raw_mercator_name);
+
end %end download_raw_data
diff --git a/Oforc_OGCM/get_file_python_mercator.m b/Oforc_OGCM/get_file_python_mercator.m
index 4320a577d1b9c49c84d543eda371e3afd694eb2d..6c27e53231df036d6d32bdc3161ba4a5d90d37b4 100644
--- a/Oforc_OGCM/get_file_python_mercator.m
+++ b/Oforc_OGCM/get_file_python_mercator.m
@@ -1,7 +1,8 @@
-function outname = get_file_python_mercator(pathmotu,mercator_type,...
- motu_url,service_id,product_id, ...
- vars, ...
- geom,date,info,outname)
+function outname = get_file_python_mercator(pathCMC,...
+ mercator_type,product_id,...
+ vars,...
+ geom,date,...
+ info,outname)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Extract a subgrid from MERCATOR to get a CROCO forcing
@@ -33,115 +34,31 @@ function outname = get_file_python_mercator(pathmotu,mercator_type,...
%
% Updated 12-Feb-2016 by P. Marchesiello
% Updated 26-Nov-2020 by G. Cambon
+% Updated 22-Feb-2024 by G. Cambon
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Usage:
-%=======
-% get_file_python_mercator(pathMotu,mercator_type,...
-% motu_url,service_id,product_id, ...
-% {'varname1' 'varname2'}, ...
-% [lonmin lonmax latmin latmax, depthmin depthmax], ...
-% {'startdate' 'enddate'},{user password},
-% 'outname.nc');
-%
-% For GLORYS 12 reanalysis
-% ========================
-% motu_url='http://my.cmems-du.eu/motu-web/Motu';
-% service_id='GLOBAL_REANALYSIS_PHY_001_030-TDS';
-% product_id='global-reanalysis-phy-001-030-daily';
-%
-% For Mercator 1/12 forecast
-% ==========================
-% motu_url='http://nrt.cmems-du.eu/motu-web/Motu';
-% service_id='GLOBAL_ANALYSIS_FORECAST_PHY_001_024-TDS';
-% product_id='global-analysis-forecast-phy-001-024';
-%
-% For Met-Office 1/4 forecast
-% ===========================
-% motu_url='http://nrtcmems.mercator-ocean.fr/motu-web/Motu';
-% service_id='GLOBAL_ANALYSIS_FORECAST_PHYS_001_015-TDS';
-% product_id='global-analysis-forecast-phys-001-015';
-%
-% Check http://marine.copernicus.eu/web/34-products-and-services-faq.php
-% for changes in the command line: addresses, file names, variable names ...
-%
-% Currently needs Python2.7.X (with X >= 10) or Python 3.X (with >=4)
-% and motuclient.py v.1.8.X (where "X" is equal or higher to "4")
-% from https://marine.copernicus.eu/services-portfolio/technical-faq/#motu-client
-%
-% python
-% ======
-% Type:
-% python --version
-% It should return: "Python 2.7.10+" or "Python 3.4+".
-%
-% motuclient
-% ==========
-% 1- Use croco's motuclient (pathMotu => Forecast_tools) or ...
-% 2- Install your own motuclient:
-% To update/upgrade and get the latest version of motuclient
-% from a previous version (<= v1.8.3), type in the following:
-% $ python -m pip install --upgrade motuclient
-% Otherwise (if there is no previous installation of motuclient),
-% type in the following:
-% $ python -m pip install motuclient
-% It should install and display motuclient package v1.8.4 (Oct. 2019).
-% To make sure, display the version:
-% $ python -m motuclient --version
-% If it does not return: "motuclient-python v1.8.X" ("X" >= "4"),
-% then type in the following:
-% $ python -m pip install motuclient==1.8.4
-%=======
+% Informations :
+% => https://help.marine.copernicus.eu/en/collections/4060068-copernicus-marine-toolbox
+% => https://help.marine.copernicus.eu/en/articles/8206131-how-to-download-a-copernicus-marine-toolbox-request-in-matlab
%
eval(['! rm ',outname])
+command = {sprintf('%s',pathCMC)
+ sprintf(' subset')
+ sprintf(' --username %s --password %s',info{1},info{2})
+ sprintf(' -i %s',product_id)
+ sprintf(' -t %s -T %s',['"',date{1},'"'],['"',date{2},'"'])
+ sprintf(' -x %f -X %f',geom(1),geom(2))
+ sprintf(' -y %f -Y %f',geom(3),geom(4))
+ sprintf(' -z %f -Z %f',geom(5),geom(6))
+ sprintf(' -o ./')
+ sprintf(' -f %s',outname)
+ sprintf(' --force-download')};
-if isempty(pathmotu), % use your own motuclient
-
- command = {'!python -m motuclient'
- sprintf(' -u %s -p %s',info{1},info{2})
- sprintf(' -m %s',motu_url)
- sprintf(' -s %s',service_id)
- sprintf(' -d %s',product_id)
- sprintf(' -t %s -T %s',date{1},date{2})
- sprintf(' -x %f -X %f',geom(1),geom(2))
- sprintf(' -y %f -Y %f',geom(3),geom(4))
- sprintf(' -z %f -Z %f',geom(5),geom(6))
- sprintf(' -o ./')
- sprintf(' --out-name %s',outname)};
-
-else % use croco's motuclient
-
- command = {'!'
- sprintf('%s',pathmotu)
- 'motuclient-python/motuclient.py'
- sprintf(' -u %s -p %s',info{1},info{2})
- sprintf(' -m %s',motu_url)
- sprintf(' -s %s',service_id)
- sprintf(' -d %s',product_id)
- sprintf(' -t %s -T %s',date{1},date{2})
- sprintf(' -x %f -X %f',geom(1),geom(2))
- sprintf(' -y %f -Y %f',geom(3),geom(4))
- sprintf(' -z %f -Z %f',geom(5),geom(6))
- sprintf(' -o ./')
- sprintf(' --out-name %s',outname)};
-end
-
for k =1:length(vars)
- command{end+1}=sprintf(' -v %s ',vars{k});
+ command{end+1}=sprintf(' -v %s',vars{k});
end
-%===
-% If using a proxy server, uncomment the following line and replace by your
-% proxy url and port server. Beware that a SPACE is needed between
-% "--proxy-server=" and "the proxy-server-name" !
-%
-%command{end+2}=sprintf('--proxy-server= http://your_proxy_server:your_proxy_port');
-%===
disp([command{:}])
-eval([command{:}])
+system([command{:}]);
-%
-% example of a python motuclient commands
-%
-%python -m motuclient --motu http://my.cmems-du.eu/motu-web/Motu --service-id GLOBAL_REANALYSIS_PHY_001_030-TDS --product-id global-reanalysis-phy-001-030-daily --longitude-min 15 --longitude-max 20 --latitude-min -40 --latitude-max -25 --date-min "2004-12-31 12:00:00" --date-max "2005-02-01 12:00:00" --depth-min 0.493 --depth-max 5727.918 --variable so --variable thetao --variable uo --variable vo --variable zos --out-dir <OUTPUT_DIRECTORY> --out-name <OUTPUT_FILENAME> --user <USERNAME> --pwd <PASSWORD>
end
diff --git a/Oforc_OGCM/make_OGCM_mercator.m b/Oforc_OGCM/make_OGCM_mercator.m
index 1619d15f8212497dd5b874344065d7cfd2e9b291..85cea3fe6a195258d9f1badbe4894e97f49ba6c6 100644
--- a/Oforc_OGCM/make_OGCM_mercator.m
+++ b/Oforc_OGCM/make_OGCM_mercator.m
@@ -111,8 +111,8 @@ if Download_data==1
for M=mo_min:mo_max
%motu_url_reana
thedatemonth=['Y',num2str(Y),'M',num2str(M)];
- download_mercator_python(pathMotu,user,password,1,...
- motu_url_reana,service_id_reana,product_id_reana, ...
+ download_mercator_python(pathCMC,user,password,1,...
+ product_id_reana, ...
Y,M, ...
lonmin,lonmax,latmin,latmax,hmax, ...
OGCM_dir,OGCM_prefix,thedatemonth,Yorig)
diff --git a/crocotools_param.m b/crocotools_param.m
index 3cd89505822980b050118dcedf15b45c769dea2c..9e18b3a99efebd0194f409a2d2973ebd610b77d2 100644
--- a/crocotools_param.m
+++ b/crocotools_param.m
@@ -433,9 +433,16 @@ if strcmp(OGCM,'mercator')
% ========================
% For GLORYS 12 reanalysis extraction + download using python motuclient
% ========================
+ % To download CMEMS data: set login/password (http://marine.copernicus.eu) and path to copernicusmarine executable
+ % see Oforc_OGCM/Copernicus_Marine_Toolbox_installation.md
+
+ pathCMC = '/path/to/home/copernicusmarine';
+ %
+ user = 'XXXX';
+ password = 'XXXX';
+ %
raw_mercator_name = 'mercator';
kdata='MONTHLY' ; % or DAILY
- motu_url_reana='http://my.cmems-du.eu/motu-web/Motu';
service_id_reana='GLOBAL_MULTIYEAR_PHY_001_030-TDS';
if strcmp(kdata,'DAILY')
product_id_reana='cmems_mod_glo_phy_my_0.083_P1D-m';
@@ -498,30 +505,19 @@ add_tides_fcst = 1; % 1: add tides
%
% MERCATOR case:
% =============
-% To download CMEMS data: set login/password (http://marine.copernicus.eu)
-% and path to croco's motuclient python package;
-% or set pathMotu='' (empty) to use your own motuclient
%
% Various sets of data are proposed in the Copernicus web site
%
if strcmp(OGCM,'mercator')
- user = 'XXX';
- password = 'XXX';
-
- pathMotu =[CROCOTOOLS_dir,'Forecast_tools/'];
mercator_type=1; % 1 --> 1/12 deg Mercator forecast
% 2 --> high-resolution regional grids
if mercator_type==1
- motu_url_fcst='http://nrt.cmems-du.eu/motu-web/Motu';
- service_id_fcst='GLOBAL_ANALYSISFORECAST_PHY_001_024-TDS';
product_id_fcst={'cmems_mod_glo_phy_anfc_0.083deg_P1D-m', ...
'cmems_mod_glo_phy-cur_anfc_0.083deg_P1D-m', ...
'cmems_mod_glo_phy-thetao_anfc_0.083deg_P1D-m', ...
'cmems_mod_glo_phy-so_anfc_0.083deg_P1D-m'};
elseif mercator_type==2
- motu_url_fcst='http://nrt.cmems-du.eu/motu-web/Motu';
- service_id_fcst='MEDSEA_ANALYSISFORECAST_PHY_006_013-TDS';
product_id_fcst={'cmems_mod_med_phy-ssh_anfc_4.2km_P1D-m', ...
'cmems_mod_med_phy-cur_anfc_4.2km_P1D-m', ...
'cmems_mod_med_phy-tem_anfc_4.2km_P1D-m', ...