1 |
|
2 |
|
3 |
==========to download this package and test run it as follows======== |
4 |
|
5 |
1) Download and setup m_map and gcmfaces toolbox |
6 |
|
7 |
2) Get MITprof toolbox |
8 |
setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
9 |
cvs login |
10 |
( enter the CVS password: "cvsanon" ) |
11 |
cvs co -d (profilesMatlabProcessing) MITgcm_contrib/gael/profilesMatlabProcessing |
12 |
|
13 |
3) Download climatological fields, and put them in the directory (profilesMatlabProcessing)/climatology |
14 |
or set directory names in profiles_prep_load_fields.m in (profilesMatlabProcessing)/profiles_process_main_v2 |
15 |
to indicate where climatological fields must be found. |
16 |
|
17 |
4) Initialize and test the configuration |
18 |
cd profilesMatlabProcessing |
19 |
matlab -nojvm |
20 |
( type profiles_process_init) |
21 |
|
22 |
=======then you can process your own data as follow (see <= signs)=== |
23 |
|
24 |
profiles_prep_select.m |
25 |
- Encapsulate the data set information in the 'dataset' structure |
26 |
- Directories, included variables, etc. are specified here. |
27 |
=> You append specs for your data set here. Included examples: argo and wod05. |
28 |
|
29 |
profiles_prep_load_fields.m |
30 |
- Load T-S atlases and variance atlas that will be used to detect |
31 |
outliers and to attribute least-squares weights (for ECCO purpose) |
32 |
=> You need to get gcmfaces package and binary files for this part. |
33 |
|
34 |
profiles_read_wod05.m/profiles_read_argo.m |
35 |
- Encapsulate the that loads the data set in memory. |
36 |
- The data set file format specifics are handled here. |
37 |
- The calling sequence must be [varargout]=profile_read_wod05(dataset,nf,m); |
38 |
where 'dataset' is the aformentioned structure, nf is the current file |
39 |
index (relative to dataset.fileInList). If m is set to 0, then varargout |
40 |
must be the number of profiles/stations. If m is a profile index, then varargout |
41 |
must be the profileCur structure that contains the current profile/station. |
42 |
=> You need to make one of those routines for your data set. The two |
43 |
profiles_read_wod05.m/profiles_read_argo.m provide you the blueprint. |
44 |
|
45 |
profiles_prep_convert.m |
46 |
- does conversions of pressure->depth, and insituT->potT, if needed. |
47 |
- profiles_prep_select.m specifies needed conversions in 'dataset'. |
48 |
- also switch [0 360] lon -> [-180 +180] lon, if needed. |
49 |
profiles_prep_locate.m |
50 |
- get the grid point index of the profiles (on the atlases grid). |
51 |
<= I need to upload gcmfaces_bindata |
52 |
|
53 |
profiles_prep_interp.m |
54 |
- interpolate to standard depth levels, as spec. in profiles_prep_select.m. |
55 |
|
56 |
profiles_prep_weights.m |
57 |
- compute instrumental + representation error profile. |
58 |
- instrumental error estimates (when available) are |
59 |
specified at read time (see profiles_read_argo.m). |
60 |
|
61 |
profiles_prep_tests.m |
62 |
- quality control/outlier detection. |
63 |
|
64 |
profiles_prep_write_nc.m |
65 |
store/write results to temporary .mat files, then final .nc file. |