#!/bin/csh

set tmin=$1
set tmax=$2
set GRILLE_DIR=$3

# ncep2gcm.x 48 32 1 4 $MAIN_DIR/NCEP/$date $NCEP_DIR/$date

#####################################################################
#####################################################################

set GRILLE_FI=$GRILLE_DIR/grilles_gcm.nc


if ( ! -f $GRILLE_FI ) then
echo creer le fichier $GRILLE_FI avec grilles_gcm_netcdf.e
exit
endif

echo OK0 $GRILLE_FI
set iip1=`ncdump -h $GRILLE_FI | grep lonu | head -1 | awk ' { print $3 } '`
set jjm=`ncdump -h $GRILLE_FI | grep latv | head -1 | awk ' { print $3 } '`
echo OK1
set jjp1=`expr $jjm + 1`
set iim=`expr $iip1 - 1`

echo iim=$iim
echo jjm=$jjm

\rm ta.nc ua.nc va.nc sst.nc u.nc v.nc T.nc ts.nc
echo LS AVANT LE DEBUT

###################################################################3
# Choix de la periode temporelle
###################################################################3
set t0="l=$tmin"

set t1tn="l=${tmin}:${tmax}"


if ( 1 == 0 ) then

set varu=ua
set varv=va
set vart=ta
set varq=hur
set varp=ps

set fileu="$NCEP_DIR/ua.nc"
set filev="$NCEP_DIR/va.nc"
set filet="$NCEP_DIR/ta.nc"
set fileq="$NCEP_DIR/hur.nc"
set filep="$NCEP_DIR/ps.nc"


else # ERA40

set ANA_DIR='http://dods-ipsl.lmd.polytechnique.fr/ipsl-bin/nph-dods/ERA40/NETCDF/4xdaily/'
set varu=u
set varv=v
set vart=T
set varq=r
set varp=msl

set an=1999
set mois=01
set fileu="$ANA_DIR/AN_PL/$an/$varu.$an$mois.aph.nc"
set filev="$ANA_DIR/AN_PL/$an/$varv.$an$mois.aph.nc"
#set filet="$ANA_DIR/AN_PL/$an/$vart.$an$mois.aph.nc"
set filet="$ANA_DIR/AN_PL/$an/t.$an$mois.aph.nc"
set fileq="$ANA_DIR/AN_PL/$an/$varq.$an$mois.aph.nc"
set filep="$ANA_DIR/AN_SF/$an/$varp.$an$mois.ash.nc"

endif

#####################################################################
#####################################################################
echo OK0
###################################################################3
# Lien avec les cihiers netcdf contenant les donnees NCEP
###################################################################3


###################################################################3
# script ferret pour interpolation
###################################################################3

cat <<eod> ! tmp.jnl
set memory/size=50
use "$GRILLE_DIR/grilles_gcm.nc"
use "$fileu"
use "$filev"
use "$filet"
use "$fileq"
use "$filep"

let uwnd=$varu
let vwnd=$varv
let air='$vart'
let rh=$varq
let sp=$varp

define axis/t=1-jan-1994:31-dec-1994:6/units=hours thour
define grid/like=uwnd[d=2]/x=grille_u[d=1]/y=grille_u[d=1] grille_u
define grid/like=vwnd[d=3]/x=grille_v[d=1]/y=grille_v[d=1] grille_v
define grid/like=air[d=4]/x=grille_v[d=1]/y=grille_u[d=1] grille_T

save/clobber/file=T.nc air[d=4,g=grille_T,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
repeat/$t1tn save/file="T.nc"/append air[d=4,g=grille_T,i=1:$iip1,j=1:$jjp1,gt=thour@asn]

save/clobber/file=u.nc uwnd[d=2,g=grille_u,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
repeat/$t1tn save/file="u.nc"/append uwnd[d=2,g=grille_u,i=1:$iip1,j=1:$jjp1,gt=thour@asn]

save/clobber/file=v.nc vwnd[d=3,g=grille_v,i=1:$iip1,j=1:$jjm,$t0,gt=thour@asn]
repeat/$t1tn save/file="v.nc"/append vwnd[d=3,g=grille_v,i=1:$iip1,j=1:$jjm,gt=thour@asn]

save/clobber/file=hur.nc rh[d=5,g=grille_T,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
repeat/$t1tn save/file="hur.nc"/append rh[d=5,g=grille_T,i=1:$iip1,j=1:$jjp1,gt=thour@asn]

save/clobber/file=ps.nc sp[d=6,g=grille_T,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
repeat/$t1tn save/file="ps.nc"/append sp[d=6,g=grille_T,i=1:$iip1,j=1:$jjp1,gt=thour@asn]


eod

ferret <<eod
go tmp.jnl
quit
eod
