graph_dislin.f90 2 KB
subroutine graph_dislin
use parameters
use dislin

IMPLICIT NONE
double precision, allocatable :: axe(:),Ep(:,:)
real ::par1,par2,par3,par4,par5,par6,par7,par8,par9,par10,par11,par12

allocate(axe(nbin))
allocate(Ep(nbin,nfreq))

axe(1)=0
do i=1,nbin-1
axe(i+1)=i*dx/1000
end do


!_____________________________________maximum floe size________________

      !CALL SETPAG('DA4P')
      !CALL METAFL('CONS')
      !CALL SCRMOD('REVERS')
      !CALL DISINI()
      !CALL PAGERA()
      !CALL SETRGB (0.5,0.2,0.4)
      !CALL GRAF (real(axe(1),4),real(axe(nbin),4),real(0,4),real(10*dx/1000,4),real(0,4),real(maxval(Dmax),4),real(0,4),real(50,4))
      !call NAME('x axis [km]','X')
      !call NAME('Dmax [m]','Y')
      !CALL TITLIN ('Maximum Floe size in the MIZ',1)
      !call TITLE

     !   call COLOR('GRAY')
    !  call GRID(10,5)

   !    call COLOR('RED')
  !    CALL CURVE (real(axe,4),real(Dmax,4),nbin)
    
 !     call ENDGRF
!call disfin

!________________________SPECTRUM_______________________________________

do i=1,nbin
do ii=1,nfreq

Ep(i,ii)=real(E(15,i,ii),4)

end do
end do

par5=real(omega(1),4)
par6=real(omega(nfreq),4)
par7=real(omega(1),4)
par8=real(1,4)
par1=real(axe(1),4)
par2=real(axe(nbin),4)
par3=real(0,4)
par4=real(10*dx/1000,4)
par9=real(0.,4)
par10=real(maxval(Ei),4)
par11=real(0.,4)
par12=real(maxval(Ei)/10,4)

!write(*,*)par1,par2,par3,par4,par5,par6,par7,par8,par9,par10,par11,par12



CALL SETPAG('DA4P')
      CALL METAFL('CONS')
      !CALL SCRMOD('REVERS')
  
    CALL DISINI()
      CALL PAGERA()
CALL VIEW3D(-5.,-5.,4.,'ABS')

do n=1,nsteps

do i=1,nbin
do ii=1,nfreq

Ep(i,ii)=real(E(n,i,ii),4)

end do
end do
call ERASE
call COLOR('WHITE')
CALL GRAF3D (par1,par2,par3,par4,par5,par6,par7,par8,par9,par10,par11,par12)
           CALL NAME('X-axis [km]','X')
      CALL NAME('Frequency','Y')
      CALL NAME('Energy','Z')
        call COLOR('BLUE')
        CALL SURMAT (real(Ep,4),nbin,nfreq,1,1)
call ENDGRF

        call system('sleep 0.1')
        

end do
call disfin


end subroutine graph_dislin