I have this code, which I read into it some input and ouput, now I want to average indivivually tbo1,tbo2,tbo3 and tbo4.My problem I dont seem to be getting it right,I thing Im getting it wrong in the write statment,I aslo want to write the output to be in colums and with heading tbo1...
Dear users
I have this short code which is suppose to read 2 in put files i.e. bn_input and HISTORY file. I have placed a short version of a history file. My bn_input contains number of configurations starting from 1 to 10 (NCONFIGS) and NFRACT which can be 0 or 1. Now the problem when I read...
Here I have included bits of the code,In the first part I assign a unique number to the atoms, then combine them, then use the combination at a later part of the code. Now I want to assign unique numbers to Ca and Mg and use them as shown below in the code
[! Assign Si = 3**0 = 1, Al = 3**1 =...
Hi microm
I atoms which are not bonded,I put them in a box,then I run an MD simulation,they react and form bonds, I know the likes bonds. I have already observed them using my other code. I also abserved the formation of triads Si-Al-Mg forming in my system. Now there is a competion between Mg...
Hi gsal
Yes its combination of 3 atoms because im interested in bonding angles, bond lengths, and other short range order properties which takes place between 2 to 3 atoms. I raised the issue of chain because microm suggestion gave an idea that I can also look at how long the chains are, which...
just for interest is it possible to work out the combination in advance because of systems with varrying number of atoms one of them is 2377 atoms = N"
Sorry for this!
I meant system size made up of 2377 atoms,with 5 atom types.
Thanks a lot guys
With these Ideas I can also be able to work out how long my Si and Al polymer chains. I like the most mikrom Idea.But one more question just for interest is it possible to work out the combination in advance because of systems with varrying number of atoms one of them is 2377...
Thanks Salgerman
Can you illustrate with an example so that I can see your approach for example demonstrate using this statement.If you are interested in combining any 3 atoms at a time, you can assign them numbers that are powers of 3; with the limitation that every atom can only appear 2...
Hi forum users
I need a set of numbers that can differentiate uniquely between the combinations.I choose powers of 3. Starting with 1, that gives multiples of 1, 2 and 3 before I reach the next number, which is 3, with it having multiples of 3, 6, 9, etc. If I write out the combinations...
Thanks a lot XWB
I this think this is an elegant way of doing it,I assume atname(k),atname(nat) are general statements I thing I should replace them with atname(1) = 'Si' as an example and shell k replace it with shell 1 ,correct me if im wrong. Then that will take care of bonding between atom...
Dear Forum users
I have this piece of code which calculate coordination numbers around a specific atom. The code below works well to calculate the coordination numbers around Si and Al atoms. Now I want to increase number of atoms ,for example at Ca and Mg.then from there I want to work out...
Hi every one
I have this code I have co developed,now I want to add a subroutine which calculated or do binning for different range of calculated angles. The code works fine but I struggling to add a subroutine which will do binning for this code. To make this easier I would like to do do...
The original code read a certain number of time origins
NTIMOR = NSTEP / IOR = 1000 / 4 = 250
and after every time origin you have to read a certain number of time steps
Code:
IOR - 1 = 4 - 1 = 3
In practice the code reads 4 steps at each iteration (the first step being a time origin) and...
Hi Macrom
I have been trying to follow your suggestion but Im getting lost: The loop should be changed as I have done below?
[code]
program read_data
implicit none
integer stepno, arrsize
real steptime, datum, datum_1,datum_2,datum_3,datum_4,datum_5,datum_6, &
datum_7,datum_8,i,k
real...
The recent data is the one I get when I run the code i.e the output from the code above. Then I broke the data into pieces to demonstrate what I want my data to be.
Hi Mikrom
I want to add the label at the top left conner of the data as shown below:
[]
STEP1
10000.100 0.62234730 -3.6641850 2.4883280 -3.6641850 1.8993200 0.29834500 2.4883280 0.29834500 3.0102229
10000.200 -1.7910060 0.15255409...
This site uses cookies to help personalise content, tailor your experience and to keep you logged in if you register.
By continuing to use this site, you are consenting to our use of cookies.