--- MITgcm/model/src/calc_phi_hyd.F	2001/08/13 23:27:53	1.15
+++ MITgcm/model/src/calc_phi_hyd.F	2001/09/26 18:09:14	1.16
@@ -1,16 +1,21 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_phi_hyd.F,v 1.15 2001/08/13 23:27:53 heimbach Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_phi_hyd.F,v 1.16 2001/09/26 18:09:14 cnh Exp $
 C $Name:  $
 
 #include "CPP_OPTIONS.h"
 
+CBOP
+C     !ROUTINE: CALC_PHI_HYD
+C     !INTERFACE:
       SUBROUTINE CALC_PHI_HYD(
      I                         bi, bj, iMin, iMax, jMin, jMax, K,
      I                         theta, salt,
      U                         phiHyd,
      I                         myThid)
-C     /==========================================================\
+C     !DESCRIPTION: \bv
+C     *==========================================================*
 C     | SUBROUTINE CALC_PHI_HYD                                  |
 C     | o Integrate the hydrostatic relation to find the Hydros. | 
+C     *==========================================================*
 C     |    Potential (ocean: Pressure/rho ; atmos = geopotential)|
 C     | On entry:                                                |
 C     |   theta,salt    are the current thermodynamics quantities|
@@ -31,7 +36,9 @@
 C     | Atmosphere:                                              |
 C     |   Integr_GeoPot allows to select one integration method  |
 C     |    (see the list below)                                  |
-C     \==========================================================/
+C     *==========================================================*
+C     \ev
+C     !USES:
       IMPLICIT NONE
 C     == Global variables ==
 #include "SIZE.h"
@@ -43,6 +50,7 @@
 #include "tamc_keys.h"
 #endif /* ALLOW_AUTODIFF_TAMC */
 
+C     !INPUT/OUTPUT PARAMETERS:
 C     == Routine arguments ==
       INTEGER bi,bj,iMin,iMax,jMin,jMax,K
       _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
@@ -52,12 +60,14 @@
       
 #ifdef INCLUDE_PHIHYD_CALCULATION_CODE
 
+C     !LOCAL VARIABLES:
 C     == Local variables ==
       INTEGER i,j, Kp1
       _RL zero, one, half
       _RL alphaRho(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RL dRloc,dRlocKp1
       _RL ddPI, ddPIm, ddPIp, ratioRp, ratioRm
+CEOP
 
       zero = 0. _d 0
       one  = 1. _d 0