1 |
for tl=1:6, |
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 |
end |
5 |
X=cos(lat).*sin(lon);Y=-cos(lat).*cos(lon);Z=sin(lat);clear lat lon |
6 |
|
7 |
lon0=-90 *pi/180; |
8 |
lat0= 0 *pi/180; |
9 |
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
10 |
ro=0.3; |
11 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
12 |
s=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
13 |
wrda('T.init',s,1,'real*8','b'); |
14 |
|
15 |
lon0= 0 *pi/180; |
16 |
lat0=35 *pi/180; |
17 |
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
18 |
ro=0.3; |
19 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
20 |
s=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
21 |
wrda('S.init',s,1,'real*8','b'); |
22 |
|
23 |
lon0=0 *pi/180; |
24 |
lat0=0 *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 |
h=1e4*(1+cos( pi*min(R/ro,1+0*R) ))/2; %cosine bell |
29 |
%h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone |
30 |
%wrda('mountain_eq.init',h,1,'real*8','b'); |
31 |
|
32 |
lon0=0 *pi/180; |
33 |
lat0=30 *pi/180; |
34 |
xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0); |
35 |
ro=0.3; |
36 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
37 |
h=1e4*(1+cos( pi*min(R/ro,1+0*R) ))/2; %cosine bell |
38 |
%h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone |
39 |
%wrda('mountain_30.init',h,1,'real*8','b'); |