addressing comments

Author: PrajaktaPMohite

cookbooks/continental_compression_with_imposed_faults/doc/continental_compression_with_imposed_faults.md

@@ -49,7 +49,8 @@ The position, dip angle and direction, thickness (5 km), and composition (plasti
  Strain rate (left) and density (right) distributions with temperature contours after 35 Myr of deformation.
 ```
 
-Deformation preferentially localizes along the two pre-existing shear zones forming a grabben-like geometry ({numref}`fig:strain_rate_and_density_0_myr`), which reflects their optimal orientation and position within the model domain to accommodate the imposed convergence. After 35 Myr of convergence ({numref}`fig:strain_rate_and_density_35_myr`) deformation remains strongly localized along these faults,  with crustal shortening and thickening occurring between them forming approximately symmetric orogenic wedge. Some degree of asymmetry develops on the right side of the wedge, where the fault farthest to the right remains active in the upper crust and connects to the fault bounding the wedge near the brittle-ductile transition zone.
+Deformation preferentially localizes along the two pre-existing shear zones forming a grabben-like geometry (Figure 2), which reflects their optimal orientation and position within the model domain to accommodate the imposed convergence. While additional smaller-structures develop after 35 Myr of convergence (Figure 3), deformation remains strongly localized along these faults with crustal shortening and thickening occurring between them.
+
 Given the nonlinearity of the rheology and governing equations, minor variations in fault strength, geometry, lithospheric structure, and boundary velocities may lead to significant variations in the spatiotemporal evolution of deformation. This cookbook provides a flexible framework for exploring the effects of these parameters, and application to hypothesis-driven questions such as fault reactivation, inversion of rift basins, or the partitioning of strain in complex orogens.
 
-These results demonstrate the potential key role of pre-existing faults in guiding the evolution of lithospheric deformation. However, given the nonlinearity of the rheology and governing equations, minor variations in fault strength, geometry, lithospheric structure, and boundary velocities may lead to significant variations in the spatiotemporal evolution of deformation. Similarly, changing the numerical resolution is likely to also affect the results, as the brittle shear band width by introduced by the plastic damper (1e21 Pa s) is still likely not fully resolved at the maximum resolution of 1e21 Pa s. Furthermore, varying degrees of minor variations in the model results will likely occur if a stricter nonlinear solver tolerance is selected. This cookbook provides a flexible framework for exploring the effects of these parameters, and application to hypothesis-driven questions such as fault reactivation, inversion of rift basins, or the partitioning of strain in complex orogens.
+These results demonstrate the potential key role of pre-existing faults in guiding the evolution of lithospheric deformation. However, given the nonlinearity of the rheology and governing equations, minor variations in fault strength, geometry, lithospheric structure, and boundary velocities may lead to significant variations in the spatiotemporal evolution of deformation. Similarly, changing the numerical resolution is likely to also affect the results, as the brittle shear band width by introduced by the plastic damper (1e21 Pa s) is still likely not fully resolved at the maximum resolution of 1e21 Pa s. Furthermore, varying degrees of minor variations in the model results will likely occur if a stricter nonlinear solver tolerance is selected. This cookbook provides a flexible framework for exploring the effects of these parameters, and application to hypothesis-driven questions such as fault reactivation, inversion of rift basins, or the partitioning of strain in complex orogens.