#!/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_ucur_*.nc ${dir_out}/temp_vcur_*.nc ${dir_out}/mean_ucur_allperiod.nc ${dir_out}/mean_vcur_allperiod.nc 

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