Skip to content
Snippets Groups Projects
Commit e657b75d authored by Dave Whipp's avatar Dave Whipp
Browse files

Reverted to old version of PDF calcualtion; cleaned up commented-out code

parent e48adc6e
Branches
Tags
No related merge requests found
Show whitespace changes
Inline Side-by-side
source/module_pdf_functions.f90
+ 7
− 11
View file @ e657b75d
...@@ -49,21 +49,17 @@ ...@@ -49,21 +49,17 @@
enddo enddo
psum=0. psum=0.
agecnt=0 agecnt=0
do i=1,lc do i=1,lc ! Loop over number of sample ages
do j=1,eratesc(i) do k=1,num+1 ! Loop over age range in PDF
do k=1,num+1
p(k)=(1./(alpha*ageu(i)*sqrt(2.*pi)))*exp(-0.5*((n(k)-age(i))/ &! Fill probability array p(k)=(1./(alpha*ageu(i)*sqrt(2.*pi)))*exp(-0.5*((n(k)-age(i))/ &! Fill probability array
(alpha*ageu(i)))**2.) (alpha*ageu(i)))**2.)
psum(k)=psum(k)+p(k) ! Fill sum array to check area under array curve psum(k)=psum(k)+p(k) ! Fill sum array to check area under array curve
enddo enddo
agecnt=agecnt+1 agecnt=agecnt+1
enddo enddo
enddo
sum=0. sum=0.
do i=1,num+1 do i=1,num+1 ! Loop over age range in PDF
!pdf(i)=(psum(i)/real(lc))*dx ! Scale PDF array to normalize area under PDF curve
!pdf(i)=(psum(i)/real(lc)) ! Scale PDF array to normalize area under PDF curve
pdf(i)=(psum(i)/real(agecnt)) ! Scale PDF array to normalize area under PDF curve pdf(i)=(psum(i)/real(agecnt)) ! Scale PDF array to normalize area under PDF curve
sum=sum+pdf(i)*dx ! Calculate area under curve sum=sum+pdf(i)*dx ! Calculate area under curve
enddo enddo
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment