[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