ice_scaling Module Reference

Functions/Subroutines
subroutine	scale_fluxes
subroutine	scale_hist_fluxes

Function/Subroutine Documentation

scale_fluxes()

subroutine ice_scaling::scale_fluxes

(

)

Definition at line 70 of file ice_scaling.f90.

!
! !REVISION HISTORY:
!
! author: C.M.Bitz ?
!
! !USES:
!
      use ice_albedo
!
! !INPUT/OUTPUT PARAMETERS:
!
!EOP
!
      real (kind=dbl_kind) :: ar   ! 1/aice

      integer (kind=int_kind) :: i, j

      do j=jlo,jhi
       do i=ilo,ihi
          if (tmask(i,j) .and. aice(i,j) > c0i) then
             ar = c1i/aice(i,j)
!             strairxT(i,j) = strairxT(i,j) * ar
!             strairyT(i,j) = strairyT(i,j) * ar
!     the dynamic is not considered

             fsens(i,j) = fsens(i,j) * ar
             flat(i,j) = flat(i,j) * ar
             fswabs(i,j) = fswabs(i,j) * ar
             flwout(i,j) = flwout(i,j) * ar
             evap(i,j) = evap(i,j) * ar
             tref(i,j) = tref(i,j) * ar  
             qref(i,j) = qref(i,j) * ar  
             fresh(i,j) = fresh(i,j) * ar
             fsalt(i,j) = fsalt(i,j) * ar
             fhnet(i,j) = fhnet(i,j) * ar
             fswthru(i,j) = fswthru(i,j) * ar
          else  ! zero out fluxes 
             strairxt(i,j) = c0i
             strairyt(i,j) = c0i
!    dynamic part using fvcom code
             fsens(i,j) = c0i
             flat(i,j) = c0i
             fswabs(i,j) = c0i
             flwout(i,j) = -stefan_boltzmann*(tffresh - tf(i,j))**4 
                                ! to make upward longwave over ocean
                                ! reasonable for history file
             evap(i,j) = c0i
             tref(i,j) = tair(i,j)
             qref(i,j) = qa(i,j)
             fresh(i,j) = c0i
             fsalt(i,j) = c0i
             fhnet(i,j) = c0i
             fswthru(i,j) = c0i    
             alvdf(i,j) = c0i ! zero out albedo where ice is absent
             alidf(i,j) = c0i
             alvdr(i,j) = c0i
             alidr(i,j) = c0i
          endif 
       enddo  
      enddo

scale_hist_fluxes()

subroutine ice_scaling::scale_hist_fluxes

(

)

Definition at line 142 of file ice_scaling.f90.

!
! !DESCRIPTION:
!
!  Divide ice fluxes by ice area used by the coupler before writing out
!  diagnostics. aice\_init is the ice area saved from coupling.  This
!  makes the fluxes written to the history file consistent with those
!  sent to the coupler.
!
! !REVISION HISTORY:
!
! author: C.M.Bitz ?
!
! !USES:
!
! !INPUT/OUTPUT PARAMETERS:
!
!EOP
!
      real (kind=dbl_kind) :: ar  ! 1/aice

      integer (kind=int_kind) :: i, j

      do j=jlo,jhi
       do i=ilo,ihi
          if (tmask(i,j) .and. aice_init(i,j) > c0i) then
             ar = c1i/aice_init(i,j)
             fresh_hist(i,j)  = fresh_hist(i,j) * ar
             fsalt_hist(i,j)  = fsalt_hist(i,j) * ar
             fhnet_hist(i,j)  = fhnet_hist(i,j) * ar
             fswthru_hist(i,j)  = fswthru_hist(i,j) * ar
          else
             fresh_hist(i,j)  = c0i
             fsalt_hist(i,j)  = c0i
             fhnet_hist(i,j)  = c0i
             fswthru_hist(i,j)  = c0i
          endif 
       enddo  
      enddo