-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathcrocotools_param.m
127 lines (114 loc) · 3.29 KB
/
crocotools_param.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
RUN_dir='/data2/matlab/CROCOv1p3/croco/PACIFICO_CENTRAL/';
FORC_DATA_DIR = [getenv('RUN_dir') './DATA/'];
CROCO_config='PACIFICO_CENTRAL';
CROCO_files_dir=[RUN_dir,'CROCO_FILES/'];
DATADIR=['/data2/matlab/CROCOv1p2p1/croco_tools/'];
topofile = [DATADIR,'Topo/etopo2.nc'];
blkname=[CROCO_files_dir,'croco_blk.nc'];
grdname=[CROCO_files_dir,'croco_grd.nc'];
frcname=[CROCO_files_dir,'croco_frc.nc'];
clmname=[CROCO_files_dir,'croco_clm.nc'];
ininame=[CROCO_files_dir,'croco_ini.nc'];
oaname=[CROCO_files_dir,'croco_oa.nc'];
blk_prefix=[CROCO_files_dir,'croco_blk']; % bulk file name
coads_time=(15:30:345);
coads_cycle=360;
woa_time=(15:30:345);
woa_cycle=360;
level=0;
Roa=0;
makeplot=0;
makeclim=0;
makeini=1;
makeoa=1;
makebry=1;
makeZbry=1;
insitu2pot=1;
bryname=[CROCO_files_dir,'croco_bry.nc'];
Zbryname=[CROCO_files_dir,'croco_Zbry.nc'];
dl=1/37;
lonmin=-120.00;
lonmax=-114.00;
latmin=7.50;
latmax=13.00;
theta_s=5;
theta_b=2;
hc=300;
hmin=20;
N=100;
tini=0;
vtransform =2;
obc=[1 1 1 1];
zref=-1000;
interp_method = 'cubic';
coads_dir='/data2/matlab/CROCOv1p2p1/croco_tools/COADS05/';
climato_dir='/data2/matlab/CROCOv1p2p1/croco_tools/WOA2009/';
CROCO_title='Pacific';
Yorig = 2000; % reference time for vector time
% in croco initial and forcing files
%
Ymin = 2022; % first forcing year
Ymax = 2022; % last forcing year
Mmin = 3; % first forcing month
Mmax = 3; % last forcing month
%
Dmin = 1; % Day of initialization
Hmin = 0; % Hour of initialization
Min_min = 0; % Minute of initialization
Smin = 0; % Second of initialization
%
SPIN_Long = 0; % SPIN-UP duration in Years
%
Mth_format = '%02d'; % Number of digits for month on input files
%
itolap_a=2;
itolap_p=2;
Download_data=1; % Habilitar la descarga de datos
hmax=5000;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 5 - Parameters for tidal forcing
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% TPXO file name (TPXO6 or TPXO7)
%
tidename=[DATADIR,'TPXO7/TPXO7.nc'];
%
% Self-Attraction and Loading GOT99.2 file name
%
sal_tides=1;
salname=[DATADIR,'GOT99.2/GOT99_SAL.nc'];
%
% Number of tides component to process
%
Ntides=10;
%
% Chose order from the rank in the TPXO file :
% "M2 S2 N2 K2 K1 O1 P1 Q1 Mf Mm"
% " 1 2 3 4 5 6 7 8 9 10"
%
tidalrank=[1 2 3 4 5 6 7 8 9 10];
add_tides_fcst = 1; % 1: add tides
% OGCM
OGCM='hycom';
if strcmp(OGCM,'hycom')
hdays=2;
fdays=3;
user = '';
password = '';
OGCM_dir = [FORC_DATA_DIR, 'HYCOM_', CROCO_config, '/'];
bry_prefix = [CROCO_files_dir, 'croco_bry_hycom_'];
clm_prefix = [CROCO_files_dir, 'croco_clm_hycom_'];
ini_prefix = [CROCO_files_dir, 'croco_ini_hycom_'];
OGCM_prefix = 'hycom_';
end
%--------------------------------------------------
% Options for make_ERA5
%--------------------------------------------------
%
ERA5_dir = [FORC_DATA_DIR,'ERA5_',CROCO_config,'/'];
My_ERA5_dir = [FORC_DATA_DIR,'ERA5/'];
itolap_era5 = 2;
%
netcdf.setDefaultFormat('FORMAT_NETCDF4');