15 |
ro=0.3; |
ro=0.3; |
16 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
17 |
t=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
t=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
18 |
wrda('T.init',t,1,'real*8','b'); |
fid=fopen('T.init','w','b'); |
19 |
|
fwrite(fid,t(:),'real*8'); |
20 |
|
fclose(fid) |
21 |
|
|
22 |
%lon0= 0 *pi/180; |
%lon0= 0 *pi/180; |
23 |
lon0=180 *pi/180; |
lon0=180 *pi/180; |
26 |
ro=0.3; |
ro=0.3; |
27 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
28 |
s=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
s=1+(1+cos( pi*min(R/ro,1+0*R) ))/2; |
29 |
wrda('S.init',s,1,'real*8','b'); |
fid=fopen('S.init','w','b'); |
30 |
|
fwrite(fid,s(:),'real*8'); |
31 |
|
fclose(fid) |
32 |
|
|
33 |
lon0=0 *pi/180; |
lon0=0 *pi/180; |
34 |
lat0=0 *pi/180; |
lat0=0 *pi/180; |
37 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
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 |
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 |
%h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone |
40 |
%wrda('mountain_eq.init',h,1,'real*8','b'); |
%fid=fopen('mountain_eq.init','w','b'); |
41 |
|
%fwrite(fid,h(:),'real*8'); |
42 |
|
%fclose(fid) |
43 |
|
|
44 |
lon0=0 *pi/180; |
lon0=0 *pi/180; |
45 |
lat0=30 *pi/180; |
lat0=30 *pi/180; |
48 |
R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 ); |
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 |
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 |
%h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone |
51 |
%wrda('mountain_30.init',h,1,'real*8','b'); |
%fid=fopen('mountain_30.init','w','b'); |
52 |
|
%fwrite(fid,h(:),'real*8'); |
53 |
|
%fclose(fid) |