[G1,wm,allsize]=qwt2(wm_row,ict)

function [G1,wm,allsize]=qwt2(wm_row,ict)
%  QWT2.M: summary of water type definitions used by OMP2.M. 
%  Calling qwt2.m from the main OMP analysis program omp2.m will 
%  produce the source water matrix G1.
%
%  You will most likely have to edit this file for your own
%  application. We recommend that you save your edited file under
%  a different file name. See the web manual for details.
%
%
% This program is part of the OMP package from:
% Institut fuer Meereskunde
% J. Karstensen 
% Troplowitzstr. 7
% 22529 Hamburg
% Germany
%
% Amendments by Matthias Tomczak March 1999
%
% BUGS: [e-mail]
%    or [e-mail]

if nargin<1
 disp(' ')
 disp(['     Please give the rwo index of SWT definitions you want to use:'])
 disp([' e.g. qwt_step([1 2 3 6]) activates row 1, 2, 3, and 6 '])
 disp('  '),return
end

% The following lines identify the water masses defined through each water type.
% There is one water mass name for each water type definition row.
% Each name has to consist of exactly 5 letters including blanks; DO NOT CHANGE
% THE LENGTH unless you want to edit the main program omp2.m, too!
%
% NOTES:  1. You can store more than one water type definition for a water mass
%            and select which one you want to use at run time. This example file
%            demonstrates this by offering two definitions for ICW. Calling
%            qwt2([1 2 3 4]) activates rows 1, 2, 3, and 4 (AAMW and the first
%            ICW set); calling qwt2([1 2 5 6]) activates row 1, 2, 5, and 6
%            (AAMW and the second ICW set).
%         2. You can use this function to build up your own inventory of water
%            type definitions by including definitions for the Atlantic (NACW,
%            SACW, MedW etc.) and other oceans (AAIW, AABW etc.). This will
%            allow you to run OMP analysis for any ocean region by selecting the
%            appropriate rows in the function call.
%         3. Water types with identical names define water masses by property-
%            property relationships (e.g. in the thermocline); their contributions
%            will be added to all contributions of water masses with the same name.
%            For example, the call qwt2([1 2 3 4 7]) will produce three results,
%            AAMW (added from 1 and 2), ICW (added from 3 and 4) and AAIW.
%            For this to work all water types with identical names for which you
%            want contributions added in the result have to be called in an
%            uninterrupted sequence. (For example, calling qwt2([1 2 7]) will add
%            the AAMW contributions, calling qwt2([1 7 2]) will not.)
%

  wm(1:5)   = ' AAMW'; %                      first row
  wm(6:10)  = ' AAMW'; %                      second row
  wm(11:15) = '  ICW'; %upper, first set      third row
  wm(16:20) = '  ICW'; %lower, first set      forth row
  wm(21:25) = '  ICW'; %upper, second set     fifth row
  wm(26:30) = '  ICW'; %lower, second set     sixth row
  wm(31:35) = ' AAIW'; %                      seventh row
  wm(36:40) = '  IEW'; %                      eigth row

all=[...

%  The following lines define the water types. The order of parameters is
%  ptemp    sal      oxy    PO4     NO3    Si    mass   pvort
%  Note: potential vorticity is multiplied by 10*8.

  10       34.56      91    2.1    30     40     1.0    0.03;...  %1 lower AAMW
  16.4     34.55     100    1.4    19     25     1.0    1.12;...  %2 upper AAMW  
   9       34.65     260    1.1    15      5     1.0    0.03;...  %3 lower ICW, first set
  18       35.8      230    0       0      0.5   1.0    0.05;...  %4 upper ICW, first set
   9       34.72     209    1.47   20      5     1.0    0.03;...   %5 lower ICW, second set
  14.35    35.4      224    0.6     6.5    0.5   1.0    0.05;...  %6 upper ICW, second set
   4.5     34.35     210    2.2    32     35     1.0    0.30;...  %7 AAIW
   8.5     35         60    2.5    35     60     1.0    0.04;...  %8 IEW

];

G1=all(wm_row,:)';
allsize = size(all);