jueves, 3 de julio de 2025

Generar video con OpenSWPC

 cd ~/Documents/OpenSWPC/test/

mpirun -np 2  /home/mariajovera/Documents/OpenSWPC/bin/swpc_psv.x -i input.inf
 

cd /home/mariajovera/Documents/OpenSWPC/test/out 


/home/mariajovera/Documents/OpenSWPC/bin/read_snp.x -i swpc.psv.xz.v.nc -ppm -pall -mul 10

o en kathleen es:

#------------------------------------------------------------------
module load hdf/5-1.8.15/intel-2015-update2
module load netcdf/4.3.3.1/intel-2015-update2
module load netcdf-fortran/4.4.1/intel-2015-update2 

file="Mref_B_full_L2b_S1.xz.v.nc";  

/home/ucfbmve/OpenSWPC-5.3.1/bin/read_snp.x -i  Mref_B_full_L2b_S1.xz.v.nc   -ppm -pall -mul 10

 

cd /home/mariajovera/Documents/OpenSWPC/test/out/swpc/psv/xz/v2


ffmpeg -i swpc.psv.xz.v2.%6d.ppm -qscale 0 -pix_fmt yuv420p -y swpc.mp4



Plotear un fotograma en especifico con GMT


# Explorar las variables del archivo .nc 
 ncdump -h swpc.psv.xz.v.nc

# Extraer un fotograma en un tiempo específico con ncks 
 ncks -d t,100 swpc.psv.xz.v.nc vz_t100.nc
 


 
# Esto crea un nuevo archivo NetCDF vz_t100.nc que contiene:
# Vz(z, x) porque ya se extrajo el tiempo.

1. Crear paleta de colores
gmt makecpt -Cturbo -T-100/100/1 > turbo.cpt


2. Graficar con GMT

gmt grdimage vz_t100.nc?Vz -JX15c/10c -B -Cturbo.cpt -P > vz_t100.ps

psconvert vz_t100.ps -Tg -A -E300

No hay comentarios:

Publicar un comentario