.

software/ors-amarisoft/k/ssb_nr_arfcn_test.py

@@ -3,31 +3,47 @@
 
 import nrarfcn as nr
 
-def default_ssb_nr_arfcn(dl_nr_arfcn):
+khz = 1e-3  # we work in MHz
+
+def default_ssb_nr_arfcn(dl_nr_arfcn, scs_khz):
     f = nr.get_frequency(nrarfcn=dl_nr_arfcn)
     gscn = nr.get_gscn_by_frequency(f)
     if nr.get_frequency_by_gscn(gscn) > f:
         gscn -= 1
+    scs = scs_khz * khz
+    while 1:
         fg = nr.get_frequency_by_gscn(gscn)
+        if fg % scs < 1e-5*scs: # == 0 with tolerating fp
+            break
+        gscn -= 1
     fg_arfcn = nr.get_nrarfcn(fg)
     return fg_arfcn
 
-def _(dl_nr_arfcn):
-    ssb_nr_arfcn = default_ssb_nr_arfcn(dl_nr_arfcn)
+def _(dl_nr_arfcn, scs_khz):
+    ssb_nr_arfcn = default_ssb_nr_arfcn(dl_nr_arfcn, scs_khz)
     f = nr.get_frequency(dl_nr_arfcn)
     fssb = nr.get_frequency(ssb_nr_arfcn)
-    print('%d  ->  %d\t;  %g MHz\t->  %g Mhz' % (dl_nr_arfcn, ssb_nr_arfcn, f, fssb))
+    gssb = nr.get_gscn_by_frequency(fssb)
+#   print('%d  ->  %d\t;  %g MHz\t->  %g Mhz' % (dl_nr_arfcn, ssb_nr_arfcn, f, fssb))
+
+    scs = scs_khz * khz  # subcarrier spacing
+    fssb_div_scs = fssb / scs
+    fssb_mod_scs = fssb % scs
+    if abs(fssb_mod_scs) < 1e-10*scs:
+       fssb_mod_scs = 0
+    print("%d %d  ->  %d\t;  %.16g MHz\t->  %.16g Mhz\t gscn %d\t/scs %.16g  %%scs %.16g·scs" %
+        (dl_nr_arfcn, scs_khz, ssb_nr_arfcn, f, fssb, gssb, fssb_div_scs, fssb_mod_scs/scs))
 
 print('# dl_nr_arfcn  ->  ssb_nr_arfcn')
-_(633300)
-_(632628)
-_(630336)
-_(627300)
-_(624288)
-_(532000)
-_(530890)
-_(526000)
-_(524650)
-_(523020)
-_(437000)
-_(176300)
+_(633300, 30)
+_(632628, 30)
+_(630336, 30)
+_(627300, 30)
+_(624288, 30)
+_(532000, 30)
+_(530890, 30)
+_(526000, 30)
+_(524650, 30)
+_(523020, 30)
+_(437000, 30)
+_(176300, 30)