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 |
---|