rayso.f90 File Reference

Functions/Subroutines
subroutine	reftra (pioc, piaero, gasym, gaero, tauc, taua, ref, tra, epsc, dqqv, mui)
	Compute reflexion and transmission.
subroutine	rayso (ivertc, k1, kmray, heuray, imer1, albe, qqv, qqqv, qqvinf, zqq, zray, qvray, qlray, fneray, romray, preray, temray, fos, rayst, ncray)
	Compute solar fluxes for both clear and cloudy atmosphere following Lacis and Hansen (1974). The multiple diffusion is taken into account by an addition method and overlapping between water vapor and liquid water with k distribution method. Some improvements from original version concerns:
double precision function	rayuoz (zh)
	Computes ozone amount above a given altitude.

Function/Subroutine Documentation

rayso()

subroutine rayso	(	integer	ivertc,
		integer	k1,
		integer	kmray,
		double precision	heuray,
		integer	imer1,
		double precision	albe,
		double precision, dimension(kmx+1)	qqv,
		double precision, dimension(kmx+1)	qqqv,
		double precision	qqvinf,
		double precision, dimension(kmx+1)	zqq,
		double precision, dimension(kmx)	zray,
		double precision, dimension(kmx)	qvray,
		double precision, dimension(kmx)	qlray,
		double precision, dimension(kmx)	fneray,
		double precision, dimension(kmx)	romray,
		double precision, dimension(kmx)	preray,
		double precision, dimension(kmx)	temray,
		double precision	fos,
		double precision, dimension(kmx)	rayst,
		double precision, dimension(kmx)	ncray )

Compute solar fluxes for both clear and cloudy atmosphere following Lacis and Hansen (1974). The multiple diffusion is taken into account by an addition method and overlapping between water vapor and liquid water with k distribution method. Some improvements from original version concerns:

• introduction of cloud fraction with hazardous recovering
• introduction of aerosol diffusion in the same way as for cloud droplets but with specific optical properties for aerosols.

  Parameters

  [in]	ivertc	index of vertical profile
  [in]	k1	index for ground level
  [in]	kmray	vertical levels number for radiation
  [in]	heuray	Universal time (Hour)
  [in]	imer1	sea index
  [in]	albe	albedo
  [in]	qqv	optical depth for water vapor (0,z)
  [in]	qqqv	idem for intermediate levels
  [in]	qqvinf	idem qqv but for altitude above 11000m
  [in]	zqq	vertical levels
  [in]	zray	altitude (physical mesh)
  [in]	qvray	specific humidity for water vapor
  [in]	qlray	specific humidity for liquid water
  [in]	fneray	cloud fraction
  [in]	romray	air density
  [in]	preray	pressure
  [in]	temray	temperature
  [out]	fos	global downward solar flux at the ground
  [out]	rayst	flux divergence of solar radiation

rayuoz()

double precision function rayuoz

(

double precision, intent(in)

zh

)

Computes ozone amount above a given altitude.

Parameters

[in]

zh

altitude

reftra()

subroutine reftra	(	double precision, intent(in)	pioc,
		double precision, intent(in)	piaero,
		double precision, intent(in)	gasym,
		double precision, intent(in)	gaero,
		double precision, intent(in)	tauc,
		double precision, intent(in)	taua,
		double precision, intent(inout)	ref,
		double precision, intent(inout)	tra,
		double precision, intent(in)	epsc,
		double precision, intent(in)	dqqv,
		double precision	mui )

Compute reflexion and transmission.

Parameters

[in]	pioc	Albedo of simple diffusion for cloud (water)
[in]	piaero	Albedo of simple diffusion for aerosol
[in]	gasym	Asymmetry factor for clouds
[in]	gaero	Asymmetry factor for aerosols
[in]	tauc	Optical depth for clouds
[in]	taua	Optical depth for aersols
[out]	ref	Reflexion
[out]	tra	Transmission
[in]	epsc	clipping threshold
[in]	dqqv	Optical depth for Water vapor