ubstracting crustal bias

webgeodyn/obsdata/obsdata.py

@@ -154,6 +154,8 @@ class ObsData:
                 dec_times.append(year + month / 12)
 
             Bs, SVs = cdf_file.varget('B_CF'), cdf_file.varget('B_SV')
+            if GO:
+                Bs -= cdf_file.varget('bias_crust')
 
             for i, (B, SV, r, th, ph, time) in enumerate(zip(Bs, SVs, radii, thetas, phis, dec_times)):
                 if self.obsGroups[obsGroupName].display_r is None and not GO: