[MITgcm-support] how to write
Nataliya Stashchuk
nataliya.stashchuk at plymouth.ac.uk
Thu Oct 28 08:01:35 EDT 2010
Hi Jean-Michel,
I tried to follow your advice and found that I still have some questions.
I am interesting what value of velocity would be in point x=199500 with time.
Here is my file:
&DIAG_STATIS_PARMS
#- regional mask: band: 1 : x <= 199500 ; 2 : x == 199500
diagSt_regMaskFile='region1.bin',
nSetRegMskFile=1,
set_regMask(1)= 1,
val_regMask(1)= 2.,
#- an example just to check the agreement with MONITOR output:
stat_fields(1,1)= 'UVEL ',
stat_fname(1)= 'u2',
stat_freq(1)= 25,
stat_phase(1)= 0.,
&
1) It seems to me that the programme does not understand where my point is.
My matlab file is:
prec='real*8';
ieee='b';
nx=1024;
ny=1;
nz=85;
xc=rdmds('XC');
regMsk=ones(size(xc));
regMsk(find (199500 == xc))=2;
fid=fopen('region1.bin','w',ieee); fwrite(fid,regMsk,prec); fclose(fid);
2) My output file gives following: Average, Std.Dev, min, max, Vol. Where should be the data for my point?
3) I use 85 levels in vertical. In my output I need to know velocity at 2 or three levels. How I can organize this?
Many thanks,
Nataliya
-----Original Message-----
From: Jean-Michel Campin [mailto:jmc at ocean.mit.edu]
Sent: 27 October 2010 20:45
To: mitgcm-support at mitgcm.org
Subject: Re: [MITgcm-support] how to write
Hi Nataliya,
You can use the STATDIAGS (part of diagnostics pkg) with a region
that cover only 1 grid point and, setting the stat_freq equal to
the time-step, will produce output at 1 grid point every time-step.
This is not well documented in the manual, but there is an example
using STATDIAGS in: verification/aim.5l_cs/input.thSI
can check data.diagnostics there and matlab script mk3regions_mask.m
Thanks,
Jean-Michel
On Wed, Oct 27, 2010 at 07:28:20PM +0100, Nataliya Stashchuk wrote:
> Dear mitgcm users,
>
> Can you give me an advise how to arrange an output file which should contain a physical value, i.e. velocity, for instance, in a single specific point, but written at every temporal step, i.e. with delT frequency? I need this time series for a spectral analysis.
>
> Thank you,
>
> Nataliya
> _______________________________________________
> 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