[MITgcm-support] linear EOS

[Jean-Michel Campin]

Hi David,

The code is right (it's why I did not answer earlier), 
the comments were (at least) confusing. I've just change
them.

Jean-Michel

On Thu, Jan 19, 2006 at 06:28:42PM -0500, dwang at ldeo.columbia.edu wrote:
> Hello all,
> 
> This seems a little bit confusing to me. In model/src/find_rho.F the code reads
> as below. Although it states that the pertubation is calculated with respect to
> the surface level, refTemp/refSalt still get their numbers from the profiles of
> tRef/sRef which are specified in data, in stead of single numbers at the
> surface level. Did I get anything wrong here? The code is the latest
> checkpoint58.
> 
> Thanks,
>  David
> -------------------------------------------------------------------------
>       IF (equationOfState.EQ.'LINEAR') THEN
> 
> C ***NOTE***
> C In the linear EOS, to make the static stability calculation meaningful
> C we alway calculate the perturbation with respect to the surface level.
> C **********
>        refTemp=tRef(kRef)
>        refSalt=sRef(kRef)
> 
>        dRho = rhoNil-rhoConst
> 
>        DO j=jMin,jMax
>         DO i=iMin,iMax
>          rholoc(i,j)=rhoNil*(
>      &     sBeta*(sFld(i,j,k,bi,bj)-refSalt)
>      &   -tAlpha*(tFld(i,j,k,bi,bj)-refTemp) )
>      &        + dRho
>         ENDDO
>        ENDDO
> --------------------------------------------------------------------------
> 
> ----------------------------------------------------------------
> This message was sent using IMP, the Internet Messaging Program.
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support