[MITgcm-support] Problem with input data

Matthew Mazloff mmazloff at MIT.EDU
Wed May 19 14:04:42 EDT 2004


I am running the MITgcm on a mac.  With Matlab, using gendata.m,  I 
create input files such as

fid=fopen('Ulinear.rect','w','b'); fwrite(fid,U,'float64'); fclose(fid);

I am using readbinaryprec = 64
All fields input fine, except for my the temperature field.  This field 
has a chunk of the data in the middle of the file set equal to zero.
I have tried other machine formats ('s', 'l', 'a', 'c'), and precisions, 
('double'), but the problem  has only been made worse.  The temperature 
field has a larger amplitude than the velocity fields.  But this is the 
only difference.  The place where the field is set equal to zero is the 
same each time.  I have checked that Matlab is writing the field 
correctly using fread, and it is fine.

Any ideas?
Thanks
Matt




More information about the MITgcm-support mailing list