REAL,parameter     :: r_aspect_thermals=2.
      REAL,parameter     :: l_mix_thermals=30.
      INTEGER,parameter  :: w2di_thermals=0