/[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.25 by jmc, Sat Feb 8 02:09:20 2003 UTC revision 1.26 by jmc, Sun Feb 9 02:00:50 2003 UTC
# Line 128  C       P(z=eta) = P(atmospheric_loading Line 128  C       P(z=eta) = P(atmospheric_loading
128          IF (k.EQ.1) THEN          IF (k.EQ.1) THEN
129            DO j=jMin,jMax            DO j=jMin,jMax
130              DO i=iMin,iMax              DO i=iMin,iMax
131                phiHyd(i,j,k) = phi0surf(i,j,bi,bj)  c             phiHyd(i,j,k) = phi0surf(i,j,bi,bj)
132  c             phiHyd(i,j,k) = 0.                phiHyd(i,j,k) = 0.
133              ENDDO              ENDDO
134            ENDDO            ENDDO
135          ENDIF          ENDIF
# Line 227  CADJ GENERAL Line 227  CADJ GENERAL
227          IF (k.EQ.1) THEN          IF (k.EQ.1) THEN
228            DO j=jMin,jMax            DO j=jMin,jMax
229              DO i=iMin,iMax              DO i=iMin,iMax
230                phiHyd(i,j,k) = phi0surf(i,j,bi,bj)  c             phiHyd(i,j,k) = phi0surf(i,j,bi,bj)
231  c             phiHyd(i,j,k) = 0.                phiHyd(i,j,k) = 0.
232              ENDDO              ENDDO
233            ENDDO            ENDDO
234          ENDIF          ENDIF
# Line 331  C--------------------------------------- Line 331  C---------------------------------------
331       &                  -((rC(K)/atm_Po)**atm_kappa) )       &                  -((rC(K)/atm_Po)**atm_kappa) )
332            DO j=jMin,jMax            DO j=jMin,jMax
333             DO i=iMin,iMax             DO i=iMin,iMax
334  c            phiHyd(i,j,K)=  c            phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+
335               phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+               phiHyd(i,j,K)=
336       &          ddPIp*maskC(i,j,K,bi,bj)       &          ddPIp*maskC(i,j,K,bi,bj)
337       &               *(tFld(I,J,K,bi,bj)-tRef(K))       &               *(tFld(I,J,K,bi,bj)-tRef(K))
338             ENDDO             ENDDO
# Line 372  C--------- Line 372  C---------
372       &                  -((rC(K)/atm_Po)**atm_kappa) )       &                  -((rC(K)/atm_Po)**atm_kappa) )
373            DO j=jMin,jMax            DO j=jMin,jMax
374             DO i=iMin,iMax             DO i=iMin,iMax
375  c            phiHyd(i,j,K)=  c            phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+
376               phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+               phiHyd(i,j,K)=
377       &          ddPIp*_hFacC(I,J, K ,bi,bj)       &          ddPIp*_hFacC(I,J, K ,bi,bj)
378       &               *(tFld(I,J, K ,bi,bj)-tRef( K ))       &               *(tFld(I,J, K ,bi,bj)-tRef( K ))
379             ENDDO             ENDDO
# Line 411  C--------- Line 411  C---------
411       &                  -((rC(Kp1)/atm_Po)**atm_kappa) )         &                  -((rC(Kp1)/atm_Po)**atm_kappa) )  
412            DO j=jMin,jMax            DO j=jMin,jMax
413             DO i=iMin,iMax             DO i=iMin,iMax
414  c            phiHyd(i,j,K)=  c            phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+
415               phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+               phiHyd(i,j,K)=
416       &        ( ddPIm*max(zero, _hFacC(i,j,K,bi,bj)-half)       &        ( ddPIm*max(zero, _hFacC(i,j,K,bi,bj)-half)
417       &         +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)-half) )       &         +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)-half) )
418       &               *(tFld(i,j, K ,bi,bj)-tRef( K ))       &               *(tFld(i,j, K ,bi,bj)-tRef( K ))
# Line 457  C--------- Line 457  C---------
457       &                  -((rC(Kp1)/atm_Po)**atm_kappa) )         &                  -((rC(Kp1)/atm_Po)**atm_kappa) )  
458            DO j=jMin,jMax            DO j=jMin,jMax
459             DO i=iMin,iMax             DO i=iMin,iMax
460  c            phiHyd(i,j,K)=  c            phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+
461               phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+               phiHyd(i,j,K)=
462       &        ( ddPIm*max(zero,(_hFacC(i,j,K,bi,bj)-one)*ratioRm+half)       &        ( ddPIm*max(zero,(_hFacC(i,j,K,bi,bj)-one)*ratioRm+half)
463       &         +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)*ratioRp     -half) )       &         +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)*ratioRp     -half) )
464       &               *(tFld(i,j, K ,bi,bj)-tRef( K ))       &               *(tFld(i,j, K ,bi,bj)-tRef( K ))
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.26
  ViewVC Help
Powered by ViewVC 1.1.22