include file: SWDIFFAC 3
      REAL DIFFAC                      ! Diffusivity factor                SWDIFFAC.2      
      PARAMETER ( DIFFAC = 1.66 )                                          SWDIFFAC.3