add montecarlo and combine ddw and sdw in twogap_ddw.f90