double precision z_dip(0:40) 
      common/z_dip/z_dip