#!/bin/bash
#Script to perform a mean on a series of directories using cdo's 
# v1 (Feb2021)-Enver Ramirez
# v2 (Oct2026)-Rosio Camayo

dir_orig=/p/projetos/monan_goc/emanuel.giarolla/monan006/outputs🏁😎/history/original
dir_out=/p/projetos/monan_goc/rosio.camayo/mom6/mean_exp006_2005-2023
#mkdir -p tmp

rm -f ${dir_out}/temp_saz_ucur_*.nc ${dir_out}/temp_saz_vcur_*.nc ${dir_out}/mean_ucur_sazonal.nc ${dir_out}/mean_vcur_sazonal.nc 

for year in `seq 2005 2023`
do
 for month in `seq -f "%02g" 6 8 `
# for month in `seq -f "%02g" 12 -1 9
 do
    cdo selvar,u ${dir_orig}/ucur_${year}_${month}.nc ${dir_out}/temp_saz_ucur_${year}_${month}.nc
    cdo selvar,v ${dir_orig}/vcur_${year}_${month}.nc ${dir_out}/temp_saz_vcur_${year}_${month}.nc
#  echo ${dir_orig}/ucur_${year}_${month}.nc
 done
done
#done >> lista
# y=`cat lista`
 echo cdo ensmean ${dir_out}/temp_saz_ucur_*.nc ${dir_out}/mean_ucur_sazonal.nc
 echo cdo ensmean ${dir_out}/temp_saz_vcur_*.nc ${dir_out}/mean_vcur_sazonal.nc
 cdo ensmean ${dir_out}/temp_saz_ucur_*.nc ${dir_out}/mean_ucur_sazonal.nc
 cdo ensmean ${dir_out}/temp_saz_vcur_*.nc ${dir_out}/mean_vcur_sazonal.nc
#rm -f ${dir_out}/temp_saz_ucur_*.nc ${dir_out}/temp_saz_vcur_*.nc 
# rm lista