[MITgcm-support] Problem with input data

Matthew Mazloff mmazloff at MIT.EDU
Thu May 20 11:24:32 EDT 2004


Hi Chris,

I have only tried it on a mac, but I don't believe this is the problem.  
The input data is read fine.  I don't know much about the masks, but I 
now believe the temperature mask may be the problem.  I see that the dry 
area of my model has the reference temperature after several time 
steps.  So initially all dry area has zero temperature.  All but 3/4 of 
one wet layer is good, but there is that one layer that is set to zero 
in the first correction step.  After 1 time step this wet area acts fine 
again (But the damage is done to the initial conditions)  and then at 
some point, im not sure where yet, the all the dry area except for the 
same layer that was messed up before, takes on the reference 
temperature! This is then maintained throughout the run.  All other 
fields are zero in the dry region.  I have attached output with a single 
timestep.  I don't know if it will help. Any ideas as to what I should 
try next?

Thanks so much,
Matt



Chris Hill wrote:

>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
>>
>>    
>>
>
>_______________________________________________
>MITgcm-support mailing list
>MITgcm-support at mitgcm.org
>http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
>  
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20040520/c8ba3026/attachment.htm>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: STDOUT.txt
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20040520/c8ba3026/attachment.txt>


More information about the MITgcm-support mailing list