[MITgcm-support] linear EOS

[dwang at ldeo.columbia.edu]

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.