Skip to content
Snippets Groups Projects
Commit 24ace58e authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann
Browse files

back to poke initalization

parent 71e00e53
No related branches found
No related tags found
No related merge requests found
SUBROUTINE inital
USE basic
USE model, ONLY : CO
USE prec_const
......@@ -39,31 +39,31 @@ SUBROUTINE init_profiles
CALL RANDOM_SEED(PUT=iseedarr)
CALL set_updatetlevel(1)
DO ikr=ikrs,ikre
DO ikr=ikrs,ikre
DO ikz=ikzs,ikze
DO ip=ips_e,ipe_e
DO ij=ijs_e,ije_e
CALL RANDOM_NUMBER(noise)
moments_e( ip,ij, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
!moments_e( ip,ij, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
END DO
END DO
DO ip=ips_i,ipe_i
DO ij=ijs_i,ije_i
CALL RANDOM_NUMBER(noise)
moments_i( ip,ij, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
!moments_i( ip,ij, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
END DO
END DO
! Poke initialization on only Ne00 and Ni00
!moments_e( 1,1, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
!moments_i( 1,1, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
moments_e( 1,1, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
moments_i( 1,1, ikr,ikz, :) = initback_moments + initnoise_moments*(noise-0.5_dp)
END DO
END DO
CALL poisson ! To set phi
END SUBROUTINE init_profiles
......@@ -77,42 +77,58 @@ SUBROUTINE load_FC_mat ! Works only for a partiular file for now with P,J <= 15,
USE futils, ONLY : openf, getarr, closef
IMPLICIT NONE
INTEGER :: fid, ip,ij, ip2,ij2
INTEGER :: ip,ij, ip2,ij2
INTEGER :: fid1, fid2, fid3, fid4
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!! Electron matrices !!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
WRITE(*,*) 'Load self electron collision matrix..'
! get the self electron colision matrix
CALL openf(selfmat_file,fid, 'r', 'D');
CALL getarr(fid, '/Caapj/Ceepj', Ceepj) ! get array (moli format)
CALL closef(fid)
CALL openf(selfmat_file,fid1, 'r', 'D');
CALL getarr(fid1, '/Caapj/Ceepj', Ceepj) ! get array (moli format)
CALL closef(fid1)
! get the Test and Back field electron ion collision matrix
WRITE(*,*) 'Load test + field electron collision matrix..'
CALL openf(eimat_file,fid, 'r', 'D');
CALL getarr(fid, '/Ceipj/CeipjT', CeipjT)
CALL getarr(fid, '/Ceipj/CeipjF', CeipjF)
CALL closef(fid)
CALL openf(eimat_file,fid2, 'r', 'D');
CALL getarr(fid2, '/Ceipj/CeipjT', CeipjT)
CALL getarr(fid2, '/Ceipj/CeipjF', CeipjF)
CALL closef(fid2)
write(25,*) eimat_file
write(25,*) 'Ceepj(3,5)'
write(25,*) Ceepj(3,5)
write(25,*) 'CeipjT(12,12)'
write(25,*) CeipjT(12,12)
write(25,*) 'CeipjF(12,2)'
write(25,*) CeipjF(12,2)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!! Ion matrices !!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
WRITE(*,*) 'Load self ion collision matrix..'
! get the self electron colision matrix
CALL openf(selfmat_file, fid, 'r', 'D');
CALL openf(selfmat_file, fid3, 'r', 'D');
IF ( (pmaxe .EQ. pmaxi) .AND. (jmaxe .EQ. jmaxi) ) THEN ! if same degrees ion and electron matrices are alike so load Ceepj
CALL getarr(fid, '/Caapj/Ceepj', Ciipj) ! get array (moli format)
CALL getarr(fid3, '/Caapj/Ceepj', Ciipj) ! get array (moli format)
ELSE
CALL getarr(fid, '/Caapj/Ciipj', Ciipj) ! get array (moli format)
CALL getarr(fid3, '/Caapj/Ciipj', Ciipj) ! get array (moli format)
ENDIF
CALL closef(fid)
CALL closef(fid3)
! get the Test and Back field ion electron collision matrix
WRITE(*,*) 'Load test + field ion collision matrix..'
CALL openf(iemat_file, fid, 'r', 'D');
CALL getarr(fid, '/Ciepj/CiepjT', CiepjT)
CALL getarr(fid, '/Ciepj/CiepjF', CiepjF)
CALL closef(fid)
END SUBROUTINE load_FC_mat
\ No newline at end of file
CALL openf(iemat_file, fid4, 'r', 'D');
CALL getarr(fid4, '/Ciepj/CiepjT', CiepjT)
CALL getarr(fid4, '/Ciepj/CiepjF', CiepjF)
CALL closef(fid4)
write(25,*) iemat_file
write(25,*) 'Ciipj(3,5)'
write(25,*) Ciipj(3,5)
write(25,*) 'CiepjT(10,10)'
write(25,*) CiepjT(10,10)
write(25,*) 'CiepjF(10,2)'
write(25,*) CiepjF(10,2)
END SUBROUTINE load_FC_mat
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment