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

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

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

revision 1.7 by jmc, Tue Mar 16 00:08:27 2010 UTC revision 1.8 by jmc, Fri Jun 10 20:55:51 2011 UTC
# Line 133  C  -- P coordinate : Phi(R_low) is simpl Line 133  C  -- P coordinate : Phi(R_low) is simpl
133    
134  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
135  c      IF ( select_rStar.GE.2 .AND. nonlinFreeSurf.GE.4 ) THEN  c      IF ( select_rStar.GE.2 .AND. nonlinFreeSurf.GE.4 ) THEN
136         IF ( select_rStar.GT.0 .AND. nonlinFreeSurf.GE.4 ) THEN         IF ( select_rStar.GE.1 .AND. nonlinFreeSurf.GE.4 ) THEN
137          DO j=jMin,jMax  C-    Integral of b.dr = rStarFac * Integral of b.dr* :
138           DO i=iMin,iMax          IF ( buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN
139    C-     Consistent with Phi'= Integr[ theta'.dPi ] :
140             DO j=jMin,jMax
141              DO i=iMin,iMax
142             phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)             phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)
143       &                          *rStarFacC(i,j,bi,bj)       &                          * rStarFacC(i,j,bi,bj)**atm_kappa
144              ENDDO
145           ENDDO           ENDDO
146          ENDDO          ELSE
147             DO j=jMin,jMax
148              DO i=iMin,iMax
149               phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)
150         &                          * rStarFacC(i,j,bi,bj)
151              ENDDO
152             ENDDO
153            ENDIF
154         ENDIF         ENDIF
155  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
156    
157          DO j=jMin,jMax         DO j=jMin,jMax
158           DO i=iMin,iMax          DO i=iMin,iMax
159             phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)             phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)
160       &            + Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)       &            + Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)
161       &            + phi0surf(i,j,bi,bj)       &            + phi0surf(i,j,bi,bj)
          ENDDO  
162          ENDDO          ENDDO
163           ENDDO
164    
165  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
166  C  -- end if k=Nr.  C  -- end if k=Nr.
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22