[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