| 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. |
Parent Directory
| 
Revision Log
| 
Revision Graph