/[MITgcm]/MITgcm/verification/advect_cs/input/gendata.m
ViewVC logotype

Diff of /MITgcm/verification/advect_cs/input/gendata.m

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

revision 1.2 by jmc, Wed Jul 27 01:02:54 2005 UTC revision 1.3 by dfer, Tue Mar 30 13:59:46 2010 UTC
# Line 15  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*c Line 15  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*c
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;
# Line 24  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*c Line 26  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*c
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;
# Line 33  ro=0.3; Line 37  ro=0.3;
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;
# Line 42  ro=0.3; Line 48  ro=0.3;
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)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22