[MITgcm-support] FW: another problem

Alistair Adcroft adcroft at MIT.EDU
Sat Oct 11 17:29:03 EDT 2003


Mike,

Given the correct sign of the input topography, it could be that the input
files are not "Big Endian" and r4/r8 is mis-matched or that they are correct
but the compile options are wrong and expecting "Little endian". As a rule
we work in big endian. The readbinaryprec parameter (check it's name in
set_defaults.F) controls whether the model reads real*4 or real*8.

I'm cc'ing this to mitgcm-support at mitgcm.org since I'm out of touch for a
week. Jean-Michel knows the R_low businesss best since he wrote it.

A.
--
Dr Alistair Adcroft            http://www.mit.edu/~adcroft
MIT Climate Modeling Initiative        tel: (617) 253-5938
EAPS 54-1523,  77 Massachusetts Ave,  Cambridge,  MA,  USA

-----Original Message-----
From: spall at whoi.edu [mailto:spall at whoi.edu] On Behalf Of Mike Spall
Sent: Thursday, October 09, 2003 12:02 PM
To: Adcroft, Alistair
Subject: another problem


Hi Alistair,

Sorry to say that I am still having a problem getting the
model running.  The first sign of trouble seems to be in
the R_low array:

(PID.TID 0000.0001)  MDSREADFIELD: opening global file: topog.bin  rlow 1
-5.4765625  rlow 2 -5.4765625  rlow 3 -5.4765625 (PID.TID 0000.0001) //
=======================================================
(PID.TID 0000.0001) // Field Model R_low (ini_masks_etc) at
iteration          1
(PID.TID 0000.0001) // CMIN =          1.000000000000000E+32
(PID.TID 0000.0001) // CMAX =         -1.000000000000000E+32
(PID.TID 0000.0001) // CINT =          0.000000000000000E+00
(PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX):
-abcdefghijklmnopqrstuvwxyz+
(PID.TID 0000.0001) //                  0.0: .
(PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):(  -2:  34:   1)
(PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):(  34:  -2:  -1)
(PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):(   1:   1:   1)
(PID.TID 0000.0001) //
=======================================================


R_low seems to have reasonable values when read in (my print statements) but
is set to 0.0 for reasons I have not been able to figure out. This then goes
on to set the initial temperature to zeros as well. Any ideas?

At least once I get going, I pretty much leave you alone!

Mike

-- 
==================================================

Michael A. Spall
Department of Physical Oceanography
W.H.O.I. MS #21
360 Woods Hole Road
Woods Hole, MA   02543

mspall at whoi.edu
(508) 289-3342 (office)
(508) 457-2181 (fax)

=================================================





More information about the MITgcm-support mailing list