high_res_cube/matlab/readme.txt

bathymetry file ETOPO2_510x6x510_ver5_filled.bin
for cs510 integrations was generated using mk_bathy3.m
followed by mk_bathy5.m

Levitus initial T and S and monthly SSS files for
cs510 integrations were generated using mk_input.m

These routines are hard-coded for running in
nireas:/home/dimitri/mitgcm/setups/cube/510
Contact menemenlis@jpl.nasa.gov for questions, etc.

=========================================================================
=========================================================================

Files contributed by Jean-Michel on 3/27/2007 and used for computation of
barotropic streamfunction:

psiLine_N2S_cs510.mat
mk_psiLine_CS.m
use_psiLine.m
v_integrate.m
psi_baro.ps.gz

>>>>>>>>>>>>>>

Hi Michael,

I put the *.m files that I was using in:
http://mitgcm.org/~jmc/
(Dimitris: might put them on CVS server later on, 
 may be in MITgcm_contrib/high_res/matlab/something ?)

"mk_psiLine_CS.m" is not really different from 
MITgcm/utils/matlab/cs_grid/bk_line (except kgr=0)

I added the psiLine_N2S_cs510.mat file (from matlab 7, 
don't think it will work with matlab 6).

v_intergrate.m (was part of use_psiLine.m) read 12 
monthly mean + integrate vertical to get the 
annual mean, vertically integrated transport.

and use_psiLine.m is just computing the stream-function 
+ making a plot.

I've just checked-in grph_CSz.m so that it works also for
CS-510 grid.

Jean-Michel

On Tue, Mar 27, 2007 at 06:37:12PM -0400, Michael Schodlok wrote:
> > 
> > Hi Jean-Michel,
> > 
> > thanks for checking this out, just got in and am still a bit
> > puzzled why mk_psiLine_cs.m didn't work here.
> > The *.mat file and modified m-files would be great.
> > Do you want put them on the ECCO2 or  the
> > MITgcm webpage?
> > 
> > cheers
> > 
> > Michael
> > 
>> > >
>> > >I tried "mk_psiLine_CS.m" on cs-510 grid, and it seems
>> > >to work (appart from a stupid plot which you can turn off
>> > >by changing kgr=3; to kgr=0; on the 1rst line).
>> > >I can put the *.mat file (psiLine_N2S_cs510.mat) somewhere 
>> > >if you want to take it.
>> > >
>> > >After I tried to compute a stream-function, using:
>> > >utils/matlab/cs_grid/bk_line/use_psiLine.m
>> > >(had to change the loading of UVELMASS & VVELMASS,
>> > > that was taking too much memory), and it works.
>> > >But the plot at the end is not right (a problem in
>> > >grph_CSz.m that I will look at later).
>> > >Again, if you want me to put those modified script somewhere, 
>> > >just tell me.
>> > >
>> > >I attached here a postscript (but gzip) of the plot, 
>> > >it's average over 1 year.
>> > >
>> > >Jean-Michel
>> > >
>> > >On Tue, Mar 27, 2007 at 09:47:25AM -0400, Jean-Michel Campin wrote:
>> > >  
>>> > >>Hi Michael,
>>> > >>
>>> > >>I will have a look. I remember that I checked the Meridional
>>> > >>stream-function, but no idea regarding this one.
>>> > >>Jean-Michel
>>> > >>
>>> > >>PS: I cc to others who might be interested.
>>> > >>
>>> > >>On Mon, Mar 26, 2007 at 06:07:25PM -0400, Michael Schodlok wrote:
>>> > >>    
>>>> > >>>Hi Jean-Michel
>>>> > >>>
>>>> > >>>I tried to use work out the streamfunction for cube510
>>>> > >>>using calcHorizPsiCube.m and run into all sorts of
>>>> > >>>problems. Well, the main problem I guess is using
>>>> > >>>the mk_psiLine_CS.m routine which seems to be
>>>> > >>>for cube32. I get reshape errors when I use it in
>>>> > >>>split_Z_cub.m. We didn't save the meta files, so
>>>> > >>>I tried to work around that but am a bit lost what
>>>> > >>>routine does what etc. Do you have any suggestions
>>>> > >>>how else to calculate streamfunction in cube?
>>>> > >>>
>>>> > >>>thanks
>>>> > >>>
>>>> > >>>Michael
1	bathymetry file ETOPO2_510x6x510_ver5_filled.bin
2	for cs510 integrations was generated using mk_bathy3.m
3	followed by mk_bathy5.m
4
5	Levitus initial T and S and monthly SSS files for
6	cs510 integrations were generated using mk_input.m
7
8	These routines are hard-coded for running in
9	nireas:/home/dimitri/mitgcm/setups/cube/510
10	Contact menemenlis@jpl.nasa.gov for questions, etc.
11
12	=========================================================================
13	=========================================================================
14
15	Files contributed by Jean-Michel on 3/27/2007 and used for computation of
16	barotropic streamfunction:
17
18	psiLine_N2S_cs510.mat
19	mk_psiLine_CS.m
20	use_psiLine.m
21	v_integrate.m
22	psi_baro.ps.gz
23
24	>>>>>>>>>>>>>>
25
26	Hi Michael,
27
28	I put the *.m files that I was using in:
29	http://mitgcm.org/~jmc/
30	(Dimitris: might put them on CVS server later on,
31	may be in MITgcm_contrib/high_res/matlab/something ?)
32
33	"mk_psiLine_CS.m" is not really different from
34	MITgcm/utils/matlab/cs_grid/bk_line (except kgr=0)
35
36	I added the psiLine_N2S_cs510.mat file (from matlab 7,
37	don't think it will work with matlab 6).
38
39	v_intergrate.m (was part of use_psiLine.m) read 12
40	monthly mean + integrate vertical to get the
41	annual mean, vertically integrated transport.
42
43	and use_psiLine.m is just computing the stream-function
44	+ making a plot.
45
46	I've just checked-in grph_CSz.m so that it works also for
47	CS-510 grid.
48
49	Jean-Michel
50
51	On Tue, Mar 27, 2007 at 06:37:12PM -0400, Michael Schodlok wrote:
52	> >
53	> > Hi Jean-Michel,
54	> >
55	> > thanks for checking this out, just got in and am still a bit
56	> > puzzled why mk_psiLine_cs.m didn't work here.
57	> > The *.mat file and modified m-files would be great.
58	> > Do you want put them on the ECCO2 or the
59	> > MITgcm webpage?
60	> >
61	> > cheers
62	> >
63	> > Michael
64	> >
65	>> > >
66	>> > >I tried "mk_psiLine_CS.m" on cs-510 grid, and it seems
67	>> > >to work (appart from a stupid plot which you can turn off
68	>> > >by changing kgr=3; to kgr=0; on the 1rst line).
69	>> > >I can put the *.mat file (psiLine_N2S_cs510.mat) somewhere
70	>> > >if you want to take it.
71	>> > >
72	>> > >After I tried to compute a stream-function, using:
73	>> > >utils/matlab/cs_grid/bk_line/use_psiLine.m
74	>> > >(had to change the loading of UVELMASS & VVELMASS,
75	>> > > that was taking too much memory), and it works.
76	>> > >But the plot at the end is not right (a problem in
77	>> > >grph_CSz.m that I will look at later).
78	>> > >Again, if you want me to put those modified script somewhere,
79	>> > >just tell me.
80	>> > >
81	>> > >I attached here a postscript (but gzip) of the plot,
82	>> > >it's average over 1 year.
83	>> > >
84	>> > >Jean-Michel
85	>> > >
86	>> > >On Tue, Mar 27, 2007 at 09:47:25AM -0400, Jean-Michel Campin wrote:
87	>> > >
88	>>> > >>Hi Michael,
89	>>> > >>
90	>>> > >>I will have a look. I remember that I checked the Meridional
91	>>> > >>stream-function, but no idea regarding this one.
92	>>> > >>Jean-Michel
93	>>> > >>
94	>>> > >>PS: I cc to others who might be interested.
95	>>> > >>
96	>>> > >>On Mon, Mar 26, 2007 at 06:07:25PM -0400, Michael Schodlok wrote:
97	>>> > >>
98	>>>> > >>>Hi Jean-Michel
99	>>>> > >>>
100	>>>> > >>>I tried to use work out the streamfunction for cube510
101	>>>> > >>>using calcHorizPsiCube.m and run into all sorts of
102	>>>> > >>>problems. Well, the main problem I guess is using
103	>>>> > >>>the mk_psiLine_CS.m routine which seems to be
104	>>>> > >>>for cube32. I get reshape errors when I use it in
105	>>>> > >>>split_Z_cub.m. We didn't save the meta files, so
106	>>>> > >>>I tried to work around that but am a bit lost what
107	>>>> > >>>routine does what etc. Do you have any suggestions
108	>>>> > >>>how else to calculate streamfunction in cube?
109	>>>> > >>>
110	>>>> > >>>thanks
111	>>>> > >>>
112	>>>> > >>>Michael