icing.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	icing (NOW)

Function/Subroutine Documentation

icing()

subroutine icing

(

type(time), intent(in)

NOW

)

Definition at line 41 of file icing.f90.

  USE mod_utils
  USE mod_force
  USE all_vars
  USE mod_time
  IMPLICIT NONE
  TYPE(TIME), INTENT(IN) :: NOW
  REAL(SP), PARAMETER :: Tfreeze = -1.75_sp 
  REAL(SP), DIMENSION(0:MT) :: ICING_WNDSPD

  IF(dbg_set(dbg_log)) write(ipt,*) "Start Icing Update"

  ! GET THE FORCING DATA
  CALL update_icing(now,icing_satmp,icing_wndx,icing_wndy)
  
  icing_wndspd = sqrt(icing_wndy**2 + icing_wndx**2)
  icing_wndspd(0) = 0.0_sp
  

  IF(dbg_set(dbg_io))THEN
     WRITE(ipt,*) "min/max(SAT)",minval(icing_satmp(1:mt)),maxval(icing_satmp(1:mt))
     WRITE(ipt,*) "min/max(WNDSPD)",minval(icing_wndspd(1:mt)),maxval(icing_wndspd(1:mt))
     WRITE(ipt,*) "min/max(T1(:,1))",minval(t1(1:mt,1)),maxval(t1(1:mt,1))
  END IF
  
  icing_0kts = (tfreeze - icing_satmp) /&
       (1.0_sp + 0.4_sp *(t1(:,1) -tfreeze))
  
  WHERE (icing_0kts < 0.0_sp) icing_0kts = 0.0_sp

  icing_10kts = (icing_wndspd + 5.0_sp) * icing_0kts

  icing_0kts = icing_0kts * icing_wndspd
  

Here is the call graph for this function:

Here is the caller graph for this function: