/[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.2 by dimitri, Fri May 30 22:13:42 2008 UTC revision 1.4 by dimitri, Fri May 30 23:18:08 2008 UTC
# Line 531  CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bi Line 531  CADJ STORE Kwz(:,:,k,bi,bj) = comlev1_bi
531  #else  #else
532           Kgm_tmp = GM_isopycK           Kgm_tmp = GM_isopycK
533  #endif  #endif
534  #ifdef ALLOW_KAPGM_CONTROL  #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL))
535       &           + GM_skewflx*kapgm(i,j,k,bi,bj)       &           + GM_skewflx*kapgm(i,j,k,bi,bj)
536  #else  #else
537       &           + GM_skewflx*GM_background_K       &           + GM_skewflx*GM_background_K
# Line 715  CADJ STORE taperFct(:,:)     = comlev1_b Line 715  CADJ STORE taperFct(:,:)     = comlev1_b
715  #else  #else
716       &     ( GM_isopycK       &     ( GM_isopycK
717  #endif  #endif
718  #ifdef ALLOW_KAPGM_CONTROL  #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL))
719       &     - GM_skewflx*kapgm(i,j,k,bi,bj)       &     - GM_skewflx*kapgm(i,j,k,bi,bj)
720  #else  #else
721       &     - GM_skewflx*GM_background_K       &     - GM_skewflx*GM_background_K
# Line 776  C-- end 2nd  loop on vertical level inde Line 776  C-- end 2nd  loop on vertical level inde
776    
777  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
778  C-- 3rd  k loop : compute Tensor Coeff. at V point  C-- 3rd  k loop : compute Tensor Coeff. at V point
779    
780  #ifdef ALLOW_KPP  #ifdef ALLOW_KPP
781        IF ( useKPP ) THEN        IF ( useKPP ) THEN
782         DO j=2-Oly,sNy+Oly         DO j=2-Oly,sNy+Oly
# Line 927  CADJ STORE taperFct(:,:)     = comlev1_b Line 928  CADJ STORE taperFct(:,:)     = comlev1_b
928  #else  #else
929       &     ( GM_isopycK       &     ( GM_isopycK
930  #endif  #endif
931  #ifdef ALLOW_KAPGM_CONTROL  #if (defined (ALLOW_AUTODIFF) && defined (ALLOW_KAPGM_CONTROL))
932       &     - GM_skewflx*kapgm(i,j,k,bi,bj)       &     - GM_skewflx*kapgm(i,j,k,bi,bj)
933  #else  #else
934       &     - GM_skewflx*GM_background_K       &     - GM_skewflx*GM_background_K
# Line 1009  CBFK Add the submesoscale contribution, Line 1010  CBFK Add the submesoscale contribution,
1010  CBFK Add in the mu vertical structure factor  CBFK Add in the mu vertical structure factor
1011          DO j=1-Oly+1,sNy+Oly-1          DO j=1-Oly+1,sNy+Oly-1
1012           DO i=1-Olx+1,sNx+Olx-1           DO i=1-Olx+1,sNx+Olx-1
 #ifdef ALLOW_KPP  
           hml=KPPhbl(i,j,bi,bj)  
 #else  
1013            hml=hMixLayer(i,j,bi,bj)            hml=hMixLayer(i,j,bi,bj)
 #endif  
1014            IF (hml.gt.0 _d 0) THEN            IF (hml.gt.0 _d 0) THEN
1015             recip_hml=1 _d 0/hml             recip_hml=1 _d 0/hml
1016            ELSE            ELSE
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22