Fortran code for sorting a data array


This program can be used to sort a set of data array either in ascending or descending order.

program sort_data
implicit none
integer, parameter :: array_length = 100
integer :: i, count
double precision, dimension(1:array_length) :: data
double precision :: temp
logical :: not_sorted
open(unit=1, file='data.input', action='read', status='old')
read(1,*) (data(i), i=1, array_length)
not_sorted = .TRUE.
do while (not_sorted)
     count = 0
     do i = 1, state - 1
          temp = data(i)
! if one needs to get ascending order need to change here .LT. to .GT.
         if (data(i) .LT. data(i + 1)) then
            data(i) = data(i+1)
            data(i+1) = temp
            count = count
         else
            count = count + 1
         end if
     end do
    if (count == (state - 1)) not_sorted = .FALSE.
end do
end program sort_data

Comments

Post a Comment

Popular posts from this blog

VMD High Resolution Image and Movie

Image
A nice wiki can be found here . The followings are my take home notes from various blogs and wikies. High resolution image Here are some steps which I usually follow to make high quality image from VMD. Set the visualization screen (vmdscene) according to requirement Render the vmdscene  using tachyon mode. File -> Render -> Tachyon Change the default Render command to have something like "/tachyon_path/tachyon_exe" -aasamples 12 %s -format TGA -res 1024 1024 -o %s.tga (The resolution can be tuned by changing the value for -res option (which is 1024 1024 in this case). Enable Ambient occlusion lighting in the Display -> Display Settings window by selecting "on" in "Amb. Occl" and "Shadows" choosers. If AO effect is being used, the choice of  material also affect the final visual effects. Use  "Diffuse", "AOshiny", "AOChalky", or "AOEdgy" materials. If th
Read more

Inset Plot in GRACE (Xmgrace)

Image
Grace is an opensource data-analysis/plotting software. Grace runs on practically any version of Unix-like OS. As well, it has been successfully ported to Win9*/NT/2000/XP (some minor functionality may be missing, though). The inset plot can be generated by following way: Give command in terminal: xmgrace dataset1.dat , (or import data by Data -> Import -> ASCII in other platforms) will create the plot of dataset1.dat in the main window (call it Graph1). Now for creating inset of another dataset let say dataset2.dat in the same graph please follow the following steps. Data -> Import -> ASCII to open the " Grace:Read Sets " dialogue box.  Graph 1 created in step 1 , is labelled G0 . In the " Read to Graph " section the " (+) G0 " graph is selected by default.  Right-click on the " (+) G0 " to open the pop-up menu (see the image).  Select " Create new ". It will create ' (+) G1 ' empty graph. 
Read more