[MITgcm-support] Re: eosType!

[Ed Hill]

Hi Van Thinh,

The following is the body of your email with the ~43MB of attachments
removed since our email list has a 1MB upper limit.  In future, please
put all your large images/attachments on a web server and send _only_
URLs referring to them.  If you have trouble finding a web server for
your large files then please contact us and we can help.

thanks,
Ed

===

Dear MITgcm-supporters!

In order to simulate internal waves in a lake (salinity S=0), I used
"eosType='UNESCO') for the equation of state. And I used the diagnostics
pkg to output the density, however I met some following problems:

1. I really don't know the ouputed results for density are based on
which formula, because when I used the formula from Gill's book based on
UNESCO EOS to convert Temp. (T) output from the MITgcm model to the
density (rho) for the pure water

(rho=999.842594+6.793952*1.E-2*t-9.095290*1.E-3*t^2+1.001685*1.E-4*t^3
    -1.120083*1.E-6*t^4+6.536332*1.E-9*t^5)

Then I got a very diffenrent result in comparison with the density
output from the MITgcm model. However, at the beginning to run the
model, I used this formula to check with density output from the model
at t=0 (RHOAnoma.0000000000.*), then I got the same result.

I attach herewith 2 graphics of rho obtained from directly the model
(rhoa24000.eps) and another obtained from converting of Temp to density
by above formula (24000_in.eps).

If anyone could know please help me to explain the reason why there
are different when I used the same UNESCO EOS?

Thanks so much!

Van Thinh


On Tue, 19 Jul 2005, Dimitris Menemenlis wrote:

> Hi Van Thinh, I do not know exact cause of your error messages.
> If I were in your situation I would try starting from scratch
> with a new, clean copy of MITgcm, modify package.cong, then
> genmake2, make depend, and make.
>
> Dimitris
>

-- 
Edward H. Hill III, PhD
office:  MIT Dept. of EAPS;  Rm 54-1424;  77 Massachusetts Ave.
             Cambridge, MA 02139-4307
emails:  eh3 at mit.edu                ed at eh3.com
URLs:    http://web.mit.edu/eh3/    http://eh3.com/
phone:   617-253-0098
fax:     617-253-4464