Changes between Version 14 and Version 15 of fortran


Ignore:
Timestamp:
Oct 18, 2016, 9:51:21 PM (4 years ago)
Author:
ivasileska
Comment:

Move out the comments

Legend:

Unmodified
Added
Removed
Modified
  • fortran

    v14 v15  
    646646
    647647== 13. Runge-Kutta method for harmonic oscillations ==
     648Using method Runge-Kuta obtain the differential equations for oscillator, the same like in program 12 but the difference is that here we are using onother function for du(1) which is velocity. Part 2
    648649{{{
    649650#!fortran
    650 
    651 !Using method Runge-Kuta obtain the differential equations for oscillator, the same like in
    652 !program 12 but the difference is that here we are using onother function for du(1) which
    653 !is velocity
    654 !Part 2
    655 
    656651external fct, out
    657652real aux(8,2)
     
    866861
    867862== 14. Runge-Kutta method for harmonic oscillations ==
     863The same like in program 12 and 12 using Runge Kutta obtain the differential equations but with another example and using !another function for du(1).Part 3
    868864{{{
    869865#!fortran
    870 !The same like in program 12 and 12 using Runge Kutta obtain the differential equations but
    871 !with another example and using !another function for du(1)
    872 !Part 3
     866
    873867external fct, out
    874868real aux(8,2)
     
    10701064}}}
    10711065== 15. Electron acceleration in uniform electric field ==
     1066The method Runge-Kuta is not used only for obtaining the differential equation for oscilloscope. The same method can be used in different tasks where the differential equations need be solve. Using method Runge-Kuta, it is easy to find the electron acceleration in electric filed.
     1067
    10721068{{{
    10731069#!fortran
    1074 !The method Runge-Kuta is not used only for obtaining the differential equation for
    1075 !oscilloscope. The same method can be used in different tasks where the differential
    1076 !equations need be solve. Using method Runge-Kuta, it is easy to find the electron
    1077 !acceleration in electric filed
    1078 
    10791070external fct, out
    10801071real aux(8,4)
     
    13031294
    13041295== 16. Determination the conditions of electron capture in SGA mode and the bunch time ==
     1296Also the method Runge Kuta is used for solving more complex problems in plasma physics, like SGA (synchrotron gyro-magnetic auto resonant) which is a self sustaining ECR plasma in magnetic field where the time is rising up. First the method Runge-Kuta is used for auto resonant condition and after that using the
     1297equations of plasma physics for the SGA mode it are determined the conditions.
    13051298{{{
    13061299#!fortran
    1307 !Also the method Runge Kuta is used for solving more complex problems in plasma physics, like
    1308 !SGA (synchrotron gyro-magnetic auto resonant) which is a self sustaining ECR plasma in
    1309 !magnetic field where the time is rising up.
    1310 !First the method Runge-Kuta is used for auto resonant condition and after that using the
    1311 !equations of plasma physics for the SGA mode it are determined the conditions.
    1312 
    13131300external fct, out
    13141301real aux(8,2)
     
    15461533
    15471534== 17. Determination the dependencies between the electron momentum and time, electron momentum and coordinate==
     1535Using the same condition like in program 16. Using method Runge-Kuta for determination the conditions of the electron motion.
    15481536{{{
    15491537#!fortran
    1550 
    1551 !  Using the same condition like in program 16
    1552 !  Using method Runge-Kuta for determination the conditions of the electron motion.
    15531538external fct, out
    15541539real aux(8,4)
     
    17871772
    17881773== 18. Electron acceleration in a uniform electric field using speed light ==
     1774The same introductions like in programs 16 and 17 but the difference is that here we have dependence of speed light. To  get the relative conditions for the motion. Using metod Runge-Kuta
    17891775{{{
    17901776#!fortran
    1791 !The same introductions like in programs 16 and 17 but the difference is that here we have
    1792 !dependence of speed light. To  get the relative conditions for the motion.
    1793 ! Using metod Runge-Kuta
    1794 
    17951777external fct, out
    17961778real aux(8,4)
     
    20222004
    20232005== 19. Charged particle motion in magnetic mirror trap under ECR ==
     2006One of the most important thing in plasma physics is to get the motion equations of the particles. To get that equations it is used Boris method. This method is divided in for sections: using the momentum of the particles, their rotations, using the electrical field and at the end to determinate the coordinates.
     2007Using this method in this program is analyzed the ECR(electron cyclotron resonance) phenomenon in parabolic approximation of the magnetic field (mirror trap) for different input values.
    20242008{{{
    20252009#!fortran
    2026 !One of the most important thing in plasma physics is to get the motion equations of the
    2027 !particles. To get that equations it is used Boris method. This method is divided in for
    2028 !sections: using the momentum of the particles, their rotations, using the electrical field
    2029 !and at the end to determinate the coordinates.
    2030 !Using this method in this program is analyzed the ECR(electron cyclotron resonance)
    2031 !phenomenon in parabolic approximation of the magnetic
    2032 !field (mirror trap) for different input values
    2033 
    20342010real :: ez = 4.8e-10, c=3.0e10, me=9.1e-28, f=2.4e09, pi=3.14159265
    20352011
     
    21942170
    21952171== 20. Solving nonlinear equation ==
     2172To find a solutions for solving nonlinear equation first need to find the roots in the initial interval. There are a lot of methods for how can these roots be obtained. In this program we are using method Newton-Raphson. This method is more rapidly convergent. To find the roots in the interval we are choosing a number which will serve as an initial approximation.
    21962173{{{
    21972174#!fortran
    2198 !To find a solutions for solving nonlinear equation first need to find the roots in the
    2199 !initial interval. There are a lot of methods for how can these roots be obtained. In this
    2200 !program we are using method Newton-Raphson. This method is more rapidly convergent. To find
    2201 !the roots in the interval we are choosing a number which will serve as an initial
    2202 !approximation.
    22032175    program heat1
    22042176
     
    22692241
    22702242== 21. Poisson equations solver ==
     2243To obtain the one dimensional Poisson equation is used the swap method and after the solving the Poisson equation can be determinate the electric filed of the electron layer. Get the plots of the dependencies between the density and the electric filed with coordinate z
    22712244{{{
    22722245#!fortran
    2273 !To obtain the one dimensional Poisson equation is used the swap method and after the
    2274 !solving the Poisson equation can be determinate the electric filed of the electron layer.
    2275 !Get the plots of the dependencies between the density and the electric filed with   
    2276 !coordinate z
    22772246      program poisn1
    22782247real d,n0
     
    23742343
    23752344== 22. Distribution of Langmuir waves in a uniform plasma  ==
     23451D plasma simulation (demo version) using numerical method created by particle cell method calculate the wavelength and the velocity of wave propagation
    23762346{{{
    23772347#!fortran
    2378 !
    2379 ! 1D plasma simulation (demo version)
    2380 ! using numerical method created by particle cell method
    2381 ! calculate the wavelength and the velocity of wave propogation
    2382 
    23832348 program plasma_1D
    23842349   real n,m0
     
    25532518
    25542519
    2555 == 23. Proton acceleration by the filed of electron leyer ==
     2520== 23. Proton acceleration by the filed of electron layer ==
     2521Using method Runge-Kutta find the max possible acceleration of the proton (maximum gradient of the magnetic filed). Plot the graphics dependencies between the power and time, coordinate difference and time
    25562522{{{
    25572523#!fortran
    2558 
    2559 
    2560 !Using method Runge-Kutta find the max possible acceleration of the proton (maximum
    2561 !gradient of the magnetic filed). Plot the graphics dependencies between the power and
    2562 !time, coordinate difference and time 
    2563 
    2564 
    25652524external fct, out
    25662525common /a/ g0, dt, k, pp2b, gme, gmi, dbdz, pi, rl, d
     
    28292788
    28302789== 24. Electron motion in uniform electric filed ==
     2790To integrate the ordinary differential equations is used method Euler. This method is applicable to a system of first order differential equations. In this program this method will help to find the velocity of the electron in uniform electric filed. 
    28312791{{{
    28322792#!fortran
    2833 !To integrate the ordinary differential equations is used method Euler. This method is
    2834 !applicable to a system of first order differential equations. In this program this method
    2835 !will help to find the velocity of the electron in uniform electric filed. 
    28362793real v,x,dt,t,dv,Fm
    28372794open(2,file='y2.dat', status='unknown')