anchor.diff:  	end function longwave
   	
   	function cloudfrac(rh,qc,nx)
  +! 		Calculate the cloud fraction based on cloud water content (qc) and relative humidity
  +! 		based on equations from Xu and Randal (1996)
   		implicit none
   		real,dimension(nx)::cloudfrac
   		real,intent(in),dimension(nx)::rh,qc
  
positive.code: 		es=611.2*exp(17.67*(t-273.15)/(t-29.65))
  		! finally return relative humidity
  		relative_humidity= e/es
  	end function relative_humidity
  	
  	function shortwave(cloud_cover, day_frac, solar_elevation,nx)
  		implicit none
  		real, dimension(nx) :: shortwave
  		real, intent(in), dimension(nx) :: day_frac, cloud_cover, solar_elevation
  
negative.code: 	end function shortwave
  	
  	function longwave(T_air, cloud_cover,nx)
  ! 		compute longwave down at the surface based on air temperature and cloud fraction
  ! 		based on Idso and Jackson (1969)
  		implicit none
  		real, dimension(nx) :: longwave
  		real, intent(in), dimension(nx) :: T_air,cloud_cover
  