ice_ocean Module Reference

Functions/Subroutines
subroutine	mixed_layer

Variables
logical(kind=log_kind)	oceanmixed_ice

Function/Subroutine Documentation

mixed_layer()

subroutine ice_ocean::mixed_layer

(

)

Definition at line 83 of file ice_ocean.f90.

!
! !USES:
! 
      use ice_flux
      use ice_grid, only: tmask
      use ice_atmo
      use ice_state
      use ice_albedo
!
! !INPUT/OUTPUT PARAMETERS:
!
!EOP
!
      integer (kind=int_kind) :: &
         i, j            ! horizontal indices

      real (kind=dbl_kind), dimension(ilo:ihi,jlo:jhi) ::&
         delt    & ! potential temperature difference   (K)
      ,  delq    & ! specific humidity difference   (kg/kg)
      ,  dummy1, dummy2, dummy3, dummy4  ! dummy arrays

      real (kind=dbl_kind) :: &
         tsfk     &! surface temperature (K)
      ,  fsh      &! sensible heat flux  (W/m^2)
      ,  flh      &! latent heat flux    (W/m^2)
      ,  swabs    &! surface absorbed shortwave heat flux (W/m^2)
      ,  flwup    &! long-wave upward heat flux  (W/m^2)
      ,  ft       &! fraction reduction of positive qdp
      ,  dtcprho  ! dt/cp_ocn*rhow

      dtcprho = dtice/(cp_ocn*rhow)

      call atmo_boundary_layer (1, 'ocn', sst, &
           dummy1, dummy2, dummy3, dummy4, delt,  delq)

      do j = jlo,jhi
      do i = ilo,ihi
       if (tmask(i,j)) then
!       if (hmix(i,j) > puny) then

         ! ocean surface temperature in Kelvin
         tsfk = sst(i,j) + tffresh

         ! shortwave radiative flux
         swabs = (c1i - albocn) * fsw(i,j)

         ! longwave radiative flux
         flwup  = -stefan_boltzmann * tsfk**4

         ! downward latent and sensible heat fluxes
         flh = lhcoef(i,j) * delq(i,j)
         fsh = shcoef(i,j) * delt(i,j)

         ! compute sst change due to exchange with atm/ice above
         ! Note: fhnet, fswthru are added in ice_therm_vertical.F
           sst(i,j) = sst(i,j) + (c1i-aice(i,j))*dtcprho/hmix(i,j) &
                         *(fsh + flh + flwup + flw(i,j) + swabs)

         ! adjust qdp if cooling of mixed layer would occur when sst le Tf
         if( sst(i,j) <= tf(i,j) .and. qdp(i,j) > c0i ) qdp(i,j) = c0i

         ! computed T change due to exchange with deep layers:
           sst(i,j) = sst(i,j) - qdp(i,j)*dtcprho/hmix(i,j)

         ! compute potential to freeze or melt ice
         frzmlt(i,j) = (tf(i,j)-sst(i,j))/dtcprho*hmix(i,j)
         frzmlt(i,j) = min(max(frzmlt(i,j),-c1000),c1000)

         ! if sst is below freezing, reset sst to Tf
         if (sst(i,j) <= tf(i,j)) sst(i,j) = tf(i,j)

       else

         frzmlt(i,j) = c0i
         sst(i,j) = tf(i,j)

!       endif                    ! hmix > puny
       endif                    ! tmask
      enddo                     ! i
      enddo                     ! j

Here is the call graph for this function:

Variable Documentation

oceanmixed_ice

logical (kind=log_kind) ice_ocean::oceanmixed_ice

Definition at line 53 of file ice_ocean.f90.

      logical (kind=log_kind) :: &
        oceanmixed_ice           ! if true, use ocean mixed layer