[MITgcm-support] different time period OBCS

Chun-Yan Zhou c.zhou at dundee.ac.uk
Thu May 3 18:03:30 EDT 2012


Hi Michael,

 Yes, actually I provide Eta as boundary condition, not velocity. the OBEeta file time space is 300s, and the file size is OK, but the OBEt file with 300s is much much larger because the Nr=40 in my case.

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
-----------------------------------------------------------------------
if the velocity period is 300s and the file OB[WESN][uv].bin (forcing files for the velocities) is a certain
size then the temperature and salinity forcing files will be exactly the same size.





________________________________________
From: mitgcm-support-request at mitgcm.org [mitgcm-support-request at mitgcm.org]
Sent: 02 May 2012 14:05
To: mitgcm-support at mitgcm.org
Subject: MITgcm-support Digest, Vol 107, Issue 3

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. Re: Initial Wind on Cubed Sphere (Jean-Michel Campin)
   2. SMF package status (Gianmaria Sannino)
   3. Re: Definition of spatial boundary in EXF package (Martin Losch)
   4. Re: SMF package status (Chris Hill)
   5. Re: SMF package status (Gianmaria Sannino)


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

Message: 1
Date: Tue, 1 May 2012 13:33:50 -0400
From: Jean-Michel Campin <jmc at ocean.mit.edu>
To: mitgcm-support at mitgcm.org
Subject: Re: [MITgcm-support] Initial Wind on Cubed Sphere
Message-ID: <20120501173350.GB13383 at ocean.mit.edu>
Content-Type: text/plain; charset=us-ascii

Chris,

This is a good question (I did not look to this script for some time,
need to remember how things are done) and should be documented at the
top ofuvLatLon2cube.m.
"spv" is an optional argument (but it's unfortunate that it's before cosalpha,sinalpha,
which you may want to specify without providing "spv") for an exclusion-value
(or "special value") that could be found in the input uFld,vFld over dead region
(e.g., continent when processing oceanic velocity).
And when "spv" is specified, the interpolation discard these "dead regions"
and, at the end, put back "spv" over these regions.
So if you specify a value for "spv" that is not present anywhere neither in
uFld nor in vFld, then it will have no effect on the result.

Cheers,
Jean-Michel

On Tue, May 01, 2012 at 04:28:43PM +0100, Chris Watkins wrote:
> Hi Jean-Michel,
>
> Thank you for such a quick response. I was looking at
> utils/matlab/cs_grid/uvLatLon2cube.m. It asks for spv as one of its input.
> Can you tell me what spv is?
>
> Thanks again,
>
> Chris
>
> On 1 May 2012 16:11, Jean-Michel Campin <jmc at ocean.mit.edu> wrote:
>
> > Hi Chris,
> >
> > Depending on the wind field, you may try one of the 2 methods:
> > 1) general method (similar to what is done within matlab script
> >  utils/matlab/cs_grid/uvLatLon2cube.m): your can define your wind field
> >  on a lat-lon grid and then use uvLatLon2cube.m to interpolate + rotate
> >  to CS-grid direction + put the uVel,vVel on to C-grid location;
> > Or, you can skip the interpolation stage and define your wind field
> >  on the CS-grid, at the grid-cell center location and then just do
> >  the rotation + the move to C-grid location.
> > 2) specific method for non-divergent wind-field (this is a special case
> >  which is not uncommon): From the analytical expression of the wind,
> >  you can derive an analytical expression of the horizontal stream-function.
> >  Then you evaluate the stream-function on the CS-grid at the grid-cell
> > corner
> >  (position XG,YG) and from there you can compute directly the wind
> >  at the C-grid location (this ensure that uVel,vVel field is non
> > divergent).
> >  There are few matlab scripts that uses this 2nd method (e.g., in
> >  MITgcm_contrib/dyncore_ASP/mat_scripts ) and also some fortran code
> >  (e.g., verification/advect_cs/code/ini_vel.F or
> >   verification/code/solid-body.cs-32x32x1/code/ini_vel.F).
> >
> > Cheers,
> > Jean-Michel
> >
> > On Tue, May 01, 2012 at 02:34:45PM +0100, Chris Watkins wrote:
> > > I am trying to set up an initial wind field using the file ini_vel.f on a
> > > cubed sphere. I know analytically what the meridional and zonal winds
> > > should be on a latitude-longitude grid. Can anyone help me understand
> > what
> > > adjustments I need to make to these values to make them into the uVel and
> > > vVel variables.
> > >
> > > I originally set uVel to the zonal wind and vVel to the meridional wind
> > but
> > > the output for U & V  for timestep 0 does not match.
> > >
> > > --
> > > Chris Watkins
> >
> > > _______________________________________________
> > > MITgcm-support mailing list
> > > MITgcm-support at mitgcm.org
> > > http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
>
>
>
> --
> Chris Watkins
>
> G+ GPlus.to/ChrisWatkins <http://gplus.to/ChrisWatkins>

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




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

Message: 2
Date: Wed, 02 May 2012 11:06:37 +0200
From: Gianmaria Sannino <gianmaria.sannino at enea.it>
To: <mitgcm-support at mitgcm.org>
Subject: [MITgcm-support] SMF package status
Message-ID: <CBC6C5BD.A166%gianmaria.sannino at enea.it>
Content-Type: text/plain; charset="us-ascii"

Dear MITgcm support / users,

I would like to know if someone is using ESMF in the 'recent' MITgcm
version.
As far as I know ESMF has been implemented some years ago, but it is not yet
available as official package.

I would be really grateful for some info about the ESMF pkg status. Looking
forward to hearing from you.
Ciao
Gianmaria
--
Gianmaria Sannino, Ph.D
ENEA, Energy and Environment Modeling Unit
Climate and Impact Modeling Lab
CR Casaccia, Bldg C59, Sp. 118, Via Anguillarese, 301
00123 Santa Maria di Galeria - Rome, Italy
Phone: +39 06 3048 6799 (fax 4264)
web: http://utmea.enea.it/people/sannino
email: gianmaria.sannino at enea.it
skype: gianmaria.sannino





==================================================
 MESSAGGIO ISTITUZIONALE
==================================================

DESTINA IL 5X1000 ALL'ENEA
Il Nostro CODICE FISCALE e' 01320740580

L'ENEA utilizzera' il ricavato per il progetto
"interattivambiente giovani":
giovani laureati ed ENEA per lo sviluppo sostenibile l'EFFICIENZA ENERGETICA ed il territorio:
"soluzioni in cerca di problemi"

Per maggiori informazioni: www.enea.it
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120502/7e5a2a12/attachment-0001.htm>

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

Message: 3
Date: Wed, 02 May 2012 11:11:28 +0200
From: Martin Losch <Martin.Losch at awi.de>
To: mitgcm-support at mitgcm.org
Subject: Re: [MITgcm-support] Definition of spatial boundary in EXF
        package
Message-ID: <83173605-8043-430A-89F8-698F28EB2869 at awi.de>
Content-Type: text/plain; CHARSET=US-ASCII

Hi Nicolas,

I think, if your first forcing grid point is a (0,-90) then your second set of parameters is appropriate.

I understand that the documentation is misleading: What it means is, that field_lon/lat0 is defined at the center of the (forcing) grid cell rather than the corner (we are all caught up in the Finite-Volume C-grid world and the description somehow carries over from the model grid, but is not quite right here).

Martin

On May 1, 2012, at 12:24 PM, Nicolas Bruneau wrote:

> Dear MITgcm support / users,
>
> I just started to work with the MITgcm model using the EXF package and I would like to have some clarification about the boundary definitions. From the documentation, I can read:
> field_lon0    <Mail Attachment.png>   starting longitude of input
> field_lon_inc <Mail Attachment.png>   increment in longitude of input
> field_lat0    <Mail Attachment.png>   starting latitude of input
> field_lat_inc <Mail Attachment.png>   increment in latitude of input
> field_nlon    <Mail Attachment.png>   number of grid points in longitude of input
> field_nlat    <Mail Attachment.png>   number of grid points in longitude of input
> If my external forcing is, for example, a 1.5-degree resolution grid starting at 0 degree in longitude and -90 in latitude, do I need to define the spatial grid in data.exf as:
>  atemp_lon0    = .75  (i.e. 0 + dX / 2)
>  atemp_lon_inc = 1.5
>  atemp_lat0    = -89.25
>  atemp_lat_inc = 120*1.5
>  atemp_nlon    = 240
>  atemp_nlat    = 121
>
> or do I define it as:
>  atemp_lon0    = 0
>  atemp_lon_inc = 1.5
>  atemp_lat0    = -90
>  atemp_lat_inc = 120*1.5
>  atemp_nlon    = 240
>  atemp_nlat    = 121
>
> When I'm looking the subroutine exf_interp.F, it looks that the "x_in" variable represents the coordinates of the input and in such a case, I would think the second definition (atemp_lon0 = 0) is the good one but following the documentation, it looks that the first one is the good one. So I imagine, I didn't get something.
>
> I would be really grateful for some help. Looking forward to hearing from you.
>
> --
> <baslogo1000.gif>     Dr Nicolas Bruneau
> British Antarctic Survey
> High Cross, Madingley Road, Cambridge, CB3 0ET, UK
> Email: nicbru at bas.ac.uk
> Phone: +44 (0) 1223 22 13 27
> --
> This message (and any attachments) is for the recipient only. NERC
> is subject to the Freedom of Information Act 2000 and the contents
> of this email and any reply you make may be disclosed by NERC unless
> it is exempt from release under the Act. Any material supplied to
> NERC may be stored in an electronic records management system.
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




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

Message: 4
Date: Wed, 2 May 2012 06:29:02 -0400
From: Chris Hill <cnh at mit.edu>
To: mitgcm-support at mitgcm.org
Subject: Re: [MITgcm-support] SMF package status
Message-ID:
        <CAOLSJZXKEh3qhwWHsDa4iYDMdJT9piKGoP=zOc9L_Jyvv+5yAg at mail.gmail.com>
Content-Type: text/plain; charset=windows-1252

Ciao Gianmaria,

 Thats exciting.
 Do you know what you want to use ESMF with?

 There is a working example of using it in the contrib area

 http://mitgcm.org/viewvc/MITgcm/MITgcm_contrib/PRM/

 I was hoping for two things, before it would be more integrated
   (1) for ESMF interfaces to stabilize in ESMF _and_
   (2) for ESMF to get to point where it could be installed easily

 Item (1) took about a decade, but did happen a while back.
 Still waiting on item (2).

Chris

On Wed, May 2, 2012 at 5:06 AM, Gianmaria Sannino
<gianmaria.sannino at enea.it> wrote:
> Dear?MITgcm support / users,
>
> I would like to?know if?someone is using?ESMF in the 'recent' MITgcm
> version.
> As far as I?know?ESMF?has been?implemented?some years?ago,?but it is not
> yet?available as?official package.
>
> I would be really grateful for some info about the ESMF pkg status. Looking
> forward to hearing from you.
> Ciao
> Gianmaria
> --
> Gianmaria Sannino, Ph.D
> ENEA, Energy and Environment Modeling Unit
> Climate and?Impact Modeling Lab
> CR Casaccia,?Bldg C59, Sp. 118,?Via Anguillarese, 301
> 00123 Santa Maria di Galeria - Rome, Italy
> Phone: +39 06 3048 6799 (fax 4264)
> web:?http://utmea.enea.it/people/sannino
> email:?gianmaria.sannino at enea.it
> skype: gianmaria.sannino
>
>
> ________________________________
> MESSAGGIO ISTITUZIONALE
> ________________________________
>
> Destina il 5 x 1000 all?ENEA
> Il nostro codice fiscale ? 01320740580
>
> L'ENEA utilizzer? il ricavato per il progetto
> "interattivambiente giovani":
> giovani laureati ed ENEA per lo sviluppo sostenibile l'EFFICIENZA ENERGETICA
> ed il territorio:
> "soluzioni in cerca di problemi"
>
> Per maggiori informazioni: www.enea.it
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>



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

Message: 5
Date: Wed, 02 May 2012 15:04:59 +0200
From: Gianmaria Sannino <gianmaria.sannino at enea.it>
To: "Mitgcm-Support at Mitgcm. Org" <mitgcm-support at mitgcm.org>
Subject: Re: [MITgcm-support] SMF package status
Message-ID: <CBC6F4C4.A311%gianmaria.sannino at enea.it>
Content-Type: text/plain;       charset="ISO-8859-1"

Hi Chris,

As you probably know, we've implemented a regional climate model for the
Mediterranean region.
We coupled MITgcm with RegCM (Atmo) via OASIS3.0
Here is a reference to our first work
An atmosphere?ocean regional climate model for the Mediterranean area:
assessment of a present climate simulation
http://www.springerlink.com/content/4k24648p5633n8w2/

After that, we developed NH high resolution gibraltar model
Numerical modelling of three-dimensional stratified tidal flow over
Camarinal Sill, Strait of Gibraltar
http://www.agu.org/pubs/crossref/pip/2011JC007093.shtml

Here you can find some very nice video of the simulation on youtube:

http://www.youtube.com/watch?v=_OU84SrS9sE
http://youtu.be/VBrCwnULcfc?hd=1


We transferred all these experiences to a new version of the Mediterranean
model.
A 1/16? with a special zoom (1/200?) in the region of Gibraltar. The model
is forced laterally and internally by tides.
These are some movies showing tides evolution in our Mediterranean model:
http://www.youtube.com/watch?v=UZpV_BU_qbg
http://www.youtube.com/watch?feature=endscreen&v=RRHqcqFCxZU&NR=1


Now we would like to couple this new Med. version to a higher resolution
atmospheric model (RegCM and NMMB).

Unfortunately the OASIS4 parallel version is no more developed, and so
we're looking for a different coupler
ables to manage larger grid without affecting scalability.
It seems that the new ESMF release (5.3.0) could be the right tool.

Ciao
Gianmaria



--
Gianmaria Sannino, Ph.D
ENEA, Energy and Environment Modeling Unit
Climate and Impact Modeling Lab
CR Casaccia, Bldg C59, Sp. 118, Via Anguillarese, 301
00123 Santa Maria di Galeria - Rome, Italy
Phone: +39 06 3048 6799 (fax 4264)
web: http://utmea.enea.it/people/sannino
email: gianmaria.sannino at enea.it
skype: gianmaria.sannino









On 02/05/12 12:29, "Chris Hill" <cnh at mit.edu> wrote:

>Ciao Gianmaria,
>
> Thats exciting.
> Do you know what you want to use ESMF with?
>
> There is a working example of using it in the contrib area
>
> http://mitgcm.org/viewvc/MITgcm/MITgcm_contrib/PRM/
>
> I was hoping for two things, before it would be more integrated
>   (1) for ESMF interfaces to stabilize in ESMF _and_
>   (2) for ESMF to get to point where it could be installed easily
>
> Item (1) took about a decade, but did happen a while back.
> Still waiting on item (2).
>
>Chris
>
>On Wed, May 2, 2012 at 5:06 AM, Gianmaria Sannino
><gianmaria.sannino at enea.it> wrote:
>> Dear MITgcm support / users,
>>
>> I would like to know if someone is using ESMF in the 'recent' MITgcm
>> version.
>> As far as I know ESMF has been implemented some years ago, but it is not
>> yet available as official package.
>>
>> I would be really grateful for some info about the ESMF pkg status.
>>Looking
>> forward to hearing from you.
>> Ciao
>> Gianmaria
>> --
>> Gianmaria Sannino, Ph.D
>> ENEA, Energy and Environment Modeling Unit
>> Climate and Impact Modeling Lab
>> CR Casaccia, Bldg C59, Sp. 118, Via Anguillarese, 301
>> 00123 Santa Maria di Galeria - Rome, Italy
>> Phone: +39 06 3048 6799 (fax 4264)
>> web: http://utmea.enea.it/people/sannino
>> email: gianmaria.sannino at enea.it
>> skype: gianmaria.sannino
>>
>>
>> ________________________________
>> MESSAGGIO ISTITUZIONALE
>> ________________________________
>>
>> Destina il 5 x 1000 all?ENEA
>> Il nostro codice fiscale ? 01320740580
>>
>> L'ENEA utilizzer? il ricavato per il progetto
>> "interattivambiente giovani":
>> giovani laureati ed ENEA per lo sviluppo sostenibile l'EFFICIENZA
>>ENERGETICA
>> ed il territorio:
>> "soluzioni in cerca di problemi"
>>
>> Per maggiori informazioni: www.enea.it
>>
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>
>_______________________________________________
>MITgcm-support mailing list
>MITgcm-support at mitgcm.org
>http://mitgcm.org/mailman/listinfo/mitgcm-support




==================================================
 MESSAGGIO ISTITUZIONALE
==================================================

DESTINA IL 5X1000 ALL'ENEA
Il Nostro CODICE FISCALE e' 01320740580

L'ENEA utilizzera' il ricavato per il progetto
"interattivambiente giovani":
giovani laureati ed ENEA per lo sviluppo sostenibile l'EFFICIENZA ENERGETICA ed il territorio:
"soluzioni in cerca di problemi"

Per maggiori informazioni: www.enea.it




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

_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
http://mitgcm.org/mailman/listinfo/mitgcm-support


End of MITgcm-support Digest, Vol 107, Issue 3
**********************************************


The University of Dundee is a registered Scottish Charity, No: SC015096




More information about the MITgcm-support mailing list