source: TOOLS/MOSAIX/Utils/CreateORCA025grid.bash @ 3665

Last change on this file since 3665 was 3665, checked in by omamce, 6 years ago

O.M. : mise a jour des properties
  • Property svn:keywords set to Date Revision HeadURL Author Id
File size: 5.9 KB
Line 
1#!/bin/bash
2set -e
3set -x
4
5
6InFile=/ccc/work/cont003/igcmg/igcmg/IGCM/OCE/NEMO/ORCA025_LIM3_PISCES/v3.6_stable/mesh_mask_eORCA025.nc
7
8
9ncks -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
10ncrename -h -v glamt,nav_lon_grid_T -v gphit,nav_lat_grid_T ORCA025_coordinates_mask.nc
11ncrename -h -v glamu,nav_lon_grid_U -v gphiu,nav_lat_grid_U ORCA025_coordinates_mask.nc
12ncrename -h -v glamv,nav_lon_grid_V -v gphiv,nav_lat_grid_V ORCA025_coordinates_mask.nc
13ncrename -h -v glamf,nav_lon_grid_F -v gphif,nav_lat_grid_F ORCA025_coordinates_mask.nc
14
15ncatted -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
57ncwa -O -h -a t ORCA025_coordinates_mask.nc tmp_ORCA025_coordinates_mask.nc ; mv tmp_ORCA025_coordinates_mask.nc ORCA025_coordinates_mask.nc
58
59cat <<EOF > add_mask.nco
60area_grid_T [y,x] = e1t (:,:) * e2t (:,:)    ;
61area_grid_U [y,x] = e1u (:,:) * e2u (:,:)    ;
62area_grid_V [y,x] = e1v (:,:) * e2v (:,:)    ;
63area_grid_F [y,x] = e1f (:,:) * e2f (:,:)    ;
64mask_T      [y,x] = tmask     (0,:,:) * 1.0d ;
65mask_U      [y,x] = umask     (0,:,:) * 1.0d ;
66mask_V      [y,x] = vmask     (0,:,:) * 1.0d ;
67mask_F      [y,x] = fmask     (0,:,:) * 1.0d ;
68maskutil_T  [y,x] = tmaskutil   (:,:) * 1.0d ;
69maskutil_U  [y,x] = umaskutil   (:,:) * 1.0d ;
70maskutil_V  [y,x] = fmaskutil   (:,:) * 1.0d ;
71maskutil_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 ;
80EOF
81
82ncap2 -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
84ncks -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
89ncatted -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
Note: See TracBrowser for help on using the repository browser.