conv_over.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	conv_over

Function/Subroutine Documentation

conv_over()

subroutine conv_over

(

)

Definition at line 45 of file conv_over.f90.

!==============================================================================|
   USE all_vars
   USE mod_utils
   IMPLICIT NONE
   REAL(SP):: AVE_T,AVE_S,AVE_R 
   INTEGER :: I,K,KK,J1,J2,J3
!==============================================================================|

   IF (dbg_set(dbg_sbr)) WRITE(ipt,*) "Start: conv_over"
 

!--APPROXIMATE CONVECTIVE OVERTURNING------------------------------------------!

   DO i=1,m
     DO k=kbm1,2,-1
       DO kk=k-1,1,-1
         IF(rho1(i,k) < rho1(i,kk)) THEN
           ave_t = sum(  t1(i,kk:k))/float(k-kk+1) 
           ave_s = sum(  s1(i,kk:k))/float(k-kk+1) 
           ave_r = sum(rho1(i,kk:k))/float(k-kk+1) 
           t1(i,kk:k)   = ave_t
           s1(i,kk:k)   = ave_s
           rho1(i,kk:k) = ave_r
         END IF
       END DO
     END DO
   END DO
       
!-----RECALCULATE ELEMENT-BASED VALUES OF SALINITY/TEMP/DENSITY----------------!

   DO i=1,n
     j1=nv(i,1) ; j2 = nv(i,2) ; j3 = nv(i,3)
     DO k=1,kbm1
       t(i,k)  = one_third*(  t1(j1,k)+  t1(j2,k)+  t1(j3,k))
       s(i,k)  = one_third*(  s1(j1,k)+  s1(j2,k)+  s1(j3,k))
     END DO
   END DO   
          
   IF (dbg_set(dbg_sbr)) WRITE(ipt,*) "End: conv_over"
       
   RETURN

Here is the call graph for this function:

Here is the caller graph for this function: