Find peak in the historgram data file - FORTRAN

[gaijin1]

Hello guys,

I have a file with X Y values (separated by space)

I want to scan the data to find the PEAK in the file. (It will probably be a lowest Y point...)

How do I do this is Fortran?

Thank you,

Alexander

 

[minger]

Seems like...

OPEN(11,file='inputfile.dat',form='formatted')
DO n=1,nPts
  READ(11,*) x,y
END DO
peak_val = MINVAL(y)
peak_loc = MINLOC(y)

Can you elaborate a little more on this "peak" that you need to find? Do you need to take derivatives?

 

[quicknote]

Hello Alexander,

In order to find the peak in your histogram data file using FORTRAN, you will need to use a combination of data reading and analysis techniques. Here are some steps you can follow:

1. First, you will need to read in the data from your file using a READ statement in FORTRAN. You can use a DO loop to read in each line of data and store the X and Y values in separate arrays.

2. Once you have all the data stored, you can use a DO loop again to iterate through the Y values and find the lowest point. You can do this by setting a variable to the first Y value and then comparing it to the rest of the Y values in the array. If a Y value is lower than the current lowest value, update the variable to that value.

3. After the loop is complete, the variable will hold the lowest Y value, which should correspond to the peak in your histogram data.

4. You can then use a PRINT statement to output the X and Y values for the peak to the screen or to a new file.

I hope this helps. Let me know if you have any further questions or if you need clarification on any of the steps.

Best,