conv_over.f90

Go to the documentation of this file.

!/===========================================================================/
! Copyright (c) 2007, The University of Massachusetts Dartmouth 
! Produced at the School of Marine Science & Technology 
! Marine Ecosystem Dynamics Modeling group
! All rights reserved.
!
! FVCOM has been developed by the joint UMASSD-WHOI research team. For 
! details of authorship and attribution of credit please see the FVCOM
! technical manual or contact the MEDM group.
!
! 
! This file is part of FVCOM. For details, see http://fvcom.smast.umassd.edu 
! The full copyright notice is contained in the file COPYRIGHT located in the 
! root directory of the FVCOM code. This original header must be maintained
! in all distributed versions.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
! AND ANY EXPRESS OR  IMPLIED WARRANTIES, INCLUDING,  BUT NOT  LIMITED TO,
! THE IMPLIED WARRANTIES OF MERCHANTABILITY AND  FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED.  
!
!/---------------------------------------------------------------------------/
! CVS VERSION INFORMATION
! $Id$
! $Name$
! $Revision$
!/===========================================================================/

!==============================================================================|
!  ADJUST THE VERTICAL WATER COLUMN WHEN UNSTABLE                              |
!==============================================================================|

   SUBROUTINE conv_over           

!==============================================================================|
   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
   END SUBROUTINE conv_over       
!==============================================================================|