| 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758 |
- set -evx
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- # This script interpolates and rotates the 10m zonal and meridional wind
- # perturbations.
- #
- # History : Virginie Guemas - Initial version - 2012
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- # Arguments
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- year='2014' # Year for which to interpolate and rotate
- dir='perturb1_era-dfs_mon' # Directory from where to take the perturbations
- # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
- for ((jmemb=1;jmemb<=4;jmemb++)) ; do
- cat > scrip_use_in <<EOF
- &remap_inputs
- remap_wgt = '/cfu/pub/scripts/interpolation/weights/rmp_ERAint_to_ORCA1t_unmasked_bicubic.nc'
- infile = '${dir}/u10_fc${jmemb}_${year}.nc'
- invertlat = TRUE
- var = 'u10'
- fromregular = TRUE
- outfile = '${dir}/u10_fc${jmemb}_eraint_${year}_orca1t_unrotated.nc'
- /
- EOF
- cp /cfu/pub/scripts/interpolation/scrip_use scrip_use
- ./scrip_use
- cat > scrip_use_in <<EOF
- &remap_inputs
- remap_wgt = '/cfu/pub/scripts/interpolation/weights/rmp_ERAint_to_ORCA1t_unmasked_bicubic.nc'
- infile = '${dir}/v10_fc${jmemb}_${year}.nc'
- invertlat = TRUE
- var = 'v10'
- fromregular = TRUE
- outfile = '${dir}/v10_fc${jmemb}_eraint_${year}_orca1t_unrotated.nc'
- /
- EOF
- ./scrip_use
- cat > namelist_rotateUVorca <<EOF
- &nam_rotUV
- Ufilein = '${dir}/u10_fc${jmemb}_eraint_${year}_orca1t_unrotated.nc'
- Uvarin = 'u10'
- Vfilein = '${dir}/v10_fc${jmemb}_eraint_${year}_orca1t_unrotated.nc'
- Vvarin = 'v10'
- meshmask = '/cfu/autosubmit/con_files/mesh_mask_nemo.N3.2_O1L42.nc'
- Ufileout = '${dir}/u10_fc${jmemb}_eraint_${year}_orca1t.nc'
- Vfileout = '${dir}/v10_fc${jmemb}_eraint_${year}_orca1t.nc'
- /
- EOF
- cp /cfu/pub/scripts/interpolation/rotateUVorca_orca1 rotateUVorca
- ./rotateUVorca
- rm -f ${dir}/u10_fc${jmemb}_eraint_${year}_orca1t_unrotated.nc ${dir}/v10_fc${jmemb}_eraint_${year}_orca1t_unrotated.nc
- done
- rm -f scrip_use scrip_use_in rotateUVorca namelist_rotateUVorca
|
