OVERVIEW

                 ICE SURFACE ALBEDO MODULE
     Routines for computing the surface albedo over ice.

OTHER MODULES USED

     utilities_mod

PUBLIC INTERFACE

use ice_albedo_mod [, only: ice_albedo_init, ice_albedo ]

PUBLIC ROUTINES

call ice_albedo_init ( t_freeze )

Input

   t_freeze   Temperature in degrees Kelvin at which sea water freezes.
                [real]

------------------------------------------------------------

call ice_albedo ( ice, thickness, temp, albedo, i1, j1 )

INPUT

   ice         ice mask (all partitions should be the same)
                 [logical, dimension(:,:,:)]

   thickness   thickness of ice in meters
                 [real, dimension(:,:,:)]

   temp        ice surface temperature in deg Kelvin
                 [real, dimension(:,:,:)]

INPUT/OUTPUT

   albedo      albedo of ice (including the no ice partition)
                 [real, dimension(:,:,:)]

OPTIONAL INPUT

   i1, j1


NOTES

  The first two dimensions are for longitude and latitude, respectively.
  The third dimension is for the number of partitions.
  The first partition is assumed to be open ocean.

NAMELIST

&ice_albedo_nml

real :: crit_thickness       = 1.00
real :: t_range              = 10.0

real :: min_ice_alb          = 0.50
real :: max_ice_alb          = 0.80

real :: const_alb            = 0.65

CHANGE HISTORY

ERROR MESSAGES

Fatal Error in ice_albedo

    initialization not called

KNOWN BUGS

     None.

NOTES

     None.

FUTURE PLANS

     None.