[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