BISMG:SarahS/jules model description

JULES-glacier model description
Long wave down radiation is put into elevation tiles (in control.F90) using the equation

$$LW down tile= \sigma T^{4} $$ where T has been adjusted for elevation.

Snowfall is put onto elevation tiles (in control.F90). The wet bulb temperature is calculated using the elevated temperature humidity. If the wet bulb temp is < 0oC then rain is converted into snow. snow = snowfall + rainfall (which is now snow)

The albedo of snow and ice is changed in tile_albedo.F90. For ice the code came from Jeff Ridley. If the tile temp is 1 deg warmer than the melt temp and snowdepth is less than 0.1m and we have an ice type then adjust the albedo from the default values (which are set in the .jin file)