/[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.1 by adcroft, Fri Sep 28 02:30:10 2001 UTC revision 1.3 by dfer, Tue Mar 30 13:59:46 2010 UTC
# Line 1  Line 1 
1  for tl=1:6,  for tl=1:6,
2  lon(:,tl,:)=rdda(sprintf('LONC.%3.3i.bin',tl),[32 32],1,'real*8','b')*pi/180;  %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;  %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  end  end
9  X=cos(lat).*sin(lon);Y=-cos(lat).*cos(lon);Z=sin(lat);clear lat lon  X=cos(lat).*sin(lon);Y=-cos(lat).*cos(lon);Z=sin(lat);clear lat lon
10    
11  lon0=-90  *pi/180;  %lon0=-90  *pi/180;
12    lon0=+90  *pi/180;
13  lat0=  0  *pi/180;  lat0=  0  *pi/180;
14  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);
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  s=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',s,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;
24  lat0=35   *pi/180;  lat0=35   *pi/180;
25  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);  xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);
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 27  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 36  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.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22