From 8fb811dd961317e97598fdc8220a72b227b83c2c Mon Sep 17 00:00:00 2001 From: Guillaume Raffy Date: Mon, 5 Jun 2023 14:35:51 +0200 Subject: [PATCH] fixed np_f=119 because there are 119 lines in fort. files --- cs_adiabats_2.f | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/cs_adiabats_2.f b/cs_adiabats_2.f index 865690e..1e6fa86 100644 --- a/cs_adiabats_2.f +++ b/cs_adiabats_2.f @@ -10,7 +10,7 @@ implicit none - integer, parameter :: np_f=102 !number of points for function, input + integer, parameter :: np_f=119 !number of points for function, input ! integer, parameter :: nps=10000 !number of points we work with in the splint subroutine. integer, parameter :: nlevtot=176 ! number of levels, some appear multiple times integer, parameter :: npar=2 ! number of parities from molscat @@ -73,7 +73,9 @@ allocate(func(np_f,nf)) ! allocate(proba(nf)) + write(6,*) 'idx_p1 = ', idx_p1 do i=np_f,1,-1 ! read in reverse to start from Rmax + write(6,*) 'i = ', i read(idx_p1,*) R(i),(func(i,j),j=1,nf) ! write(6,*) 'R(', i, ') = ', R(i) ! write(6,*) 'func( ', i,',nf=', nf,') = ', func(i,nf)