Status of Week 3, 2024

[PaulWessel]

Hi @Esteban82 and @gd-a. Happy new year. Just to summarise what I have been doing the last few days/weeks:

1. The GMT Team finally released GMT 6.5 which we needed for our animation paper - so it held progress back a bit.
2. GMT 6.5 also comes with much new data sets, at least one (Moon) is featured in a figure.
3. I will use the word "figure" for either a static PNG or a MP4 movie. I hope the AGU journal we will publish in (I will push for G^3) will finally allow us to have movies on the fly as Fig (X) as we will absolutely not settle for retarded zip supplements etc. If AGU cannot do this we will go to Elsevier.
4. I have renamed the repo to gmt-2024-animation since we are now in 2024. You may need to clone this one as the gmt-2023-animation is gone and I had some trouble doing this. All is well for me now.
5. All scripts (makes no difference if figs or movies) will be named WED-AFig?.sh. I have done this for all my figs but I have left yours as is.
6. Your "homework" is among other things to rename your scripts (and products) accordingly. Have a look at any of mine and you will see a common comment block with authors, fig number, purpose etc. Please put that at the top of your script. We will add maybe a zenodo URL to these at some point; maybe that is how we distribute some data.
7. I have made sure mine are all using SI units (plus points); no more inches.
8. I have eliminated sub-directories for output. All products land in the main dir and all *.png/mp4 are excluded by .gitignore.
9. The Makefile can make and remove stuff; try make help for overview.
10. I will add a science animation to this as well. @gd-a has that but @Esteban82 only has his "toy" plots. Since we are pushing this article for improving and simplify animation in science we should each have one such movie. @federico does nice seismology ones - might that be one extra from you? I will do some plate flexure or seamount related.
11. The science animations will require data sets that may or may not be online somewhere and are too large to stick in this repo. Our repo should allow anyone to clone it, run make all, and basically recreate the entire paper.
12. The biggest and ambitious scripts are probably too large to print in the journal but could be git links to the repo.
13. Make a plan how you will make data available to users who gets the repo. I may know more about this once I work on my science movie and realise I need to somehow share tons of Gb...

OK gents, please get started. There is of course more typing to be done in the paper but finalising the scripts is probably top priority.

[PaulWessel]

Because it is a private repo I cannot assign both of you to do a review etc so I will just ping you via @gd-a and @Esteban82 in the issues. Please let me know that you are aware of this post even if you have no time to act the coming week.

[Esteban82]

I had an animation where I calculated the upward continuation of the Free Air Anomaly for a range of altitudes. It shows and helps to choose the best parameters to apply to the data.

[gd-a]

As far as I'm concerned, the precipitation movie is done. I can review the code to fit your requirements but there's not much left I believe. Concerning the seismicity one, bugs in grdview and -perspective are bothering me and, with my new position, I don't have so much time to find a workaround.

[PaulWessel]

All good, Guillaume. I have bug troubles too so shall see if I even can wrap up another science movie. Yep, please do what polish you can on your precipitation, rename, comment etc and I will have a look. Federico, feel free to add you other movie.

[PaulWessel]

We will go from simple to complicated so the preciptation movie will culminate the paper. I will add one on seamount densities, and Federico can add the gravity or seismicity, depending. I think I will add a data directory for files used by scripts. TBD.

[gd-a]

the preciptation movie will culminate the paper.

Isn't the sea-ice concentration one better ? IMO, the precipitation is a little boring.

[PaulWessel]

Oh sure, up to you to pick it I think. Sea-ice is great.

[gd-a]

I think it's up to you really. Your the main author :) People really liked the heat map too. It shouldn't be too difficult to adjust it to be a cumulative animation... (with seismicity ?)

[gd-a]

8. I have eliminated sub-directories for output. All products land in the main dir and all *.png/mp4 are excluded by .gitignore.

Are you sure about that ? No "TMP" directory ?

13. Make a plan how you will make data available to users who gets the repo. I may know more about this once I work on my science movie and realise I need to somehow share tons of Gb...

You can always provide a "download data" script. Mine points to Copernicus (ECMWF)... Though, it might takes more than just a few hours to download it all...

[PaulWessel]

Do you need a special TMP dir? Can't use /tmp ? Yes, if it is most practical with a download script then the makefile can run that if files are not locally placed.

[PaulWessel]

I see your code commit, e,g

he following lines of code aggregate several components :
#   * On the go :
#       A simple routine to declare the working directory
#
#   * INCLUDE.SH :
#       Generated file that retains all the constants the movie needs. It is overly dense to
#       highlight some of the many customization GMT offers
#
#   * PRE.SH :
#       Generated fi

etc. You want me to add that to the new repo?

[gd-a]

Sorry, I deleted the comment because GitHub directly transformed my commit number into a link where I was able to retrieve the content. Precipitation (completed) is in gd-a branch Seismicity (static, unfinished) is in gd-a branch

If you find half an hour in you schedule, you can run "WED-A_Fig_S.sh" to see where my problem is. I think I will modify it to add the heat map aforementioned

[PaulWessel]

OK running it now - will let you know how it goes.

[PaulWessel]

Here is the full output from bash -xv

bash -xv WED-A_Fig_S.sh
#!/usr/bin/bash 

#   #   #   #   #   #   #
#   A GMT SCRIPT :
#   - modules : 
#       movie, events, grdview, plot3d, coast, basemap, math, info
#
#   - acknowledgment :
#       Institut de Physique du Globe de Paris [IPGP] / Research Collection
#       for more info, visit https://dataverse.ipgp.fr/dataset.xhtml?persistentId=doi:10.18715/IPGP.2022.l3ylhaal
#   
#   - data :    * .csv              :
#                   event timestamp, location and identifier
#               * @earth_relief_    :
#                   remote data from IGPP Global Earth Relief
#
#   - objective : provide an 3D animation of subduction's zone seismicity
#

# Hypocenters data
DATA_file="DATA/2014-2019_catalog.csv"
+ DATA_file=DATA/2014-2019_catalog.csv

# # Topography data
res_topo="01s"
+ res_topo=01s
gmt grdmath -RGP @earth_synbath_${res_topo}_g 0 GT 0 NAN @earth_synbath_${res_topo}_g MUL = topo_island.nc
+ gmt grdmath -RGP @earth_synbath_01s_g 0 GT 0 NAN @earth_synbath_01s_g MUL = topo_island.nc
grdblend [NOTICE]: Remote data courtesy of GMT data server oceania [http://oceania.generic-mapping-tools.org]

grdblend [NOTICE]: Earth Relief at 1x1 arc seconds tiles provided by SRTMGL1 (land only) [NASA/USGS].
grdblend [NOTICE]:   -> Download 1x1 degree grid tile (earth_relief_01s_g): N15W062
grdblend [NOTICE]:   -> Download 1x1 degree grid tile (earth_relief_01s_g): N16W062
grdblend [NOTICE]: SRTM15 Earth Relief v2.5.5 original at 15x15 arc seconds [Tozer et al., 2019].
grdblend [NOTICE]:   -> Download 10x10 degree grid tile (earth_relief_15s_p): N10W070

# Area of Interest
SOUFRIERE="-61.663406/16.043829"
+ SOUFRIERE=-61.663406/16.043829

max_dist_volc=12
+ max_dist_volc=12
west_lim=$(gmt math -Q $(echo ${SOUFRIERE} | awk -F[/] '{print $1}') ${max_dist_volc}  KM2DEG SUB -Vq =)
+++ echo -61.663406/16.043829
+++ awk '-F[/]' '{print $1}'
++ gmt math -Q -61.663406 12 KM2DEG SUB -Vq =
+ west_lim=-61.7713244711
east_lim=$(gmt math -Q $(echo ${SOUFRIERE} | awk -F[/] '{print $1}') ${max_dist_volc}  KM2DEG ADD -Vq =)
+++ echo -61.663406/16.043829
+++ awk '-F[/]' '{print $1}'
++ gmt math -Q -61.663406 12 KM2DEG ADD -Vq =
+ east_lim=-61.5554875289
south_lim=$(gmt math -Q $(echo ${SOUFRIERE} | awk -F[/] '{print $2}') ${max_dist_volc}  KM2DEG SUB -Vq  =)
+++ echo -61.663406/16.043829
+++ awk '-F[/]' '{print $2}'
++ gmt math -Q 16.043829 12 KM2DEG SUB -Vq =
+ south_lim=15.9359105289
north_lim=$(gmt math -Q $(echo ${SOUFRIERE} | awk -F[/] '{print $2}') ${max_dist_volc}  KM2DEG ADD -Vq  =)
+++ echo -61.663406/16.043829
+++ awk '-F[/]' '{print $2}'
++ gmt math -Q 16.043829 12 KM2DEG ADD -Vq =
+ north_lim=16.1517474711

DOMAIN_region=$(gmt mapproject -RGP,MQ+r1 -WR)
++ gmt mapproject -RGP,MQ+r1 -WR
+ DOMAIN_region=-R298/300/14/17
DOMAIN_volc="-R${west_lim}/${east_lim}/${south_lim}/${north_lim}"
+ DOMAIN_volc=-R-61.7713244711/-61.5554875289/15.9359105289/16.1517474711

# Figure's dimension
width_volc=$(gmt math -Q 9 2 SQRT MUL =) # diagonal from 9cm edge
++ gmt math -Q 9 2 SQRT MUL =
+ width_volc=12.7279220614
posx_volc=2
+ posx_volc=2
posy_volc=2
+ posy_volc=2

width_region=5
+ width_region=5
posx_region=$(gmt math -Q ${width_volc} ${posx_volc} ADD 9 ADD =)
++ gmt math -Q 12.7279220614 2 ADD 9 ADD =
+ posx_region=23.7279220614
posy_region=3.5
+ posy_region=3.5

width_profile=${width_region}
+ width_profile=5
posx_profile=${posx_region} # $(gmt math -Q ${posx_region} 1.5 SUB =)
+ posx_profile=23.7279220614
posy_profile=$(gmt math -Q ${posy_region} $(gmt mapproject ${DOMAIN_region} -JM${width_region} -Wh) ADD 2 ADD =)
+++ gmt mapproject -R298/300/14/17 -JM5 -Wh
++ gmt math -Q 3.5 7.73569317523 ADD 2 ADD =
+ posy_profile=13.2356931752

# 3D view
azimuth=135
+ azimuth=135
elevation=20
+ elevation=20
PERSPECTIVE="${azimuth}/${elevation}"
+ PERSPECTIVE=135/20

max_depth_data=$(gmt math -C3 -o3 ${DATA_file} -1 MUL LOWER -Sf =)
++ gmt math -C3 -o3 DATA/2014-2019_catalog.csv -1 MUL LOWER -Sf =
gmtmath [ERROR]: DATA/2014-2019_catalog.csv is not a number, operator or file name
+ max_depth_data=
max_depth_3d=$(gmt math -Q $max_depth_data 1000 MUL =)
++ gmt math -Q 1000 MUL =
gmtmath [ERROR]: Operation "MUL" requires 2 operands
+ max_depth_3d=
# max_depth_coupe=-137000 #gmt histogram $DATA_file -i3 -Glightblue -W -T1 -Baf -png test
max_alt=$(gmt grdinfo topo_island.nc -Mf | grep "v_max" | awk '{print $5}')
++ gmt grdinfo topo_island.nc -Mf
++ grep v_max
++ awk '{print $5}'
+ max_alt=1456

# Topo/bathy profiles
cat <<- EOF > track.txt
-61.6525093993956 16.02593225832284
-61.67133858104974 16.0570823829057
EOF
+ cat
gmt grdtrack ${DOMAIN_region} track.txt -G@earth_synbath_${res_topo}_g -C200k/0.1/0.25+v -Sm+sstack.txt -Vq > table.txt
+ gmt grdtrack -R298/300/14/17 track.txt -G@earth_synbath_01s_g -C200k/0.1/0.25+v -Sm+sstack.txt -Vq
grdblend [NOTICE]: Remote data courtesy of GMT data server oceania [http://oceania.generic-mapping-tools.org]

grdblend [NOTICE]: Earth Relief at 1x1 arc seconds tiles provided by SRTMGL1 (land only) [NASA/USGS].
grdblend [NOTICE]:   -> Download 1x1 degree grid tile (earth_relief_01s_g): N14W062
gmt math stack.txt 0 LE 0 NAN stack.txt MUL = sea_area.txt
+ gmt math stack.txt 0 LE 0 NAN stack.txt MUL = sea_area.txt

# # # # # # # # # 
gmt begin test png
+ gmt begin test png
    gmt makecpt -Coleron -T0/${max_alt} -H > relief.cpt
+ gmt makecpt -Coleron -T0/1456 -H
    arrow_lon=$(gmt math -Q ${azimuth} -1 MUL =)
++ gmt math -Q 135 -1 MUL =
+ arrow_lon=-135

    # 3D view : hypocenters
    gmt plot3d ${DOMAIN_volc}/${max_depth_3d}/0 -JM${width_volc} -JZ12c  \
        ${DATA_file} -i2,1,3+d-0.001 -Se${azimuth}/0.5/$(gmt math -Q 0.5 ${elevation} SIND MUL =) -Gred \
        -BwsneZ3+b+zlightbrown -Bz50000+l"depth (m)" \
        -Xa${posx_volc} -Ya${posy_volc} \
        -p${PERSPECTIVE}
++ gmt math -Q 0.5 20 SIND MUL =
+ gmt plot3d -R-61.7713244711/-61.5554875289/15.9359105289/16.1517474711//0 -JM12.7279220614 -JZ12c DATA/2014-2019_catalog.csv -i2,1,3+d-0.001 -Se135/0.5/0.171010071663 -Gred -BwsneZ3+b+zlightbrown '-Bz50000+ldepth (m)' -Xa2 -Ya2 -p135/20
plot3d [ERROR]: Cannot find file DATA/2014-2019_catalog.csv

    # 3D view : surface
    gmt coast ${DOMAIN_volc} -JM${width_volc} \
    -Wthick,black -Df -Slightblue -t75 \
    -Xa${posx_volc} -Ya$(gmt math -Q ${posy_volc} 12 ${elevation} COSD MUL ADD =) \
    -p
++ gmt math -Q 2 12 20 COSD MUL ADD =
+ gmt coast -R-61.7713244711/-61.5554875289/15.9359105289/16.1517474711 -JM12.7279220614 -Wthick,black -Df -Slightblue -t75 -Xa2 -Ya13.2763114494 -p

    gmt grdview ${DOMAIN_volc}/0/${max_alt} -JZ1c \
        topo_island.nc -Crelief.cpt -Qc -I+d \
        -Bafg \
        -Xa${posx_volc} -Ya$(gmt math -Q ${posy_volc} 12 ${elevation} COSD MUL ADD =) \
        -p
++ gmt math -Q 2 12 20 COSD MUL ADD =
+ gmt grdview -R-61.7713244711/-61.5554875289/15.9359105289/16.1517474711/0/1456 -JZ1c topo_island.nc -Crelief.cpt -Qc -I+d -Bafg -Xa2 -Ya13.2763114494 -p

# # # 
# # # 
# # # 
    # Minimap : coast
    gmt coast ${DOMAIN_region} -JM${width_region}c -W -Slightblue -Glightbrown \
        -TdjBL+f1+l \
        -Ba0.5f0.5g0.5 \
        -Xa${posx_region} -Ya${posy_region} \
        --FORMAT_GEO_MAP="ddd.xxx" --FONT_ANNOT_PRIMARY=4p
+ gmt coast -R298/300/14/17 -JM5c -W -Slightblue -Glightbrown -TdjBL+f1+l -Ba0.5f0.5g0.5 -Xa23.7279220614 -Ya3.5 --FORMAT_GEO_MAP=ddd.xxx --FONT_ANNOT_PRIMARY=4p

    gmt basemap -Bwesn -p90 \
        -LjBR+w$(gmt math -Q 1 DEG2KM RINT -Vq =)k+o-0.15c/$(gmt math -Q -${width_region} 0.3 ADD =)c+u+c+f \
        -Xa${posx_region} -Ya${posy_region} \
        --FONT_ANNOT_PRIMARY=4p --MAP_SCALE_HEIGHT=3p 
++ gmt math -Q 1 DEG2KM RINT -Vq =
++ gmt math -Q -5 0.3 ADD =
+ gmt basemap -Bwesn -p90 -LjBR+w111k+o-0.15c/-4.7c+u+c+f -Xa23.7279220614 -Ya3.5 --FONT_ANNOT_PRIMARY=4p --MAP_SCALE_HEIGHT=3p

    # Minimap : hypocenters (within radius)
    gmt select ${DATA_file} -i2,1 -fg -C${SOUFRIERE}+d50k > subset.txt
+ gmt select DATA/2014-2019_catalog.csv -i2,1 -fg -C-61.663406/16.043829+d50k
gmtselect [ERROR]: Cannot find file DATA/2014-2019_catalog.csv
gmtselect [ERROR]: Option -C: Cannot read file -61.663406/16.043829!
    gmt plot ${DATA_file} -i2,1 -Sc0.1c -Glightgray  -Xa${posx_region} -Ya${posy_region} -t80
+ gmt plot DATA/2014-2019_catalog.csv -i2,1 -Sc0.1c -Glightgray -Xa23.7279220614 -Ya3.5 -t80
plot [ERROR]: Cannot find file DATA/2014-2019_catalog.csv
    gmt plot subset.txt -Sc0.1c -Gblack -Xa${posx_region} -Ya${posy_region} -t50
+ gmt plot subset.txt -Sc0.1c -Gblack -Xa23.7279220614 -Ya3.5 -t50

    # Minimap : swath
    gmt plot table.txt -Wfaint,black -Xa${posx_region} -Ya${posy_region}
+ gmt plot table.txt -Wfaint,black -Xa23.7279220614 -Ya3.5

    # # Minimap : dome
    # echo $(echo ${SOUFRIERE} | awk -F[/] '{printf "%f %f", $1,$2}') | gmt plot -Sx0.25c -Wred -Xa${posx_region} -Ya${posy_region} # Marker of La Soufriere Volcano

    # Minimap : rotating arrow
    gmt plot -Rg -JG180/90/3c -SV0.4c+ea+h0+a60+gred -Wfatter,red \
    -Xa$(gmt math -Q ${posx_region} 0.66 SUB =) -Ya$(gmt math -Q ${posy_region} 3.75 ADD =) <<- EOF
    ${arrow_lon} 5 0 0.75c
    EOF
++ gmt math -Q 23.7279220614 0.66 SUB =
++ gmt math -Q 3.5 3.75 ADD =
+ gmt plot -Rg -JG180/90/3c -SV0.4c+ea+h0+a60+gred -Wfatter,red -Xa23.0679220614 -Ya7.25
    # gmt basemap -Bafg -Xa$(gmt math -Q ${posx_region} 0.66 SUB =) -Ya$(gmt math -Q ${posy_region} 3.75 ADD =) # Frame for "polar" projection

# # # 
# # # 
# # # 
    # Profile : hypocenters
    gmt plot subset.txt  -R-40/100/-10000/3000 -JX${width_profile}/${width_region} \
        -Xa${posx_profile} -Ya${posy_profile}
+ gmt plot subset.txt -R-40/100/-10000/3000 -JX5/5 -Xa23.7279220614 -Ya13.2356931752
plot [WARNING]: File subset.txt is empty!

    # Profile : sea_area
    gmt plot sea_area.txt -Glightblue -Sb1q+b0 -Xa${posx_profile} -Ya${posy_profile}
+ gmt plot sea_area.txt -Glightblue -Sb1q+b0 -Xa23.7279220614 -Ya13.2356931752
    gmt plot -W -Xa${posx_profile} -Ya${posy_profile} <<- EOF
    -40 0
    100 0
    EOF
+ gmt plot -W -Xa23.7279220614 -Ya13.2356931752

    # Profile : topo thickness
    gmt plot -Wthick stack.txt -i0,1,5,6 -L+b -Glightgray -Xa${posx_profile} -Ya${posy_profile}
+ gmt plot -Wthick stack.txt -i0,1,5,6 -L+b -Glightgray -Xa23.7279220614 -Ya13.2356931752
    echo "WSW" | gmt text -F+cTL+fHelvetica-Bold -DJ0.1 -Xa${posx_profile} -Ya${posy_profile}
+ echo WSW
+ gmt text -F+cTL+fHelvetica-Bold -DJ0.1 -Xa23.7279220614 -Ya13.2356931752
    echo "ENE" | gmt text -F+cTR+fHelvetica-Bold -DJ0.1 -Xa${posx_profile} -Ya${posy_profile}
+ echo ENE
+ gmt text -F+cTR+fHelvetica-Bold -DJ0.1 -Xa23.7279220614 -Ya13.2356931752

    gmt basemap  -Bxafg1000+l"Distance from dome (km)" -Byaf+l"Depth (m)" -BWSne -Xa${posx_profile} -Ya${posy_profile}
+ gmt basemap '-Bxafg1000+lDistance from dome (km)' '-Byaf+lDepth (m)' -BWSne -Xa23.7279220614 -Ya13.2356931752

gmt end show
+ gmt end show

[gd-a]

So what I'd like is :

1. Control the rotation of the left part with [-p+wlon0/lat0] (https://github.com/GenericMappingTools/gmt/issues/7066)
2. avoid the whitespace between the topography and the "sea area" (coast not aligned with grdview)
3. make grdview more transparent (I tried different options and none worked, haven't tried again since 6.5 upgrade)
4. project the hypocenters on the top right figure (profile from track)

That being said, I'm not sure the left part can be satisfying enough without fixing (1), (2) and (3)... So if it's too complicated, maybe we should go for another approach.

[PaulWessel]

Will get back to you on this but may not happen right away - other stuff is happening...

[Esteban82]

Hi @PaulWessel. In this link I had post it my animation of the upward continuation of the free air anomaly. I think I could change the progression of the altitudes. What do you think?

https://forum.generic-mapping-tools.org/t/animation-of-upward-continuation-of-free-air-anomaly/838

[PaulWessel]

Had a hard time seeing much movements in the grid (just the spectrum lines). Not sure if it is the exciting ... what other things might you have?

[gd-a]

Update, I tried :

cat << 'EOF' > main.sh
SOUFRIERE="-61.663406/16.043829"

gmt begin
gmt grdview -R-62/-61/15.75/16.75/-300/2000 @earth_relief_01s -JM10c -JZ1c -Bafg -I+d -Qc -C -p${MOVIE_FRAME}/45+w${SOUFRIERE} -t50
echo "-61.663406 16.043829" | gmt plot -Sc0.3 -Gred -p
gmt end
EOF

gmt movie main.sh -Ntest -C720p -T0/180/10 -D5 -Fgif

And it works fine apparently. I close the issue

EDIT : except when I run my script, it doesn't work... at all (png = 15MB) Waiting for your insight to spend time on it again.

[PaulWessel]

Question. I am working on my density grid movie and find stuff like this in the movie call:

# The awkward -C settings was selected to give HD or UHD widths and
# even integer heights.
-C16.5517241379cx6.03448275862cx116

The idea here is that width * 116 gives integer pixels:

gmt math -Q 16.5517241379 116 MUL =
1920
gmt math -Q 6.03448275862 116 MUL =
700

It would be cleaner if I was able to give -C1920x700x116 to instead internally compute the dimensions 16.5517241379 cm times 6.03448275862 cm. How would we implement this? We dont know if the dots per unit is dpi or dpi. I think we need a modifier to -C, maybe

-C1920x700x116[+c|i]

where the modifier sets the units of the pixels per unit value 116. Did not make GMT 6.5 so if we use it we hopefully will have it in 6.6 or 6.5.1 before then. Is this OK with us as well as @seisman and @joa-quim ?

[Esteban82]

Had a hard time seeing much movements in the grid (just the spectrum lines). Not sure if it is the exciting ... what other things might you have?

Here (https://youtube.com/shorts/HKOIEnuea3M?si=hsSoxu68zr0ppWCm) I made another version with a different progression of the altitudes (not linear as in the other example). Now I think it is more exciting.

gmt math -T0/1040000/240+n -o1 T = tmp_z --FORMAT_FLOAT_OUT=%0.0f

[Esteban82]

... all *.png/mp4 are excluded by .gitignore.

This produce that also the file IndianaJones_Logo.png is excluded (which is used in the title slide). Paul, was this a mistake? Should I add the logo?

[PaulWessel]

Sorry, yes add back what is needed.

[Esteban82]

It would be cleaner if I was able to give -C1920x700x116 to instead internally compute the dimensions 16.5517241379 cm times 6.03448275862 cm. How would we implement this? We dont know if the dots per unit is dpi or dpi. I think we need a modifier to -C,

Since it is animation, I think it would be better to say the values in pixels (-C1920x700) with the pixel density (dpi/dpc).

[PaulWessel]

Today I have asked the Editor of G^3 if we live in 2024 and may have animation in an article. They better.

[Esteban82]

• avoid the whitespace between the topography and the "sea area" (coast not aligned with grdview)

I think it is not possible to do this with gmt coast because the data in GSHHG don't math the coastline of the relief grid. I am looking for an alternative (but I found an issue).

[PaulWessel]

Let me know if there are any scripts you wish me to look at for possible polishing or simplification. I have added a paragraph to the paper about my density model so I think I am close to being finished with text. May need to improve Conclusion and ensure we cit what should be cited in the References. Please describe your science animations and what it is about etc. Would be nice if we could submit this paper by mid-February, no?

[gd-a]

Well, you can have a look at the precipitation's one (WED-A_Fig_P) under gd-a branch

[PaulWessel]

I think the order of figures and movies are close to being set so please label your movies Fig_#.mp4 (I have indicated what # is in the title comments) and please use those names in the actual scripts aka well. We are not calling figure "IndianaJones .." etc since in articles they are all Fig*. Please insert the PNG version of your movies into the document and add a figure captions like I did.

[PaulWessel]

Well, you can have a look at the precipitation's one (WED-A_Fig_P) under gd-a branch

Great, weekend is coming!

[PaulWessel]

I made a few changes to WED-A_Fig_S.sh you can review;

1 You create a region that is 24x24 km centred on La Soufrière volcano, but 12 km N-S differ in degrees fro 12 km -E-W at the latitude of 16.043829. I did it exactly (div 12 km by cos latitude). You can approve it if that is what you mean.

Avoiding some long awk statements by storing lon and lat separately, then join them to a point for different uses.

You have a lot of commented out things (minimal) and I did not want to just delete it. Please clean up those section so the only code left is what is needed.dddd

[gd-a]

your movies Fig_#.mp4

So... Fig_7 for seismicity (in progress), Fig_8 for precipitation (waiting for review) and Fig_9 for sea-ice (to dust) ?

[PaulWessel]

Yes, carve out some paragraphs that sets the stage. unlikely that we will publish the scripts in the text - that should only happen for the smallest of cases. Place caption and ideally the PNG of the movies there for now.

[Esteban82]

I have the impressions that the manuscript might have too many animations: 2 from Paul. 3 from me (and maybe a 4th one). 2 (or 3) from Guillaume.

How many animations should we have? Maybe we could remove the simple Indiana Jones animation. What do you think?

I am asking because the manuscript could be too long and maybe hard to read.

A related issue is that there are fees for overlength articles in AGU (Excess Page Fee; https://www.agu.org/Publish-with-AGU/Publish/Author-Resources/Publication-fees ).

[PaulWessel]

Yes, thought the same. I think 2 is OK for me. I did think of reducing Indiana to just one somehow. Let us aim for 2 each. You could just show the few extra lines to do sound and mention the use of narration as well. Excess fee is a laugh - I got the money.

[Esteban82]

Ok, I think 2 for each is enough.

I will look how to merge both Indiana's animation.

[Esteban82]

Today I have asked the Editor of G^3 if we live in 2024 and may have animation in an article.

Paul, do you have any news about the journal?

Edit: I asked mainly to see the guidelines of the journal.

[PaulWessel]

No news yet. They are awfully slow.

[PaulWessel]

Good news from AGU/Wiley: https://authorservices.wiley.com/author-resources/Journal-Authors/Prepare/manuscript-preparation-guidelines.html/embedded-rich-media.html

But, we may need to apply som restrictions on file sizes and use more compression to pass their limits, Anyway, we should assume submission to G^3 as soon as we can, and organise the repo to have as much as possible to lighten the article load.

[Esteban82]

Great.

I am watching the restrictions

File specifications

Maximum file size 300 MB (note that the combined manuscript files for a submission, including video, audio, tables, figures, and text must not exceed 350 MB)

This means that the total size of all files should do not exceed 350 MB? Or is it for each video?

In any case, I think we should submit low resolutions animations. I suggest to post the high resolution animations on youtube and to add a link in the manuscript for each animation.

Also, we definetely must include only one Indiana Jones animation.

[Esteban82]

These are sizes of the followings videos: WED-A_Fig_4.mp4 248K WED-A_Fig_6.mp4 5.6M WED-A_Fig_7.mp4 4.3M WED-A_Fig_8.mp4 5.7M WED-A_Fig_9.mp4 11M Total: 27.9 MB

[Esteban82]

In any case, I think we should submit low resolutions animations. I suggest to post the high resolution animations on youtube and to add a link in the manuscript for each animation.

Also, we definetely must include only one Indiana Jones animation.

Sorry. I thought it said KB. Our movies are ok.

[PaulWessel]

It says no file can exceed 300 Mb and that the sum of all files must be <= 350 Mb.

[Esteban82]

It says no file can exceed 300 Mb

It says this but the system doesn't allow me to upload a file bigger than 100 Mb. I already send an email asking.