/[MITgcm]/MITgcm/model/src/calc_phi_hyd.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_phi_hyd.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.39 by heimbach, Fri Feb 13 21:56:48 2009 UTC revision 1.40 by jmc, Tue Mar 16 00:08:27 2010 UTC
# Line 250  C  --  end if integr_GeoPot = ... Line 250  C  --  end if integr_GeoPot = ...
250  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
251        ELSEIF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN        ELSEIF ( buoyancyRelation .EQ. 'OCEANICP' ) THEN
252  C       This is the hydrostatic pressure calculation for the Ocean  C       This is the hydrostatic pressure calculation for the Ocean
253  C       which uses the FIND_RHO() routine to calculate density  C       which uses the FIND_RHO() routine to calculate density before
254  C       before integrating (1/rho)'*dp over the current layer/interface  C       integrating (1/rho)_prime*dp over the current layer/interface
255  #ifdef      ALLOW_AUTODIFF_TAMC  #ifdef      ALLOW_AUTODIFF_TAMC
256  CADJ GENERAL  CADJ GENERAL
257  #endif      /* ALLOW_AUTODIFF_TAMC */  #endif      /* ALLOW_AUTODIFF_TAMC */
# Line 283  CADJ &     kind = isbyte Line 283  CADJ &     kind = isbyte
283            ENDDO            ENDDO
284          ENDIF          ENDIF
285    
286  C--     Calculate specific volume anomaly : alpha' = 1/rho - alpha_Cst  C--     Calculate specific volume anomaly : alpha_prime = 1/rho - alpha_Cst
287          DO j=jMin,jMax          DO j=jMin,jMax
288            DO i=iMin,iMax            DO i=iMin,iMax
289              locAlpha=alphaRho(i,j)+rhoConst              locAlpha=alphaRho(i,j)+rhoConst
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22