Add get_china_mask function to calculate mask

Currently only province mask is supported.