[MITgcm-support] rbcs advice

Jean-Michel Campin jmc at ocean.mit.edu
Wed Mar 5 14:05:14 EST 2014


Hi Jody,
We will see what we can do about this and will let you know
of any update.
Cheers,
Jean-Michel

On Tue, Mar 04, 2014 at 07:48:36PM -0800, Jody Klymak wrote:
> Hi again,
> 
> OK, to answer my own question, no, there is no equivalent of rbcsForcingOffset for obcs.  
> 
> Can I suggest in obcs_fields_load.F there be an offset variable, and that this be added to OBCS_PARAMS.h and obcs_readparms.F?  
> 
>       CALL GET_PERIODIC_INTERVAL(
>      O                  iRecP, iRec0, iRec1, bWght, aWght,
>      I                  externForcingCycle, externForcingPeriod,
>      I                  deltaTclock, myTime-obcsForcingOffset, myThid )
> 
> If you want to use a snapshot, then, it is just a matter of setting 
> 
> obcsForcingOffset=-externForcingPeriod/2
> 
> Or, one could imagine this being called "externForcingOffset" and being used in "external_fields_load.F" as well....
> 
> Cheers,   Jody
> 
> 
> 
> On Mar 4, 2014, at  17:52 PM, Jody Klymak <jklymak at uvic.ca> wrote:
> 
> > One other quick question that I know has come up before but I can't seem to find...
> > 
> > I have snapshots at t=0,3600,7200, and I want to use them for my obcs.  However, I seem to remember that externalForcingPeriod=3600 assumes the first file entry is at t=1800, the second at 5400s etc.  Is there an offset to change this so that the first is t=0, etc, or do I have to hack get_periodic_interval.F?  
> > 
> > Thanks,   Jody
> > 
> > 
> > On Mar 4, 2014, at  9:18 AM, Jody Klymak <jklymak at uvic.ca> wrote:
> > 
> >> Hi Dimitris,
> >> 
> >> I wrote you offline as well, but...
> >> 
> >> On Mar 4, 2014, at  9:06 AM, Dimitris Menemenlis <dmenemenlis at gmail.com> wrote:
> >> 
> >>> Hi Jody, I think yes.  What grid spacing are you running your regional domain?
> >>> And what resolution (space and time) is your open boundary conditions?
> >> 
> >> Right now I have hourly (well 3720 s to be an even divisor of an M2 tidal cycle) output on a 5- to 1-km grid, with 1 km near my region of interest, and 5 km at the boundaries.  The forcing of this model was an analytical source function (a couple of internal tide line-sources outside the domain), and the model was spun up for 12 tidal cycles allowing reflections from a rough continental slope to set up. Probably accurate out to mode 10 or so.  
> >> 
> >> I am trying to telescope this down to about 100x200 m in a 50x100 km region to look at the non-linear internal waves on a continental slope.  I nest this inside a 400x300 km region using telescoping resolution, and will initialize and force with the larger 1-km grid. 
> >> 
> >> All the telescoping is to help make the BCs not too important, but it would be nice to make them as close as possible to reduce reflections. 
> >> 
> >> Thanks!   Jody
> >> 
> >>> 
> >>> Michael Schodlok has a lot of experience going from
> >>> CS510 (18 km) to 1 km domains, then adding tides,
> >>> and possibly obcs sponge to get rid of some edge artifacts.
> >>> 
> >>> But with sufficiently high-resolution (in time) obcs,
> >>> I am guessing that you will not need sponging at all.
> >>> At least that's the theory.
> >>> Chris Hill was experimenting with "replaying" bits
> >>> of the 1/48th using obcs.  I don't know how far he got.
> >>> 
> >>> Let me know if you want to access the 1/48th output
> >>> and we can figure out best way to proceed.
> >>> 
> >>> Cheers, Dimitris
> >>> 
> >>> On Mar 4, 2014, at 8:50 AM, Jody Klymak <jklymak at uvic.ca> wrote:
> >>> 
> >>>> Do you think its true that I could just force w/ obcs and do away w/ the sponge? I've used sponges for so long I hadn't thought of giving this a try, but maybe it is all I need to do, and I can ditch rbcs.  Certainly my forcing should have all the modes that could conceivably reach the boundaries by the time my integration is done.  
> >>>> 
> >>>> Second, 1/48th degree global with tides!  That would be very interesting!
> >>> 
> >>> 
> >>> _______________________________________________
> >>> MITgcm-support mailing list
> >>> MITgcm-support at mitgcm.org
> >>> http://mitgcm.org/mailman/listinfo/mitgcm-support
> >> 
> >> --
> >> Jody Klymak    
> >> http://web.uvic.ca/~jklymak/
> > 
> > 
> > _______________________________________________
> > 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



More information about the MITgcm-support mailing list