#!/bin/bash

umask 022

export PATH=/cray_home/local/python/python_3.7/bin:${PATH}
export PYTHONPATH=/cray_home/local/python/python_3.7/bin
export GDAL_DATA=/cray_home/local/python/python_3.7/share/gdal

export GRADSB=${MHOME}/produtos/WEB/opengrads-2.2.1.oga.1/Contents
export GADDIR=/cray_home/local/grads/grads-2.1.0.oga.1/bin/Resources/SupportData
export GADLIB=/cray_home/local/grads/grads-2.1.0.oga.1/bin/Linux/Versions/2.1.0.oga.1/x86_64/libs
export GASCRP="${MHOME}/produtos/WEB/opengrads-2.2.1.oga.1/Contents/Resources/Scripts ${MHOME}/produtos/WEB/scripts/grads"
export PATH=${GRADSB}/ImageMagick-7.0.8-45/utilities/:${GRADSB}:$PATH

module switch PrgEnv-cray PrgEnv-gnu

export LD_LIBRARY_PATH=/usr/lib/:${LD_LIBRARY_PATH}
export LIBPATH=/usr/lib/
export MAGICK_CONFIGURE_PATH=${GRADSB}/ImageMagick-7.0.8-45/config

i=0
for arq in $(ls /lustre_xc50/ioper/data/BRAMS/ams_15km/2022090800G/model/POSPROCESS/profile*ctl /lustre_xc50/ioper/data/BRAMS/ams_15km/2022090900/model/POSPROCESS/profile*ctl ); do
cat <<EOF> aot_anim.gs
'open ${arq}'
'set display color white'
'c'

'set grads off'

'q time'
say result
fct=subwrd(result, 3)

'set lat -23.267'
'set lon -47.408'
'set z 1 32'
*'set xaxis 0 350'
'set xint 20'
*'set vrange 0 350'

'd pm25'
'draw title BRAMS 15km CPTEC/INPE Mat. Part. (dp<2.5um)(ug/m3)\'fct'(Inicializacao: 00Z08SEP2022)'
'printim perfil_mp25_${i}.png'
'quit'

EOF

grads -bpc aot_anim.gs
let i=$i+1
sleep 01

done

magick -limit memory 1 -limit map 1 $(ls -ltr aot_550_*.png | awk '{print $9}') ./final_pm25_perfil_animation.gif
exit 0