From 300353ae5f91150e8e129595a7fa925fb89c7707 Mon Sep 17 00:00:00 2001
From: Filipe-Eduard <48655670+Filipe-Eduard@users.noreply.github.com>
Date: Sat, 1 Oct 2022 13:33:28 -0300
Subject: [PATCH] Update principal.f90
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Modifica o output do campo magnético (H) e de sua derivada temporal (dH/dt) para que ambos sejam multiplicados pelo parâmetro de Langevin (alpha). Modifica o diagrama de bifurcação de omega para que a frequência varrida seja correta.
---
 principal.f90 | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/principal.f90 b/principal.f90
index e5b700c..bca10a7 100644
--- a/principal.f90
+++ b/principal.f90
@@ -308,7 +308,7 @@ subroutine principal
 campo(k) = campo(k-1) + (dt_inicial/6.0)*(g1 + 2.0*g2 + 2.0*g3 + g4)
 tempototal(k) = tempototal(k-1) + dt_inicial 
 
-write(300*rea,510) tempototal(k),campo(k),y(k)
+write(300*rea,510) tempototal(k),alpha*campo(k), alpha*y(k)
 
 end do
 end if
@@ -368,8 +368,8 @@ subroutine principal
 if(contfreqinteiro1.ne.contfreqinteiro2) then
 multiplofreq=multiplofreq+1 
 end if
-campo(k)=sin((freqcampo+(bifmax-freqcampo)*multiplofreq)*tempototal(k)) 
-y(k)=(freqcampo + (bifmax-freqcampo)*multiplofreq)*cos((freqcampo + (bifmax-freqcampo)*multiplofreq)*tempototal(k))
+campo(k)=sin((freqcampo+(bifmax-freqcampo)*multiplofreq/nfreq)*tempototal(k)) 
+y(k)=(freqcampo + (bifmax-freqcampo)*multiplofreq/nfreq)*cos((freqcampo + (bifmax-freqcampo)*multiplofreq/nfreq)*tempototal(k))
 else
 tempototal(k) = tempototal(k-1) + dt_inicial 
 campo(k) = sin(freqcampo*tempototal(k)) 
@@ -976,7 +976,7 @@ subroutine principal
 derivada2=(magtempo(2,k)-magtempo(2,k-1))/dt_inicial
 derivada3=(magtempo(3,k)-magtempo(3,k-1))/dt_inicial
 
-write(5*rea,2024) campo(k),y(k), magtempo(1,k),magtempo(2,k),magtempo(3,k), derivada1, derivada2, derivada3, k*dt_inicial
+write(5*rea,2024) alpha*campo(k),alpha*y(k), magtempo(1,k),magtempo(2,k),magtempo(3,k), derivada1, derivada2, derivada3, k*dt_inicial
 
 contfreqinteiro1= ((k-1)*dt_inicial)/intervalo
 contfreqinteiro2= (k*dt_inicial)/intervalo
@@ -989,7 +989,7 @@ subroutine principal
 
 ! Escrever aqui na unidade certa um arquivo mag_tempo para cada frequencia
  
-write(400*rea+multiplofreq+1,2024) campo(k),y(k), magtempo(1,k),magtempo(2,k),magtempo(3,k), derivada1, derivada2, derivada3, k*dt_inicial
+write(400*rea+multiplofreq+1,2024) alpha*campo(k), alpha*y(k), magtempo(1,k),magtempo(2,k),magtempo(3,k), derivada1, derivada2, derivada3, k*dt_inicial
 
 end if