/[MITgcm]/MITgcm_contrib/mlosch/interp_llc/check_HxHy.m
ViewVC logotype

Contents of /MITgcm_contrib/mlosch/interp_llc/check_HxHy.m

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


Revision 1.1 - (show annotations) (download)
Thu May 3 21:07:20 2007 UTC (18 years, 2 months ago) by mlosch
Branch: MAIN
CVS Tags: HEAD
initial checkin of topography and hydrography interpolation scripts for
the llc-grid, based on old matlab scripts by Alistair Adcroft
Let's hope, they are useful.

1 function []=check_HxHy(H2,Hx2,Hy2)
2 % Check that Hx2 and Hy2 don't violate basic rules
3 % This should *not* ever happen
4 %return; % spead up process
5 [Hxcc,Hycc]=gen_hxhy(H2);
6 ic=find(Hx2<Hxcc-0.001);
7 if prod(size(ic)) ~= 0
8 %[Hx2(ic)/1e3 Hxcc(ic)/1e3 Hx2(ic)-Hxcc(ic)]
9 error('The X cross-sectional area turned out inconsistently')
10 end
11 jc=find(Hy2<Hycc-0.001);
12 if prod(size(jc)) ~= 0
13 %[Hy2(jc)/1e3 Hycc(jc)/1e3 Hy2(jc)-Hycc(jc)]
14 error('The Y cross-sectional area turned out inconsistently')
15 end

  ViewVC Help
Powered by ViewVC 1.1.22