1 |
adcroft |
1.1 |
for tl=1:6, |
2 |
jmc |
1.2 |
%lon(:,tl,:)=rdda(sprintf('LONC.%3.3i.bin',tl),[32 32],1,'real*8','b')*pi/180; |
3 |
|
|
%lat(:,tl,:)=rdda(sprintf('LATC.%3.3i.bin',tl),[32 32],1,'real*8','b')*pi/180; |
4 |
|
|
xx=rdda(sprintf('grid_cs32.face%3.3i.bin',tl),[33 33],1,'real*8','b'); |
5 |
|
|
lon(:,tl,:)=xx(1:32,1:32)*pi/180; |
6 |
|
|
xx=rdda(sprintf('grid_cs32.face%3.3i.bin',tl),[33 33],2,'real*8','b'); |
7 |
|
|
lat(:,tl,:)=xx(1:32,1:32)*pi/180; |
8 |
adcroft |
1.1 |
end |
9 |
|
|
X=cos(lat).*sin(lon);Y=-cos(lat).*cos(lon);Z=sin(lat);clear lat lon |
10 |
|
|
|
11 |
jmc |
1.2 |
%lon0=-90 *pi/180; |
12 |
|
|
lon0=+90 *pi/180; |
13 |
adcroft |
1.1 |
lat0= 0 *pi/180; |
14 |
|
|
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
15 |
|
|
ro=0.3; |
16 |
|
|
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
17 |
jmc |
1.2 |
t=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
18 |
dfer |
1.3 |
fid=fopen('T.init','w','b'); |
19 |
|
|
fwrite(fid,t(:),'real*8'); |
20 |
|
|
fclose(fid) |
21 |
adcroft |
1.1 |
|
22 |
jmc |
1.2 |
%lon0= 0 *pi/180; |
23 |
|
|
lon0=180 *pi/180; |
24 |
adcroft |
1.1 |
lat0=35 *pi/180; |
25 |
|
|
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
26 |
|
|
ro=0.3; |
27 |
|
|
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
28 |
|
|
s=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
29 |
dfer |
1.3 |
fid=fopen('S.init','w','b'); |
30 |
|
|
fwrite(fid,s(:),'real*8'); |
31 |
|
|
fclose(fid) |
32 |
adcroft |
1.1 |
|
33 |
|
|
lon0=0 *pi/180; |
34 |
|
|
lat0=0 *pi/180; |
35 |
|
|
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
36 |
|
|
ro=0.3; |
37 |
|
|
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
38 |
|
|
h=1e4*(1+cos( pi*min(R/ro,1+0*R) ))/2; %cosine bell |
39 |
|
|
%h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone |
40 |
dfer |
1.3 |
%fid=fopen('mountain_eq.init','w','b'); |
41 |
|
|
%fwrite(fid,h(:),'real*8'); |
42 |
|
|
%fclose(fid) |
43 |
adcroft |
1.1 |
|
44 |
|
|
lon0=0 *pi/180; |
45 |
|
|
lat0=30 *pi/180; |
46 |
|
|
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
47 |
|
|
ro=0.3; |
48 |
|
|
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
49 |
|
|
h=1e4*(1+cos( pi*min(R/ro,1+0*R) ))/2; %cosine bell |
50 |
|
|
%h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone |
51 |
dfer |
1.3 |
%fid=fopen('mountain_30.init','w','b'); |
52 |
|
|
%fwrite(fid,h(:),'real*8'); |
53 |
|
|
%fclose(fid) |