ncep2gcm.x
ncep2gcm.x — text/x-csh, 3 ko (3583 bytes)
Contenu du fichier
#!/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