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