Commit 5554f0d0 authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

fix pour Alens

parent ae3084b2
......@@ -2,19 +2,19 @@ BEGIN{
print "engine=pico"
Alens=0
N_nu=0
M_ncdm=0
m_ncdm=0
}
#############
NF==3 && $2 == "Alens"{
Alens=$3
print "fix Alens cosmo ",Alens
$2 == "Alens" {
Alens=$4
}
##############
NF==3 && $2 == "m_ncdm"{
/class/ && $2 == "m_ncdm"{
m_ncdm=$3
print "fix m_ncdm cosmo ",m_ncdm
print "fix m_ncdm cosm ",m_ncdm
}
#############
$2 == "N_eff"{
if (NF==3) {
N_nu=$3+1
......
......@@ -18,10 +18,10 @@ print,+"# free params="+strn(n_elements(names))
print,names
; prof2
;ichi=sz[1]-12
ichi=sz[1]-12
;Mnu
ichi=sz[1]-11
;ichi=sz[1]-11
;
imu=ichi-1
;imu=ichi-2 ; Alens class
......@@ -44,7 +44,7 @@ for iu=0,Nu-1 do begin &$
w=where(m eq mu[iu]) &$
sub=t[w,*] &$
c=min(sub[*,ichi],imin) &$
print,mu[iu],double(c) &$
;print,mu[iu],double(c) &$
bf[iu,*]=transpose(sub[imin,*]) &$
endfor
......@@ -77,4 +77,6 @@ printf,lun,strjoin(strtrim(string(line),2)," ") &$
endfor
free_lun,lun
print,"best fit***************"
print,bf[imin,*]
print,"CHI2=",bf[imin,ichi]
......@@ -11,15 +11,14 @@
var=tau_reio
val=(0.105 0.110 .115)
#val=(0.025 0.035 0.045 0.055 0.065 0.075 0.085 0.095)
#val=(0.01 0.02 0.025 0.03 0.035 0.04 0.045 0.05 0.055 0.06 0.065 0.07 0.075 0.08 0.085 0.09 0.095 0.10 0.12)
#N=20
#xmin=0.02
#xmax=0.10
#for ((i=0;i<N;i++)) ; do
#val[$i]=$(echo | awk -v xmin=$xmin -v xmax=$xmax -v N=$N -v i=$i 'BEGIN{step=(xmax-xmin)/(N-1)} {printf("%f",xmin+(i)*step)}')
#done
N=40
xmin=0.05
xmax=0.25
for ((i=0;i<N;i++)) ; do
val[$i]=$(echo | awk -v xmin=$xmin -v xmax=$xmax -v N=$N -v i=$i 'BEGIN{step=(xmax-xmin)/(N-1)} {printf("%f",xmin+(i)*step)}')
done
file=$(readlink -e $1)
......@@ -134,7 +133,7 @@ cp -f bestfit_class "$OUTDIR"/best_fitClass\${SGE_TASK_ID}
EOBATCH
#####################################################################################################################################
qsub -P P_planck_prod -t 1-5 -pe openmpi_8 8 -q pa_long camelrun
qsub -P P_planck_prod -t 6-10 -pe multicores 8 -q mc_long camelrun
#qsub -P P_planck_prod -t 1-10 -pe openmpi_8 8 -q pa_long camelrun
qsub -P P_planck_prod -t 1-10 -pe multicores 8 -q mc_long camelrun
#####################################################################################################################################
done
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment