Horizontal Grid (CPP coordinate)

Implementing SCHISM model as part of NOAA's Integrated Water Modeling Projects

Overview

In addition to hgrid.ll (lon/lat), some pre-processing scripts need the horizontal grid on a projected coordinate system. The hgrid.cpp described here is based on CPP projection (see a description of the coordinate system here).


Scripts

A Fortran script is provide in Git to perform the projection from lon/lat to CPP:

$schism_git_dir/src/Utility/Pre-Processing/NWM/Grid_manipulation/cpp.f90

, which requires the user input of the origin of the system in lon/lat. This is also provided in:

$schism_git_dir/src/Utility/Pre-Processing/NWM/Grid_manipulation/cpp.in

To generate hgrid.cpp from an existing hgrid.ll, you can "cp -rL" the "Grid_manipulation" folder to your run dir; "cd Grid_manipulation"; then do "./auto.pl" (which in turn runs the fortran script as "./cpp < cpp.in" and do some housekeeping. Ignore the warning messages on hgrid.utm.*, which is not used for the 2D setup.

The "Grid_manipulation" folder also contains small utilities that reproduce some frequently used functions of the Acetool.