/[MITgcm]/MITgcm/pkg/kpp/kpp_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_calc.F

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

revision 1.4 by heimbach, Mon Nov 13 16:37:02 2000 UTC revision 1.8 by heimbach, Mon May 14 21:42:07 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "KPP_OPTIONS.h"  #include "KPP_OPTIONS.h"
5    
# Line 102  c     to OLx=1, OLy=1. Line 103  c     to OLx=1, OLy=1.
103  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
104  #include "tamc.h"  #include "tamc.h"
105  #include "tamc_keys.h"  #include "tamc_keys.h"
       INTEGER    isbyte  
       PARAMETER( isbyte = 4 )  
106  #else /* ALLOW_AUTODIFF_TAMC */  #else /* ALLOW_AUTODIFF_TAMC */
107        integer ikey        integer ikey
108  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
# Line 340  c initialize arrays to zero Line 339  c initialize arrays to zero
339       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)       &         vddiff(I,J,1,2) * SurfaceTendencyS(i,j,bi,bj)
340       &       ) *       &       ) *
341       &       delZ(1) / work2(I,J)       &       delZ(1) / work2(I,J)
342          bosol(I,J) = - gravity * vddiff(I,J,1,1) * Qsw(i,j,bi,bj) *          bosol(I,J) = gravity * vddiff(I,J,1,1) * Qsw(i,j,bi,bj) *
343         &       recip_Cp*recip_rhoNil*recip_dRf(1) *
344       &       delZ(1) / work2(I,J)       &       delZ(1) / work2(I,J)
345         END DO         END DO
346        END DO        END DO
# Line 678  C     the bottom of the mixing layer. Line 678  C     the bottom of the mixing layer.
678    
679        RETURN        RETURN
680        END        END
681    
682          subroutine KPP_CALC_DUMMY(
683         I     bi, bj, myTime, myThid )
684    C     /==========================================================\
685    C     | SUBROUTINE KPP_CALC_DUMMY                                |
686    C     | o Compute all KPP fields defined in KPP.h                |
687    C     | o Dummy routine for TAMC
688    C     |==========================================================|
689    C     | This subroutine serves as an interface between MITGCMUV  |
690    C     | code and NCOM 1-D routines in kpp_routines.F             |
691    C     \==========================================================/
692          IMPLICIT NONE
693    
694    #include "SIZE.h"
695    #include "EEPARAMS.h"
696    #include "PARAMS.h"
697    #include "KPP.h"
698    #include "KPP_PARAMS.h"
699    #include "GRID.h"
700    
701    c Routine arguments
702    c     bi, bj - array indices on which to apply calculations
703    c     myTime - Current time in simulation
704    
705          INTEGER bi, bj
706          INTEGER myThid
707          _RL     myTime
708    
709    #ifdef ALLOW_KPP
710    
711    c Local constants
712          integer i, j, k
713    
714          DO j=1-OLy,sNy+OLy
715             DO i=1-OLx,sNx+OLx
716                KPPhbl (i,j,bi,bj) = 1.0
717                KPPfrac(i,j,bi,bj) = 0.0
718                DO k = 1,Nr
719                   KPPghat   (i,j,k,bi,bj) = 0.0
720                   KPPviscAz (i,j,k,bi,bj) = viscAz
721                   KPPdiffKzT(i,j,k,bi,bj) = diffKzT
722                   KPPdiffKzS(i,j,k,bi,bj) = diffKzS
723                ENDDO
724             ENDDO
725          ENDDO
726          
727    #endif
728          RETURN
729          END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22