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

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

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

revision 1.140 by jmc, Wed Jul 30 03:31:35 2014 UTC revision 1.141 by jmc, Wed Aug 6 23:17:11 2014 UTC
# Line 127  C     i,j,k         :: loop indices Line 127  C     i,j,k         :: loop indices
127        LOGICAL  DIAGNOSTICS_IS_ON        LOGICAL  DIAGNOSTICS_IS_ON
128        EXTERNAL DIAGNOSTICS_IS_ON        EXTERNAL DIAGNOSTICS_IS_ON
129  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
130    #ifdef ALLOW_AUTODIFF
131          _RL thetaRef
132    #endif /* ALLOW_AUTODIFF */
133  CEOP  CEOP
134    
135  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
# Line 674  C     dry grid boxes of rhoInSitu.  See Line 676  C     dry grid boxes of rhoInSitu.  See
676          ELSE          ELSE
677  C-        fluid is not water:  C-        fluid is not water:
678            DO k=1,Nr            DO k=1,Nr
679               IF ( select_rStar.GE.1 .OR. selectSigmaCoord.GE.1 ) THEN
680    C-    isothermal (theta=const) reference state
681                 thetaRef = thetaConst
682               ELSE
683    C-    horizontally uniform (tRef) reference state
684                 thetaRef = tRef(k)
685               ENDIF
686             DO j=1-OLy,sNy+OLy             DO j=1-OLy,sNy+OLy
687              DO i=1-OLx,sNx+OLx              DO i=1-OLx,sNx+OLx
688                rhoInSitu(i,j,k,bi,bj) = 0.               rhoInSitu(i,j,k,bi,bj) =
689         &         ( theta(i,j,k,bi,bj)
690         &              *( salt(i,j,k,bi,bj)*atm_Rq + oneRL )
691         &         - thetaRef )*maskC(i,j,k,bi,bj)
692              ENDDO              ENDDO
693             ENDDO             ENDDO
694            ENDDO            ENDDO
Legend:
Removed from v.1.140  
changed lines
  Added in v.1.141
  ViewVC Help
Powered by ViewVC 1.1.22