[MITgcm-devel] MNC "tile assembly" : assemble per-file tiles into a global grid
Martin Losch
mlosch at awi-bremerhaven.de
Tue Dec 14 11:29:15 EST 2004
Hi Ed,
I figured you'd like some feedback on your mnc_assembly.m.
I tried your mnc_assembly script and FAILED (o:
1. I think you should replace all occurances of "&&" and "||" with "&"
and "|" if you want it to work with Matlab6.1 (don't laugh). After
making those changes, I could run the test_assem_test.m script
successfully, that is without obvious errors.
2. In your documentation (help mnc_assembly) your call your second
argument a "cell" array, but in fact it is a "struct", as can be seen
from test_assem_test.m (a cell array doesn't work, it took me some time
to figure that out).
3. I can now run everything on my output (20CPUs, 2D vertical slab,
much like plume_one_slope), but the resulting all.00000.nc contains
only the values of one tile/CPU (some tile in the middle). I cannot
figure out why, so I put now a 140MB gzipped archive on
sea.mit.edu:/net/sea2/scratch/mlosch/for_ed.tgz for you to grab and see
what's going, if you find the time. I tried this:
>> vars=struct([]);
>> vars(1).name='Temp';
>>vars(2).name='W';
>>[n,m]=mnc_assembly('./state.0000.%06d.nc',vnames)
Found 20 files matching the pattern: "./state.0000.%06d.nc"
Looking ...
Copying ...
n=0
m=0
>>
Martin
On Dec 10, 2004, at 8:56 PM, Ed Hill wrote:
>
> Hi folks,
>
> A first working version of a MatLAB script that does "MNC tile
> assembly"
> for both the "exch1" and "exch2" tile geometries is now available. You
> can get it along with two example data sets at:
>
> $ cvs co development/edhill/mnc_assembly
>
> and run it using:
>
> $ cd development/edhill/mnc_assembly
> $ matlab
>>> mnc_assem_test
>
> There are numerous things that should (and will!) be done to improve
> the
> script:
>
> - respect the 2GB limit
> - allow for selection of arbitrary T and Z levels
> - better handling of attributes and coord-vars
>
> but this is just a first working (Beta?) version.
>
> Ed
>
> --
> Edward H. Hill III, PhD
> office: MIT Dept. of EAPS; Rm 54-1424; 77 Massachusetts Ave.
> Cambridge, MA 02139-4307
> emails: eh3 at mit.edu ed at eh3.com
> URLs: http://web.mit.edu/eh3/ http://eh3.com/
> phone: 617-253-0098
> fax: 617-253-4464
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-devel
More information about the MITgcm-devel
mailing list