Fix rotation of input winds: Cassini to earth relative
commit1da64c73937c56cb9f0fdfe4f1b03152ebbba934
authorDave Gill <gill@ucar.edu>
Fri, 1 Dec 2017 19:10:22 +0000 (1 12:10 -0700)
committerMichael Duda <duda@ucar.edu>
Tue, 25 Sep 2018 21:29:21 +0000 (25 15:29 -0600)
tree7e5a965ad098d4acde830c786e3107e6205f833c
parent3ab4522e777aefbf4c4b4cb0938849f695a04d5d
Fix rotation of input winds: Cassini to earth relative

Modify the existing met_to_map routine to use the correct rotation
angle, alpha. Previously, the rotation angle was computed with the lat,lon
for the d(lat)/d(y) and d(lon)/d(y) from the interpolated data on the
WRF domain. The correct lat and lon for these d/dy values is to use the
lat lon from the original input data.

The d(lat)/d(y) and d(lon)/d(y) computations use the existing latlon_to_ij
and the inverse ij_to_latlon in the module_map_utils.F file. This module is
now required to be "used" in rotate_winds_module.F, and the associated
Makefile for metgrid is also modified to include this new dependency.
metgrid/src/Makefile
metgrid/src/rotate_winds_module.F