Custom Query (964 matches)

I have several Matlab grid manipulation scripts that may be useful to some Users.

• cartesian_grid.m: Given a larger donor grid (D), this function extracts a rotated Cartesian grid centered at (Xc, Yc) and rotated by angle theta. The larger donor grid is used during plotting to display the location of the rotated grid.

• spherical_grid.m: Given a larger donor grid (D), this function extracts a rotated Spherical grid centered at (Tlon, Tlat) and rotated by angle alpha.

• grid_rotated.m: Given a larger Grid NetCDF file (Ginp), this function extracts and creates a rotated subdomain (Cartesian or Spherical) Grid NetCDF file (Gout) centered at (Xmid,Ymid) and rotated by angle theta.

• roms_metrics.m: Computes the several ROMS grid metrics from the (lon,lat) coordinates.

I highly recommend using roms_metrics when building the ROMS GRID NetCDF file to have consistent grid variables. The function grid_metrics is deprecated. Use the robust roms_metrics instead. The grid_metrics function will be removed in the future.

Notice that the rotated grids are experimental. They were designed to have rotated nested grids. Such capability is not available in the ROMS numerical kernel.

This update includes changes to few routines associated with nesting. The changes are primarily due to land/sea masking (MASKING and WET_DRY). Many thanks to John Warner for the corrections.

• Updated nesting.F, main3d.F, u2dbc_im.F, and v3dbc_im.F. Also, updated the TLM, RPM, and ADM versions.

• Corrected adjoint bugs in ad_u2dbc_im.F in implementing the western edge Shchepetkin boundary condition formulation.

• Updated the documentation and reporting of 4D-Var input parameters in s4dvar.in and read_asspar.F, respectively, for the background quality control of input observations (BGQC). The threshold values P_bgqc and S_bgqc are the squared number of standard deviations used in the quality control for rejecting observations. The information was incomplete, missing the fact that the user needs to enter a squared number value. Many thanks to Andy Moore for bringing this to my attention.

• Corrected typos in v2dbc_im.F. We need to have the isVbar parameter instead of isUbar in the LBC conditinal statement. This bug was introduced ​src:ticket:880. Similar correction were required for tl_v2dbc_im.F, rp_v2dbc_im.F, and ad_v2dbc_im.F. Many thanks to John Warner for bringing this issue to my attention.

• Updated module interpolate.F to include an interpolation object for the ROMS-JEDI interface.