How to calculate Percentiles of Multidimensional metrix for extreme precipitation value?

Dear all,

I have three dimension data sets (i.e stations,year, day). I want to calculate the 90 and 95th percentile to make the threshold for extreme precipitaion value for individual staion peryear. Also my data sets has NaNs in them. Not sure if that affects the commands. please help, thankyou very much. attached data (rainfall) in format of 302(staions)*13 (year)*365(days).

clc
clear
for s=1:302;%%%% number fo stations
    for  y=1:13; %%%% number fo year
        for d=1:365; %%%% number fo day
            percentile(s,y,d)= prctile(rainfall,[90 95],'all')
        end
    end