summaryrefslogtreecommitdiff
path: root/fortran/navy_four_levels/FourLevelPulseProp_Double.f95
Commit message (Collapse)AuthorAgeFilesLines
* zpts -> 200 to match xmds defaultsEugeniy Mikhailov2011-08-261-1/+1
| | | | waiting time is manageable about 477 seconds
* Prints progress percentEugeniy Mikhailov2011-08-261-0/+1
| | | | To keep bored operator from killing the long calculation task
* Increased maximum number of nodes along zEugeniy Mikhailov2011-08-261-1/+1
|
* Looks like this numbers match xmds codeEugeniy Mikhailov2011-08-261-6/+12
|
* Small clarifying commentsEugeniy Mikhailov2011-08-261-1/+2
|
* Added pulse descriptionEugeniy Mikhailov2011-08-261-5/+8
|
* delc --> del3. The coupling field is just a field 3 nowEugeniy Mikhailov2011-08-261-6/+6
|
* Omc -> Om3 to match xmds code, Rabi frequencies match xmdsEugeniy Mikhailov2011-08-261-24/+35
| | | | Whatch for factor of 2 in Rabi. NEED TO DOUBLE CHECK!
* M_PI -> piEugeniy Mikhailov2011-08-261-1/+1
|
* Added total decay rates of upper levels 3,4Eugeniy Mikhailov2011-08-261-7/+15
| | | | | Decay rates calculated via branching ratios according to xmds code. Decays match xmds code with in factor 2. NEED TO CHECK IT !
* Names g1, ..., g8 --> gc1, ... ,gc8 to avoid confusion with decay ratesEugeniy Mikhailov2011-08-261-11/+11
|
* Comment about W variables as decay ratesEugeniy Mikhailov2011-08-251-0/+1
|
* Remove unused comment and common blockEugeniy Mikhailov2011-08-251-4/+0
|
* Moved forward the definition for the peak Rabi frequenciesEugeniy Mikhailov2011-08-231-6/+20
|
* 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
|
* 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
|
* Small format fixEugeniy Mikhailov2011-08-081-1/+1
|
* Get rid of unused variables, labels, format stringsEugeniy Mikhailov2011-08-081-11/+3
|
* 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-121-0/+708
and edited makefile accordingly. Added t variable to the output file Movie4level.dat for plotting purposes.