changed yinit to fixed length for nesi

laynep · Sep 25, 2012 · 9689156 · 9689156
1 parent 7305dfd
commit 9689156
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 4 deletions.
diff --git a/d_hybrid_initialconditions.f90 b/d_hybrid_initialconditions.f90
@@ -211,7 +211,6 @@ end function metric
 		Y1(x) = sgn*sqrt(2_dp*rho_kinetic - (Y1(y)*Y1(y)))
 		!Exit condition.
 		if (metric(y0(2:5),y1(2:5)) .le. tol ) exit do1
-    print*,"bad ic"
 	end do do1
 
 	!Set initial conditions.

diff --git a/hybrid_subroutines.f90 b/hybrid_subroutines.f90
@@ -375,7 +375,7 @@ subroutine ic_zoom_init(y0, yref, iccounter, toler)
 	implicit none
 
 	real(dp), dimension(:), intent(inout) :: y0, yref
-  real(dp), dimension(size(yref)) :: yinit
+  real(dp), dimension(5) :: yinit
 	real(dp), intent(inout) :: toler
 	integer, intent(inout) :: iccounter
   integer :: u

diff --git a/parameters_hybrid.txt b/parameters_hybrid.txt
@@ -17,7 +17,7 @@
 	mu = 500D3,
 	nu = .015D3,
 	lambda = 6.8D-3,
-	energy_scale = 1D-1/
+	energy_scale = 1D-2/
 &ics
 	points=10,
 	IC=6,
@@ -39,7 +39,7 @@
 
 !NOTE:   N,phi,psi,phi_dot,psi_dot
 &zoom
-	yinit=0D0, 4.7719688234190150, 168.42510851795384, 3.43932060249110307E-003, 1.37154090993970388E-002,
+	yinit=0D0, 1.5094201026659038, 5.8994452241335926, 5.42316087791986994E-005, 1.13126401897043315E-004,
   toler=1D-2
 /