[CIG-SEISMO] How to get GMT files of movie and shakemap?

[Pieyre Le Loher]

Dear Hu Yuanxin,

Did you run ./configure again before compiling ? If not you should because
I also modified Makefile.in.

Best regards.

Pieyre Le Loher


> Dear Pieyre,
>
> Thank you very much.
>
> I have tested the new code "create_movie_shakemap_AVS_DX_GMT.f90". The
> code has not been passed by compiler. My compiler is
> gFortran-4.1.2-44. The following is error information resulted from
> "make create_movie_shakemap_AVS_DX_GMT.f90" and "gfortran -o
> xcreate_movie_shakemap_AVS_DX_GMT
> create_movie_shakemap_AVS_DX_GMT.f90", respectively.
> -----------make create_movie_shakemap_AVS_DX_GMT.f90---------------------
> make create_movie_shakemap_AVS_DX_GMT
> gfortran  -std=gnu -fimplicit-none -frange-check -O3 -pedantic
> -pedantic-errors -Waliasing -Wampersand -Wline-truncation -Wsurprising
> -Wunderflow -fno-trapping-math -c -o
> obj/create_movie_shakemap_AVS_DX_GMT.o
> create_movie_shakemap_AVS_DX_GMT.f90
> gfortran  -std=gnu -fimplicit-none -frange-check -O3 -pedantic
> -pedantic-errors -Waliasing -Wampersand -Wline-truncation -Wsurprising
> -Wunderflow -fno-trapping-math -o xcreate_movie_shakemap_AVS_DX_GMT
> obj/create_movie_shakemap_AVS_DX_GMT.o obj/read_parameter_file.o
> obj/read_value_parameters.o obj/get_value_parameters.o
> obj/param_reader.o
> obj/create_movie_shakemap_AVS_DX_GMT.o: In function `MAIN__':
> create_movie_shakemap_AVS_DX_GMT.f90:(.text+0x6c74): undefined
> reference to `utm_geo_'
> collect2: ld returned 1 exit status
> make: *** [xcreate_movie_shakemap_AVS_DX_GMT] Error 1
> -------------------------------------------------------------------------------------------------------
>
> -----------------------------gfortran
> -o-------------------------------------------------------------
> gfortran -o xcreate_movie_shakemap_AVS_DX_GMT
> create_movie_shakemap_AVS_DX_GMT.f90
> /tmp/ccc1DeP2.o: In function `MAIN__':
> create_movie_shakemap_AVS_DX_GMT.f90:(.text+0x3f8): undefined
> reference to `read_parameter_file_'
> create_movie_shakemap_AVS_DX_GMT.f90:(.text+0x41f): undefined
> reference to `get_value_string_'
> create_movie_shakemap_AVS_DX_GMT.f90:(.text+0x5e38): undefined
> reference to `utm_geo_'
> collect2: ld returned 1 exit status
> --------------------------------------------------------------------------------------------------------
>
> I will test the new two *.pl in /UTILS/Visualization/GMT/ in two days.
>
> Best regards,
>
> Hu Yuanxin
>
> 2010/8/27 Pieyre Le Loher <Pieyre.Le_Loher at inria.fr>:
>> Dear Hu Yuanxin,
>>
>> I just committed a few modifications in SPECFEM3D so you should run 'svn
>> update'. I modified the path '/opt/seismo-util/lib/perl' in the two Perl
>> scripts you mentioned. Basing on SPECFEM3D-1.4.3, I also modified
>> "create_movie_shakemap_AVS_DX_GMT.f90" in SPECFEM-2.0 to support the GMT
>> format in output.
>>
>> Please tell me if everything works well with you dataset.
>> Thanks.
>>
>> Best regards.
>>
>> Pieyre
>>
>>
>>
>>
>>> Dear Pieyre,
>>>
>>> Thank you for your quick replies.
>>>
>>> I am reading the code "create_movie_shakemap_AVS_DX_GMT.f90" for
>>> getting the output. I found that the code in SPECFEM3D_SESAME could
>>> not get GMT file. Because I will use shakemap in ARCGIS, I must get
>>> GMT *.xyz files since openDX usage is restricted and only allows
>>> viewing of files. Can I replace the code with old code
>>> "create_movie_shakemap_AVS_DX_GMT.f90" in SPECFEM3D-1.4.3 because that
>>> the old one can get GMT files?
>>>
>>> Can I use the script plot_shakemap.pl and movie2gif.pl in
>>> /UTILS/Visualization/GMT/ to convert the binary movie and shakemap
>>> data into GMT xyz files? But I don't know what is the lib (
>>> '/opt/seismo-util/lib/perl' ) in the two script.
>>>
>>> Best regards,
>>>
>>> Hu Yuanxin
>>>
>>
>>
>>
>