gaus

Computes gaussian latitudes and weights.

Prototype

	function gaus (
		nlat [1] : integer   
	)

	return_val [2*nlat,2] :  double

Arguments

A scalar integer equal to the number of latitude points PER hemisphere.

Description

A two dimensional array is returned that contains both the gaussian latitudes (first dimension) and the gaussian weights (second dimension). By default the returned values are of type double.

See Also

There are numerous other latitude/longitude functions.

latGau, latGauWgt, latGlobeF, latGlobeFo, lonGlobeF, lonGlobeFo

Examples

Example 1

Generate gaussian latitudes and weights (type double) for a specified number of latitudes:

  nlat = 64  ; for globe

  gau_info = gaus(nlat/2)   ; divide by 2 to get "per hemisphere"
  glat     = gau_info(:,0)  ; gaussian latitudes ( 1st dimension of gau_info)
  gwgt     = gau_info(:,1)  ; gaussian weights   ( 2nd dimension of gau_info)

Use doubletofloat to generate gaussian latitudes and weights (type float) for a specified number of latitudes:

  nlat     = 94     ; for globe
; conversion is embedded
  gau_info = doubletofloat(gaus(nlat/2))  
  glat     = gau_info(:,0)                ; glat ==> float
  gwgt     = gau_info(:,1)                ; gwgt ==> float