[MITgcm-support] [Fwd: Auto-discard notification]

[Ed Hill]

Hi Jane,

I subscribed you to the MITgcm-support list and am forwarding your
message.

Ed

-----Forwarded Message-----
> From: jane at ocean.fsu.edu
> To: MITgcm-support at mitgcm.org
> Cc: dewar at ocean.fsu.edu
> Subject: periodic boundary conditions in X, walls in Y
> Date: Thu, 23 Sep 2004 17:29:39 -0400
> 
> i'm trying to set up a channel where there are perodic boundary conditions
> in X and walls in Y.  i doesn't seem to work.  any suggestions would be
> helpful.
> ---------------------------------------------------------------------------
> data.pkg
> ---------------------------------------------------------------------------
> # Packages
>  &PACKAGES
>  useOBCS=.FALSE.,
>  &
> ---------------------------------------------------------------------------
> data
> ---------------------------------------------------------------------------
> # ====================
> # | Model parameters |
> # ====================
> #
> # Continuous equation parameters
>  &PARM01
>  tRef= 50*0,
>  sRef= 50*35.,
>  viscAz=1.E-4,
>  viscAh=2.E1,
>  no_slip_sides=.FALSE.,
>  no_slip_bottom=.FALSE.,
> 
>  viscA4=0.E12,
>  diffKhT=2.E1,
>  diffKzT=1.E-5,
>  diffKhS=2.E1,
>  diffKzS=1.E-5,
>  f0=1.e-4,
>  beta=0.E-11,
>  tAlpha=2.E-4,
>  sBeta =0.E-4,
>  gravity=9.81,
>  gBaro=9.81,
>  rigidLid=.FALSE.,
>  implicitFreeSurface=.TRUE.,
>  eosType='LINEAR',
>  hFacMin=0.2,
>  nonHydrostatic=.FALSE.,
>  readBinaryPrec=64,
>  globalFiles=.TRUE.,
>  &
> 
> # Elliptic solver parameters
>  &PARM02
>  cg2dMaxIters=1000,
>  cg2dTargetResidual=1.E-13,
>  cg3dMaxIters=400,
>  cg3dTargetResidual=1.E-13,
>  &
> 
> # Time stepping parameters
>  &PARM03
>  niter0=0,
>  nTimeSteps=6048,
>  deltaT=200.0,
>  abEps=0.1,
>  pChkptFreq=604800.0,
>  chkptFreq=0.0,
>  dumpFreq=86400.,
>  monitorFreq=86400.,
>  &
> 
> # Gridding parameters
>  &PARM04
>  usingCartesianGrid=.TRUE.,
>  usingSphericalPolarGrid=.FALSE.,
>  delX=200*5.e3,
>  delY=50*5.e3,
>  delZ=50*30,
>  &
> 
> # Input datasets
>  &PARM05
>  bathyFile='topog.channel',
>  hydrogThetaFile='theta.bin',
>  &
> 
> ---------------------------------------------------------------------------
> gendata.m
> ---------------------------------------------------------------------------
> 
> % This is a matlab script that generates the input data
> 
> % Dimensions of grid
> nx=200;
> ny=50;
> nz=50;
> % Nominal depth of model (meters)
> H=1500;
> % Scale of bump (m)
> L=25e3;
> % Height of bump (m)
> %dh=0.90*H;
> % Horizontal resolution (m)
> dx=5e3;
> % Rotation
> f=1e-4;
> % Stratification
> N=1.5 * f*L/H;
> 
> % Gravity
> g=9.81;
> % E.O.S.
> alpha=2.e-4;
> 
> z=[15:30:1500];
> dz=30;
> % dz=H/nz;
> sprintf('delZ = %d * %7.6g',nz,dz)
> 
> x=(1:nx)*dx;x=x-mean(x);x=x+nx/3*dx;
> y=(1:ny)*dx;y=y-mean(y);
> [Y,X]=meshgrid(y,x);
> % z=-dz/2:-dz:-H;
> %z=-(cumsum(dz)-dz/2);
> 
> % Temperature profile
> %Tz=N^2/(g*alpha)
> %phi=sprintf(' %8.6g',Tz*z+sum(Tz*dz)/2);
> %Tref=str2num(phi);
> %[sprintf('Tref =') sprintf(' %8.6g,',Tz*z+sum(Tz*dz)/2)]
> 
> ieee='b';
> accuracy='real*8';
> 
> % Create three dimensional temperature profile
> t3d=zeros(nx,ny,nz);
> r0=20e3;
> r2=X.^2 + Y.^2;
> z0=300;
> for k=1:nz
>  t3d(:,:,k)=10-10*z(k)/H;
>  for i=1:nx
>  for j=1:ny
>   if Y(i,j) >= 0
>    if z(k) >=240 
>      zp=z(k)-240;
>      fz=(abs(zp)*exp(-abs(zp)/z0))/z0;
>    else
>      fz=0;
>    end
>   else
>    if z(k) < 1260
>      zp=z(k)-1260;
>      fz=(abs(zp)*exp(-abs(zp)/z0))/z0;
>    else
>      fz=0;
>    end
>   end
> %  Tp=fz.*(r0.*Y(i,j)./r0^2).*exp(-r2(i,j)./r0^2).*4.*2.7*2.7;
>   zp=z(k)-750;
>   Tp=3.*exp(-zp*zp/300^2);
>   if r2(i,j) > 2*r0^2
>     rp=r2(i,j)-2*r0^2;
>     Tp=Tp*exp(-rp^2/r0^2);
>   end
>   t3d(i,j,k)=t3d(i,j,k)+Tp;
> % t3d(i,j,k)=Tp;
>  end
>  end
> end
> fid=fopen('theta.bin','w',ieee); fwrite(fid,t3d,accuracy); fclose(fid);
> 
> % Simple channel
> h=zeros(nx,ny);
> %h(:,1)=0;
> h(:,2:ny-1)=-H;
> %h(:,ny)=0;
> fid=fopen('topog.channel','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
> 
> 
> ----------------------------------------------------------------
> This message was sent using IMP, the Internet Messaging Program.
-- 
Edward H. Hill III, PhD
office:  MIT Dept. of EAPS;  Rm 54-1424;  77 Massachusetts Ave.
             Cambridge, MA 02139-4307
emails:  eh3 at mit.edu                ed at eh3.com
URLs:    http://web.mit.edu/eh3/    http://eh3.com/
phone:   617-253-0098
fax:     617-253-4464