[3621] | 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 | |
---|