[MITgcm-support] single profile output

[Jean-Michel Campin]

Hi Li,

On Mon, Apr 06, 2009 at 06:52:48PM -0700, q li wrote:
> I have two practical problem, maybe easy for everyone.
> 
> First, I only need output of T (Temperature). How can I close the output for
> others, such as U, V, W?

You can disable all snap-shot output (by setting dumpFreq=0., in file "data")
and then use the diagnostics package, which allow to select which variable
you want. If you specify a positive frequency, it will be time-averaged,
with a negative frequency (+ pick a "timePhase") it's just instantaneous
output:
#  frequency(n):< 0 : write snap-shot output every |frequency| seconds
#               > 0 : write time-average output every frequency seconds

> 
> Second, I also need temperature output at one single point every 1 minute,
> in order to compare with our observations. It looks MITgcm deosn't give you
> too much flexibility to do that. I also cannot find how to revise this
> numerical model from the manual.

A way to do this would be to define a region-mask which corresponds to just
one grid point, and then to use the "stat-diags" (part of the diagnostics
package, 2nd namelist in file "data.diagnostics"). You can have a look at
 MITgcm/verification/global_ocean.cs32x15/input.thsice/data.diagnostics
I know that it has been done for single or multiple profile output.

Thanks,
Jean-Michel

> 
> I would appreciate if anyone could take a minute to answer it. It must be
> too simple for the experts.
> 
> Li
> 
> 
> 
>       
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support