rtsig() reads and transforms sigma file. More...

Functions/Subroutines
subroutine	getomega (jcap, nc, km, idvc, idvm, idrt, idsl, nvcoord, vcoord, lonb, latb, ijn, j1, j2, jc, sd, sps, psi, ti, ui, vi, wi, pm, pd)

subroutine	modstuff (km, idvc, idsl, nvcoord, vcoord, ps, psx, psy, d, u, v, pi, pm, om)
	modstuff() computes model coordinate dependent functions. More...

subroutine	modstuff2 (im, ix, km, idvc, idsl, nvcoord, vcoord, ps, psx, psy, d, u, v, pi, pm, om, me)
	modstuff2() computes model coordinate dependent functions. More...

subroutine	rtsig (lusig, head, k1, k2, kgds, ijo, levs, ntrac, jcap, lnt2, me, h, p, px, py, t, u, v, d, trc, iret)

subroutine	trssc (jcap, nc, km, ntrac, idvc, idvm, idsl, nvcoord, vcoord, cpi, idrt, lonb, latb, ijl, ijn, j1, j2, jc, chgq0, szs, sps, st, sd, sz, sq, gfszs, gfsps, gfsp, gfsdp, gfst, gfsu, gfsv, gfsq, gfsw)
	trssc() transforms sigma spectral fields to grid. More...

Detailed Description

rtsig() reads and transforms sigma file.

This subprogram reads a sigma file and transforms the fields to a designated global grid. Add Iredells subroutine to read sigma files.

Parameters

[out]	lusig	integer(sigio_intkind) sigma file unit number.
[out]	head	type(sigio_head) sigma file header.
[out]	k1	integer first model level to return.
[out]	k2	integer last model level to return.
[out]	kgds	integer (200) GDS to which to transform.
[out]	ijo	integer dimension of output fields.
[out]	levs	integer number of total vertical levels.
[out]	ntrac	integer number of output tracers.
[out]	jcap	integer number of waves.
[out]	lnt2	integer (jcap+1)*(jcap+2).
[out]	h	real (ijo) surface orography (m).
[out]	p	real (ijo) surface pressure (Pa).
[out]	px	real (ijo) log surface pressure x-gradient (1/m).
[out]	py	real (ijo) log surface pressure y-gradient (1/m).
[out]	t	real (ijo,k1:k2) temperature (K).
[out]	tx	real (ijo,k1:k2) virtual temperature x-gradient (K/m).
[out]	ty	real (ijo,k1:k2) virtual temperature y-gradient (K/m).
[out]	u	real (ijo,k1:k2) x-component wind (m/s).
[out]	v	real (ijo,k1:k2) y-component wind (m/s).
[out]	d	real (ijo,k1:k2) wind divergence (1/s).
[out]	trc	real (ijo,k1:k2,ntrac) tracers. 1 = specific humidity (kg/kg) 2 = Ozone mixing ratio (kg/kg) 3 = cloud condensate mixing ratio (kg/kg) atomic oxyge, oxygen etc
[out]	iret	Integer return code.

Program History Log



Date Programmer Comments



1999-10-18 Mark Iredell      Initial


2013-04-19 Jun  Wang         Add option to get tmp and ps(in pascal) from enthalpy and ps(cb) option


2013-05-06 Shrinivas Moorthi Initialize midea to 0


2013-05-07 Shrinivas Moorthi Remove mo3, mct, midea and define io3, ict etc correctly and get correct cloud condensate.


2013-08-02 Shrinivas Moorthi Rewrote the whole routine to read the sigma file differently and to read all tracers. Added sptezj for two 2d fields


2014-02-20 Shrinivas Moorthi Modified conversion from spectral to grid taking advantage of threding in SP library. This really speeds up the code. Also threaded loop for Temperature from Tv


Author
Mark Iredell np23 
Date
1999-10-18

Definition in file GFSPOSTSIG.F.

Function/Subroutine Documentation

subroutine modstuff	(	integer, intent(in)	km,
		integer, intent(in)	idvc,
		integer, intent(in)	idsl,
		integer, intent(in)	nvcoord,
		real, dimension(km+1,nvcoord), intent(in)	vcoord,
		real, intent(in)	ps,
		real, intent(in)	psx,
		real, intent(in)	psy,
		real, dimension(km), intent(in)	d,
		real, dimension(km), intent(in)	u,
		real, dimension(km), intent(in)	v,
		real*8, dimension(km+1), intent(out)	pi,
		real*8, dimension(km), intent(out)	pm,
		real, dimension(km), intent(out)	om
	)

modstuff() computes model coordinate dependent functions.

This subprogram computes fields which depend on the model coordinate such as pressure thickness and vertical velocity.

Parameters

[in]	km	integer number of levels.
[in]	idvc	integer vertical coordinate id (1 for sigma and 2 for hybrid).
[in]	idsl	integer type of sigma structure (1 for phillips or 2 for mean).
[in]	nvcoord	integer number of vertical coordinates.
[in]	vcoord	real (km+1,nvcoord) vertical coordinates.
[in]	ps	real surface pressure (Pa).
[in]	psx	real log surface pressure x-gradient (1/m).
[in]	psy	real log surface pressure y-gradient (1/m).
[in]	d	real (km) wind divergence (1/s).
[in]	u	real (km) x-component wind (m/s).
[in]	v	real (km) y-component wind (m/s).
[out]	pi	real (km+1) interface pressure (Pa).
[out]	pm	real (km) mid-layer pressure (Pa).
[out]	om	real (km) vertical velocity (Pa/s).

Program History Log

Date	Programmer	Comments
1999-10-18	Mark Iredell	Initial
2013-04-19	Jun Wang	Add option to get pi by using 8 byte real computation

Author: Mark Iredell np23

Date: 1999-10-18

Definition at line 255 of file GFSPOSTSIG.F.

subroutine modstuff2	(	integer, intent(in)	im,
		integer, intent(in)	ix,
		integer, intent(in)	km,
		integer, intent(in)	idvc,
		integer, intent(in)	idsl,
		integer, intent(in)	nvcoord,
		real, dimension(km+1,nvcoord), intent(in)	vcoord,
		real, dimension(ix), intent(in)	ps,
		real, dimension(ix), intent(in)	psx,
		real, dimension(ix), intent(in)	psy,
		real, dimension(ix,km), intent(in)	d,
		real, dimension(ix,km), intent(in)	u,
		real, dimension(ix,km), intent(in)	v,
		real*8, dimension(ix,km+1), intent(out)	pi,
		real*8, dimension(ix,km), intent(out)	pm,
		real, dimension(ix,km), intent(out)	om,
		integer, intent(in)	me
	)

modstuff2() computes model coordinate dependent functions.

This subprogram computes fields which depend on the model coordinate such as pressure thickness and vertical velocity.

Parameters

[in]	im	integer inner computational domain.
[in]	ix	integer maximum inner dimension.
[in]	km	integer number of levels.
[in]	idvc	integer vertical coordinate id (1 for sigma and 2 for hybrid).
[in]	idsl	integer type of sigma structure (1 for phillips or 2 for mean).
[in]	nvcoord	integer number of vertical coordinates.
[in]	vcoord	real (km+1,nvcoord) vertical coordinates.
[in]	ps	real surface pressure (Pa).
[in]	psx	real log surface pressure x-gradient (1/m).
[in]	psy	real log surface pressure y-gradient (1/m).
[in]	d	real (km) wind divergence (1/s).
[in]	u	real (km) x-component wind (m/s).
[in]	v	real (km) y-component wind (m/s).
[out]	pi	real (km+1) interface pressure (Pa).
[out]	pm	real (km) mid-layer pressure (Pa).
[out]	om	real (km) vertical velocity (Pa/s).

Program History Log

Date	Programmer	Comments
1999-10-18	Mark Iredell	Initial
2013-04-19	Jun Wang	Add option to get pi by using 8 byte real computation
2013-08-13	Shrinivas Moorthi	Modified to include im points and thread

Author: Mark Iredell np23

Date: 1999-10-18

Definition at line 333 of file GFSPOSTSIG.F.

Referenced by INITPOST_GFS_NEMS_MPIIO().

subroutine trssc	(	integer, intent(in)	jcap,
		integer, intent(in)	nc,
		integer, intent(in)	km,
		integer, intent(in)	ntrac,
		integer, intent(in)	idvc,
		integer, intent(in)	idvm,
		integer, intent(in)	idsl,
		integer, intent(in)	nvcoord,
		real*8, dimension(km+1,nvcoord), intent(in)	vcoord,
		real, dimension(0:ntrac), intent(in)	cpi,
		integer, intent(in)	idrt,
		integer, intent(in)	lonb,
		integer, intent(in)	latb,
		integer, intent(in)	ijl,
		integer, intent(in)	ijn,
		integer, intent(in)	j1,
		integer, intent(in)	j2,
		integer, intent(in)	jc,
		integer, intent(in)	chgq0,
		real, dimension(nc), intent(in)	szs,
		real, dimension(nc), intent(in)	sps,
		real, dimension(nc,km), intent(in)	st,
		real, dimension(nc,km), intent(in)	sd,
		real, dimension(nc,km), intent(in)	sz,
		real, dimension(nc,km*ntrac), intent(in)	sq,
		real, dimension(ijn), intent(inout)	gfszs,
		real, dimension(ijn), intent(inout)	gfsps,
		real, dimension(ijn,km), intent(inout)	gfsp,
		real, dimension(ijn,km), intent(inout)	gfsdp,
		real, dimension(ijn,km), intent(inout)	gfst,
		real, dimension(ijn,km), intent(inout)	gfsu,
		real, dimension(ijn,km), intent(inout)	gfsv,
		real, dimension(ijn,km*ntrac), intent(inout)	gfsq,
		real, dimension(ijn,km), intent(inout)	gfsw
	)

trssc() transforms sigma spectral fields to grid.

Transforms sigma spectral fields to grid and converts log surface pressure to surface pressure and virtual temperature to temperature.

Parameters

[in]	jcap	integer spectral truncation.
[in]	nc	integer first dimension (nc>=(jcap+1)*(jcap+2)).
[in]	km	integer number of levels.
[in]	ntrac	integer number of tracers.
[in]	idvm	integer mass variable id.
[in]	idrt	integer data representation type.
[in]	lonb	integer number of longitudes.
[in]	latb	integer number of latitudes.
[in]	ijl	integer horizontal dimension.
[in]	j1	integer first latitude.
[in]	j2	integer last latitude.
[in]	jc	integer number of cpus.
[in]	szs	real (nc) orography.
[in]	sps	real (nc) log surface pressure.
[in]	st	real (nc,levs) virtual temperature.
[in]	sd	real (nc,levs) divergence.
[in]	sz	real (nc,levs) vorticity.
[in]	sq	real (nc,levs*ntrac) tracers.
[out]	zs	real (ijl) orography.
[out]	ps	real (ijl) surface pressure.
[out]	t	real (ijl,km) temperature.
[out]	u	real (ijl,km) zonal wind.
[out]	v	real (ijl,km) meridional wind.
[out]	q	real (ijl,km*ntrac) tracers.

Program History Log

Date	Programmer	Comments
1999-10-18	Mark Iredell	Initial

Author: Mark Iredell w/nmc23

Date: 1992-10-31

Definition at line 444 of file GFSPOSTSIG.F.

Functions/Subroutines

Detailed Description

Program History Log

Function/Subroutine Documentation

Program History Log

Program History Log

Program History Log