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

Contents of /MITgcm/verification/exp5/input/gendata.m

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


Revision 1.3 - (show annotations) (download)
Thu Apr 24 19:28:08 2008 UTC (13 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: HEAD
Changes since 1.2: +0 -0 lines
FILE REMOVED
remove this experiment (replaced by tutorial_deep_convection)

1 % This is a matlab script that generates the input data
2
3 % Dimensions of grid
4 nx=64;
5 ny=64;
6 nz=20;
7 % Nominal depth of model (meters)
8 H=2000;
9 % Size of domain
10 Lx=32e3;
11 % Scale of disk (m)
12 L=8e3;
13 % Horizontal resolution (m)
14 dx=Lx/nx;
15 % Rotation
16 f=1e-4;
17 % Stratification
18 N=0. * f*L/H;
19 % Flux
20 Qo=1;
21
22 % Gravity
23 g=9.81;
24 % E.O.S.
25 alpha=2.e-4;
26
27 Tz=N^2/(g*alpha)
28
29 dz=H/nz;
30 sprintf('delZ = %d * %7.6g,',nz,dz)
31
32 x=(1:nx)*dx;x=x-mean(x);
33 y=(1:ny)*dx;y=y-mean(y);
34 z=-dz/2:-dz:-H;
35
36 [Y,X]=meshgrid(y,x);
37
38 % Temperature profile
39 [sprintf('Tref =') sprintf(' %8.6g,',Tz*z-mean(Tz*z))]
40
41 % Gaussian bump
42 Q=Qo*(1+0.01*rand([nx,ny]));
43 r=sqrt(X.^2+Y.^2);
44 %% Q( find(r>L) )=0;
45 fid=fopen('Qnet.circle','w','b'); fwrite(fid,Q,'real*8'); fclose(fid);

  ViewVC Help
Powered by ViewVC 1.1.22