Difference between revisions of "easygrid"

From WikiROMS
Jump to navigationJump to search
Line 121: Line 121:
{{note}} '''NOTE:''' The default settings in EASYGRID are for the [[FJORD_TIDAL_CASE|Fjord tidal case]], hence you can just run EASYGRID (using the included <span class="red">Fjord_bathy.mat</span> and <span class="red">Fjord_coast.mat</span> files) to generate the grid and initial-conditions files required for the [[FJORD_TIDAL_CASE|Fjord tidal case]] tutorial.
{{note}} '''NOTE:''' The default settings in EASYGRID are for the [[FJORD_TIDAL_CASE|Fjord tidal case]], hence you can just run EASYGRID (using the included <span class="red">Fjord_bathy.mat</span> and <span class="red">Fjord_coast.mat</span> files) to generate the grid and initial-conditions files required for the [[FJORD_TIDAL_CASE|Fjord tidal case]] tutorial.


  <span style="color:green; font-size:95%; line-height:99%;">%*******************************************************************************************************************************
  <span style="color:green; font-size:95%; line-height:99%;">%****************************************************************************************************************************************************
  % USER SETTINGS ****************************************************************************************************************
  % USER SETTINGS *************************************************************************************************************************************
  %*******************************************************************************************************************************
  %****************************************************************************************************************************************************
%
% Switches -------------------------
      <font color="black">save_grid = 1;</font>  % Save GRID in a NetCDF file ( yes = 1; no = 0 )
      <font color="black">save_init = 1;</font>  % Create (and save) INITIAL CONDITIONS (from grid) ( yes = 1; no = 0 )
  <font color="black">screen_output = 1;</font>  % Displays (on the screen) many parameters that need to be copy-pasted in the ocean.in file ( yes = 1; no = 0 )
      <font color="black">plot_grid = 1;</font>  % Plot grid ( yes = 1; no = 0 )
    <font color="black">smooth_grid = 1;</font>  % Smooth bathymetry ( yes = 1; no = 0 ) using H. Arango's smth_bath.m , which applies a Shapiro filter to the bathymetry
  %
  %
% -------------------------------------------------------------------------
% SWITCHES ( ON = 1; OFF = 0 ) --------------------------------------------
% -------------------------------------------------------------------------
    <span class="black">create_grid = 1;</span>  % Create GRID. Turn OFF to work with a previously created grid (i.e. grid variables existing on Workspace)
      <span class="black">plot_grid = 1;</span>  % Plot grid
    <span class="black">smooth_grid = 1;</span>  % Smooth bathymetry using H. Arango's smth_bath.m , which applies a Shapiro filter to the bathymetry
      <span class="black">edit_mask = 0;</span>  % Edit rho-mask using interactive plot. Use this to manually change sea-pixels into land-pixels and vice-versa
  <span class="black">screen_output = 1;</span>  % Displays (on screen) many parameters that need to be copy-pasted in the ocean.in file   
      <span class="black">save_grid = 1;</span>  % Save GRID in a NetCDF file
      <span class="black">save_init = 1;</span>  % Create (and save) INITIAL CONDITIONS (from grid)
                        % ON = 1; OFF = 0
%   
  % -------------------------------------------------------------------------
  % -------------------------------------------------------------------------
  % GRID SETTINGS -----------------------------------------------------------
  % GRID SETTINGS -----------------------------------------------------------
Line 143: Line 148:
         <span class="black">Y = 2550;</span>            % Length of domain (meters)
         <span class="black">Y = 2550;</span>            % Length of domain (meters)
  <span class="black">rotangle = -43;</span>            % Angle (degrees) to rotate the grid conterclock-wise
  <span class="black">rotangle = -43;</span>            % Angle (degrees) to rotate the grid conterclock-wise
     <span class="black">resol = 30;</span>              % Cell width and height (i.e. Resolution)in meters. Grid cells are forced to be (almost) square.
     <span class="black">resol = 100;</span>              % Cell width and height (i.e. Resolution)in meters. Grid cells are forced to be (almost) square.
         <span class="black">N = 10;</span>              % Number of vertical levels
         <span class="black">N = 10;</span>              % Number of vertical levels
  %  
  %  
  %     
  %     
  % Bathymetry -------------- % Bathymetry for ROMS is measured positive downwards (zeros are not allowed) see: https://www.myroms.org/wiki/index.php/Grid_Generation#Bathymetry
  % Bathymetry -------------- % Bathymetry for ROMS is measured positive downwards (zeros are not allowed) see: https://www.myroms.org/wiki/index.php/Grid_Generation#Bathymetry
                             % To allow variations in surface elevation (eg. tides) while keeping all depths positive,
                             % To allow variations in surface elevation (eg. tides) while keeping all depths positive,
                             % an arbitrary offset (see minh below) is added to the depth vector.
                             % an arbitrary offset (see minh below) is added to the depth vector.
  %  
  %    
       <span class="black">hh = nan;</span>            % Analytical Depth (meters) used to create a uniform-depth grid. If using a bathymetry file, leave hh = nan;
       <span class="black">hh = nan;</span>            % Analytical Depth (meters) used to create a uniform-depth grid. If using a bathymetry file, leave hh = nan;
     <span class="black">minh = 4;</span>              % Arbitrary depth offset in meters (see above). minh should be a little more than the maximum expected tidal variation.
     <span class="black">minh = 4;</span>              % Arbitrary depth offset in meters (see above). minh should be a little more than the maximum expected tidal variation.
     <span class="black">Bathy =</span> <span class="purple">'Test_bathy.mat'</span><span class="black">;</span>% Bathymetry file. If using the analytical depth above (i.e. hh ~= nan), then Bathy will not be used.
     <span class="black">Bathy =</span> <span class="purple">'Fjord_bathy.mat'</span>;% Bathymetry file. If using the analytical depth above (i.e. hh ~= nan), then Bathy will not be used.
                             % The bathymetry file should be a .mat file containing 3 vectors (xbathy, ybathy and zbathy). where xbathy = Longitude,  
                             % The bathymetry file should be a .mat file containing 3 vectors (xbathy, ybathy and zbathy). where xbathy = Longitude,  
                             % ybathy = Latitude and zbathy = depth (positive, in meters). xbathy and ybathy are in decimal degrees See: http://en.wikipedia.org/wiki/Decimal_degrees
                             % ybathy = Latitude and zbathy = depth (positive, in meters). xbathy and ybathy are in decimal degrees See: http://en.wikipedia.org/wiki/Decimal_degrees
Line 166: Line 171:
  %
  %
  % Coastline ----------------------
  % Coastline ----------------------
     <span class="black">Coast = load(</span><span class="purple">'Test_coast.mat'</span><span class="black">);</span> % If there isn't a coastline file... comment-out this line: e.g. %Coast = load('test_coast.mat');
     <span class="black">Coast = load(</span><span class="purple">'Fjord_coast.mat'</span><span class="black">);</span> % If there isn't a coastline file... comment-out this line: e.g. %Coast = load('Fjord_coast.mat');
                                    % The coastline is only used for plotting. The coastline .mat file should contain 2 vectors named "lat" and "lon"
                                    % The coastline is only used for plotting. The coastline .mat file should contain 2 vectors named "lat" and "lon"
  %                                     
  %                                     
  %       
  %       
Line 173: Line 178:
  % OUTPUT: File naming and NetCDF descriptors ------------------------------
  % OUTPUT: File naming and NetCDF descriptors ------------------------------
  % -------------------------------------------------------------------------
  % -------------------------------------------------------------------------
  <span class="black">Grid_filename =</span> <span class="purple">'Fjord'</span>;     % Filename for grid and initial conditions files (don't include extension).  
  <span class="black">Grid_filename = <span class="purple">'Fjord'</span>;</span>   % Filename for grid and initial conditions files (don't include extension).  
                              % "_grd.nc" is added to grid file and "_ini.nc" is added to initial conditions file
                                % "_grd.nc" is added to grid file and "_ini.nc" is added to initial conditions file
  <span class="black">Descrip_grd  =</span> <span class="purple">'Test grid'</span>;              %Description for grid .nc file
  <span class="black">Descrip_grd  = <span class="purple">'Test grid'</span>;</span>               %Description for grid .nc file
  <span class="black">Descrip_ini  =</span> <span class="purple">'Test initial conditions'</span>; %Description for initial conditions .nc file
  <span class="black">Descrip_ini  = <span class="purple">'Test initial conditions'</span>;</span> %Description for initial conditions .nc file
  <span class="black">Author        =</span> <span class="purple">'John Smith'</span>;
  <span class="black">Author        = <span class="purple">'John Smith'</span>;</span>
  <span class="black">Computer      =</span> <span class="purple">'My Computer'</span>;
  <span class="black">Computer      = <span class="purple">'My Computer'</span>;</span>
  %
  %
  % -------------------------------------------------------------------------
  % -------------------------------------------------------------------------
Line 184: Line 189:
  % -------------------------------------------------------------------------
  % -------------------------------------------------------------------------
  <span class="blue">if</span> <span class="black">save_init == 1;</span> % See "Switches section" (above) to turn this ON or OFF
  <span class="blue">if</span> <span class="black">save_init == 1;</span> % See "Switches section" (above) to turn this ON or OFF
     % Initial conditions will be constant throught the grid domain
     % Initial conditions will be constant throughout the grid domain
     %--------------------------------------------------------------------------
     %--------------------------------------------------------------------------
     <span class="black">NH4          = 0.1;</span>    % Ammonium concentration (millimole_NH4 meter-3)
     <span class="black">NH4          = 0.1;</span>    % Ammonium concentration (millimole_NH4 meter-3)
Line 211: Line 216:
  %
  %
  %
  %
  %  
  %****************************************************************************************************************************************************
%*******************************************************************************************************************************
  % END OF USER SETTINGS ******************************************************************************************************************************
  % END OF USER SETTINGS *********************************************************************************************************
  %****************************************************************************************************************************************************</span>
  %*******************************************************************************************************************************</span>
<br>
<br>
<br>
<br>
===Geographical/Grid parameters===
===Geographical/Grid parameters===
After you got the bathymetry and coastline of your study region, now it is time to place a grid on it... In the <span class="green">USER SETTINGS</span> section, you have to specify <span class="red">lat lon X Y</span> and <span class="red">rotangle</span>. This is done in an iterative manner, where first you input your best-guesses, then you plot the grid. From the plot you can adjust your geographical/grid parameters... plot... adjust... plot... adjust... and so on.  
After you got the bathymetry and coastline of your study region, now it is time to place a grid on it... In the <span class="green">USER SETTINGS</span> section, you have to specify <span class="red">lat lon X Y</span> and <span class="red">rotangle</span>. This is done in an iterative manner, where first you input your best-guesses, then you plot the grid. From the plot you can adjust your geographical/grid parameters... plot... adjust... plot... adjust... and so on.  

Revision as of 01:04, 25 April 2008

EASYGRID
Warning This WikiROMS article is currently under HEAVY construction. This message will be erased when active construction is finished (in 1 or 2 days).
April 21, 2008

EASYGRID is a simple matlab script to make ROMS grids and initialization files. It is an austere grid-making alternative designed as a beginner's tool, yet capable of making grids for realistic applications. Since most of EASYGRID is in one file (easygrid.m), it is easier to follow and understand the steps required to generate a simple grid (i.e. useful as a learning or teaching tool) and also reduces the required steps to get the software up and running (i.e. downloading and compiling dependencies).


This WikiROMS page is mainly a tutorial, where I provide test bathymetry and coastline files so you can experiment with EASYGRID. Also, this page is where I will post news and updates.

QUICK LINKS:
Hosting page at Google-Code
Discussion Forum (link soon)
NEWS:
24 April 2008: Version 0 is out
Warning WARNING:
EASYGRID v1 is in BETA stage and could have problems. Report HERE (link soon) not only bugs, but also report successful use. Once the bug:success ratio decreases... I will remove this sign.



Download EASYGRID

  • Download EASYGRID (This download also includes, bathymetry, coastline and other files required for this tutorial)
  • Extract (i.e. unzip or unpack) the file in a location where it can stay indefinitely.
  • Add the location where you placed EASYGRID to your Matlab search path (see #5 as an example).
  • Dependencies: EASYGRID uses MEXNC, SNCTOOLS and ROMS Matlab tool-kit (only a few scripts like spheric_dist.m , smth_bath.m , pcolorjw.m).
Note For a tutorial on how to download/install these dependencies, CLICK HERE.



Getting bathymetry

EASYGRID needs a bathymetry .mat file containing 3 vectors (xbathy, ybathy and zbathy), where:

  • xbathy = longitude
  • ybathy = latitude
  • zbathy = depth (positive, in meters... For more information about bathymetry for ROMS, click HERE)
Note NOTE: Vectors xbathy and ybathy are in decimal degrees.


There easiest way to get bathymetry data (that I know) is to use Rich Signell's read_srtm30plus.m Matlab script. Instant bathymetry in 1 line of Matlab code! Click HERE for instructions to get and use read_srtm30plus.m.


Fjord latlon.jpg

The most difficult way to get bathymetry data (that I know) is to scan a chart and then digitize it using your mouse and a digitizing software like Didger (sorry, I am not aware of any good digitizing free software).


For this tutorial, you will need Fjord_bathy.mat . You probably already have this file since it is included in EASYGRID_v0.rar (file from the "Download" section). I created Fjord_bathy.mat by merging two bathymetry files... one high-resolution digitized from a chart and another coarser-resolution using read_srtm30plus.m (see Figure).


Note If you want to reproduce this plot as an exercise, copy-paste the following code on Matlab's Command Window:

load Fjord_bathy.mat
plot(xbathy,ybathy,'.');
xlabel('Longitude');
ylabel('Latitude');
axis square;



Getting coastline

In EASYGRID, the coastline is only used for plotting and it therefore not indispensable. However, it is handy to have a coastline for reference, especially when editing the land/sea mask. For EASYGRID, the coastline should be 2 vectors named "lat" and "lon" (both in units of decimal degrees) contained in a .mat file.


Fjord coast.jpg

An easy way to get coastline data is at NOAA's Coastline Extractor. Remember to select Matlab as format option (not default). SEAGRID's WikiROMS page has excellent information on how to use the Coastline Extractor, and tools to correct and manipulate coastline data.


For this tutorial, you will need Fjord_coast.mat . You probably already have this file since it is included in EASYGRID_v0.rar (file from the "Download" section). I created Fjord_coast.mat from the digitized chart I use to make the bathymetry file (see Figure).


Note If you want to reproduce this plot as an exercise, copy-paste the following code on Matlab's Command Window:

load Fjord_coast.mat
plot(lon,lat);
xlabel('Longitude');
ylabel('Latitude');
axis square;




Making your grid

You probably will need to run EASYGRID several times in a iterative manner to create, smooth and mask-edit your grid. You can turn ON or OFF different functionalities of EASYGRID (like the mask-editing and bathymetry-smoothing routines) using the SWITCHES located in the USER SETTINGS section, which is the first part of script (after the help lines). Below a snippet from EASYGRID's code that contains the SWITCHES section.

%****************************************************************************************************************************************************
% USER SETTINGS *************************************************************************************************************************************
%****************************************************************************************************************************************************
%
% -------------------------------------------------------------------------
% SWITCHES ( ON = 1; OFF = 0 ) --------------------------------------------
% -------------------------------------------------------------------------
   create_grid = 1;   % Create GRID. Turn OFF to work with a previously created grid (i.e. grid variables existing on Workspace)
     plot_grid = 1;   % Plot grid
   smooth_grid = 1;   % Smooth bathymetry using H. Arango's smth_bath.m , which applies a Shapiro filter to the bathymetry
     edit_mask = 0;   % Edit rho-mask using interactive plot. Use this to manually change sea-pixels into land-pixels and vice-versa 
 screen_output = 1;   % Displays (on screen) many parameters that need to be copy-pasted in the ocean.in file    
     save_grid = 1;   % Save GRID in a NetCDF file
     save_init = 1;   % Create (and save) INITIAL CONDITIONS (from grid)




A typical procedure to create a grid will follow (more or less) these steps:

  1. TURN ON create_grid and plot_grid and run EASYGRID.
  2. Tune-up Geographical and Grid parameters and run EASYGRID again.
  3. Repeat step 2 above until satisfied with the grid location, size and resolution.
  4. TURN OFF create_grid and proceed to smooth bathymetry by turning ON smooth_grid and re-running EASYGRID.
  5. Tune-up smoothing parameters and repeat the step above as necessary... when finished, turn OFF smooth_grid.
  6. Edit land/sea mask by turning ON edit_mask and re-running EASYGRID.
  7. Tune-up the land/sea mask as many times as necessary... when done, turn OFF edit_mask.
  8. When satisfied with the grid, bathymetry and mask... save it all to a GRID NetCDF file by turnning ON save_grid and re-running EASYGRID one more time.



Setting up USER SETTINGS

The first thing to do is to edit EASYGRID's USER SETTINGS. To make things easier for you, there are detailed instructions and explanations beside each variable entry line in the USER SETTINGS. In fact, most of the "user-manual" for EASYGRID is included within itself. Therefore I will simply copy-paste below a snippet from EASYGRID's code that contains the USER SETTINGS/instructions section.

Note NOTE: The default settings in EASYGRID are for the Fjord tidal case, hence you can just run EASYGRID (using the included Fjord_bathy.mat and Fjord_coast.mat files) to generate the grid and initial-conditions files required for the Fjord tidal case tutorial.

%****************************************************************************************************************************************************
% USER SETTINGS *************************************************************************************************************************************
%****************************************************************************************************************************************************
%
% -------------------------------------------------------------------------
% SWITCHES ( ON = 1; OFF = 0 ) --------------------------------------------
% -------------------------------------------------------------------------
    create_grid = 1;   % Create GRID. Turn OFF to work with a previously created grid (i.e. grid variables existing on Workspace)
      plot_grid = 1;   % Plot grid
    smooth_grid = 1;   % Smooth bathymetry using H. Arango's smth_bath.m , which applies a Shapiro filter to the bathymetry
      edit_mask = 0;   % Edit rho-mask using interactive plot. Use this to manually change sea-pixels into land-pixels and vice-versa 
  screen_output = 1;   % Displays (on screen) many parameters that need to be copy-pasted in the ocean.in file    
      save_grid = 1;   % Save GRID in a NetCDF file
      save_init = 1;   % Create (and save) INITIAL CONDITIONS (from grid)
                       % ON = 1; OFF = 0
%    
% -------------------------------------------------------------------------
% GRID SETTINGS -----------------------------------------------------------
% -------------------------------------------------------------------------
%
% Geographical and Grid parameters --------
     lat =  44.8125;        % Latitude  (degrees) of the bottom-left corner of the grid.
     lon = -62.8855;        % Longitude (degrees) of the bottom-left corner of the grid. 
%
       X = 9100;            % Width of domain (meters)
       Y = 2550;            % Length of domain (meters)
rotangle = -43;             % Angle (degrees) to rotate the grid conterclock-wise
   resol = 100;              % Cell width and height (i.e. Resolution)in meters. Grid cells are forced to be (almost) square.
       N = 10;              % Number of vertical levels
%    
%     
% Bathymetry -------------- % Bathymetry for ROMS is measured positive downwards (zeros are not allowed) see: https://www.myroms.org/wiki/index.php/Grid_Generation#Bathymetry
                            % To allow variations in surface elevation (eg. tides) while keeping all depths positive,
                            % an arbitrary offset (see minh below) is added to the depth vector.
%      
      hh = nan;             % Analytical Depth (meters) used to create a uniform-depth grid. If using a bathymetry file, leave hh = nan;
    minh = 4;               % Arbitrary depth offset in meters (see above). minh should be a little more than the maximum expected tidal variation.
   Bathy = 'Fjord_bathy.mat';% Bathymetry file. If using the analytical depth above (i.e. hh ~= nan), then Bathy will not be used.
                            % The bathymetry file should be a .mat file containing 3 vectors (xbathy, ybathy and zbathy). where xbathy = Longitude, 
                            % ybathy = Latitude and zbathy = depth (positive, in meters). xbathy and ybathy are in decimal degrees See: http://en.wikipedia.org/wiki/Decimal_degrees
       %Bathymetry smoothing routine...  See "Switches section" (above) to turn this ON or OFF
       if smooth_grid == 1;
           order = 2;       % Order of Shapiro filter (2,4,8)... default: 2
            rlim = 0.35;    % Maximum r-factor allowed (0.35)... default: 0.35
           npass = 50;      % Maximum number of passes.......... default: 50
       end
%---------------------------------
%                      
%
% Coastline ----------------------
   Coast = load('Fjord_coast.mat'); % If there isn't a coastline file... comment-out this line: e.g. %Coast = load('Fjord_coast.mat');
                                    % The coastline is only used for plotting. The coastline .mat file should contain 2 vectors named "lat" and "lon"
%                                    
%      
% -------------------------------------------------------------------------
% OUTPUT: File naming and NetCDF descriptors ------------------------------
% -------------------------------------------------------------------------
Grid_filename = 'Fjord'; 	   % Filename for grid and initial conditions files (don't include extension). 
                               % "_grd.nc" is added to grid file and "_ini.nc" is added to initial conditions file
Descrip_grd   = 'Test grid';               %Description for grid .nc file
Descrip_ini   = 'Test initial conditions'; %Description for initial conditions .nc file
Author        = 'John Smith';
Computer      = 'My Computer';
%
% -------------------------------------------------------------------------
% INITIAL CONDITIONS ------------------------------------------------------
% -------------------------------------------------------------------------
if save_init == 1; % See "Switches section" (above) to turn this ON or OFF
    % Initial conditions will be constant throughout the grid domain
    %--------------------------------------------------------------------------
    NH4          = 0.1;     % Ammonium concentration (millimole_NH4 meter-3)
    NO3          = 10;      % Nitrate concentration (millimole_N03 meter-3)
    chlorophyll1 = 0.3;     % Chlorophyll concentration in small phytoplankyon (milligrams_chlorophyll meter-3)
    chlorophyll2 = 0.3;     % Chlorophyll concentration in large phytoplankyon (milligrams_chlorophyll meter-3)
    detritus1    = 0.03;    % Small detritus concentration (millimole_nitrogen meter-3)
    detritus2    = 0.03;    % Large detritus concentration (millimole_nitrogen meter-3)
    detritusC1   = 1;       % Small detritus carbon concentration (millimole_carbon meter-3)
    detritusC2   = 0.2;     % Large detritus carbon concentration (millimole_carbon meter-3)
    phyto1       = 0.02;    % Small phytoplankton concentration (millimole_nitrogen meter-3)
    phyto2       = 0.02;    % Large phytoplankton concentration (millimole_nitrogen meter-3)
    phytoC1      = 0.2;     % Small phytoplankton carbon concentration (millimole_carbon meter-3)
    phytoC2      = 0.1;     % Small phytoplankton carbon concentration (millimole_carbon meter-3)
    salt         = 30;      % Salinity (PSU)
    temp         = 9;       % Potential temperature (Celsius)
    u            = 0;       % u-momentum component (meter second-1)
    ubar         = 0;       % Vertically integrated u-momentum component (meter second-1)
    v            = 0;       % v-momentum component (meter second-1)
    vbar         = 0;       % Vertically integrated v-momentum component (meter second-1)
    zeta         = 0;       % Free-surface (meters)
    zooplankton  = 0.01;    % Zooplankton concentration "millimole_nitrogen meter-3"
    zooplanktonC = 0.5;     % Zooplankton carbon concentration "millimole_carbon meter-3"
    %--------------------------------------------------------------------------
end
%
%
%****************************************************************************************************************************************************
% END OF USER SETTINGS ******************************************************************************************************************************
%****************************************************************************************************************************************************



Geographical/Grid parameters

After you got the bathymetry and coastline of your study region, now it is time to place a grid on it... In the USER SETTINGS section, you have to specify lat lon X Y and rotangle. This is done in an iterative manner, where first you input your best-guesses, then you plot the grid. From the plot you can adjust your geographical/grid parameters... plot... adjust... plot... adjust... and so on.

Note You may want to turn off some switches (see below) to speed up the plotting process.

% Switches -------------------------
     save_grid = 0;   % Save GRID in a NetCDF file ( yes = 1; no = 0 )
     save_init = 0;   % Create (and save) INITIAL CONDITIONS (from grid) ( yes = 1; no = 0 )
 screen_output = 0;   % Displays (on the screen) many parameters that need to be copy-pasted in the ocean.in file ( yes = 1; no = 0 )
     plot_grid = 1;   % Plot grid ( yes = 1; no = 0 )
   smooth_grid = 0;   % Smooth bathymetry ( yes = 1; no = 0 ) using H. Arango's smth_bath.m , which applies a Shapiro filter to the bathymetry


Fjord NOTsmooth.jpg
This is the resulting grid (without smoothing) if you run EASYGRID with the default settings.
NoteIf you want to see the 4 different grids that ROMS uses, zoom-in (in the grid plot), and then copy-paste the following in Matlab's Command Window:
plot(lon_rho,lat_rho,'r.')
plot(lon_u,lat_u,'bv')
plot(lon_v,lat_v,'g^')
plot(lon_psi,lat_psi,'m+')
  • The rho-grid represents the grid centers (red dots)
  • The u-grid represents the grid East-West sides (blue triangles)
  • The v-grid represents the grid North-South sides (green triangles)
  • The psi-grid represents the grid corners (purple crosses)
Fjord grid.jpg




Bathymetry smoothing

ROMS doesn't like rough bathymetry with abrupt changes in topography. Therefore, after you got your grid where you wanted it to be, the next step is to smooth the bathymetry. Turn ON the smooth_grid variable in the %Switches section. Usually simply enabling smoothing (with the default settings) will do a pretty good job. But if unsatisfied, you can iteratively tweak the parameters in the %Bathymetry somoothing section of the USER SETTINGS until obtaining the desired smoothness. Below is a comparison of grids without and with smoothing (default settings).

Fjord NOTsmoothNsmooth.jpg




Initial Conditions

ROMS initialization NetCDF files contain information of the spatial distribution of one (or multiple) variable(s) that ROMS uses at the START of the simulation. The value of each initialization variable is specified for each corresponding cell in the ROMS grid... hence, inside a typical initial-conditions file, there will be one three-dimensional matrix (matching grid dimensions) for each specified variable. EASYGRID can generate Initial-Conditions NetCDF files for ROMS, however, only in a rudimentary manner... the user can only choose a single value for each initial-condition variable... EASYGRID makes each variable CONSTANT over the entire grid domain. This may be only useful in some simple cases and to do "quick-tests".


To change the value of a initial-condition variable, simply go to the INITIAL CONDITIONS section in the USER SETTINGS and change the value of the variable of interest.


If you need to add another initial-condition variable... (1) add in the INITIAL CONDITIONS section of the USER SETTINGS:

new_variable = 0; % Variable description (units)


and (2) add the following snippet to the end of the easygrid.m file (substitute the parts in red)

%--------------------------------------------------------------------------
       dims                    = { 'time'; 's_rho'; 'eta_rho'; 'xi_rho'};
       nc{ 'newvar'}           = ncdouble(dims);
       nc{ 'newvar'}(:,:,:)    = ones(N,Mp,Lp).* 'new_variable';
       nc{ 'newvar'}.time      = 'ocean_time' ;


Warning WARNING: The dimensions of the new initial-conditions variable should be the same as the grid that will use it. For example, the u momentum initial-condition matrix should have the same dimensions as the u-grid. In a biological setting however, most of the newly-specified variables are likely to be used by the rho-grid.

Parameter output on screen

If screen_output is turned ON (in the USER SETTINGS section), EASYGRID will output some parameters to screen. Make sure you write this numbers (Lm, MM, N and DT) since you will need to enter them in the ocean.in file prior to running ROMS.

|  COPY-PASTE the following parameters into your ocean.in file
|  ----------------------------------------------------------------------------------------------
|
|
|     Lm == 90         ! Number of I-direction INTERIOR RHO-points
|     Mm == 24         ! Number of J-direction INTERIOR RHO-points
|      N == 10         ! Number of vertical levels
|
|
|  Make sure the Baroclinic time-step (DT) in your ocean.in file is less than: 3.3882 seconds
|  ----------------------------------------------------------------------------------------------
  • The parameters Lm, Mm and N are the basically the dimensions of the RHO-grid.
  • The parameter DT is the baroclinic time-step and it should be less than SQRT( (dx^2+dy^2)/(g*h(x,y) ). EASYGRID solves this equation and tells you the result (3.3882 sec in the case above). You can read a bit more about this in this forum post (middle of page).



Saving .nc files

Once you got your grid in place, your bathymetry smoothed and your initial-conditions adjusted... it is time to save the NetCDF files for ROMS. To do this, simply turn on ALL the Switches in the USER SETTING section and run EASYGRID. Besides the plot, this time it will also save the grid and initial conditions .nc files. EASYGRID will also output some parameters to screen. WarningMake sure you write down this numbers (Lm, MM, N and DT), since you will need to enter them in the ocean.in file prior to running ROMS.



Editing Masks

After you created your grid.nc file, you will likely need to edit the mask. That is, you will need to change some land-pixes of the mask into sea-pixels... and vice versa. To edit your grid's mask, you can use editmask.m , which is a GUI script included in the ROMS Matlab tool-kit (inside the rmask directory).

Warning The current version of editmask.m (from the ROMS toolkit) doesn't work for me. I included in the EASYGRID_v1.rar file (that you probably already downloaded) a revised version of editmask.m that works for me.



Fjord editmask1.jpg

Here are step-by-step instructions to use editmask.m:

  • Start the GUI by typing in Matlab's Command Window:
editmask
  • Select the NetCDF grid file to edit (Fjord_grd.nc in this tutorial)
  • Select the .mat coastline file that corresponds to the grid file (Fjord_coast.mat in this tutorial)
  • Select "Set Land" ("Edit Mode" section on the right) and mask 1-cell bays like the examples circled in A and B.
  • Proceed to mask disconnected lakes, like one in example circled in E.
  • Although tiny islands apparently don't cause troubles... you may want to select "Set Sea" and unmask "islands" that you consider unrealistic or unnecessary, like the ones circled by C and D.
  • Once you finished editing the mask, CLICK SAVE. This will update your NetCDF grid file for ROMS.
  • Close the GUI.







Note This concludes this tutorial. However, you should consider reviewing the Tidal Fjord TUTORIAL, which uses the grid and initial-condition files generated here to setup a simple tidal realistic application.