[MITgcm-support] External forcing

[Martin Losch]

Elja,
has anyone answered you yet?

As far as I know there is technology available to read 3d fields that  
you could use for your purpose (e.g. see the rbcs-pkg where 3D fields  
are used for restoring), but I don't think that you can just turn on  
some already coded mechanism to import and apply 3D+time forcing.

Do you have global problems in mind? Then it's probably easiest to  
compute the tidal potential online and apply it as a pressure force  
(add it to the phi0surf in external_forcing_surf.F). Alternatively  
you could (ab-)use the surface pressure field for applying  
equilibrium tides (there is an input field (pload and ploadfile) for  
that; have a look at PARAMS.h, FFIELDS.h).
For regional applications you can use the obcs pkg to for the  
boundaries with tidal velocities.

In both case there are no LSA (loading and self attraction) terms  
available, yet, (but I've heard rumors that they might be coming soon  
in this theater).

Martin


On 17 Mar 2008, at 17:59, Elja Huibregtse wrote:
> Dear all,
>
> In order to simulate tidal forcing, I would like to apply time  
> dependent periodic external forcing to MITgcm. So far, I found the  
> option to use (x,y,t)-files. This data only holds for the surface,  
> I assume? It is, however, necessary to apply the tidal forcing to  
> all layers.
>
> Is it, therefore, possible to apply external forcing datasets,  
> ordered in a (x,y,z,t) fashion?
>
> Thank you.
> Best regards,
>
> Elja
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support