%
%  Ed Hill
%  Mon Aug  2 23:43:41 EDT 2004
%

matlab -nojvm -nodisplay

  clear all ; close all

  %  Get the 1-yr list of days
  load ilist
  deltatT = 1200;                % model time step size (s)
  tavefreq = 259200;             % averaging period (s)
  startDate = datenum(1992,1,1); % model integration starting date
  for iy = 1992:2004
    fid = fopen(sprintf('ilist_%d.txt',iy),'w');
    for it = 1:length(ilist)
      iter = ilist(it);
      dv = datevec(startDate + (iter*deltatT-tavefreq/2)/60/60/24);
      newyear = dv(1);
      if newyear == iy
        a = sprintf('%10.10d %d',iter,dv(1));
        fprintf(fid,'%s\n',a);
      end
    end
    fclose(fid);
  end
  

  clear all ; close all
  nlay = 50;
  ne = 510;
  
  Spat = 'Stave/data/Stave.%10.10d.data';
  Tpat = 'Ttave/data/Ttave.%10.10d.data';
  Upat = '_Veltave/data/uVeltave.%10.10d.data';
  Vpat = '_Veltave/data/vVeltave.%10.10d.data';
  
  load ilist

  calcUVBave(Spat, Tpat, Upat, Vpat, ilist, nlay, ne)

  
  
#=======  split into 4 nodes (3 each)  =======
cat ilist_199[2-4].txt | awk '{print $1}' > ilist_92--94
cat ilist_199[5-7].txt | awk '{print $1}' > ilist_95--97
cat ilist_199[89].txt ilist_2000.txt | awk '{print $1}' > ilist_98--00
cat ilist_200[1-3].txt | awk '{print $1}' > ilist_01--03
wc -l ilist_??--??
scp ilist_??--?? ta54:

#=======  split into 6 nodes (2 each)  =======
cat ilist_1992.txt ilist_1993.txt  | awk '{print $1}' > ilist_40
cat ilist_1994.txt ilist_1995.txt  | awk '{print $1}' > ilist_41
cat ilist_1996.txt ilist_1997.txt  | awk '{print $1}' > ilist_42
cat ilist_1998.txt ilist_1999.txt  | awk '{print $1}' > ilist_43
cat ilist_2000.txt ilist_2001.txt  | awk '{print $1}' > ilist_44
cat ilist_2002.txt ilist_2003.txt  | awk '{print $1}' > ilist_45
wc -l ilist_4[0-5]
scp ilist_4[0-5] ta54:


/scratch/edhill/matlab_r13/bin/matlab -nojvm -nodisplay

clear all ; close all
  nlay = 50;
  ne = 510;
  
  Spat = 'c5l/Stave/data/Stave.%10.10d.data';
  Tpat = 'c5l/Ttave/data/Ttave.%10.10d.data';
  Upat = 'c5l/_Veltave/data/uVeltave.%10.10d.data';
  Vpat = 'c5l/_Veltave/data/vVeltave.%10.10d.data';
  load ilist

  !/sbin/ifconfig eth0 | grep 'inet addr:' | sed -e 's|\.| |g' | awk '{print $5}' > myip
  load myip
  ips = [ myip  201 200 200 200 ]

  calcUVBave_ttcp(ips, Spat, Tpat, Upat, Vpat, ilist, nlay, ne)