1 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
---|
2 | % |
---|
3 | % Build a ROMS initial file from Levitus Data |
---|
4 | % |
---|
5 | % Extrapole and interpole temperature and salinity from a |
---|
6 | % Climatology to get initial conditions for |
---|
7 | % ROMS (initial netcdf files) . |
---|
8 | % Get the velocities and sea surface elevation via a |
---|
9 | % geostrophic computation. |
---|
10 | % |
---|
11 | % Data input format (netcdf): |
---|
12 | % temperature(T, Z, Y, X) |
---|
13 | % T : time [Months] |
---|
14 | % Z : Depth [m] |
---|
15 | % Y : Latitude [degree north] |
---|
16 | % X : Longitude [degree east] |
---|
17 | % |
---|
18 | % Data source : IRI/LDEO Climate Data Library (World Ocean Atlas 1998) |
---|
19 | % http://ingrid.ldgo.columbia.edu/ |
---|
20 | % http://iridl.ldeo.columbia.edu/SOURCES/.NOAA/.NODC/.WOA98/ |
---|
21 | % |
---|
22 | % Pierrick Penven, IRD, 2005. |
---|
23 | % Olivier Aumont, IRD, 2006. |
---|
24 | % Patricio Marchesiello, IRD 2007 |
---|
25 | % Christophe Eugene Raoul Menkes, IRD 2007 |
---|
26 | % |
---|
27 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
---|
28 | clear all |
---|
29 | close all |
---|
30 | %%%%%%%%%%%%%%%%%%%%% USERS DEFINED VARIABLES %%%%%%%%%%%%%%%%%%%%%%%% |
---|
31 | % |
---|
32 | % Title |
---|
33 | % |
---|
34 | romstools_param |
---|
35 | % |
---|
36 | % Data climatologies file names: |
---|
37 | % |
---|
38 | no3_month_data = [woapisces_dir,'no3_seas.cdf']; |
---|
39 | no3_ann_data = [woapisces_dir,'no3_ann.cdf']; |
---|
40 | po4_month_data = [woapisces_dir,'po4_seas.cdf']; |
---|
41 | po4_ann_data = [woapisces_dir,'po4_ann.cdf']; |
---|
42 | o2_month_data = [woapisces_dir,'o2_seas.cdf']; |
---|
43 | o2_ann_data = [woapisces_dir,'o2_ann.cdf']; |
---|
44 | sio3_month_data = [woapisces_dir,'sio3_seas.cdf']; |
---|
45 | sio3_ann_data = [woapisces_dir,'sio3_ann.cdf']; |
---|
46 | dic_month_data = [woapisces_dir,'dic_seas.cdf']; |
---|
47 | dic_ann_data = [woapisces_dir,'dic_ann.cdf']; |
---|
48 | talk_month_data = [woapisces_dir,'talk_seas.cdf']; |
---|
49 | talk_ann_data = [woapisces_dir,'talk_ann.cdf']; |
---|
50 | doc_month_data = [woapisces_dir,'doc_seas.cdf']; |
---|
51 | doc_ann_data = [woapisces_dir,'doc_ann.cdf']; |
---|
52 | fer_month_data = [woapisces_dir,'fer_seas.cdf']; |
---|
53 | fer_ann_data = [woapisces_dir,'fer_ann.cdf']; |
---|
54 | dust_month_data = [woapisces_dir,'dust_seas.cdf']; |
---|
55 | dust_ann_data = [woapisces_dir,'dust_ann.cdf']; |
---|
56 | % |
---|
57 | % |
---|
58 | %%%%%%%%%%%%%%%%%%% END USERS DEFINED VARIABLES %%%%%%%%%%%%%%%%%%%%%%% |
---|
59 | % |
---|
60 | % Title |
---|
61 | % |
---|
62 | if makeini |
---|
63 | disp(' ') |
---|
64 | disp([' Adding PISCES data into initial file: ',ininame]) |
---|
65 | % |
---|
66 | % Initial file |
---|
67 | % |
---|
68 | add_ini_pisces(ininame,'write'); |
---|
69 | % |
---|
70 | % Horizontal and vertical interp/extrapolations |
---|
71 | % |
---|
72 | disp(' ') |
---|
73 | disp(' Interpolations / extrapolations') |
---|
74 | disp(' ') |
---|
75 | disp('nitrate ...') |
---|
76 | ext_tracers_ini(ininame,grdname,no3_month_data,no3_ann_data,... |
---|
77 | 'nitrate','NO3','r',tini); |
---|
78 | disp(' ') |
---|
79 | disp('phosphate ...') |
---|
80 | ext_tracers_ini(ininame,grdname,po4_month_data,po4_ann_data,... |
---|
81 | 'phosphate','PO4','r',tini); |
---|
82 | disp(' ') |
---|
83 | disp('Silicate ...') |
---|
84 | ext_tracers_ini(ininame,grdname,sio3_month_data,sio3_ann_data,... |
---|
85 | 'silicate','Si','r',tini); |
---|
86 | disp(' ') |
---|
87 | disp(' Oxygen ...') |
---|
88 | ext_tracers_ini(ininame,grdname,o2_month_data,o2_ann_data,... |
---|
89 | 'oxygen','O2','r',tini); |
---|
90 | disp(' ') |
---|
91 | disp('Dissolved Inorganic Carbon ...') |
---|
92 | ext_tracers_ini(ininame,grdname,dic_month_data,dic_ann_data,... |
---|
93 | 'dic','DIC','r',tini); |
---|
94 | disp(' ') |
---|
95 | disp('Total alkalinity ...') |
---|
96 | ext_tracers_ini(ininame,grdname,talk_month_data,talk_ann_data,... |
---|
97 | 'talk','TALK','r',tini); |
---|
98 | disp(' ') |
---|
99 | disp('Dissolved Organic Carbon...') |
---|
100 | ext_tracers_ini(ininame,grdname,doc_month_data,doc_ann_data,... |
---|
101 | 'doc','DOC','r',tini); |
---|
102 | disp(' ') |
---|
103 | disp('Iron ...') |
---|
104 | ext_tracers_ini(ininame,grdname,fer_month_data,fer_ann_data,... |
---|
105 | 'fer','FER','r',tini); |
---|
106 | % |
---|
107 | % Make a few plots |
---|
108 | % |
---|
109 | disp(' ') |
---|
110 | disp(' Make a few plots...') |
---|
111 | test_clim(ininame,grdname,'NO3',1,coastfileplot) |
---|
112 | figure |
---|
113 | test_clim(ininame,grdname,'PO4',1,coastfileplot) |
---|
114 | figure |
---|
115 | test_clim(ininame,grdname,'Si',1,coastfileplot) |
---|
116 | figure |
---|
117 | test_clim(ininame,grdname,'O2',1,coastfileplot) |
---|
118 | figure |
---|
119 | test_clim(ininame,grdname,'DIC',1,coastfileplot) |
---|
120 | figure |
---|
121 | test_clim(ininame,grdname,'TALK',1,coastfileplot) |
---|
122 | figure |
---|
123 | test_clim(ininame,grdname,'DOC',1,coastfileplot) |
---|
124 | figure |
---|
125 | test_clim(ininame,grdname,'FER',1,coastfileplot) |
---|
126 | |
---|
127 | end %if makeini |
---|
128 | % |
---|
129 | % End |
---|
130 | % |
---|
131 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
---|