Add new common function for get initial condition layers
Add additional function for testing TAVAMP
Update the soil splitting function
1. Could be used for initial condition layers
2. Add handling for fixed two top layer (SLLB[0]=5, SLLB[1]=15)