real,parameter     :: r_aspect_thermals=2.
      real,parameter     :: l_mix_thermals=30.
      integer,parameter  :: w2di_thermals=0