1.13
* intialize some parameters like thb in the main routine: shape
* change stop into warning when CIC is invalid
* a timer is added 
* replace Inf in CIC with a real number 

1.14
* replace graphics:::plot with graphics::plot in the function plotshape