Calculating Return Levels for Extreme Precipitation

27 views (last 30 days)
hello!
I have a matlab script that calculates the return levels for the return periods 2, 10, 25, 50, 100years using GEV distribution. However, when i plot it, it shows a blank plot. I would be so very grateful is i could get some help with this. I have also attached my data below. Thank you!!
This is the script i used:
returnT = [2;10;25;50;100];
returnP = 1-(1./returnT);
parmhatP=gevfit(prcp);
w1 = gevinv(returnP,parmhatP(1),parmhatP(2),parmhatP(3));
returnpV = 1:10:100;
returnTp = 1-(1./returnpV);
figure();
returnLevel_P = gevinv(returnTp,parmhatP(1),parmhatP(2),parmhatP(3));
plot(returnpV,returnLevel_P,'k',returnT,w1,'k*');hold on
ylim([150 500]);
title('Precipitation');
xlabel('Return period [years]');
ylabel('Return level [mm/mo]');
grid on

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!