[MITgcm-support] query ?

[Martin Losch]

Hi,
if you really want to use the POLY3 equation of state, then there are 
fortran programs in
${ROOT}/utils/knudsen2
(knudsen2.f, and unesco.f). both compute the coefficents of a third 
order polynomial fit to the "real" equation of state (which is the 
knudsen formula in one case and the unesco formula in the other). 
You'll have to edit the files for your configuration (change #levels, 
deltaZ etc.)

However, I recommend, that you use a different equation of state, e.g. 
JMD95Z, JMD95P, MDJWF, all of which are "full" nonlinear equations of 
state. the first two are based on the unesco formula with adjusted 
coefficients to compute the density from salinity and potential 
temperature (ref: Jackett+McDougall 1995, JTECH or JPO, i forget 
which). the third one is a different polynomial according to McDougall 
etal (2003, JTECH). JMDP95P and MDJWF use pressure levels to compute 
density, JMD95Z simply uses z-levels (p=-rho0*g*z).

Martin

On Mar 2, 2005, at 11:54 AM, Amitabh Mitra wrote:

>   
>
>  Dear Sir,
>
>  In the verification experiment
>
>  global_ocean.90x40x15 how do I create
>
>  the input file POLY3.COEFFS when the
>
>  no. of vertical levels is changed from
>
>  15 to 17(Nr=17).
>
>  I downloaded the required data and applied
>
>  in the experiment. After neccesssary modifications
>
>  in the file SIZE.h and input/data as I run the
>
>  model the following error is displayed.
>  ---------------------------------------------------------------
>  SGI_2100 33% mpirun -np 1 ./mitgcmuv
>  ini_eos: attempt to read POLY3.COEFFS failed
>            because bad # of levels in data
>  STOP Bad data in POLY3.COEFFS
>  STOP
>  MPI: MPI_COMM_WORLD rank 0 has terminated without calling 
> MPI_Finaliz()
>  SGI_2100 34%
>  ---------------------------------------------------------------
>  Can someone help me solve my problem ?
>
>  Thanking You,
>
>  Amitabh Mitra
>
>
>  <inbox.gif>_______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-support