# User defined
diret='.'
outpath='./SeaIce'
year_ini=2000
year_end=2008
# end: User defined


source $MODULESHOME/init/sh
module load cdo/1.8.2_cmor2
#module load cdo

# Escolher os anos inicial e final
#set -x
for i in `seq $year_ini $year_end`
do
  for x in `ls $diret/ice_month_${i}_??.nc`; do

# Adicionar informações do grid_BESM_cmip6a_T.nc ao arquivo que se quer interpolar
  cp $x tmp.nc
  ncatted -h -a coordinates,,c,c,"GEOLON GEOLAT" tmp.nc

## *** interpolando 1 (n48 = Gaussian 192x96)
#  cdo remapbil,n48 tmp.nc ${outpath}/T062_${x##*/}

## *** interpolando 2 (grade regular 0.5 x 0.5)
##        - Definindo a grade de saída
#cat > mygrid.txt << EOF
#gridtype = lonlat
#xsize    = 720
#ysize    = 360
#xfirst   = -279.5 
#xinc     = 0.5
#yfirst   = -89.5
#yinc     = 0.5
#EOF
#
## Finalmente interpolando
#  cdo remapbil,mygrid.txt tmp.nc ${outpath}/remap720x360_${x##*/}
#rm mygrid.txt

# Grade especial
echo  cdo --no_warnings remapbil,GridIce.txt tmp.nc ${outpath}/remap_${x##*/}
      cdo --no_warnings remapbil,GridIce.txt tmp.nc ${outpath}/remap_${x##*/}
  ncatted -O -a calendar,Time,o,c,"julian" ${outpath}/remap_${x##*/}

rm tmp.nc 
  done
done