[MITgcm-support] sponge layer

[Stephanie Dutkiewicz]

David -
I recently checked in a pkg/rbcs which essentially does what
I think you are after, but only for ptracers. The packages reads
in a three-d field that will set the relaxation timescale (so this 
timescale can be spatially variable). I set this pkg up to read in 
temperature and salinity files to relax to...but never added the
code to do so
(didn't want to offend any coding guru's by implementing such
  a crude sponge layer for T/S when the OBCS were there...), but
if there is a need for such a thing, it will be easy to do...adding
something like rbcs_add_tendency for t and s.

Anyone object?

steph
ps. I include a sample data.rbcs so you can get an idea of the
switches.

On Wed, 14 Dec 2005, David Wang wrote:

> Hello...
>
> I'm going to implement "sponge layers" near artificial domain boundaries 
> within which model temperature may be restored to some temperature 
> climatology (say Levitus94) at different relaxation timescales. I've already 
> checked out the sponge layer option in the obcs pkg, where the relaxation 
> temperature is the linear combination of model temperature and the prescribed 
> temperature at the boundary. This is not really what I want. If I need to 
> write some extra code to achieve my purpose (which is not uncommon), could 
> anyone provide any guidelines... or did I miss anything?
>
> btw, I did notice the 3-dim relaxation code emerged in the latest checkpoint, 
> but its relaxation timescale is set to be spatically constant...
>
> Thanks, David
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
# Relax Boundary Conditions parameters
 &RBCS_PARM01
 tauRelax = 2592000.,
 relaxBoundaryFile = 'mask1.bin',
 relaxTFile = '',
 relaxSFile = '',
 rbcsIniter = 0,
 rbcsForcingPeriod = 2592000.,
 rbcsForcingCycle = 15552000.,
 useRBCtemp = .FALSE.,
 useRBCsalt = .FALSE.,
 useRBCptracers = .TRUE.,
 &end

 &RBCS_PARM02
 useRBCptrnum(1)=.TRUE.
 useRBCptrnum(2)=.FALSE.
 relaxPtracerFile(1) = 'ptr1_3.bin',
 relaxPtracerFile(2) = ' ',
 &end