2D_example/input/genibcao.m

% This is a matlab script that generates the input data
% variable x resolution

% Initialize...
clear
if 0
load IBCAO_2km.mat
LON=LON-90;
aaa=find(LON>180);
LON(aaa)=LON(aaa)-360;
XC=[2:1451];YC=[1451:2900];
LON=LON(XC,YC);
LAT=LAT(XC,YC);
IBCAO=IBCAO(XC,YC);
end

prec='real*8';
ieee='b';

% Dimensions of grid
nx=200;
ny=1;
nz=100;
% Nominal depth of model (meters)
H=500.0;

% Horizontal resolution (m)
% dx
iswitch=80;
dx=zeros(nx,1);
for i=1:nx
  dx(i) = 35 + 16*tanh((iswitch-i)/50);
  xG(i)=sum(dx(1:i));
end
Lx=sum(dx(:))
sum(dx(51:end))
plot(dx)
%dy = Lx/nx
writebin('dx.bin',dx,1,prec);


break
% dlon
re = 6378137.0d0;
slat=70.4;
dlon2dx=re*2*pi*cosd(slat)/360;
dx2dlon=1./dlon2dx;
dlon=zeros(nx,1);
lonG=zeros(nx+1,1);
lonC=zeros(nx,1);
lonG(nx+1)=-50.6;
for i=nx:-1:1
  dlon(i)=dx(i).*dx2dlon;
  lonG(i)=lonG(i+1)-dlon(i);
  lonC(i)=0.5*lonG(i+1)+0.5*lonG(i); 
end
latC(1)=70.42;
dlat=dlon./8;
for i=1:nx-1
latC(i+1)=latC(i)-dlat(i);
end
latC=latC';
Llon=sum(dlon(:));
 writebin('dlon.bin',dlon,1,prec);
break
% Interpolate bathymetry to my gridC
%[tri,w]=interp_TS(lonC,latC,LON,LAT);
% z = IBCAO(:).'; 
% ibcao_store_2d = sum(z(tri) .* w,2);
% ibcao_store_2d(nx)=0.0;
% writebin('bathy_store_2D.ibcao',ibcao_jks_2d,1,prec);
 bathy=-500*ones(200,1);
 bathy(nx)=0;
 writebin('bathy_store_2D.500',bathy,1,prec);

% ICEFRONT depth and circum files
frontdepth=zeros(nx,ny);
frontdepth(nx-1)=bathy(nx-1);
frontlength=zeros(nx,ny);
frontlength(nx-1)=1./dx(nx-1);
writebin('frontdepth.bin',frontdepth,1,prec);
writebin('frontlength.bin',frontlength,1,prec);

1	% This is a matlab script that generates the input data
2	% variable x resolution
3
4	% Initialize...
5	clear
6	if 0
7	load IBCAO_2km.mat
8	LON=LON-90;
9	aaa=find(LON>180);
10	LON(aaa)=LON(aaa)-360;
11	XC=[2:1451];YC=[1451:2900];
12	LON=LON(XC,YC);
13	LAT=LAT(XC,YC);
14	IBCAO=IBCAO(XC,YC);
15	end
16
17	prec='real*8';
18	ieee='b';
19
20	% Dimensions of grid
21	nx=200;
22	ny=1;
23	nz=100;
24	% Nominal depth of model (meters)
25	H=500.0;
26
27	% Horizontal resolution (m)
28	% dx
29	iswitch=80;
30	dx=zeros(nx,1);
31	for i=1:nx
32	dx(i) = 35 + 16*tanh((iswitch-i)/50);
33	xG(i)=sum(dx(1:i));
34	end
35	Lx=sum(dx(:))
36	sum(dx(51:end))
37	plot(dx)
38	%dy = Lx/nx
39	writebin('dx.bin',dx,1,prec);
40
41
42	break
43	% dlon
44	re = 6378137.0d0;
45	slat=70.4;
46	dlon2dx=re2pi*cosd(slat)/360;
47	dx2dlon=1./dlon2dx;
48	dlon=zeros(nx,1);
49	lonG=zeros(nx+1,1);
50	lonC=zeros(nx,1);
51	lonG(nx+1)=-50.6;
52	for i=nx:-1:1
53	dlon(i)=dx(i).*dx2dlon;
54	lonG(i)=lonG(i+1)-dlon(i);
55	lonC(i)=0.5lonG(i+1)+0.5lonG(i);
56	end
57	latC(1)=70.42;
58	dlat=dlon./8;
59	for i=1:nx-1
60	latC(i+1)=latC(i)-dlat(i);
61	end
62	latC=latC';
63	Llon=sum(dlon(:));
64	writebin('dlon.bin',dlon,1,prec);
65	break
66	% Interpolate bathymetry to my gridC
67	%[tri,w]=interp_TS(lonC,latC,LON,LAT);
68	% z = IBCAO(:).';
69	% ibcao_store_2d = sum(z(tri) .* w,2);
70	% ibcao_store_2d(nx)=0.0;
71	% writebin('bathy_store_2D.ibcao',ibcao_jks_2d,1,prec);
72	bathy=-500*ones(200,1);
73	bathy(nx)=0;
74	writebin('bathy_store_2D.500',bathy,1,prec);
75
76	% ICEFRONT depth and circum files
77	frontdepth=zeros(nx,ny);
78	frontdepth(nx-1)=bathy(nx-1);
79	frontlength=zeros(nx,ny);
80	frontlength(nx-1)=1./dx(nx-1);
81	writebin('frontdepth.bin',frontdepth,1,prec);
82	writebin('frontlength.bin',frontlength,1,prec);
83