[CIG-MC] citcoms

[yangan at mail.ustc.edu.cn]

hello sir:
         I met a problem when i use citcoms.py without pyre.In a model when i change Solver=cgrad to =multigrid.
The error information is this:
CitcomSFull: Stokes_flow_Incomp.c:304: solve_Ahat_p_fhat_CG: Assertion `r1dotz1 != 0.0' failed.
how could I solve it?
     
this is my inputfile.            Thank you for your time
# specify output files
datadir="/home/hjs/yangan/data/data94"
datafile="case9"

output_format=ascii
output_optional=surf,botm,horiz_avg
output_ll_max=20

# specify input files
# Important: when file_vbcs is 1, make sure topvbc is also turned on.
file_vbcs=0
vel_bound_file="velbc/velbc_full"

coor=0
coor_file=""

mat_control=0
mat_file=""

lith_age=0
lith_age_file=""

tracer=0
tracer_file=""


# processors, mesh grid information
nproc_surf=12
nprocx=1
nprocy=1
nprocz=1
nodex=33
nodey=33
nodez=33
mgunitx=2
mgunity=2
mgunitz=2
levels=5


# the calculation region definition
radius_inner=0.55
radius_outer=1.0


# restart info
restart=0
post_p=0
datadir_old="."
datafile_old=""
solution_cycles_init=0
zero_elapsed_time=1

# timestepping information
minstep=1
maxstep=11000
maxtotstep=11000
storage_spacing=1000
checkpointFrequency=1000
cpu_limits_in_seconds=360000000

# Initial conditions
tic_method=3
num_perturbations=1
perturbmag=0.01
perturbl=3
perturbm=2
perturblayer=5


# Boundary Conditions
topvbc=0
     topvbxval=0.0
     topvbyval=0.0
botvbc=0
     botvbxval=0.0
     botvbyval=0.0

toptbc=1
     toptbcval=0.0
bottbc=1
     bottbcval=1.0

lith_age_time=0
lith_age_depth=0.031400
mantle_temp=1.000000
temperature_bound_adj=0
depth_bound_adj=0.157000
width_bound_adj=0.087270


# miscellaneous information
stokes_flow_only=0
inputdiffusicity=1.0
rayleigh=7.681756e04
Q0=0


# required information
Problem=convection
Geometry=sphere
Solver=multigrid
node_assemble=1

down_heavy=2
up_heavy=2
vlowstep=20
vhighstep=2
accuracy=0.0001

# Depth Information
# for material above z_lith,    mat = 1;
# for material above z_410,     mat = 2;
# for material above z_lmantle, mat = 3;
# for material below z_lmantle, mat = 4;
z_lith=0.0007848
z_410=0.03924                   #also depth of 410-phase change
z_lmantle=0.06435               #also depth of 670-phase change
z_cmb=0.439                     #     depth of another phase change


# Viscosity Information
Viscosity=system
rheol=1
visc_smooth_method=3

VISC_UPDATE=off
num_mat=4

TDEPV=off
viscE=2.99573,2.99573,2.99573,2.99573
viscT=.5,.5,.5,.5
visc0=1,1,1,1

SDEPV=off
sdepv_expt=1,1,1,1
sdepv_misfit=0.020

low_visc_channel=off
low_visc_wedge=off
lv_min_radius=0.9764
lv_max_radius=0.9921
lv_channel_thickness=0.0047
lv_reduction=0.5

VMIN=on
visc_min=0.01
VMAX=on
visc_max=1000


# Tracer information
tracer=off

tracer_ic_method=0
tracers_per_element=10
tracer_file="tracer.dat"

tracer_flavors=0
ic_method_for_flavors=0
z_interface=0.700000

regular_grid_deltheta=1.0
regular_grid_delphi=1.0

chemical_buoyancy=off
buoy_type=1
buoyancy_ratio=1.0
reset_initial_composition=off


# Phase Change Information
Ra_410=0.0
clapeyron410=0.0235
transT410=0.78
width410=0.0058

Ra_670=0.0
clapeyron670=-0.0235
transT670=0.875
width670=0.0058

Ra_cmb=0.0
clapeyroncmb=-0.0235
transTcmb=0.875
widthcmb=0.0058


# Dimensional Information
radius=6371e3
density=3340.0
thermdiff=1.0e-6
gravacc=9.81
thermexp=3.0e-5
refvisc=1e21
cp=1200
density_above=1030.0
density_below=6600.0