/[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.2 by jmc, Wed Jul 27 01:02:54 2005 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');  wrda('T.init',t,1,'real*8','b');
19    
20  lon0=  0  *pi/180;  %lon0=  0  *pi/180;
21    lon0=180  *pi/180;
22  lat0=35   *pi/180;  lat0=35   *pi/180;
23  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);
24  ro=0.3;  ro=0.3;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22