--- MITgcm/model/src/calc_phi_hyd.F	2009/02/13 21:56:48	1.39
+++ MITgcm/model/src/calc_phi_hyd.F	2010/03/16 00:08:27	1.40
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_phi_hyd.F,v 1.39 2009/02/13 21:56:48 heimbach Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_phi_hyd.F,v 1.40 2010/03/16 00:08:27 jmc Exp $
 C $Name:  $
 
 #include "PACKAGES_CONFIG.h"
@@ -250,8 +250,8 @@
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
       ELSEIF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
 C       This is the hydrostatic pressure calculation for the Ocean
-C       which uses the FIND_RHO() routine to calculate density
-C       before integrating (1/rho)'*dp over the current layer/interface
+C       which uses the FIND_RHO() routine to calculate density before
+C       integrating (1/rho)_prime*dp over the current layer/interface
 #ifdef      ALLOW_AUTODIFF_TAMC
 CADJ GENERAL
 #endif      /* ALLOW_AUTODIFF_TAMC */
@@ -283,7 +283,7 @@
           ENDDO
         ENDIF
 
-C--     Calculate specific volume anomaly : alpha' = 1/rho - alpha_Cst
+C--     Calculate specific volume anomaly : alpha_prime = 1/rho - alpha_Cst
         DO j=jMin,jMax
           DO i=iMin,iMax
             locAlpha=alphaRho(i,j)+rhoConst