/[MITgcm]/MITgcm/pkg/gmredi/gmredi_calc_tensor.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_calc_tensor.F

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

revision 1.4 by adcroft, Fri Feb 2 21:36:29 2001 UTC revision 1.7 by heimbach, Tue Aug 21 15:27:19 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "GMREDI_OPTIONS.h"  #include "GMREDI_OPTIONS.h"
5    
# Line 117  C       Now we convert deltaH to a non-d Line 118  C       Now we convert deltaH to a non-d
118          if (K.eq.2) VisbeckK(i,j,bi,bj)=0.          if (K.eq.2) VisbeckK(i,j,bi,bj)=0.
119  Calt?   if (dSigmaDrReal(i,j).NE.0.) then  Calt?   if (dSigmaDrReal(i,j).NE.0.) then
120  Calt?    N2=(-Gravity*recip_Rhonil)*dSigmaDrReal(i,j)  Calt?    N2=(-Gravity*recip_Rhonil)*dSigmaDrReal(i,j)
121          if (dRdSigmaLtd(i,j).NE.0.) then          if ( dRdSigmaLtd(i,j).NE.0. .AND. Ssq.NE.0. ) then
122           N2=(-Gravity*recip_Rhonil)/dRdSigmaLtd(i,j)           N2=(-Gravity*recip_Rhonil)/dRdSigmaLtd(i,j)
123           SN=sqrt(Ssq*N2)           SN=sqrt(Ssq*N2)
124           VisbeckK(i,j,bi,bj)=VisbeckK(i,j,bi,bj)+deltaH           VisbeckK(i,j,bi,bj)=VisbeckK(i,j,bi,bj)+deltaH
# Line 131  C       Limit range that KapGM can take Line 132  C       Limit range that KapGM can take
132  #endif /* GM_VISBECK_VARIABLE_K */  #endif /* GM_VISBECK_VARIABLE_K */
133    
134    
135  #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE  #ifdef ALLOW_TIMEAVE
136  C--     Time-average  C--     Time-average
137          GM_Kwx_T(i,j,k,bi,bj)=GM_Kwx_T(i,j,k,bi,bj)          GM_Kwx_T(i,j,k,bi,bj)=GM_Kwx_T(i,j,k,bi,bj)
138       &                       +Kwx(i,j,k,bi,bj)*deltaTclock       &                       +Kwx(i,j,k,bi,bj)*deltaTclock
# Line 144  C--     Time-average Line 145  C--     Time-average
145       &  Visbeck_K_T(i,j,bi,bj)=Visbeck_K_T(i,j,bi,bj)       &  Visbeck_K_T(i,j,bi,bj)=Visbeck_K_T(i,j,bi,bj)
146       &                       +VisbeckK(i,j,bi,bj)*deltaTclock       &                       +VisbeckK(i,j,bi,bj)*deltaTclock
147  #endif  #endif
148  #endif /* INCLUDE_DIAGNOSTICS_INTERFACE_CODE */  #endif /* ALLOW_TIMEAVE */
149         ENDDO         ENDDO
150        ENDDO        ENDDO
151    
152  #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE  #ifdef ALLOW_TIMEAVE
153        GM_TimeAve(k,bi,bj)=GM_TimeAve(k,bi,bj)+deltaTclock        GM_TimeAve(k,bi,bj)=GM_TimeAve(k,bi,bj)+deltaTclock
154  #endif  #endif
155    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22