--- MITgcm/model/src/calc_phi_hyd.F 1998/08/22 17:51:07 1.2 +++ MITgcm/model/src/calc_phi_hyd.F 1998/09/05 17:52:13 1.3 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_phi_hyd.F,v 1.2 1998/08/22 17:51:07 cnh Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_phi_hyd.F,v 1.3 1998/09/05 17:52:13 cnh Exp $ #include "CPP_EEOPTIONS.h" @@ -38,9 +38,10 @@ C (This is now the actual hydrostatic pressure|height at the T/S points) DO j=jMin,jMax DO i=iMin,iMax - phiHyd(i,j,K)=phiHyd(i,j,Km1)+rhoConst*halfLayer + phiHyd(i,j,K)=phiHyd(i,j,Km1)-rhoConst*halfLayer & *0.5 _d 0*(drF(Km1)+drF(K)) & *0.5 _d 0*( buoyKM1(i,j)+buoyKP1(i,j) ) + & *rkFac ENDDO ENDDO