[MITgcm-support] Problem with input data
Chris Hill
cnh at mit.edu
Thu May 20 10:01:32 EDT 2004
Matt,
Is it only on the mac? Also, can you put a full STDOUT somewhere on the
web, or attached to a mail.
Thanks,
CHris
> -----Original Message-----
> From: mitgcm-support-bounces at mitgcm.org
> [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of
> Matthew Mazloff
> Sent: Wednesday, May 19, 2004 7:57 PM
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] Problem with input data
>
> 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
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>
More information about the MITgcm-support
mailing list