/[MITgcm]/MITgcm_contrib/high_res_cube/eddy_flux/calcUVBave.m
ViewVC logotype

Diff of /MITgcm_contrib/high_res_cube/eddy_flux/calcUVBave.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by edhill, Tue Aug 3 15:54:26 2004 UTC revision 1.3 by edhill, Thu Aug 12 15:50:38 2004 UTC
# Line 109  for il = 1:nlay Line 109  for il = 1:nlay
109        bcube(:,:,icf) = eh3densjmd95( squeeze(scube(:,:,icf)), ...        bcube(:,:,icf) = eh3densjmd95( squeeze(scube(:,:,icf)), ...
110                                       squeeze(tcube(:,:,icf)), p );                                       squeeze(tcube(:,:,icf)), p );
111        if il > 1        if il > 1
112          p = -9.998*10^2 * 9.81 * R(il-1) * 0.0001;          p = 9.998*10^2 * 9.81 * R(il-1) * 0.0001;
113          bcubem1(:,:,icf) = eh3densjmd95( scube(:,:,icf), ...          bcubem1(:,:,icf) = eh3densjmd95( scube(:,:,icf), ...
114                                           tcube(:,:,icf), p );                                           tcube(:,:,icf), p );
115        end        end
116      end      end
117            tmask   = double(tcube ~= 0.0);
118      %  Fill the (ne+1)*(ne+1) grid with T values      bcube   = bcube .* tmask;
119        bcubem1 = bcubem1 .* tmask;
120    
121        %  Fill the (ne+1)*(ne+1) grid with B,T values
122      ni = ne;    nip1 = ni + 1;      ni = ne;    nip1 = ni + 1;
123      nj = ne;    njp1 = nj + 1;      nj = ne;    njp1 = nj + 1;
124      bcubep1 = zeros( [ nip1 njp1 6 ] );      bcubep1 = zeros( [ nip1 njp1 6 ] );
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22