32 % 11: REMOVE SUPER-ENSEMBLE OUTLIERS |
32 % 11: REMOVE SUPER-ENSEMBLE OUTLIERS |
33 % 12: RE-FORM SUPER ENSEMBLES |
33 % 12: RE-FORM SUPER ENSEMBLES |
34 % 13: (RE-)LOAD SADCP DATA |
34 % 13: (RE-)LOAD SADCP DATA |
35 % 14: CALCULATE INVERSE SOLUTION |
35 % 14: CALCULATE INVERSE SOLUTION |
36 % 15: CALCULATE SHEAR SOLUTION |
36 % 15: CALCULATE SHEAR SOLUTION |
37 % 16: CALCULATE DIFFUSIVITY PROFILE |
37 % 16: PLOT RESULTS & SHOW WARNINGS |
38 % 17: PLOT RESULTS & SHOW WARNINGS |
38 % 17: SAVE OUTPUT |
39 % 18: SAVE OUTPUT |
|
40 |
39 |
41 %====================================================================== |
40 %====================================================================== |
42 % P R O C E S S _ C A S T . M |
41 % P R O C E S S _ C A S T . M |
43 % doc: Thu Jun 24 16:54:23 2004 |
42 % doc: Thu Jun 24 16:54:23 2004 |
44 % dlm: Fri Sep 26 15:57:32 2014 |
43 % dlm: Tue May 12 08:36:53 2015 |
45 % (c) 2004 A.M. Thurnherr |
44 % (c) 2004 A.M. Thurnherr |
46 % uE-Info: 83 80 NIL 0 0 72 2 2 8 NIL ofnI |
45 % uE-Info: 378 0 NIL 0 0 72 2 2 8 NIL ofnI |
47 %====================================================================== |
46 %====================================================================== |
48 |
47 |
49 % NOTES: |
48 % NOTES: |
50 % - changing this function should not be required, except to fix bugs |
49 % - changing this function should not be required, except to fix bugs |
51 % - in order to preserve state variables across load/save calls in matlab, |
50 % - in order to preserve state variables across load/save calls in matlab, |
79 % with CLIVAR P18 leg1 casts 46, 48, 49 & 50 |
78 % with CLIVAR P18 leg1 casts 46, 48, 49 & 50 |
80 % Sep 18, 2008: - BUG: p.navdata was non-existent field when no nav |
79 % Sep 18, 2008: - BUG: p.navdata was non-existent field when no nav |
81 % data were loaded |
80 % data were loaded |
82 % Apr 26, 2012: - finally removed finestructure kz code |
81 % Apr 26, 2012: - finally removed finestructure kz code |
83 % Sep 26, 2014: - added support for p.orig in [saveres.m] (patch by Dan Torres) |
82 % Sep 26, 2014: - added support for p.orig in [saveres.m] (patch by Dan Torres) |
|
83 % May 12, 2015: - finally removed entire step 16 (diffusivity) |
84 |
84 |
85 %---------------------------------------------------------------------- |
85 %---------------------------------------------------------------------- |
86 % STEP 0: EXECUTE ALWAYS |
86 % STEP 0: EXECUTE ALWAYS |
87 %---------------------------------------------------------------------- |
87 %---------------------------------------------------------------------- |
88 |
88 |
366 [di,p]=loadsadcp(f,di,p); |
366 [di,p]=loadsadcp(f,di,p); |
367 end |
367 end |
368 dino=di; |
368 dino=di; |
369 lanarrow |
369 lanarrow |
370 diary on |
370 diary on |
|
371 if existf(d,'ctdprof_p') |
|
372 dr.ctd_t=interp1q(d.ctdprof_z,d.ctdprof_t,dr.z); |
|
373 dr.ctd_s=interp1q(d.ctdprof_z,d.ctdprof_s,dr.z); |
|
374 end |
|
375 if existf(d,'ctdprof_ss') |
|
376 dr.ctd_ss=interp1q(d.ctdprof_z,d.ctdprof_ss,dr.z); |
|
377 end |
371 end |
378 end |
372 |
379 |
373 end_processing_step; |
380 end_processing_step; |
374 end % OF STEP 11: REMOVE SUPER-ENSEMBLE OUTLIERS |
381 end % OF STEP 11: REMOVE SUPER-ENSEMBLE OUTLIERS |
375 |
382 |
457 |
464 |
458 end_processing_step; |
465 end_processing_step; |
459 end % OF STEP 15: CALCULATE SHEAR SOLUTION |
466 end % OF STEP 15: CALCULATE SHEAR SOLUTION |
460 |
467 |
461 %---------------------------------------------------------------------- |
468 %---------------------------------------------------------------------- |
462 % STEP 16: CALCULATE DIFFUSIVITY PROFILE |
469 % STEP 16: PLOT RESULTS & SHOW WARNINGS |
463 %---------------------------------------------------------------------- |
|
464 |
|
465 pcs.cur_step = pcs.cur_step + 1; |
|
466 if pcs.begin_step <= pcs.cur_step |
|
467 pcs.step_name = 'CALCULATE DIFFUSIVITY PROFILE'; begin_processing_step; |
|
468 |
|
469 % |
|
470 % prepare ctd data for output |
|
471 % |
|
472 if existf(d,'ctdprof_p') |
|
473 % first save CTD data with profile |
|
474 dr.ctd_t=interp1q(d.ctdprof_z,d.ctdprof_t,dr.z); |
|
475 dr.ctd_s=interp1q(d.ctdprof_z,d.ctdprof_s,dr.z); |
|
476 end |
|
477 if existf(d,'ctdprof_ss') |
|
478 dr.ctd_ss=interp1q(d.ctdprof_z,d.ctdprof_ss,dr.z); |
|
479 end |
|
480 |
|
481 end_processing_step; |
|
482 end % OF STEP 16: CALCULATE DIFFUSIVITY PROFILE |
|
483 |
|
484 %---------------------------------------------------------------------- |
|
485 % STEP 17: PLOT RESULTS & SHOW WARNINGS |
|
486 %---------------------------------------------------------------------- |
470 %---------------------------------------------------------------------- |
487 |
471 |
488 pcs.cur_step = pcs.cur_step + 1; |
472 pcs.cur_step = pcs.cur_step + 1; |
489 if pcs.begin_step <= pcs.cur_step |
473 if pcs.begin_step <= pcs.cur_step |
490 pcs.step_name = 'PLOT RESULTS & SHOW WARNINGS'; begin_processing_step; |
474 pcs.step_name = 'PLOT RESULTS & SHOW WARNINGS'; begin_processing_step; |
527 |
511 |
528 streamer([p.name,' Figure 11']); |
512 streamer([p.name,' Figure 11']); |
529 pause(0.01) |
513 pause(0.01) |
530 |
514 |
531 end_processing_step; |
515 end_processing_step; |
532 end % OF STEP 17: PLOT RESULTS & SHOW WARNINGS |
516 end % OF STEP 16: PLOT RESULTS & SHOW WARNINGS |
533 |
517 |
534 %---------------------------------------------------------------------- |
518 %---------------------------------------------------------------------- |
535 % STEP 18: SAVE OUTPUT |
519 % STEP 17: SAVE OUTPUT |
536 %---------------------------------------------------------------------- |
520 %---------------------------------------------------------------------- |
537 |
521 |
538 pcs.cur_step = pcs.cur_step + 1; |
522 pcs.cur_step = pcs.cur_step + 1; |
539 if pcs.begin_step <= pcs.cur_step |
523 if pcs.begin_step <= pcs.cur_step |
540 pcs.step_name = 'SAVE OUTPUT'; begin_processing_step; |
524 pcs.step_name = 'SAVE OUTPUT'; begin_processing_step; |