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

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

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

revision 1.10 by jmc, Mon Feb 5 03:22:44 2007 UTC revision 1.11 by jmc, Tue Mar 16 00:08:27 2010 UTC
# Line 155  C--   Zonal & Meridional gradient of pot Line 155  C--   Zonal & Meridional gradient of pot
155  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
156        IF (select_rStar.GE.2 .AND. nonlinFreeSurf.GE.1 ) THEN        IF (select_rStar.GE.2 .AND. nonlinFreeSurf.GE.1 ) THEN
157         IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN         IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN
158  C--    z* coordinate slope term: rho'/rho0 * Grad_r(g.z)  C--    z* coordinate slope term: rho_prime/rho0 * Grad_r(g.z)
159          factorZ = gravity*recip_rhoConst*recip_rhoFacC(k)*0.5 _d 0          factorZ = gravity*recip_rhoConst*recip_rhoFacC(k)*0.5 _d 0
160          DO j=jMin,jMax          DO j=jMin,jMax
161           DO i=iMin,iMax           DO i=iMin,iMax
# Line 180  C--    z* coordinate slope term: rho'/rh Line 180  C--    z* coordinate slope term: rho'/rh
180           ENDDO           ENDDO
181          ENDDO          ENDDO
182         ELSEIF (buoyancyRelation .EQ. 'OCEANICP' ) THEN         ELSEIF (buoyancyRelation .EQ. 'OCEANICP' ) THEN
183  C--    p* coordinate slope term: alpha' * Grad_r( p )  C--    p* coordinate slope term: alpha_prime * Grad_r( p )
184          factorP = 0.5 _d 0          factorP = 0.5 _d 0
185          DO j=jMin,jMax          DO j=jMin,jMax
186           DO i=iMin+1,iMax           DO i=iMin+1,iMax
# Line 199  C--    p* coordinate slope term: alpha' Line 199  C--    p* coordinate slope term: alpha'
199           ENDDO           ENDDO
200          ENDDO          ENDDO
201         ELSEIF ( buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN         ELSEIF ( buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN
202  C--    p* coordinate slope term: alpha' * Grad_r( p )  C--    p* coordinate slope term: alpha_prime * Grad_r( p )
203          conv_theta2T = (rC(k)/atm_Po)**atm_kappa          conv_theta2T = (rC(k)/atm_Po)**atm_kappa
204          factorP = (atm_Rd/rC(k))*conv_theta2T*0.5 _d 0          factorP = (atm_Rd/rC(k))*conv_theta2T*0.5 _d 0
205          DO j=jMin,jMax          DO j=jMin,jMax
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22