| 8 |
These routines are hard-coded for running in |
These routines are hard-coded for running in |
| 9 |
nireas:/home/dimitri/mitgcm/setups/cube/510 |
nireas:/home/dimitri/mitgcm/setups/cube/510 |
| 10 |
Contact menemenlis@jpl.nasa.gov for questions, etc. |
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 |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch