1 | #!/bin/bash |
---|
2 | set -e |
---|
3 | set -x |
---|
4 | |
---|
5 | |
---|
6 | InFile=/ccc/work/cont003/igcmg/igcmg/IGCM/OCE/NEMO/ORCA025_LIM3_PISCES/v3.6_stable/mesh_mask_eORCA025.nc |
---|
7 | |
---|
8 | |
---|
9 | ncks -4 -h -O -v glamt,glamu,glamv,glamf,gphit,gphiu,gphiv,gphif,tmask,umask,vmask,fmask,tmaskutil,umaskutil,vmaskutil,fmaskutil,e1t,e2t,e1u,e1f,e2u,e1v,e2v,e2f,e2f ${InFile} ORCA025_coordinates_mask.nc |
---|
10 | ncrename -h -v glamt,nav_lon_grid_T -v gphit,nav_lat_grid_T ORCA025_coordinates_mask.nc |
---|
11 | ncrename -h -v glamu,nav_lon_grid_U -v gphiu,nav_lat_grid_U ORCA025_coordinates_mask.nc |
---|
12 | ncrename -h -v glamv,nav_lon_grid_V -v gphiv,nav_lat_grid_V ORCA025_coordinates_mask.nc |
---|
13 | ncrename -h -v glamf,nav_lon_grid_F -v gphif,nav_lat_grid_F ORCA025_coordinates_mask.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 units,nav_lon_grid_T,c,c,"degrees_east" \ |
---|
32 | -a standard_name,nav_lon_grid_T,c,c,"Longitude" \ |
---|
33 | -a long_name,nav_lon_grid_T,c,c,"Longitude" \ |
---|
34 | -a units,nav_lat_grid_T,c,c,"degrees_north" \ |
---|
35 | -a standard_name,nav_lat_grid_T,c,c,"Latitude" \ |
---|
36 | -a long_name,nav_lat_grid_T,c,c,"Latitude" \ |
---|
37 | -a units,nav_lon_grid_U,c,c,"degrees_east" \ |
---|
38 | -a standard_name,nav_lon_grid_U,c,c,"Longitude" \ |
---|
39 | -a long_name,nav_lon_grid_U,c,c,"Longitude" \ |
---|
40 | -a units,nav_lat_grid_U,c,c,"degrees_north" \ |
---|
41 | -a standard_name,nav_lat_grid_U,c,c,"Latitude" \ |
---|
42 | -a long_name,nav_lat_grid_U,c,c,"Latitude" \ |
---|
43 | -a units,nav_lon_grid_V,c,c,"degrees_east" \ |
---|
44 | -a standard_name,nav_lon_grid_V,c,c,"Longitude" \ |
---|
45 | -a long_name,nav_lon_grid_V,c,c,"Longitude" \ |
---|
46 | -a units,nav_lat_grid_V,c,c,"degrees_north" \ |
---|
47 | -a standard_name,nav_lat_grid_V,c,c,"Latitude" \ |
---|
48 | -a long_name,nav_lat_grid_V,c,c,"Latitude" \ |
---|
49 | -a units,nav_lon_grid_F,c,c,"degrees_east" \ |
---|
50 | -a standard_name,nav_lon_grid_F,c,c,"Longitude" \ |
---|
51 | -a long_name,nav_lon_grid_F,c,c,"Longitude" \ |
---|
52 | -a units,nav_lat_grid_F,c,c,"degrees_north" \ |
---|
53 | -a standard_name,nav_lat_grid_F,c,c,"Latitude" \ |
---|
54 | -a long_name,nav_lat_grid_F,c,c,"Latitude" \ |
---|
55 | ORCA025_coordinates_mask.nc |
---|
56 | |
---|
57 | ncwa -O -h -a t ORCA025_coordinates_mask.nc tmp_ORCA025_coordinates_mask.nc ; mv tmp_ORCA025_coordinates_mask.nc ORCA025_coordinates_mask.nc |
---|
58 | |
---|
59 | cat <<EOF > add_mask.nco |
---|
60 | area_grid_T [y,x] = e1t (:,:) * e2t (:,:) ; |
---|
61 | area_grid_U [y,x] = e1u (:,:) * e2u (:,:) ; |
---|
62 | area_grid_V [y,x] = e1v (:,:) * e2v (:,:) ; |
---|
63 | area_grid_F [y,x] = e1f (:,:) * e2f (:,:) ; |
---|
64 | mask_T [y,x] = tmask (0,:,:) * 1.0d ; |
---|
65 | mask_U [y,x] = umask (0,:,:) * 1.0d ; |
---|
66 | mask_V [y,x] = vmask (0,:,:) * 1.0d ; |
---|
67 | mask_F [y,x] = fmask (0,:,:) * 1.0d ; |
---|
68 | maskutil_T [y,x] = tmaskutil (:,:) * 1.0d ; |
---|
69 | maskutil_U [y,x] = umaskutil (:,:) * 1.0d ; |
---|
70 | maskutil_V [y,x] = fmaskutil (:,:) * 1.0d ; |
---|
71 | maskutil_F [y,x] = fmaskutil (:,:) * 1.0d ; |
---|
72 | //maskutil_T [:, 0] = 0.0d ; |
---|
73 | //maskutil_T [:, 1441] = 0.0d ; |
---|
74 | //maskutil_U [:, 0] = 0.0d ; |
---|
75 | //maskutil_U [:, 1441] = 0.0d ; |
---|
76 | //maskutil_V [:, 0] = 0.0d ; |
---|
77 | //maskutil_V [:, 1441] = 0.0d ; |
---|
78 | //maskutil_F [:, 0 ] = 0.0d ; |
---|
79 | //maskutil_F [:, 1441] = 0.0d ; |
---|
80 | EOF |
---|
81 | |
---|
82 | ncap2 -h -S add_mask.nco ORCA025_coordinates_mask.nc tmp_ORCA025_coordinates_mask.nc ; mv tmp_ORCA025_coordinates_mask.nc ORCA025_coordinates_mask.nc |
---|
83 | |
---|
84 | ncks -h -x --variable e1t,e2t,e1u,e2u,e1v,e2v,e1f,e2f,tmask,umask,vmask,fmask,tmaskutil,umaskutil,vmaskutil,fmaskutil ORCA025_coordinates_mask.nc tmp_ORCA025_coordinates_mask.nc ; |
---|
85 | |
---|
86 | #python creates_nemo_bounds.py |
---|
87 | #mv tmp_ORCA025_coordinates_mask.nc ORCA025_coordinates_mask.nc |
---|
88 | |
---|
89 | ncatted -h \ |
---|
90 | -a coordinates,area_grid_T,c,c,"nav_lat_grid_T nav_lon_grid_T" \ |
---|
91 | -a coordinates,area_grid_U,c,c,"nav_lat_grid_U nav_lon_grid_U" \ |
---|
92 | -a coordinates,area_grid_V,c,c,"nav_lat_grid_V nav_lon_grid_V" \ |
---|
93 | -a coordinates,area_grid_F,c,c,"nav_lat_grid_F nav_lon_grid_F" \ |
---|
94 | -a units,area_grid_T,c,c,"m^2" \ |
---|
95 | -a units,area_grid_U,c,c,"m^2" \ |
---|
96 | -a units,area_grid_V,c,c,"m^2" \ |
---|
97 | -a units,area_grid_F,c,c,"m^2" \ |
---|
98 | -a coordinates,mask_T,c,c,"nav_lat_grid_T nav_lon_grid_T" \ |
---|
99 | -a coordinates,mask_U,c,c,"nav_lat_grid_U nav_lon_grid_U" \ |
---|
100 | -a coordinates,mask_V,c,c,"nav_lat_grid_V nav_lon_grid_V" \ |
---|
101 | -a coordinates,mask_F,c,c,"nav_lat_grid_F nav_lon_grid_F" \ |
---|
102 | -a coordinates,maskutil_T,c,c,"nav_lat_grid_T nav_lon_grid_T" \ |
---|
103 | -a coordinates,maskutil_U,c,c,"nav_lat_grid_U nav_lon_grid_U" \ |
---|
104 | -a coordinates,maskutil_V,c,c,"nav_lat_grid_V nav_lon_grid_V" \ |
---|
105 | -a coordinates,maskutil_F,c,c,"nav_lat_grid_F nav_lon_grid_F" \ |
---|
106 | -a cell_measures,mask_T,c,c,"area: area_grid_T" \ |
---|
107 | -a cell_measures,mask_U,c,c,"area: area_grid_U" \ |
---|
108 | -a cell_measures,mask_V,c,c,"area: area_grid_V" \ |
---|
109 | -a cell_measures,mask_F,c,c,"area: area_grid_F" \ |
---|
110 | -a cell_measures,maskutil_T,c,c,"area: area_grid_T" \ |
---|
111 | -a cell_measures,maskutil_U,c,c,"area: area_grid_U" \ |
---|
112 | -a cell_measures,maskutil_V,c,c,"area: area_grid_V" \ |
---|
113 | -a cell_measures,maskutil_F,c,c,"area: area_grid_F" \ |
---|
114 | -a bounds,nav_lon_grid_T,c,c,"bounds_lon_grid_T" \ |
---|
115 | -a bounds,nav_lat_grid_T,c,c,"bounds_lat_grid_T" \ |
---|
116 | -a bounds,nav_lon_grid_U,c,c,"bounds_lon_grid_U" \ |
---|
117 | -a bounds,nav_lat_grid_U,c,c,"bounds_lat_grid_U" \ |
---|
118 | -a bounds,nav_lon_grid_V,c,c,"bounds_lon_grid_V" \ |
---|
119 | -a bounds,nav_lat_grid_V,c,c,"bounds_lat_grid_V" \ |
---|
120 | -a bounds,nav_lon_grid_F,c,c,"bounds_lon_grid_F" \ |
---|
121 | -a bounds,nav_lat_grid_F,c,c,"bounds_lat_grid_F" \ |
---|
122 | ORCA025_coordinates_mask.nc |
---|
123 | |
---|