1. General, what is in this folder.
Routines:
  Wind_driven_sea_level.m : main routine.
  cmocean.m               : function of colormap for plotting.
  get_baroclinic_RWH.m    : function to calculate Rossby wave height driven by wind stress curls and/or SSH at eastern boundary.
  get_wind_stress_curl.m  : function to calculate wind stress curls from wind stress.
  latlon2xy.m             : function to calculate distance between two grids, note the sign of the distance is defined positively eastward in x and northward in y.
  ra_windstr.m            : function to calculate wind stress from wind speeds.
Data:
  uwnd.10m.mon.mean.nc
  vwnd.10m.mon.mean.nc    : NECP 20CRv3, monthly mean u- & v-wind speeds at 10 meter height.
  ETOPO1_Ice_g_gmt4.grd.nc: Elevation data of ETOPO1 for creating the ocean mask.
  SODA_NA_Qiang_from_Soenker.mat: SODA SSH in North Atlantic regions.


2. Sections in the main routine
Section 1: Calculate the wind stress curls. It also setups the spatial grids (regular 1-degree resolutions as NCEP), and ocean mask based on ETOPO1.
Section 2: Get SSH data (from SODA). The SODA SSH data is downscaled into regular 1-degree grids.
Section 3: Calculate the Rossby wave height.
           Setup your calculating domains and forcing at the eastern boundary.
             Line 196-197: define the domain of  RWH outputs
               lon_min (the western boundary)  lon_max (the eastern boundary)
		       lat_min (the southmost boundary) lat_max (the northmost boundary)
		     Line 202: define the time vector for RWH outputs
		       currently from 1850 to 2013, monthly
		     Line 208: define the eastern boundary of forcing