/[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.24 by jmc, Tue Jun 20 22:55:08 2006 UTC revision 1.25 by heimbach, Wed Feb 7 00:01:15 2007 UTC
# Line 325  CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bi Line 325  CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bi
325  C-    express the Tensor in term of Diffusivity (= m**2 / s )  C-    express the Tensor in term of Diffusivity (= m**2 / s )
326        DO j=1-Oly+1,sNy+Oly-1        DO j=1-Oly+1,sNy+Oly-1
327         DO i=1-Olx+1,sNx+Olx-1         DO i=1-Olx+1,sNx+Olx-1
328    #ifdef ALLOW_KAPGM_CONTROL
329            Kgm_tmp = GM_isopycK + GM_skewflx*kapgm(i,j,k,bi,bj)
330    #else
331          Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K          Kgm_tmp = GM_isopycK + GM_skewflx*GM_background_K
332    #endif
333  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
334       &          + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx)           &          + VisbeckK(i,j,bi,bj)*(1. _d 0 + GM_skewflx)    
335  #endif  #endif
# Line 414  CADJ STORE taperFct(:,:)     = comlev1_b Line 418  CADJ STORE taperFct(:,:)     = comlev1_b
418          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly+1,sNy+Oly-1
419           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx+1,sNx+Olx-1
420            Kuz(i,j,k,bi,bj) =            Kuz(i,j,k,bi,bj) =
421    #ifdef ALLOW_KAPGM_CONTROL
422         &     ( GM_isopycK - GM_skewflx*kapgm(i,j,k,bi,bj)
423    #else
424       &     ( GM_isopycK - GM_skewflx*GM_background_K       &     ( GM_isopycK - GM_skewflx*GM_background_K
425    #endif
426  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
427       &     +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i-1,j,bi,bj))*GM_advect       &     +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i-1,j,bi,bj))*GM_advect
428  #endif  #endif
# Line 533  CADJ STORE taperFct(:,:)     = comlev1_b Line 541  CADJ STORE taperFct(:,:)     = comlev1_b
541          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly+1,sNy+Oly-1
542           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx+1,sNx+Olx-1
543            Kvz(i,j,k,bi,bj) =            Kvz(i,j,k,bi,bj) =
544    #ifdef ALLOW_KAPGM_CONTROL
545         &     ( GM_isopycK - GM_skewflx*kapgm(i,j,k,bi,bj)
546    #else
547       &     ( GM_isopycK - GM_skewflx*GM_background_K       &     ( GM_isopycK - GM_skewflx*GM_background_K
548    #endif
549  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
550       &     +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i,j-1,bi,bj))*GM_advect       &     +op5*(VisbeckK(i,j,bi,bj)+VisbeckK(i,j-1,bi,bj))*GM_advect
551  #endif  #endif
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22