Bodges

  1. Require surface Exner function in order to couple Exner function - problem determining departure point otherwise
  2. Coupling mcl needs two changes to avoid floating point exceptions, e.g.
    • In LSP_AUTOC 
      ! Avoiding doing a power to a very small number - marc 21/4/16
!   IF (qcl(i) >  0.0.AND.cfliq(i) >  0.0) THEN
   IF (qcl(i) >  1.0e-200.AND.cfliq(i) >  0.0) THEN
    • In COSP_PRECIP_MXRATIO 
      ! Avoding divisions by very small numbers - marc 23/4/15
! ...AND.(flux(i,k) /= 0.0)) THEN
...AND.(abs(flux(i,k)) > 1.0e-100)) THEN