add calc_dpexp.py (calculate different reference area)