import skrf as rf
from pylab import *
import numpy as np
from skrf.calibration import (IEEEP370_SE_NZC_2xThru, IEEEP370_MM_NZC_2xThru,
IEEEP370_SE_ZC_2xThru, IEEEP370_MM_ZC_2xThru)
# load single-ended data
se_2xthru = rf.Network('CAL_K1_SN100815_V11_RI.s4p')
# transform to mixed-modes
mm_2xthru = se_2xthru.copy()
mm_2xthru.se2gmm(p = 2)
# extrapolate to DC for time step
mm_2xthru_dc = mm_2xthru.extrapolate_to_dc(kind='cubic')
# de-embedding
z0 = 50
dm_mmnzc = IEEEP370_MM_NZC_2xThru(dummy_2xthru = se_2xthru, z0 = z0, name = 'mmnzc')