| 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061 |
- #!/bin/bash
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- # This script interpolates the erainterim forcing fields from the 1.5*1.5 grid
- # to the ORCA1 grid.
- #
- # History : Virginie Guemas - Initial version - 2012
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- # Arguments
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- year0=2012
- yearf=2014
- lstvars=( 't2' 'q2' 'precip' 'snow' 'qsw' 'qlw' )
- lstmembs=( 0 )
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- set -evx
- cp /cfu/pub/scripts/interpolation/scrip_use_extrap .
- for ((year=${year0};year<=${yearf};year++)) ; do
- for memb in ${lstmembs[@]} ; do
- for var in ${lstvars[@]} ; do
- case $var in
- 'u10'|'v10') mid='_fc$'${memb}'_' ;;
- *) mid='_' ;;
- esac
- cp /cfu/scratch/vguemas/ERAint/${var}${mid}eraint_${year}.nc ${var}_eraint_${year}.nc
- case $var in
- 'u10'|'v10') ncpdq -O -U ${var}_eraint_${year}.nc ${var}_eraint_${year}.nc ;;
- esac
- case $var in
- 't2'|'q2') nstep=1460 ; method='bicubic' ; filewgts='/cfu/pub/scripts/interpolation/weigths/rmp_ERAint_ocean_to_ORCA1t_v3.2_bicubic.nc' ;;
- 'precip'|'snow') nstep=12; method='bicubic' ; filewgts='/cfu/pub/scripts/interpolation/weigths/rmp_ERAint_ocean_to_ORCA1t_v3.2_bicubic.nc' ;;
- 'qsw'|'qlw') nstep=365; method='bicubic' ; filewgts='/cfu/pub/scripts/interpolation/weigths/rmp_ERAint_ocean_to_ORCA1t_v3.2_bicubic.nc' ;;
- esac
-
- cat > scrip_use_in <<EOF
- &remap_inputs
- remap_wgt = '${filewgts}'
- infile = '${var}_eraint_${year}.nc'
- invertlat = TRUE
- var = '${var}'
- fromregular = TRUE
- outfile = '${var}_eraint_${year}_orca1.nc'
- /
- EOF
- ./scrip_use_extrap
- rm -f ${var}_eraint_${year}.nc
- ncks -O -v ${var} ${var}_eraint_${year}_orca1.nc ${var}_eraint_${year}_orca1.nc
- done
- done
- done
- rm -f scrip_use_in scrip_use_extrap
|
