[MITgcm-support] Problem with input data

[Matthew Mazloff]

Update:

The data is read in correctly as the output shows

 Field Initial Temperature at iteration          1
(PID.TID 0000.0001) // CMIN =          4.057084608000000E-04
(PID.TID 0000.0001) // CMAX =          4.093527510783521E-04
(PID.TID 0000.0001) // CINT =          1.349737140130389E-07

However at the first checkpoint, after nothing has happened but a 
calculation of the pressure (I am running hydrostatic) the output reads

(PID.TID 0000.0001) %MON dynstat_theta_max            =   
4.0935275107835E-04
(PID.TID 0000.0001) %MON dynstat_theta_min            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   
2.0377073219221E-04

note theta_min has gone to zero!

What could be happening?

Thanks,
Matt


Matthew Mazloff wrote:

> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support