/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_fluxform.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_fluxform.F

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

revision 1.50 by jmc, Wed Feb 12 00:45:56 2014 UTC revision 1.51 by jmc, Fri Apr 4 19:53:30 2014 UTC
# Line 26  C stresses as well as internal viscous s Line 26  C stresses as well as internal viscous s
26  CEOI  CEOI
27    
28  #include "MOM_FLUXFORM_OPTIONS.h"  #include "MOM_FLUXFORM_OPTIONS.h"
29    #ifdef ALLOW_AUTODIFF
30    # include "AUTODIFF_OPTIONS.h"
31    #endif
32  #ifdef ALLOW_MOM_COMMON  #ifdef ALLOW_MOM_COMMON
33  # include "MOM_COMMON_OPTIONS.h"  # include "MOM_COMMON_OPTIONS.h"
34  #endif  #endif
# Line 59  C     == Global variables == Line 62  C     == Global variables ==
62  #ifdef ALLOW_MOM_COMMON  #ifdef ALLOW_MOM_COMMON
63  # include "MOM_VISC.h"  # include "MOM_VISC.h"
64  #endif  #endif
65  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF
66  # include "tamc.h"  # include "tamc.h"
67  # include "tamc_keys.h"  # include "tamc_keys.h"
68  # include "MOM_FLUXFORM.h"  # include "MOM_FLUXFORM.h"
# Line 462  C-- account for 3.D divergence of the fl Line 465  C-- account for 3.D divergence of the fl
465           DO j=jMin,jMax           DO j=jMin,jMax
466            DO i=iMin,iMax            DO i=iMin,iMax
467             gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)             gU(i,j,k,bi,bj) = gU(i,j,k,bi,bj)
468       &     - (rStarExpW(i,j,bi,bj) - 1. _d 0)/deltaTfreesurf       &     - (rStarExpW(i,j,bi,bj) - 1. _d 0)/deltaTFreeSurf
469       &       *uVel(i,j,k,bi,bj)       &       *uVel(i,j,k,bi,bj)
470            ENDDO            ENDDO
471           ENDDO           ENDDO
# Line 689  C-- account for 3.D divergence of the fl Line 692  C-- account for 3.D divergence of the fl
692           DO j=jMin,jMax           DO j=jMin,jMax
693            DO i=iMin,iMax            DO i=iMin,iMax
694             gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)             gV(i,j,k,bi,bj) = gV(i,j,k,bi,bj)
695       &     - (rStarExpS(i,j,bi,bj) - 1. _d 0)/deltaTfreesurf       &     - (rStarExpS(i,j,bi,bj) - 1. _d 0)/deltaTFreeSurf
696       &       *vVel(i,j,k,bi,bj)       &       *vVel(i,j,k,bi,bj)
697            ENDDO            ENDDO
698           ENDDO           ENDDO
Legend:
Removed from v.1.50  
changed lines
  Added in v.1.51
  ViewVC Help
Powered by ViewVC 1.1.22