#!/bin/bash

grid_file="ocean_static_0001_03.nc"

ncks -v geolat,geolon ocean_static_0001_03.nc grid_geolon_geolat.nc

#for infile in $(ls ocean_sfc_ts_*.nc); do
for infile in $(ls ocean_sfc_ts_2230_04.nc); do
outfile=exp010_remap720x360_tic_`basename $infile .nc`.nc

echo $infile $outfile

rm tmp.nc

cp ${infile} tmp.nc

ncks -A $grid_file tmp.nc 

ncatted -h -a coordinates,,c,c,"geolon geolat" tmp.nc

# * 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 d-remap_Rosio/${outfile}

done