[MITgcm-support] value of OBCS sponge U/Vrelaxobcsinner U/Vrelaxobcsbound

Tue Feb 14 06:36:46 EST 2012

Hell all,
I run a model with tide velocity(M2, S2, K1, O1 constituents) at the boundary, and I changing the U/Vrelaxobcsinner value will influence the Eta inside the domain a lot. I don't know the rule to give the most proper value of them. I noticed the following suggestion was given earlier but it does not work for my case. Maybe the tidal wave does not belong to mode-1 internal waves?
 Is there any suggestions or the reference about the sponge theory?
 Any help is much appreciate!
>>> Say you have a forcing frequency w.
>>>  - number of points: as long as the horizontal wavelength of mode-1 internal waves of frequency w,
>>>  - [UV]relaxobcsinner = (2*pi/w)/1000,
>>>  - [UV]relaxobcsbound = 2*pi/w.

Chunyan Zhou
Division of Civil Engineering
School of Engineering, Physics and Mathematics
College of Art, Science and Engineering
Fulton Building,G19
University of Dundee
Dundee, UK DD1 4HN
Tel. 01382 385431

________________________________________
From: mitgcm-support-request at mitgcm.org [mitgcm-support-request at mitgcm.org]
Sent: 13 February 2012 17:00
To: mitgcm-support at mitgcm.org
Subject: MITgcm-support Digest, Vol 104, Issue 12

Send MITgcm-support mailing list submissions to
        mitgcm-support at mitgcm.org

To subscribe or unsubscribe via the World Wide Web, visit
        http://mitgcm.org/mailman/listinfo/mitgcm-support
or, via email, send a message with subject or body 'help' to
        mitgcm-support-request at mitgcm.org

You can reach the person managing the list at
        mitgcm-support-owner at mitgcm.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of MITgcm-support digest..."

Today's Topics:

   1. CS Grid Generator and Convertor Matlab Scripts
      (Christopher Watkins)

----------------------------------------------------------------------

Message: 1
Date: Mon, 13 Feb 2012 16:29:17 +0000 (GMT)
From: Christopher Watkins <c.watkins at qmul.ac.uk>
To: "." <mitgcm-support at mitgcm.org>
Subject: [MITgcm-support] CS Grid Generator and Convertor Matlab
        Scripts
Message-ID:
        <1846019721.94520.1329150557623.JavaMail.root at planck.maths.qmul.ac.uk>
Content-Type: text/plain; charset=utf-8

Hi,

I'm trying to generate cubed sphere grids of different sizes for different planetary radii.

I have downloaded the generator and convertor scripts from http://mitgcm.org/viewvc/MITgcm/MITgcm_contrib/high_res_cube/

I use gengrids.m changing line 9 to for nx=1+[16 32 64] to generate CS grids with faces of 16x16, 32x32 and 64x64 squares.

I edit line 167 in gengrid_fn.m so that Rsphere is equal to the radius of the planet.

This produces three directories: cs-conf-c-16-128, cs-conf-c-32-256 and cs-conf-c-64-512.

The directories consist of the following files:

CUBE_3DGRID.mat
CUBE_HGRID.mat
DXC.bin
DXF.bin
DXG.bin
DXV.bin
DYC.bin
DYF.bin
DYG.bin
DYU.bin
grid.info
LATC.001.bin
LATC.002.bin
LATC.003.bin
LATC.004.bin
LATC.005.bin
LATC.006.bin
LATG.001.bin
LATG.002.bin
LATG.003.bin
LATG.004.bin
LATG.005.bin
LATG.006.bin
LONC.001.bin
LONC.002.bin
LONC.003.bin
LONC.004.bin
LONC.005.bin
LONC.006.bin
LONG.001.bin
LONG.002.bin
LONG.003.bin
LONG.004.bin
LONG.005.bin
LONG.006.bin
RA.bin
RAS.bin
RAW.bin
RAZ.bin
TUV.mat

I then use convertMITgrid.m in each directory to create the .bin files required to run using a GS grid. I edit line 1 to be 17, 33 or 65 (ie the number used in gengrids.m line 9 plus 1)

However, in this last step I get errors listed shown below. I've hunted through the support forum and Google but can't find what I'm doing wrong. Can anyone advise me?

Many thanks,

Chris

ans =

Error while opening file:
No such file or directory

ans =

Error while opening file:
No such file or directory

ans =