[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