[CIG-SEISMO] Fortran runtime Error:

Thu Sep 19 02:40:20 PDT 2013

yes.... it looks like that you have a updated format of
nummaterial_velocity_file
with 8 par and  not updated version of decompose..... see
http://www.geodynamics.org/wsvn/cig/seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90

On Thu, Sep 19, 2013 at 11:19 AM, <s.khan-1 at utwente.nl> wrote:

> Dear All,
>
> The nummaterial_velocity_file has the following contents:
>
> .............................................................................................................................................................................................................
> 2   1          2300.000000          2800.000000          1500.000000
>    9000.000000         90000.000000  0
>
> ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
>
> Thanks
>
> Kind regards,
> Saad Khan
> ________________________________________
> From: cig-seismo-bounces at geodynamics.org [
> cig-seismo-bounces at geodynamics.org] on behalf of emanuele casarotti [
> emanuele.casarotti at gmail.com]
> Sent: Wednesday, September 18, 2013 7:05 PM
> To: Dimitri Komatitsch
> Cc: cig-seismo at geodynamics.org
> Subject: Re: [CIG-SEISMO] Fortran runtime Error:
>
> Looking at the fortran code attached....Maybe it is the opposite...
> the qk could be in the nummaterial file but not in the decompose ...
> Not the current version of the code?
>
> Emanuele
>
> Il giorno 18/set/2013, alle ore 18:52, Dimitri Komatitsch
> <komatitsch at lma.cnrs-mrs.fr> ha scritto:
>
> >
> > Dear all,
> >
> > Could it be the Qkappa missing from nummaterial_velocity_file?
> >
> > Best regards,
> > Dimitri.
> >
> > On 09/18/2013 06:14 PM, emanuele casarotti wrote:
> >> How is your nummaterial_velocity_file?
> >>
> >> Emanuele
> >>
> >> Il giorno 18/set/2013, alle ore 17:03, "s.khan-1 at utwente.nl
> >> <mailto:s.khan-1 at utwente.nl>" <s.khan-1 at utwente.nl
> >> <mailto:s.khan-1 at utwente.nl>> ha scritto:
> >>
> >>> Dear All,
> >>>
> >>> I am facing problem with the decomposer. 'make xdecompose_mesh_SCOTCH'
> >>> worked successfully.
> >>> I have 4 processor and typing command as:
> >>>
> ........................................................................................................................................................................................................................................
> >>> khan31332 at itclx03:~/SPECFEM3D$ *./bin/xdecompose_mesh_SCOTCH 4
> >>> examples/homogeneous_halfspace/MESH/ in_out_files/DATABASES_MPI/*
> >>> total number of nodes:
> >>>   nnodes =        23273
> >>> total number of spectral elements:
> >>>   nspec =        20736
> >>> materials:
> >>>   num_mat =            1
> >>>   defined =            1 undefined =            0
> >>> *At line 326 of file decompose_mesh_SCOTCH.f90*
> >>> *Fortran runtime error: Bad integer for item 7 in list input*
> >>>
> >>>
> .........................................................................................................................................................................................................................................
> >>> I am looking into decompose_mesh_SCOTCH.f90 (the part is given below).
> >>> But I am unable to find any problem there.
> >>>
> >>>
> .........................................................................................................................................................................................
> >>> ! material definitions
> >>>       !
> >>>       ! format: note that we save the arguments in a slightly
> >>> different order in mat_prop(:,:)
> >>>       !              #(6) material_domain_id #(0) material_id  #(1)
> >>> rho #(2) vp #(3) vs #(4) Q_mu #(5) anisotropy_flag
> >>>       !
> >>>       !read(98,*) idomain_id,num_mat,rho,vp,vs,Q_mu,aniso_flag
> >>>       ! reads lines unti it reaches a defined material
> >>>       num_mat = -1
> >>>       do while( num_mat < 0 .and. ier == 0)
> >>>         read(98,'(A256)',iostat=ier) line
> >>>         read(line,*) idomain_id,num_mat
> >>>       enddo
> >>>       if( ier /= 0 ) stop 'error reading in defined materials in
> >>> nummaterial_velocity_file'
> >>>
> >>>       ! reads in defined material properties
> >>>       read(line,*) idomain_id,num_mat,rho,vp,vs,qmu,aniso_flag
> >>>
> >>>       ! checks material_id bounds
> >>>       if(num_mat < 1 .or. num_mat > count_def_mat)  stop "ERROR :
> >>> Invalid nummaterial_velocity_file file."
> >>>
> >>>       if(idomain_id == 1 .or. idomain_id == 2) then
> >>>         ! material is elastic or acoustic
> >>>
> >>>         !read(98,*) num_mat, mat_prop(1,num_mat),mat_prop(2,num_mat),&
> >>>         !
> >>> mat_prop(3,num_mat),mat_prop(4,num_mat),mat_prop(5,num_mat)
> >>>         mat_prop(1,num_mat) = rho
> >>>         mat_prop(2,num_mat) = vp
> >>>         mat_prop(3,num_mat) = vs
> >>>         mat_prop(4,num_mat) = qmu
> >>>         mat_prop(5,num_mat) = aniso_flag
> >>>         mat_prop(6,num_mat) = idomain_id
> >>>
> >>>       else
> >>>         ! material is poroelastic
> >>>         if( use_poroelastic_file .eqv. .false. ) stop 'error
> >>> poroelastic material requires nummaterial_poroelastic_file'
> >>>
> .........................................................................................................................................................................................
> >>>
> >>>
> >>>
> >>>
> >>> Thanks in advance for your reply and attention to my problem.
> >>>
> >>> Kind regards,
> >>> Saad Khan
> >>>
> >>> _______________________________________________
> >>> CIG-SEISMO mailing list
> >>> CIG-SEISMO at geodynamics.org <mailto:CIG-SEISMO at geodynamics.org>
> >>> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
> >>
> >>
> >> _______________________________________________
> >> CIG-SEISMO mailing list
> >> CIG-SEISMO at geodynamics.org
> >> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
> >
> > --
> > Dimitri Komatitsch
> > CNRS Research Director (DR CNRS), Laboratory of Mechanics and Acoustics,
> > UPR 7051, Marseille, France    http://komatitsch.free.fr
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-seismo/attachments/20130919/8570a605/attachment.html>