| 1 | #!/bin/bash |
|---|
| 2 | set -e |
|---|
| 3 | set -x |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | InDir=/ccc/work/cont003/gencmip6/bedidil/InterpolToy/DATA/GRIDS |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | ncks -4 -h -O -v icoh.lon,icoh.lat,icoh.clo,icoh.cla ${InDir}/grids.nc dynamico450_grid.nc |
|---|
| 10 | ncrename -h -v icoh.lon,lon dynamico450_grid.nc |
|---|
| 11 | ncrename -h -v icoh.lat,lat dynamico450_grid.nc |
|---|
| 12 | ncrename -h -v icoh.clo,bounds_lon dynamico450_grid.nc |
|---|
| 13 | ncrename -h -v icoh.cla,bounds_lat dynamico450_grid.nc |
|---|
| 14 | |
|---|
| 15 | ncatted -h\ |
|---|
| 16 | -a history,global,d,c,c,"" \ |
|---|
| 17 | -a Comment,global,d,c,c,"" \ |
|---|
| 18 | -a source,global,d,c,c,"" \ |
|---|
| 19 | -a group,global,d,c,c,"" \ |
|---|
| 20 | -a production,global,d,c,c,"" \ |
|---|
| 21 | -a cr,global,d,c,c,"" \ |
|---|
| 22 | -a grid,global,d,c,c,"" \ |
|---|
| 23 | -a Institution,global,d,c,c,"" \ |
|---|
| 24 | -a Model,global,d,c,c,"" \ |
|---|
| 25 | -a HOSTNAME,global,d,c,c,"" \ |
|---|
| 26 | -a LOGNAME,global,d,c,c,"" \ |
|---|
| 27 | -a name,global,d,c,c,"" \ |
|---|
| 28 | -a title,global,d,c,c,"" \ |
|---|
| 29 | -a timeStamp,global,d,c,c,"" \ |
|---|
| 30 | -a uuid,global,d,c,c,"" \ |
|---|
| 31 | -a bounds,lon,c,c,"bounds_lon" \ |
|---|
| 32 | -a bounds,lat,c,c,"bounds_lat" \ |
|---|
| 33 | -a units,lon,c,c,"degrees_east" \ |
|---|
| 34 | -a standard_name,lon,c,c,"Longitude" \ |
|---|
| 35 | -a long_name,lon,c,c,"Longitude" \ |
|---|
| 36 | -a units,lat,c,c,"degrees_north" \ |
|---|
| 37 | -a standard_name,lat,c,c,"Latitude" \ |
|---|
| 38 | -a long_name,lat,c,c,"Latitude" \ |
|---|
| 39 | dynamico450_grid.nc |
|---|
| 40 | |
|---|
| 41 | ncwa -h -a y_icoh dynamico450_grid.nc tmp_dynamico450_grid.nc ; mv tmp_dynamico450_grid.nc dynamico450_grid.nc |
|---|
| 42 | ncatted -h \ |
|---|
| 43 | -a cell_methods,lon,d,c,"" \ |
|---|
| 44 | -a cell_methods,lat,d,c,"" \ |
|---|
| 45 | -a cell_methods,bounds_lon,d,c,"" \ |
|---|
| 46 | -a cell_methods,bounds_lat,d,c,"" \ |
|---|
| 47 | dynamico450_grid.nc |
|---|
| 48 | |
|---|
| 49 | ncrename -h -d x_icoh,cell dynamico450_grid.nc |
|---|
| 50 | ncrename -h -d crn_icoh,nvertex dynamico450_grid.nc |
|---|
| 51 | |
|---|
| 52 | ncpdq -h --permute cell,nvertex dynamico450_grid.nc tmp_dynamico450_grid.nc ; mv tmp_dynamico450_grid.nc dynamico450_grid.nc |
|---|
| 53 | |
|---|
| 54 | |
|---|
| 55 | cat <<EOF > add_mask.nco |
|---|
| 56 | defdim("time_counter", 2) ; |
|---|
| 57 | defdim("axis_nbounds", 2) ; |
|---|
| 58 | time_centered [time_counter] = { 43200.0, 129600.0 } ; |
|---|
| 59 | time_centered_bounds [time_counter, axis_nbounds] = { 0.0, 86400.0, 86400.0, 172800.0 } ; |
|---|
| 60 | aire [cell] = ( 3.14159*6400000)^2/2016012.0 ; |
|---|
| 61 | OceMask [time_counter,cell] = 0.0d ; |
|---|
| 62 | EOF |
|---|
| 63 | |
|---|
| 64 | ncap2 -h -S add_mask.nco dynamico450_grid.nc tmp_dynamico450_grid.nc ; mv tmp_dynamico450_grid.nc dynamico450_grid.nc |
|---|
| 65 | |
|---|
| 66 | ncatted -h \ |
|---|
| 67 | -a coordinates,OceMask,c,c,"time_centered lat lon" \ |
|---|
| 68 | -a standard_name,time_centered,c,c,"time" \ |
|---|
| 69 | -a long_name,time_centered,c,c,"Time axis" \ |
|---|
| 70 | -a calendar,time_centered,c,c,"360_day" \ |
|---|
| 71 | -a units,time_centered,c,c,"seconds since 2000-01-01 00:00:00" \ |
|---|
| 72 | -a time_origin,time_centered,c,c,"2000-01-01 00:00:00" \ |
|---|
| 73 | -a bounds,time_centered,c,c,"time_centered_bounds" \ |
|---|
| 74 | dynamico450_grid.nc |
|---|
| 75 | |
|---|
| 76 | ncks -O --mk_rec time_counter dynamico450_grid.nc tmp_dynamico450_grid.nc ; mv tmp_dynamico450_grid.nc dynamico450_grid.nc |
|---|
| 77 | |
|---|
| 78 | ncatted -h \ |
|---|
| 79 | -a long_name,aire,c,c,"Grid area" \ |
|---|
| 80 | -a units,aire,c,c,"m^2" \ |
|---|
| 81 | -a online_operation,aire,c,c,"once" \ |
|---|
| 82 | -a coordinates,aire,c,c,"lat lon" \ |
|---|
| 83 | dynamico450_grid.nc |
|---|
