| @@ -485,6 +485,9 c between k and the magnetic field (radians) | |||
|
|
485 | 485 | sum2=sum2+fi(i)*yi(i) |
|
|
486 | 486 | end do |
|
|
487 | 487 | dl=ak**2*dlf*te/densmks |
|
|
488 | c write(*,fmt='("Before YE")') | |
|
|
489 | c write(*,*) imode | |
|
|
490 | c call exit | |
|
|
488 | 491 | |
|
|
489 | 492 | if(imode.eq.1.or.imode.eq.3) then |
|
|
490 | 493 | |
| @@ -493,6 +496,8 c between k and the magnetic field (radians) | |||
|
|
493 | 496 | phi=(omegae/ak)/(sqrt(2.0)*vte) |
|
|
494 | 497 | psie=(ven/ak)/(sqrt(2.0)*vte) |
|
|
495 | 498 | ye=y_electron(thetae,phi,psie,alpha) |
|
|
499 | write(*,fmt='("AFTER YE")') | |
|
|
500 | c call exit | |
|
|
496 | 501 | |
|
|
497 | 502 | else if(imode.eq.2) then |
|
|
498 | 503 | c |
| @@ -503,7 +508,8 c | |||
|
|
503 | 508 | alpha2=abs(pi/2.0-alpha)*180.0/pi |
|
|
504 | 509 | c write(*,*) "ye: ", ye |
|
|
505 | 510 | call collision(densmks, te, freq, alpha2, ye) |
|
|
506 | c write(*,fmt='(" geobfield: time is before earliest model.")') | |
|
|
511 | c write(*,fmt='("AFTER COLLISION")') | |
|
|
512 | c call exit | |
|
|
507 | 513 | ye=ye*omega+cmplx(0.0,1.0) |
|
|
508 | 514 | |
|
|
509 | 515 | end if |
| @@ -512,12 +518,12 c write(*,fmt='(" geobfield: time is before earliest model.")') | |||
|
|
512 | 518 | p=(cabs(ye))**2*real(sum2)+cabs(sum1+cmplx(0.0,dl))**2*real(ye) |
|
|
513 | 519 | p=p/(cabs(yed+sum1))**2 |
|
|
514 | 520 | spect1=p*2.0e0/(omega*pi) |
|
|
515 |
write(*,*) "spect1 |
|
|
|
521 | write(*,*) "spect1:",spect1 | |
|
|
516 | 522 | return |
|
|
517 | 523 | end |
|
|
518 | 524 | |
|
|
519 |
subroutine acf2(wl, tau, te1, ti1, fi1, ven1, vin1, wi1, |
|
|
|
520 | & alpha1, dens1, bfld1, acf) | |
|
|
525 | subroutine acf2(wl, tau, te1, ti1, fi1, ven1, vin1, wi1, | |
|
|
526 | & alpha1, dens1, bfld1, acf, nion1) | |
|
|
521 | 527 | c |
|
|
522 | 528 | c computes autocorrelation function for given plasma parameters |
|
|
523 | 529 | c by integrating real spectrum |
| @@ -531,33 +537,44 c | |||
|
|
531 | 537 | real pi |
|
|
532 | 538 | integer nion1 |
|
|
533 | 539 | integer wi1(nion1) |
|
|
534 | integer i,j,k | |
|
|
540 | integer i,j,k,imode | |
|
|
541 | common /mode/imode | |
|
|
535 | 542 | c |
|
|
543 | c write(*,*) "INITIAL acf:",wl,tau,te1,ti1,fi1,ven1,vin1,wi1,alpha1 | |
|
|
544 | write(*,*) "INITIAL acf:",dens1, bfld1, acf, nion1 | |
|
|
545 | c write(*,fmt='("INIT")') | |
|
|
536 | 546 | pi=4.0*atan(1.0) |
|
|
537 | 547 | c |
|
|
538 | 548 | c copy arguments to common block |
|
|
539 | 549 | c |
|
|
540 | 550 | ak=2.0*pi/wl |
|
|
541 | 551 | imode=2 |
|
|
552 | write(*,*) "imode:",imode | |
|
|
542 | 553 | |
|
|
543 | 554 | nion=nion1 |
|
|
544 | 555 | alpha=alpha1 |
|
|
545 | 556 | te=te1 |
|
|
546 | 557 | ven=ven1 |
|
|
558 | c write(*,fmt='("INIT2")') | |
|
|
547 | 559 | do i=1,nion |
|
|
560 | c write(*,fmt='("INIT2.5")') | |
|
|
548 | 561 | ti(i)=ti1(i) |
|
|
549 | 562 | fi(i)=fi1(i) |
|
|
550 | 563 | vin(i)=vin1(i) |
|
|
551 | 564 | wi(i)=wi1(i) |
|
|
552 | 565 | end do |
|
|
566 | c write(*,fmt='("INIT3")') | |
|
|
553 | 567 | dens=dens1 |
|
|
554 | 568 | bfld=bfld1 |
|
|
555 | 569 | |
|
|
556 | 570 | c write(*,*) wl,alpha1,bfld1,dens |
|
|
557 | 571 | c call exit |
|
|
558 | ||
|
|
572 | c write(*,fmt='("Before Gauss")') | |
|
|
559 | 573 | call gaussq(tau,acf) |
|
|
560 | 574 | |
|
|
575 | write(*,*) "FINAL acf:",acf | |
|
|
576 | ||
|
|
577 | c write(*,fmt='("After Gauss")') | |
|
|
561 | 578 | return |
|
|
562 | 579 | end |
|
|
563 | 580 | |
| @@ -590,7 +607,7 c write(*,*) leniw,lenw | |||
|
|
590 | 607 | nrmom=0 |
|
|
591 | 608 | ksave=0 |
|
|
592 | 609 | momcom=0 |
|
|
593 | ||
|
|
610 | write(*,*) "Before qc25f:",acf,imode | |
|
|
594 | 611 | c much faster, more robust |
|
|
595 | 612 | c write(*,*) "acf_in: ",acf |
|
|
596 | 613 | call qc25f(spect1,a,b,tau,integr,nrmom,maxp1,ksave,acf, |
| @@ -598,7 +615,7 c write(*,*) "acf_in: ",acf | |||
|
|
598 | 615 | c write(*,*) "acf_out: ",acf |
|
|
599 | 616 | c call qawf(spect1,a,tau,integr,epsabs,acf,abserr,neval, |
|
|
600 | 617 | c & ier,limlst,lst,leniw,maxp1,lenw,iwork,work) |
|
|
601 | ||
|
|
602 | ||
|
|
618 | write(*,*) "After qc25f:",acf | |
|
|
619 | c call exit | |
|
|
603 | 620 | return |
|
|
604 | 621 | end |
General Comments 0
You need to be logged in to leave comments.
Login now