…as well as lakes and wetlands outside of the landmask

with new code:

<        ! Assume wetland, glacier and/or lake when dataset landmask implies ocean

<           if (pctgla(n) < 1.e-6_r8) then
<               pctwet(n)    = 100._r8 - pctlak(n)
<               pctgla(n)    = 0._r8
<           else
<               pctwet(n)    = 100._r8 - pctgla(n) - pctlak(n)
<           end if

from

>        ! Assume wetland and/or lake when dataset landmask implies ocean

>           pctwet(n)        = 100._r8 - pctlak(n)
>           pctgla(n)        = 0._r8