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

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

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

Revision 1.2 - (hide annotations) (download)
Wed Jul 27 01:02:54 2005 UTC (18 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57s_post, checkpoint58b_post, checkpoint57y_post, checkpoint62c, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint62a, checkpoint62d, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint57t_post, checkpoint57v_post, checkpoint60, checkpoint61, checkpoint62, checkpoint58w_post, checkpoint57y_pre, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, checkpoint58r_post, checkpoint58n_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, checkpoint58k_post, checkpoint62b, checkpoint58v_post, checkpoint58l_post, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint58h_post, checkpoint58j_post, checkpoint57o_post, checkpoint61q, checkpoint61z, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61x, checkpoint61y
Changes since 1.1: +12 -6 lines 
updated (to read new grid files).
1 adcroft 1.1 for tl=1:6,
2 jmc 1.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     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 adcroft 1.1 end
9     X=cos(lat).*sin(lon);Y=-cos(lat).*cos(lon);Z=sin(lat);clear lat lon
10    
11 jmc 1.2 %lon0=-90 *pi/180;
12     lon0=+90 *pi/180;
13 adcroft 1.1 lat0= 0 *pi/180;
14     xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);
15     ro=0.3;
16     R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 );
17 jmc 1.2 t=1+(1+cos( pi*min(R/ro,1+0*R) ))/2;
18     wrda('T.init',t,1,'real*8','b');
19 adcroft 1.1
20 jmc 1.2 %lon0= 0 *pi/180;
21     lon0=180 *pi/180;
22 adcroft 1.1 lat0=35 *pi/180;
23     xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);
24     ro=0.3;
25     R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 );
26     s=1+(1+cos( pi*min(R/ro,1+0*R) ))/2;
27     wrda('S.init',s,1,'real*8','b');
28    
29     lon0=0 *pi/180;
30     lat0=0 *pi/180;
31     xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);
32     ro=0.3;
33     R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 );
34     h=1e4*(1+cos( pi*min(R/ro,1+0*R) ))/2; %cosine bell
35     %h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone
36     %wrda('mountain_eq.init',h,1,'real*8','b');
37    
38     lon0=0 *pi/180;
39     lat0=30 *pi/180;
40     xo=cos(lat0).*sin(lon0);yo=-cos(lat0).*cos(lon0);zo=sin(lat0);
41     ro=0.3;
42     R=sqrt( (X-xo).^2 + (Y-yo).^2 + (Z-zo).^2 );
43     h=1e4*(1+cos( pi*min(R/ro,1+0*R) ))/2; %cosine bell
44     %h=1e4*(1-( min(R/ro,1+0*R) ))/2; %cone
45     %wrda('mountain_30.init',h,1,'real*8','b');
  ViewVC Help
Powered by ViewVC 1.1.22