Hello to everyone.

I am a new user to R-CRISIS and even though I've managed to perform a traditional seismic hazard project (compute PGA, PGV, and PGD values via a logic tree approach and mapped them in form of contours), I have had a hard time trying to figure out disaggregation results. I know there is a manual and a different discussion here in the forum but I cannot perform correctly the disaggregation of seismic hazard even though I carefully read both the manual and the other post. 

First of all, I have input a GMPE model that estimates spectrum acceleration (Sa) in order to create Uniform Hazard Spectrums (UHS) for a few towns of my study area. As it was expected, the results were much greater than those of design spectrums (for instance Eurocode 8) which is not an issue taking into consideration that the probability of exceedance is the same for all the spectral ordinates.

So this is were seismic hazard disaggregation needs to be utilized, when I click on the menu of disaggregation (not batch disaggregation), there is a x-y-z diagram, where x is distance, y is magnitude and z is either intensity or probability of exceedance depending on the user's input in the above menu. Now here is the real problem:

1)Which seismic hazard results should I disaggregate? I mean if i go to the disaggregation menu and the GMPE model is the one which estimates Sa, I will need to select a period (spectral ordinate), which is tricky, will I select the ordinate whose Sa value is the greatest of them all? Or should I disaggregate seismic hazard results from another project for the same area which computes PGA values (which means that the GMPE model won't have spectral ordinates), and then see the couple M, R that contribute most to the hazard and assume that it will be the same for the GMPE model that estimates Sa? To justify myself, the range of M and R is the same for both GMPEs because they were proposed by the same scientists. 

2)Which value should I insert in the intensity section?

3)Is there any way to calculate the epsilon value and why there are two epsilon values in the menu? May I leave them both 1?

4)I've seen in papers that when the disaggregation is performed, the results are shown via a x-y-z diagram, where x is distance, y is magnitude and z is the % contribution to hazard. Apart from the x-axis and y-axis, I do not know how to create the z-axis. Is it something about the diagram's title that I am missing?

5)Last but not least, should I keep all the area seismic sources activated when I am performing a disaggregation of seismic hazard? Or I ought to keep only the source where my study-area's town is within it?

I want to thank everyone, first of all for the time you spend to read this post, and second for whichever answer you provide me with. It is all appreciated!

Also I am inserting an image which is a screenshot of my disaggregation menu in case you manage to find something that is inserted wrong. In that image all sources are active.

Seismic hazard disaggregation when the GMPE that estimates Sa is inserted.