lside = index(controlcard,"SIDE").gt.0
efree = index(controlcard,"EFREE").gt.0
call readi(controlcard,'NTEMP',nT,1)
- write (iout,*) "nT",nT
+ cumul_prob=nt.lt.0
+ nT=iabs(nT)
+ write (iout,*) "nT",nT," cumul_prob",cumul_prob
call multreada(controlcard,'TEMPER',beta_h,nT,300.0d0)
write (iout,*) "nT",nT
write (iout,*) 'beta_h',(beta_h(i),i=1,nT)