diff --git a/csst_mci_sim/csst_mci_sim.py b/csst_mci_sim/csst_mci_sim.py
index 58a7426d4f7d6e73c689b3043fc9dc0c7a27ed77..78bfd15c51728b084e837a8af48fa55e3b78b876 100644
--- a/csst_mci_sim/csst_mci_sim.py
+++ b/csst_mci_sim/csst_mci_sim.py
@@ -741,8 +741,8 @@ def centroid(data):
     -------
     cx: the centroid column number,  in horizontal direction definet in python image show
     cy: the centroid row number  ,   in vertical direction
-    
-    '''
+    ################# '''
+    ###
     ###
     from scipy import ndimage
     cy, cx = ndimage.center_of_mass(data)
@@ -1036,7 +1036,7 @@ class MCIsimulator():
 
         Creates lensing parameters;
         """
-        
+#############################
         self.image_g = np.zeros(
             (self.information['ysize'], self.information['xsize']), dtype=float)
         self.image_r = np.zeros(
@@ -1062,7 +1062,7 @@ class MCIsimulator():
         self.ghostMax = self.information['ghostratio']
         self.log.info('Maximum in the ghost model %e' % self.ghostMax)
         return
-   # ##############################################
+# ##############################################
 
     def readCosmicRayInformation(self):
         """
@@ -2745,7 +2745,7 @@ class MCIsimulator():
         Because the pixel-to-pixel non-uniformity effect (i.e. multiplicative) flat fielding takes place
         before CTI and other effects, the flat field file must be the same size as the pixels that see
         the sky."""
-        
+
         # generate flatfiledfile, with
         flatsigma = self.information['flatsigma']
 
@@ -3239,7 +3239,7 @@ class MCIsimulator():
         ####
 
         ofd_g = fits.PrimaryHDU()
-        ####   World coordinate system and related parameters  #####
+        # World coordinate system and related parameters  #####
         hdu_g = fits.ImageHDU(data)
 
         hdu_g.header['WCSAXES'] = (np.int16(2),    'number of WCS axes')
@@ -4340,7 +4340,7 @@ class MCIsimulator():
 
         hdu_r.header['BUNIT'] = ('ADU', 'physical unit of array values')
 
-        ###### file information  ###############
+        # ##### file information  ###############
 
         temp = t3.utcnow()
         data_time = temp.strftime("%Y-%m-%dT%H:%M:%S.%f")
@@ -4381,7 +4381,7 @@ class MCIsimulator():
 
         # ####  detector and  Filter information  #####
         hdu_r.header['FILTER'] = (self.filter_r[:6], 'filter band')
-        #####  Detector information ####
+        # ####  Detector information ####
         hdu_r.header['DETSN'] = (
             'E2V-CCD-290-0000000', 'detector serial number')
         hdu_r.header['DETNAME'] = ('blue-opt', 'detector name')
@@ -4891,7 +4891,7 @@ class MCIsimulator():
 
         hdu_i.header['BUNIT'] = ('ADU', 'physical unit of array values')
 
-        ###### file information  ###############
+        # ##### file information  ###############
 
         temp = t3.utcnow()
         data_time = temp.strftime("%Y-%m-%dT%H:%M:%S.%f")