/[MITgcm]/MITgcm/pkg/ggl90/ggl90_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/ggl90/ggl90_calc.F

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

revision 1.24 by mlosch, Mon Nov 5 08:11:45 2012 UTC revision 1.25 by mlosch, Mon Nov 5 13:11:11 2012 UTC
# Line 360  C     ... across x-faces Line 360  C     ... across x-faces
360  #ifdef ISOTROPIC_COS_SCALING  #ifdef ISOTROPIC_COS_SCALING
361       &      *CosFacU(j,bi,bj)       &      *CosFacU(j,bi,bj)
362  #endif /* ISOTROPIC_COS_SCALING */  #endif /* ISOTROPIC_COS_SCALING */
      &         *maskW(i,j,k,bi,bj) ! paranoia setting  
363           ENDDO           ENDDO
364          ENDDO          ENDDO
365  C     ... across y-faces  C     ... across y-faces
# Line 375  C     ... across y-faces Line 374  C     ... across y-faces
374  #ifdef ISOTROPIC_COS_SCALING  #ifdef ISOTROPIC_COS_SCALING
375       &      *CosFacV(j,bi,bj)       &      *CosFacV(j,bi,bj)
376  #endif /* ISOTROPIC_COS_SCALING */  #endif /* ISOTROPIC_COS_SCALING */
      &         *maskS(i,j,k,bi,bj) ! paranoia setting  
377           ENDDO           ENDDO
378          ENDDO          ENDDO
379  C     Compute divergence of fluxes  C     Compute divergence of fluxes
# Line 388  C     Compute divergence of fluxes Line 386  C     Compute divergence of fluxes
386       &      gTKE(i,j) = -recip_drC(k)*recip_rA(i,j,bi,bj)/hFac       &      gTKE(i,j) = -recip_drC(k)*recip_rA(i,j,bi,bj)/hFac
387       &         *((dfx(i+1,j)-dfx(i,j))       &         *((dfx(i+1,j)-dfx(i,j))
388       &          +(dfy(i,j+1)-dfy(i,j)) )       &          +(dfy(i,j+1)-dfy(i,j)) )
      &         *maskC(i,j,k,bi,bj) ! paranoia setting  
389           ENDDO           ENDDO
390          ENDDO          ENDDO
391  C      end if GGL90diffTKEh .eq. 0.  C      end if GGL90diffTKEh .eq. 0.
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22