CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

help with user specified probability distribution of galaxies in icosmo
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> iCosmo
View previous topic :: View next topic  
Author Message
Augusta Abrahamse



Joined: 23 Jun 2010
Posts: 1
Affiliation: UC Davis

PostPosted: June 23 2010  Reply with quote

I was wondering how to generate Cl's by specifying my own probability distribution of galaxies. (I see someone else has asked about this, but I haven't seen a response...) It looks like the code is already (somewhat) set up to do this, but I'm finding a few bugs. I don't know if someone has worked out the changes to make this run...?

To set things up I do the following:

; (1) set fiducial parameters:
fid=set_fiducial(cosmo_in={omega_m:.25d},calc_in={fit_tk:2},expt_in={sv1_n_zbin:2,sv1_zerror:0.02d})
; (2) create a survey structure:
sv=mk_survey(fid,'sv1',pz_input={z:zvals, pz:pzvals})

here zvals is a vector of z values and pzvals (for now) is a 2rowx401column matrix, (2 rows are the P(z) for each of two z bins and 401 is the length of the z vector):

This however crashes in mk_survey.pro at line 94:
"if keyword_set(pz_input) then d"

If I remove the "d" and combine line 94 with line 95 as:
"if keyword_set(pz_input) then pz=mk_pz(z,dndztype,dndzp,dndzz, n_g=n_g, z_med=z_med, z_min=z_min ,z_max=z_max ,pz_input=pz_input)"

then mk_pz will be called, but various errors keep it from running... I was wondering if someone might have debugged this portion of the code between 2008 (the date in the release) and now. Or is there another way of using this option?
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> iCosmo All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.