Module containing utilites used by the orog program. More...

Functions/Subroutines
subroutine, public	find_nearest_pole_points (i_north_pole, j_north_pole, i_south_pole, j_south_pole, im, jm, is_north_pole, is_south_pole)
	Find the point on the model grid tile closest to the north and south pole.

subroutine, public	find_poles (geolat, nx, ny, i_north_pole, j_north_pole, i_south_pole, j_south_pole)
	Find the point on the model grid tile closest to the north and south pole.

subroutine, public	get_index (imn, jmn, npts, lono, lato, delxn, jst, jen, ilist, numx)
	Determine the location of a cubed-sphere point within the high-resolution orography data.

real function, public	get_lat_angle (dy)
	Convert the 'y' direction distance of a cubed-sphere grid point to the corresponding distance in latitude.

real function, public	get_lon_angle (dx, lat_in)
	Convert the 'x' direction distance of a cubed-sphere grid point to the corresponding distance in longitude.

real function, public	get_xnsum (lon1, lat1, lon2, lat2, imn, jmn, glat, zavg, zslm, delxn)
	Count the number of high-resolution orography points that are higher than the model grid box average orography height.

subroutine, public	get_xnsum2 (lon1, lat1, lon2, lat2, imn, jmn, glat, zavg, delxn, xnsum1, xnsum2, hc)
	Count the number of high-resolution orography points that are higher than a critical value inside a model grid box (or a portion of a model grid box).

subroutine, public	get_xnsum3 (lon1, lat1, lon2, lat2, imn, jmn, glat, zavg, delxn, xnsum1, xnsum2, hc)
	Count the number of high-resolution orography points that are higher than a critical value inside a model grid box (or a portion of a model grid box).

logical function, public	inside_a_polygon (lon1, lat1, npts, lon2, lat2)
	Check if a point is inside a polygon.

subroutine, public	latlon2xyz (siz, lon, lat, x, y, z)
	Convert from latitude and longitude to x,y,z coordinates.

subroutine, public	minmax (im, jm, a, title, imax, jmax)
	Print out the maximum and minimum values of an array and optionally pass back the i/j location of the maximum.

subroutine, public	remove_isolated_pts (im, jm, slm, oro, var, var4, oa, ol)
	Remove isolated model points.

real function	spherical_angle (v1, v2, v3)
	Compute spherical angle.

real function, public	timef ()
	Get the date/time from the system clock.

subroutine, public	transpose_mask (imn, jmn, mask)
	Transpose the global landmask by flipping the poles and moving the starting longitude to Greenwich.

subroutine, public	transpose_orog (imn, jmn, glob)
	Transpose the global orography data by flipping the poles and moving the starting longitude to Greenwich.

Variables
real, parameter	deg2rad = 3.14159265358979/180.
	degrees per radians.

real, parameter	earth_radius = 6371200.
	earth radius in meters.

real, parameter	pi =3.1415926535897931
	pi.

real, parameter	rad2deg = 180./3.14159265358979
	radians per degrees.

Detailed Description

Module containing utilites used by the orog program.

Author: George Gayno NOAA/EMC

Function/Subroutine Documentation

◆ find_nearest_pole_points()

subroutine, public orog_utils::find_nearest_pole_points	(	integer, intent(in)	i_north_pole,
		integer, intent(in)	j_north_pole,
		integer, intent(in)	i_south_pole,
		integer, intent(in)	j_south_pole,
		integer, intent(in)	im,
		integer, intent(in)	jm,
		logical, dimension(im,jm), intent(out)	is_north_pole,
		logical, dimension(im,jm), intent(out)	is_south_pole )

Find the point on the model grid tile closest to the north and south pole.

Parameters

[in]	i_north_pole	'i' index of north pole. '0' if pole is outside of grid.
[in]	j_north_pole	'j' index of north pole. '0' if pole is outside of grid.
[in]	i_south_pole	'i' index of south pole. '0' if pole is outside of grid.
[in]	j_south_pole	'j' index of south pole. '0' if pole is outside of grid.
[in]	im	i-dimension of model tile
[in]	jm	j-dimension of model tile
[out]	is_north_pole	'true' for points surrounding the north pole.
[out]	is_south_pole	'true' for points surrounding the south pole.

Author: GFDL Programmer

Definition at line 484 of file orog_utils.F90.

Referenced by io_utils::read_mdl_grid_file().

◆ find_poles()

subroutine, public orog_utils::find_poles	(	real, dimension(nx+1,ny+1), intent(in)	geolat,
		integer, intent(in)	nx,
		integer, intent(in)	ny,
		integer, intent(out)	i_north_pole,
		integer, intent(out)	j_north_pole,
		integer, intent(out)	i_south_pole,
		integer, intent(out)	j_south_pole )

Find the point on the model grid tile closest to the north and south pole.

Parameters

[in]	geolat	Latitude on the supergrid.
[in]	nx	i-dimension of the supergrid.
[in]	ny	j-dimension of the supergrid.
[out]	i_north_pole	'i' index of north pole. '0' if pole is outside of grid.
[out]	j_north_pole	'j' index of north pole. '0' if pole is outside of grid.
[out]	i_south_pole	'i' index of south pole. '0' if pole is outside of grid.
[out]	j_south_pole	'j' index of south pole. '0' if pole is outside of grid.

Author: GFDL Programmer

Definition at line 409 of file orog_utils.F90.

Referenced by io_utils::read_mdl_grid_file().

◆ get_index()

subroutine, public orog_utils::get_index	(	integer, intent(in)	imn,
		integer, intent(in)	jmn,
		integer, intent(in)	npts,
		real, dimension(npts), intent(in)	lono,
		real, dimension(npts), intent(in)	lato,
		real, intent(in)	delxn,
		integer, intent(out)	jst,
		integer, intent(out)	jen,
		integer, dimension(imn), intent(out)	ilist,
		integer, intent(out)	numx )

Determine the location of a cubed-sphere point within the high-resolution orography data.

The location is described by the range of i/j indices on the high-res grid.

Parameters

[in]	imn	'i' dimension of the high-resolution orography data set.
[in]	jmn	'j' dimension of the high-resolution orography data set.
[in]	npts	Number of vertices to describe the cubed-sphere point.
[in]	lonO	The longitudes of the cubed-sphere vertices. Must range from 0 - 360.
[in]	latO	The latitudes of the cubed-sphere vertices.
[in]	delxn	Resolution of the high-resolution orography data set.
[out]	jst	Starting 'j' index on the high-resolution grid.
[out]	jen	Ending 'j' index on the high-resolution grid.
[out]	ilist	List of 'i' indices on the high-resolution grid.
[out]	numx	The number of 'i' indices on the high-resolution grid.

Author: GFDL programmer

Definition at line 689 of file orog_utils.F90.

Referenced by make_mask(), makemt2(), makeoa2(), and makepc2().

◆ get_lat_angle()

real function, public orog_utils::get_lat_angle ( real, intent(in) dy )

Convert the 'y' direction distance of a cubed-sphere grid point to the corresponding distance in latitude.

Parameters

[in] dy Distance along the 'y' direction of a cubed-sphere point in meters.

Returns: get_lat_angle Corresponding latitudinal distance in degrees.

Author: GFDL programmer

Definition at line 127 of file orog_utils.F90.

References earth_radius, get_lat_angle(), and rad2deg.

Referenced by get_lat_angle(), and makeoa2().

◆ get_lon_angle()

real function, public orog_utils::get_lon_angle	(	real, intent(in)	dx,
		real, intent(in)	lat_in )

Convert the 'x' direction distance of a cubed-sphere grid point to the corresponding distance in longitude.

Parameters

[in]	dx	Distance along the 'x' direction of a cubed-sphere grid point in meters.
[in]	lat_in	Latitude of the cubed-sphere point in degrees.

Returns: get_lon_angle Corresponding distance in longitude in degrees.

Author: GFDL programmer

Definition at line 151 of file orog_utils.F90.

References deg2rad, earth_radius, get_lon_angle(), and rad2deg.

Referenced by get_lon_angle(), and makeoa2().

◆ get_xnsum()

real function, public orog_utils::get_xnsum	(	real, intent(in)	lon1,
		real, intent(in)	lat1,
		real, intent(in)	lon2,
		real, intent(in)	lat2,
		integer, intent(in)	imn,
		integer, intent(in)	jmn,
		real, dimension(jmn), intent(in)	glat,
		integer, dimension(imn,jmn), intent(in)	zavg,
		integer, dimension(imn,jmn), intent(in)	zslm,
		real, intent(in)	delxn )

Count the number of high-resolution orography points that are higher than the model grid box average orography height.

Parameters

[in]	lon1	Longitude of corner point 1 of the model grid box.
[in]	lat1	Latitude of corner point 1 of the model grid box.
[in]	lon2	Longitude of corner point 2 of the model grid box.
[in]	lat2	Latitude of corner point 2 of the model grid box.
[in]	imn	'i' dimension of the high-resolution orography data.
[in]	jmn	'j' dimension of the high-resolution orography data.
[in]	glat	Latitude of each row of the high-resolution orography data.
[in]	zavg	The high-resolution orography.
[in]	zslm	The high-resolution land mask.
[in]	delxn	Resolution of the high-res orography data.

Returns: get_xnsum The number of high-res points above the mean orography.

Author: GFDL Programmer

Definition at line 791 of file orog_utils.F90.

References get_xnsum().

Referenced by get_xnsum(), and makeoa2().

◆ get_xnsum2()

subroutine, public orog_utils::get_xnsum2	(	real, intent(in)	lon1,
		real, intent(in)	lat1,
		real, intent(in)	lon2,
		real, intent(in)	lat2,
		integer, intent(in)	imn,
		integer, intent(in)	jmn,
		real, dimension(jmn), intent(in)	glat,
		integer, dimension(imn,jmn), intent(in)	zavg,
		real, intent(in)	delxn,
		real, intent(out)	xnsum1,
		real, intent(out)	xnsum2,
		real, intent(out)	hc )

Count the number of high-resolution orography points that are higher than a critical value inside a model grid box (or a portion of a model grid box).

The critical value is a function of the standard deviation of orography.

Parameters

[in]	lon1	Longitude of corner point 1 of the model grid box.
[in]	lat1	Latitude of corner point 1 of the model grid box.
[in]	lon2	Longitude of corner point 2 of the model grid box.
[in]	lat2	Latitude of corner point 2 of the model grid box.
[in]	imn	'i' dimension of the high-resolution orography data.
[in]	jmn	'j' dimension of the high-resolution orography data.
[in]	glat	Latitude of each row of the high-resolution orography data.
[in]	zavg	The high-resolution orography.
[in]	delxn	Resolution of the high-res orography data.
[out]	xnsum1	The number of high-resolution orography above the critical value inside a model grid box.
[out]	xnsum2	The number of high-resolution orography points inside a model grid box.
[out]	hc	Critical height.

Author: GFDL Programmer

Definition at line 902 of file orog_utils.F90.

Referenced by makeoa2().

◆ get_xnsum3()

subroutine, public orog_utils::get_xnsum3	(	real, intent(in)	lon1,
		real, intent(in)	lat1,
		real, intent(in)	lon2,
		real, intent(in)	lat2,
		integer, intent(in)	imn,
		integer, intent(in)	jmn,
		real, dimension(jmn), intent(in)	glat,
		integer, dimension(imn,jmn), intent(in)	zavg,
		real, intent(in)	delxn,
		real, intent(out)	xnsum1,
		real, intent(out)	xnsum2,
		real, intent(in)	hc )

Count the number of high-resolution orography points that are higher than a critical value inside a model grid box (or a portion of a model grid box).

Unlike routine get_xnsum2(), this routine does not compute the critical value. Rather, it is passed in.

Parameters

[in]	lon1	Longitude of corner point 1 of the model grid box.
[in]	lat1	Latitude of corner point 1 of the model grid box.
[in]	lon2	Longitude of corner point 2 of the model grid box.
[in]	lat2	Latitude of corner point 2 of the model grid box.
[in]	imn	'i' dimension of the high-resolution orography data.
[in]	jmn	'j' dimension of the high-resolution orography data.
[in]	glat	Latitude of each row of the high-resolution orography data.
[in]	zavg	The high-resolution orography.
[in]	delxn	Resolution of the high-res orography data.
[out]	xnsum1	The number of high-resolution orography above the critical value inside a model grid box.
[out]	xnsum2	The number of high-resolution orography points inside a model grid box.
[in]	hc	Critical height.

Author: GFDL Programmer

Definition at line 1004 of file orog_utils.F90.

Referenced by makeoa2().

◆ inside_a_polygon()

logical function, public orog_utils::inside_a_polygon	(	real, intent(in)	lon1,
		real, intent(in)	lat1,
		integer, intent(in)	npts,
		real, dimension(npts), intent(in)	lon2,
		real, dimension(npts), intent(in)	lat2 )

Check if a point is inside a polygon.

Parameters

[in]	lon1	Longitude of the point to check.
[in]	lat1	Latitude of the point to check.
[in]	npts	Number of polygon vertices.
[in]	lon2	Longitude of the polygon vertices.
[in]	lat2	Latitude of the polygon vertices.

Returns: inside_a_polygon When true, point is within the polygon.

Author: GFDL programmer

Definition at line 317 of file orog_utils.F90.

References inside_a_polygon(), latlon2xyz(), pi, and spherical_angle().

Referenced by inside_a_polygon(), make_mask(), makemt2(), makeoa2(), and makepc2().

◆ latlon2xyz()

subroutine, public orog_utils::latlon2xyz	(	integer, intent(in)	siz,
		real, dimension(siz), intent(in)	lon,
		real, dimension(siz), intent(in)	lat,
		real, dimension(siz), intent(out)	x,
		real, dimension(siz), intent(out)	y,
		real, dimension(siz), intent(out)	z )

Convert from latitude and longitude to x,y,z coordinates.

Parameters

[in]	siz	Number of points to convert.
[in]	lon	Longitude (radians) of points to convert.
[in]	lat	Latitude (radians) of points to convert.
[out]	x	'x' Coordinate of the converted points.
[out]	y	'y' Coordinate of the converted points.
[out]	z	'z' Coordinate of the converted points.

Author: GFDL programmer

Definition at line 100 of file orog_utils.F90.

Referenced by inside_a_polygon().

◆ minmax()

subroutine, public orog_utils::minmax	(	integer, intent(in)	im,
		integer, intent(in)	jm,
		real, dimension(im,jm), intent(in)	a,
		character(len=8), intent(in)	title,
		integer, intent(out), optional	imax,
		integer, intent(out), optional	jmax )

Print out the maximum and minimum values of an array and optionally pass back the i/j location of the maximum.

Parameters

[in]	im	The 'i' dimension of the array.
[in]	jm	The 'j' dimension of the array.
[in]	a	The array to check.
[in]	title	Name of the data to be checked.
[out]	imax	The 'i' location of the maximum.
[out]	jmax	The 'j' location of the maximum.

Author: Jordan Alpert NOAA/EMC

Definition at line 49 of file orog_utils.F90.

Referenced by tersub().

◆ remove_isolated_pts()

subroutine, public orog_utils::remove_isolated_pts	(	integer, intent(in)	im,
		integer, intent(in)	jm,
		real, dimension(im,jm), intent(inout)	slm,
		real, dimension(im,jm), intent(inout)	oro,
		real, dimension(im,jm), intent(inout)	var,
		real, dimension(im,jm), intent(inout)	var4,
		real, dimension(im,jm,4), intent(inout)	oa,
		real, dimension(im,jm,4), intent(inout)	ol )

Remove isolated model points.

Parameters

[in]	im	'i' dimension of a model grid tile.
[in]	jm	'j' dimension of a model grid tile.
[in,out]	slm	Land-mask on the model tile.
[in,out]	oro	Orography on the model tile.
[in,out]	var	Standard deviation of orography on the model tile.
[in,out]	var4	Convexity on the model tile.
[in,out]	oa	Orographic asymmetry on the model tile.
[in,out]	ol	Orographic length scale on the model tile.

Author: Jordan Alpert NOAA/EMC

Definition at line 560 of file orog_utils.F90.

Referenced by tersub().

◆ spherical_angle()

real function orog_utils::spherical_angle	(	real, dimension(3), intent(in)	v1,
		real, dimension(3), intent(in)	v2,
		real, dimension(3), intent(in)	v3 )

private

Compute spherical angle.

Parameters

[in]	v1	Vector 1.
[in]	v2	Vector 2.
[in]	v3	Vector 3.

Returns: spherical_angle Spherical Angle.

Author: GFDL programmer

Definition at line 261 of file orog_utils.F90.

References pi, and spherical_angle().

Referenced by inside_a_polygon(), and spherical_angle().

◆ timef()

real function, public orog_utils::timef

Get the date/time from the system clock.

Returns: timef

Author: Mark Iredell

Definition at line 1063 of file orog_utils.F90.

References timef().

Referenced by tersub(), and timef().

◆ transpose_mask()

subroutine, public orog_utils::transpose_mask	(	integer, intent(in)	imn,
		integer, intent(in)	jmn,
		integer(1), dimension(imn,jmn), intent(inout)	mask )

Transpose the global landmask by flipping the poles and moving the starting longitude to Greenwich.

Parameters

[in]	imn	i-dimension of landmask data.
[in]	jmn	j-dimension of landmask data.
[in,out]	mask	The global landmask data.

Author: G. Gayno

Definition at line 176 of file orog_utils.F90.

Referenced by io_utils::read_global_mask().

◆ transpose_orog()

subroutine, public orog_utils::transpose_orog	(	integer, intent(in)	imn,
		integer, intent(in)	jmn,
		integer(2), dimension(imn,jmn), intent(inout)	glob )

Transpose the global orography data by flipping the poles and moving the starting longitude to Greenwich.

Parameters

[in]	imn	i-dimension of orography data.
[in]	jmn	j-dimension of orography data.
[in,out]	glob	The global orography data.

Author: G. Gayno

Definition at line 219 of file orog_utils.F90.

Referenced by io_utils::read_global_orog().

Variable Documentation

◆ deg2rad

real, parameter orog_utils::deg2rad = 3.14159265358979/180.

private

degrees per radians.

Definition at line 17 of file orog_utils.F90.

Referenced by get_lon_angle().

◆ earth_radius

real, parameter orog_utils::earth_radius = 6371200.

private

earth radius in meters.

Definition at line 14 of file orog_utils.F90.

Referenced by get_lat_angle(), and get_lon_angle().

◆ pi

real, parameter orog_utils::pi =3.1415926535897931

private

pi.

Definition at line 15 of file orog_utils.F90.

Referenced by inside_a_polygon(), and spherical_angle().

◆ rad2deg

real, parameter orog_utils::rad2deg = 180./3.14159265358979

private

radians per degrees.

Definition at line 16 of file orog_utils.F90.

Referenced by get_lat_angle(), and get_lon_angle().

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ find_nearest_pole_points()

◆ find_poles()

◆ get_index()

◆ get_lat_angle()

◆ get_lon_angle()

◆ get_xnsum()

◆ get_xnsum2()

◆ get_xnsum3()

◆ inside_a_polygon()

◆ latlon2xyz()

◆ minmax()

◆ remove_isolated_pts()

◆ spherical_angle()

◆ timef()

◆ transpose_mask()

◆ transpose_orog()

Variable Documentation

◆ deg2rad

◆ earth_radius

◆ pi

◆ rad2deg