[CIG-SEISMO] Fortran runtime Error:

Thu Sep 19 07:37:36 PDT 2013

Dear Emanuele,

I am grateful for your response. I have decompose_mesh.F90 and decompose_mesh_SCOTCH.F90 in src/decompose_mesh_SCOTCH. I can replace decompse_mesh.F90 from the repository. But the error exists in decompose_mesh_SCOTCH.F90. Which I cannot find the repository. What and where shall I make changes?

Thanks once again.

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: Thursday, September 19, 2013 11:40 AM
To: cig-seismo at geodynamics.org
Subject: Re: [CIG-SEISMO] Fortran runtime Error:

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<mailto: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<mailto:cig-seismo-bounces at geodynamics.org> [cig-seismo-bounces at geodynamics.org<mailto:cig-seismo-bounces at geodynamics.org>] on behalf of emanuele casarotti [emanuele.casarotti at gmail.com<mailto:emanuele.casarotti at gmail.com>]
Sent: Wednesday, September 18, 2013 7:05 PM
To: Dimitri Komatitsch
Cc: cig-seismo at geodynamics.org<mailto: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<mailto: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>
>> <mailto: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>
>> <mailto: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> <mailto: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<mailto: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<mailto:CIG-SEISMO at geodynamics.org>
http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
_______________________________________________
CIG-SEISMO mailing list
CIG-SEISMO at geodynamics.org<mailto: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/143a7991/attachment-0001.html>