/[MITgcm]/MITgcm_contrib/submesoscale/code/gmredi_calc_tensor.F
ViewVC logotype

Diff of /MITgcm_contrib/submesoscale/code/gmredi_calc_tensor.F

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

revision 1.5 by zhc, Fri Mar 12 18:31:00 2010 UTC revision 1.7 by zhc, Fri Mar 19 19:29:47 2010 UTC
# Line 1085  C Psiy are on faces above V Line 1085  C Psiy are on faces above V
1085            SM_PsiY(i,j)=op5*(SM_Lf(i,j+1)+SM_Lf(i,j))*dBdyAV(i,j)            SM_PsiY(i,j)=op5*(SM_Lf(i,j+1)+SM_Lf(i,j))*dBdyAV(i,j)
1086       &                        *_maskS(i,j,k,bi,bj)       &                        *_maskS(i,j,k,bi,bj)
1087    
1088    c hzhang: clipping here, ref to Baylor paper 3 appendix A MOM
1089            hml=hMixLayer(i,j,bi,bj)
1090            IF (hml .lt. (delR(1)+delR(2)+delR(3)+delR(4))) THEN
1091                    SM_PsiX(i,j)=0 _d 0
1092                    SM_PsiY(i,j)=0 _d 0
1093            ENDIF
1094    
1095            hml=.5 * delR(k)
1096            IF ( abs(SM_PsiX(i,j)) .gt. hml ) THEN
1097                    SM_PsiX(i,j)=SIGN( hml, SM_PsiX(i,j) )
1098            ENDIF
1099            IF ( abs(SM_PsiY(i,j)) .gt. hml ) THEN
1100                    SM_PsiY(i,j)=SIGN( hml, SM_PsiY(i,j) )
1101            ENDIF
1102    c hzhang: clipping done
1103    
1104    
1105  #ifndef GM_BOLUS_ADVEC  #ifndef GM_BOLUS_ADVEC
1106  C Kwx,Kwy are on rVel Points  C Kwx,Kwy are on rVel Points
1107            Kwx(i,j,k,bi,bj)  = Kwx(i,j,k,bi,bj)            Kwx(i,j,k,bi,bj)  = Kwx(i,j,k,bi,bj)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22