[3621] | 1 | #!/bin/bash |
---|
| 2 | |
---|
| 3 | |
---|
| 4 | ### =========================================================================== |
---|
| 5 | ## |
---|
| 6 | ## Default configuration |
---|
| 7 | ## =========================================================================== |
---|
| 8 | |
---|
| 9 | ### =========================================================================== |
---|
| 10 | ### |
---|
| 11 | ### You should not change anything below this line .... |
---|
| 12 | ### |
---|
| 13 | ### =========================================================================== |
---|
| 14 | |
---|
| 15 | ## |
---|
| 16 | ## Command line parameters |
---|
| 17 | ## =========================================================================== |
---|
| 18 | set -o verbose |
---|
| 19 | set -o xtrace |
---|
| 20 | set -e |
---|
| 21 | |
---|
| 22 | |
---|
| 23 | while [[ ${1} = -* ]] ; do |
---|
| 24 | case ${1} in |
---|
| 25 | ( -- ) shift ; break ;; |
---|
| 26 | ( -v | --verbose ) set -o verbose ;; |
---|
| 27 | ( -x | --xtrace ) set -o xtrace ;; |
---|
| 28 | ( -e ) set -e ;; |
---|
| 29 | ( -V | --noverbose ) set +o verbose ;; |
---|
| 30 | ( -X | --noxtrace ) set +o xtrace ;; |
---|
| 31 | ( -E ) set +e ;; |
---|
| 32 | ( -* ) echo ${Bold}"Unknown option : ${1}"${Norm} ; return 1 ;; |
---|
| 33 | esac |
---|
| 34 | shift |
---|
| 35 | done |
---|
| 36 | |
---|
| 37 | ## |
---|
| 38 | AtmFile=${1:-/ccc/store/cont003/dsm/p25sepul/IGCM_OUT/IPSLCM5A2/PROD/piControl/CM5A2.1.pi.00/ATM/Output/MO/CM5A2.1.pi.00_40100101_40191231_1M_histmth.nc} |
---|
| 39 | |
---|
| 40 | dim_lat=$(ncdump -h ${AtmFile} | grep "lat *=" | awk '{print $3}' ) |
---|
| 41 | dim_lon=$(ncdump -h ${AtmFile} | grep "lon *=" | awk '{print $3}' ) |
---|
| 42 | (( dim_lat = dim_lat - 1 )) |
---|
| 43 | |
---|
| 44 | GridFile=lmdz${dim_lon}x${dim_lat}_grid.nc |
---|
| 45 | |
---|
| 46 | ncks --overwrite -h -d time_counter,0,1 -v lon,lat,aire,fract_oce,fract_sic ${AtmFile} ${GridFile} |
---|
| 47 | #ncwa --overwrite -h -a time_counter ${GridFile} tmp_${GridFile} ; mv tmp_${GridFile} ${GridFile} |
---|
| 48 | |
---|
| 49 | ncatted -h -a coordinates,aire,m,c,"lat lon" ${GridFile} |
---|
| 50 | ncatted -h -a coordinates,fract_oce,m,c,"lat lon" ${GridFile} |
---|
| 51 | ncatted -h -a coordinates,fract_sic,m,c,"lat lon" ${GridFile} |
---|
| 52 | |
---|
| 53 | ncks --overwrite -h --xcl -v time_counter,time_centered,time_counter_bounds ${GridFile} tmp_${GridFile} ; mv tmp_${GridFile} ${GridFile} |
---|
| 54 | |
---|
| 55 | cat <<EOF >add_fields.nco |
---|
| 56 | time_centered [time_counter] = { 43200.0, 129600.0 } ; |
---|
| 57 | time_centered_bounds [time_counter, axis_nbounds] = { 0.0, 86400.0, 86400.0, 172800.0 } ; |
---|
| 58 | fract_oce_plus_sic = fract_oce + fract_sic ; |
---|
| 59 | where (fract_oce_plus_sic < 0.00001 ) fract_oce_plus_sic=0.0 ; |
---|
| 60 | where (fract_oce_plus_sic > 0.99999 ) fract_oce_plus_sic=1.0 ; |
---|
| 61 | // Fill masked values to land values |
---|
| 62 | OceMask[time_counter,lat,lon] = 0.0 ; |
---|
| 63 | where (fract_oce_plus_sic > 0.0 ) OceMask=1.0 ; |
---|
| 64 | EOF |
---|
| 65 | |
---|
| 66 | ncap2 -h --overwrite --script-file add_fields.nco ${GridFile} tmp_${GridFile} ; mv tmp_${GridFile} ${GridFile} |
---|
| 67 | |
---|
| 68 | ncatted -h \ |
---|
| 69 | -a coordinates,OceMask,c,c,"time_centered lat lon" \ |
---|
| 70 | -a standard_name,time_centered,c,c,"time" \ |
---|
| 71 | -a long_name,time_centered,c,c,"Time axis" \ |
---|
| 72 | -a calendar,time_centered,c,c,"360_day" \ |
---|
| 73 | -a units,time_centered,c,c,"seconds since 2000-01-01 00:00:00" \ |
---|
| 74 | -a time_origin,time_centered,c,c,"2000-01-01 00:00:00" \ |
---|
| 75 | -a bounds,time_centered,c,c,"time_centered_bounds" \ |
---|
| 76 | dynamico450_grid.nc |
---|
| 77 | |
---|
| 78 | ncks -O --mk_rec time_counter ${GridFile} tmp_${GridFile} ; mv tmp_${GridFile} ${GridFile} |
---|
| 79 | |
---|
| 80 | ncatted -h \ |
---|
| 81 | -a long_name,aire,c,c,"Grid area" \ |
---|
| 82 | -a units,aire,c,c,"m^2" \ |
---|
| 83 | -a online_operation,aire,c,c,"once" \ |
---|
| 84 | -a coordinates,aire,c,c,"lat lon" \ |
---|
| 85 | ${GridFile} |
---|
| 86 | |
---|
| 87 | ### =========================================================================== |
---|
| 88 | rm -f add_fields.nco |
---|
| 89 | |
---|
| 90 | |
---|
| 91 | ### =========================================================================== |
---|
| 92 | ### |
---|
| 93 | ### That's all folk's !!! |
---|
| 94 | ### |
---|
| 95 | ### =========================================================================== |
---|
| 96 | |
---|