/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_som_adv_y.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_som_adv_y.F

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

revision 1.2 by jmc, Wed Aug 22 00:26:01 2007 UTC revision 1.3 by jmc, Tue Jan 8 19:57:34 2008 UTC
# Line 78  C  vLoc         :: volume transported (p Line 78  C  vLoc         :: volume transported (p
78        PARAMETER( two   = 2. _d 0 )        PARAMETER( two   = 2. _d 0 )
79        PARAMETER( three = 3. _d 0 )        PARAMETER( three = 3. _d 0 )
80        INTEGER i,j        INTEGER i,j
81          _RL  recip_dT
82        _RL  slpmax, s1max, s1new, s2new        _RL  slpmax, s1max, s1new, s2new
83        _RL  vLoc, alf1, alf1q, alpmn        _RL  vLoc, alf1, alf1q, alpmn
84        _RL  alfp, alpq, alp1, locTp        _RL  alfp, alpq, alp1, locTp
# Line 108  C  vLoc         :: volume transported (p Line 109  C  vLoc         :: volume transported (p
109        _RL  fn_yz(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  fn_yz(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
110  CEOP  CEOP
111    
112          recip_dT = 0.
113          IF ( deltaTloc.GT.0. _d 0 ) recip_dT = 1.0 _d 0 / deltaTloc
114    
115        IF ( limiter.EQ.1 ) THEN        IF ( limiter.EQ.1 ) THEN
116         DO j=1-OLy,sNy+OLy         DO j=1-OLy,sNy+OLy
117          DO i=1-OLx,sNx+OLx          DO i=1-OLx,sNx+OLx
# Line 175  C       use same indexing as velocity, " Line 179  C       use same indexing as velocity, "
179          fn_zz(i,j) = aln(i,j)*sm_zz(i, j )          fn_zz(i,j) = aln(i,j)*sm_zz(i, j )
180          fn_xz(i,j) = aln(i,j)*sm_xz(i, j )          fn_xz(i,j) = aln(i,j)*sm_xz(i, j )
181  C--    Save zero-order flux:  C--    Save zero-order flux:
182          vT(i,j)    = fp_o(i,j) - fn_o(i,j)          vT(i,j) = ( fp_o(i,j) - fn_o(i,j) )*recip_dT
183         ENDDO         ENDDO
184        ENDDO        ENDDO
185    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22