[CIG-MC] CitcomS: adding input parameter to pyre version

Eh Tan tan2 at geodynamics.org
Tue Jun 22 19:15:10 PDT 2010


Hi Magali,

You also need to add your input parameter on the python side, in 
CitcomS/Components/BC.py.

Eh


On 6/23/10 6:46 AM, Magali Billen wrote:
> Hello,
>
> I have a question about a parameter that can be read from the input 
> file in the standard version of CitcomS,
> but is not recognized in the pyre version. I'd like to know what 
> changes to make so that it gets read in.
>
> This is actually a more general question than just this one variable, 
> as I'd like to make some other modifications to CitcomS, which
> require creating a few new parameters that need to get read in and I'd 
> like to use nested solvers, so I need to learn how to do this for the 
> pyre version.
>
> Here are the details:
> I'd like to use the variable E->convection.blob_bc_persist which 
> exists in the non pyre version of CitcomS:
>
> [java:CitcomS/CitcomS-3.1.1/lib] billen% grep blob_bc_persist *.c
> Initial_temperature.c:     
>  input_boolean("blob_bc_persist",&(E->convection.blob_bc_persist),"off",m);
> Initial_temperature.c:        if(E->convection.blob_bc_persist){
>
> but does not have a binding
>
> [java:CitcomS/CitcomS-3.1.1/module] billen% grep blob_bc_persist *.c
>
> The other blob related variables are listed in setProperties.c
>
> [java:CitcomS/CitcomS-3.1.1/module] billen% grep blob *.c
> setProperties.c:    getFloatVectorProperty(properties, "blob_center", 
> E->convection.blob_center, 3, fp);
> setProperties.c:    if( E->convection.blob_center[0] == -999.0 && 
> E->convection.blob_center[1] == -999.0 && E->convection.blob_center[2] 
> == -999.0 ) {
> setProperties.c:        E->convection.blob_center[0] = 
> 0.5*(E->control.theta_min+E->control.theta_max);
> setProperties.c:        E->convection.blob_center[1] = 
> 0.5*(E->control.fi_min+E->control.fi_max);
> setProperties.c:        E->convection.blob_center[2] = 
> 0.5*(E->sphere.ri+E->sphere.ro);
> setProperties.c:    getFloatProperty(properties, "blob_radius", 
> E->convection.blob_radius, fp);
> setProperties.c:    getFloatProperty(properties, "blob_dT", 
> E->convection.blob_dT, fp);
>
> So, I added the following line to setProperties.c
>
> getIntProperty(properties, "blob_bc_persist", 
> E->convection.blob_bc_persist, fp);
>
> re-ran make (it did recompile/link setProperties), and then reran my 
> model, but I still get the following error:
>
> >> inplume1.cfg:39:
>  -- pyre.inventory(error)
>  -- regional.ic.blob_bc_persist <- '1'
>  -- unrecognized property 'regional.ic.blob_bc_persist'
> usage: citcoms [--<property>=<value>] 
> [--<facility>.<property>=<value>] [FILE.cfg] ...
> component 'CitcomS'
>     properties: controller, help, help-components, help-persistence, 
> help-properties, job, launcher, nodes, sc
> heduler, solver, steps, typos, weaver
>     facilities: controller,job,launcher,scheduler,solver,weaver
> For more information:
>   --help-properties: prints details about user settable properties
>   --help-components: prints details about user settable facilities and 
> components
> CitcomS: configuration error(s)
>
> Can you tell me what else I have to do to make this get read in by the 
> Pyre version?
>
> Thanks,
> Magali
>
>
> _______________________________________________
> CIG-MC mailing list
> CIG-MC at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc
>    

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-mc/attachments/20100623/20b7a988/attachment.htm 


More information about the CIG-MC mailing list