| 1 |
% |
| 2 |
% %Id: $ |
| 3 |
% |
| 4 |
% Ed Hill |
| 5 |
% |
| 6 |
% Create a netCDF file defining a global Lat-Lon grid for SCRIP. |
| 7 |
|
| 8 |
clear all |
| 9 |
close all |
| 10 |
|
| 11 |
debug_lev = 0; |
| 12 |
|
| 13 |
% |
| 14 |
% JMC wants: 128x64 |
| 15 |
% |
| 16 |
nlat = 64; |
| 17 |
nlon = 128; |
| 18 |
latcell = 180/nlat; |
| 19 |
loncell = 360/nlon; |
| 20 |
lat = linspace(-90+latcell/2, 90-latcell/2, nlat); |
| 21 |
% |
| 22 |
% Note: the APE grid needs to be [0, MAX_EAST] |
| 23 |
% lon = linspace(-180+loncell/2, 180-loncell/2, nlon); |
| 24 |
lon = linspace(0+loncell/2, 360-loncell/2, nlon); |
| 25 |
|
| 26 |
dlat_t = diff(lat); |
| 27 |
dlon_t = diff(lon); |
| 28 |
dlat = 0.5*( [dlat_t(1) ; dlat_t(:)] + ... |
| 29 |
[dlat_t(:) ; dlat_t(length(dlat_t))] ); |
| 30 |
dlon = 0.5*( [dlon_t(1) ; dlon_t(:)] + ... |
| 31 |
[dlon_t(:) ; dlon_t(length(dlon_t))] ); |
| 32 |
clear dlat_t dlon_t |
| 33 |
|
| 34 |
grid_size = nlat*nlon; |
| 35 |
grid_rank = 2; |
| 36 |
grid_corners = 4; |
| 37 |
|
| 38 |
! rm -f remap_grid_LL128x64.cdf |
| 39 |
|
| 40 |
nc = netcdf(['remap_grid_LL128x64.cdf'], 'clobber'); |
| 41 |
nc.title = 'Global 2-degree Lat-Lon grid'; |
| 42 |
nc.author = 'Ed Hill'; |
| 43 |
nc('grid_size') = grid_size; |
| 44 |
nc('grid_rank') = grid_rank; |
| 45 |
nc('grid_corners') = grid_corners; |
| 46 |
|
| 47 |
nc{'grid_dims'} = ncint('grid_rank'); |
| 48 |
nc{'grid_dims'}(:) = [ length(lon) length(lat) ]; |
| 49 |
|
| 50 |
nc{'grid_imask'} = ncint('grid_size'); |
| 51 |
nc{'grid_imask'}(:) = ones(length(grid_size),1); |
| 52 |
|
| 53 |
nc{'grid_center_lat'} = ncdouble('grid_size'); |
| 54 |
nc{'grid_center_lon'} = ncdouble('grid_size'); |
| 55 |
nc{'grid_center_lat'}.units = 'degrees'; |
| 56 |
nc{'grid_center_lon'}.units = 'degrees'; |
| 57 |
nc{'grid_corner_lat'} = ncdouble('grid_size', 'grid_corners'); |
| 58 |
nc{'grid_corner_lon'} = ncdouble('grid_size', 'grid_corners'); |
| 59 |
nc{'grid_corner_lat'}.units = 'degrees'; |
| 60 |
nc{'grid_corner_lon'}.units = 'degrees'; |
| 61 |
|
| 62 |
cen_lat = zeros(grid_size,1); |
| 63 |
cen_lon = zeros(grid_size,1); |
| 64 |
cor_lat = zeros(grid_size,4); |
| 65 |
cor_lon = zeros(grid_size,4); |
| 66 |
j = 0; |
| 67 |
for ilon = 1:length(lon) |
| 68 |
for ilat = 1:length(lat) |
| 69 |
j = j + 1; |
| 70 |
latc = lat(ilat); |
| 71 |
lonc = lon(ilon); |
| 72 |
cen_lat(j) = latc; |
| 73 |
cen_lon(j) = lonc; |
| 74 |
% Note: grid points *MUST* be ordered CCW and here |
| 75 |
% we start from the lower-left corner |
| 76 |
da = dlat(ilat)/2.0; |
| 77 |
do = dlon(ilon)/2.0; |
| 78 |
cor_lat(j,:) = [ latc-da latc-da latc+da latc+da ]; |
| 79 |
cor_lon(j,:) = [ lonc-do lonc+do lonc+do lonc-do ]; |
| 80 |
% if debug_lev > 1 |
| 81 |
% pause(1) |
| 82 |
% hold on |
| 83 |
% plot([cor_lon(j,:) cor_lon(j,1)],... |
| 84 |
% [cor_lat(j,:) cor_lat(j,1)],'o-g') |
| 85 |
% plot(cen_lon(j),cen_lat(j),'xr') |
| 86 |
% hold off |
| 87 |
% end |
| 88 |
end |
| 89 |
end |
| 90 |
nc{'grid_center_lat'}(:) = cen_lat; |
| 91 |
nc{'grid_center_lon'}(:) = cen_lon; |
| 92 |
nc{'grid_corner_lat'}(:) = cor_lat; |
| 93 |
nc{'grid_corner_lon'}(:) = cor_lon; |
| 94 |
|
| 95 |
close(nc) |
| 96 |
|
Parent Directory
| 
Revision Log
| 
Revision Graph