summaryrefslogtreecommitdiff
path: root/fortran
Commit message (Collapse)AuthorAgeFilesLines
* Added note about time execution estimateEugeniy Mikhailov2011-08-231-1/+17
|
* Output of the pulse before and after the cellEugeniy Mikhailov2011-08-231-3/+14
|
* Added output at z=0 and at the very last non skipped time pointEugeniy Mikhailov2011-08-231-5/+13
|
* Removed unused variable tptsEugeniy Mikhailov2011-08-231-1/+1
|
* Small description about Rabi frequencies assumptions addedEugeniy Mikhailov2011-08-231-1/+2
|
* Physics constants to Fortran constantsEugeniy Mikhailov2011-08-231-5/+5
|
* Time to uS at plotsEugeniy Mikhailov2011-08-231-0/+4
|
* No more then given number of step along time is outputedEugeniy Mikhailov2011-08-231-3/+9
|
* Format specifier to exponential notationEugeniy Mikhailov2011-08-231-2/+2
|
* More readable close to xmds parameter namesEugeniy Mikhailov2011-08-231-22/+25
|
* Added some new max points parametersEugeniy Mikhailov2011-08-231-13/+15
|
* Fortran code plotting addedEugeniy Mikhailov2011-08-231-0/+17
|
* Small format fixEugeniy Mikhailov2011-08-081-1/+1
|
* Remove conversion warning: too noisyEugeniy Mikhailov2011-08-081-2/+2
|
* Get rid of unused variables, labels, format stringsEugeniy Mikhailov2011-08-081-11/+3
|
* Add some helpful warning compiler optionsEugeniy Mikhailov2011-08-081-2/+12
|
* Add tmin and move pulse center to t=0Eugeniy Mikhailov2011-07-161-5/+6
|
* Make dt = dzEugeniy Mikhailov2011-07-161-5/+5
| | | | | | | * as Simon discovered dt must equal to dz (in normilized units) then some term cancels in derivatives calculations * as result tpts is not a free parameter anymore but locked to length of the cell and number of zpts
* t and z limits are given independentlyEugeniy Mikhailov2011-07-141-0/+2
| | | | | * something fishy with navy code, I decrease length by factor of 2 and delays for vacuum field is exactly the same
* Canonic EIT - control field (3) is offEugeniy Mikhailov2011-07-131-1/+1
|
* Adjust grid sizeEugeniy Mikhailov2011-07-131-1/+1
| | | | | Close to optimal grid size for this problem. Increasing the number of points by 10 in both directions almost does not change output
* Output of vacuum pulse dataEugeniy Mikhailov2011-07-131-3/+2
|
* Grid set in advance independent of pulse time and widthEugeniy Mikhailov2011-07-131-7/+12
|
* Code proper indent and beatificationEugeniy Mikhailov2011-07-131-63/+46
|
* Remove more of unused codeEugeniy Mikhailov2011-07-131-3/+0
|
* One more loop finished with 'end do'Eugeniy Mikhailov2011-07-131-24/+24
|
* Indent file to improve readabilityEugeniy Mikhailov2011-07-131-70/+71
|
* Do loop switched to F90 format: they finished with 'end do'Eugeniy Mikhailov2011-07-131-6/+6
| | | | * this improves readability
* Remove a lot of unused codeEugeniy Mikhailov2011-07-131-68/+8
|
* Unused subroutines removedEugeniy Mikhailov2011-07-121-70/+0
|
* Remove 1st space from the formated outputEugeniy Mikhailov2011-07-121-1/+1
| | | | * I see no need for this space
* Add the output of the Gaussian peak position in timeEugeniy Mikhailov2011-07-121-0/+1
|
* Calculate canonic EIT caseEugeniy Mikhailov2011-07-121-7/+7
| | | | | | | | * field 1 is weak and has Gaussian temporal profile * field 2 is strong * field 3 is off * vacuum propagated field start similar to field 1 * cell length enlarged to see pulse delay
* Output format changed. Output |Om| instead of amplitudeEugeniy Mikhailov2011-07-121-1/+1
| | | | | * this take care of complex part which seems to be disregarded with current format for write statement
* Remove unused subroutinesEugeniy Mikhailov2011-07-121-124/+0
|
* Permissions fix: sources do not need executable flagEugeniy Mikhailov2011-07-121-0/+0
|
* Replaced FourLevelPulseProp_v*_Double.f95 with FourLevelPulseProp_Double.f95 ↵Simon Rochester2011-07-123-716/+8
| | | | | | and edited makefile accordingly. Added t variable to the output file Movie4level.dat for plotting purposes.
* v2 and v3 are now the same except output filenamesEugeniy Mikhailov2011-07-102-25/+29
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | The v2 and v3 were design to be the same except difference of light fields time profiles in initial conditions. But looks like navy guys during edits drift both versions from correct realization. After bug fixes and mergers of insignificant details. I went peace by peace comparison and removed insignificant differences * print time elapsed * print the same other statements to output * same spaces in the statements so they are not shown as diffs * same number of t point * z frames number calculated the same * L is the same * density is the same * remove trailing semicolons to make lines the same * Omcpeak=0.1 in both files Then I look at crucial differences * k1=1-(W32+W31+W34)*dt -> W34 taken in account * j2=W21*dt * j8=W41*dt added to v3 to take in account decay of rho44 to rho11 * W42=1 in v2 it used to be zero * Omc is constant in both cases * rho14 in v3 had product of two rhos, now it the sum as in v2 * rho24 in v3 was calculated differently form v2: no conjg(Omc) was used * the fix is taken from v2 and seems to match our xmds code The only difference remained is output file names.
* Remove references of Ompeak. Ompeak -> Omcpeak to fix navy's typoEugeniy Mikhailov2011-07-101-2/+2
|
* b2 now initialized before first useEugeniy Mikhailov2011-07-102-2/+2
|
* output file rename to reflect main file versionEugeniy Mikhailov2011-07-101-1/+1
|
* replcae PAUSE with WRITE(*,*) to comply with Fortran 95Eugeniy Mikhailov2011-07-102-2/+2
|
* grumpy about fortranEugeniy Mikhailov2011-07-101-0/+4
|
* output file names match executable versionEugeniy Mikhailov2011-07-101-2/+2
|
* both fortran codes compile nowEugeniy Mikhailov2011-07-101-6/+11
|
* navy v3 added and modified to be compalible with fortranEugeniy Mikhailov2011-07-101-0/+704
|
* edit navy version of code to compile with gfortranEugeniy Mikhailov2011-07-082-0/+724