|
UPP
11.0.0
|
|
UPP
11.0.0
|
rtsig() reads and transforms sigma file. More...
Go to the source code of this file.
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) |
| getomega() computes omega. More... | |
| 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) |
| rtsig() reads and transforms sigma file. More... | |
| 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... | |
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.
| [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. |
| 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 |
Definition in file GFSPOSTSIG.F.
| subroutine getomega | ( | integer, intent(in) | jcap, |
| integer, intent(in) | nc, | ||
| integer, intent(in) | km, | ||
| integer, intent(in) | idvc, | ||
| integer, intent(in) | idvm, | ||
| integer, intent(in) | idrt, | ||
| integer, intent(in) | idsl, | ||
| integer, intent(in) | nvcoord, | ||
| real*8, dimension(km+1,nvcoord), intent(in) | vcoord, | ||
| integer, intent(in) | lonb, | ||
| integer, intent(in) | latb, | ||
| integer, intent(in) | ijn, | ||
| integer, intent(in) | j1, | ||
| integer, intent(in) | j2, | ||
| integer, intent(in) | jc, | ||
| real, dimension(nc,km), intent(in) | sd, | ||
| real, dimension(nc), intent(in) | sps, | ||
| real, dimension(ijn), intent(in) | psi, | ||
| real, dimension(ijn,km), intent(in) | ti, | ||
| real, dimension(ijn,km), intent(in) | ui, | ||
| real, dimension(ijn,km), intent(in) | vi, | ||
| real, dimension(ijn,km), intent(out) | wi, | ||
| real, dimension(ijn,km), intent(out) | pm, | ||
| real, dimension(ijn,km), intent(out) | pd | ||
| ) |
getomega() computes omega.
| [in] | jcap | integer spectral truncation. |
| [in] | nc | integer first dimension (nc>=(jcap+1)*(jcap+2)). |
| [in] | km | integer number of levels. |
| [in] | idvc | integer |
| [in] | idvm | integer mass variable id. |
| [in] | idrt | integer data representation type. |
| [in] | idsl | integer |
| [in] | nvcoord | integer number of vertical coordinates. |
| [in] | vcoord | real (km+1,nvcoord) vertical coordinates. |
| [in] | lonb | integer number of longitudes. |
| [in] | latb | integer number of latitudes. |
| [in] | ijn | integer |
| [in] | j1 | integer first latitude. |
| [in] | j2 | integer last latitude. |
| [in] | jc | integer number of cpus. |
| [in] | sd | real (nc,km) divergence. |
| [in] | sps | real (nc) log surface pressure. |
| [in] | psi | real (ijn) |
| [in] | ti | real (ijn,km) |
| [in] | ui | real (ijn,km) |
| [in] | vi | real (ijn,km) |
| [out] | wi | real (ijn,km) |
| [out] | pm | real (ijn,km) |
| [out] | pd | real (ijn,km) |
Definition at line 689 of file GFSPOSTSIG.F.
Referenced by trssc().
| 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.
| [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). |
| 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 |
Definition at line 297 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.
| [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). |
| [in] | me | integer |
| 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 |
Definition at line 376 of file GFSPOSTSIG.F.
Referenced by initpost_gfs_nems_mpiio().
| subroutine rtsig | ( | integer(sigio_intkind), intent(in) | lusig, |
| type(sigio_head), intent(in) | head, | ||
| integer, intent(in) | k1, | ||
| integer, intent(in) | k2, | ||
| integer, dimension(200), intent(in) | kgds, | ||
| integer, intent(in) | ijo, | ||
| integer, intent(in) | levs, | ||
| integer, intent(in) | ntrac, | ||
| integer, intent(in) | jcap, | ||
| integer, intent(in) | lnt2, | ||
| integer, intent(in) | me, | ||
| real, dimension(ijo), intent(out) | h, | ||
| real, dimension(ijo), intent(out) | p, | ||
| real, dimension(ijo), intent(out) | px, | ||
| real, dimension(ijo), intent(out) | py, | ||
| real, dimension(ijo,k1:k2), intent(out) | t, | ||
| real, dimension(ijo,k1:k2), intent(out) | u, | ||
| real, dimension(ijo,k1:k2), intent(out) | v, | ||
| real, dimension(ijo,k1:k2), intent(out) | d, | ||
| real, dimension(ijo,k1:k2,ntrac), intent(out), target | trc, | ||
| integer, intent(out) | iret | ||
| ) |
rtsig() reads and transforms sigma file.
| [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). |
| [in] | me | integer |
| [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] | 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. |
| 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 |
Definition at line 94 of file GFSPOSTSIG.F.
| 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.
| [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] | idvc | integer |
| [in] | idvm | integer mass variable id. |
| [in] | idsl | integer |
| [in] | nvcoord | integer number of vertical coordinates. |
| [in] | vcoord | real (km+1,nvcoord) vertical coordinates. |
| [in] | cpi | real |
| [in] | idrt | integer data representation type. |
| [in] | lonb | integer number of longitudes. |
| [in] | latb | integer number of latitudes. |
| [in] | ijl | integer horizontal dimension. |
| [in] | ijn | integer |
| [in] | j1 | integer first latitude. |
| [in] | j2 | integer last latitude. |
| [in] | jc | integer number of cpus. |
| [in] | chgq0 | integer |
| [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. |
| [in,out] | gfszs | real (ijn) |
| [in,out] | gfsps | real (ijn) |
| [in,out] | gfsp | real (ijn, km) |
| [in,out] | gfsdp | real (ijn, km) |
| [in,out] | gfst | real (ijn, km) |
| [in,out] | gfsu | real (ijn, km) |
| [in,out] | gfsv | real (ijn, km) |
| [in,out] | gfsq | real (ijn, km*ntrac) |
| [in,out] | gfsw | real (ijn, km) |
| Date | Programmer | Comments |
|---|---|---|
| 1999-10-18 | Mark Iredell | Initial |
Definition at line 497 of file GFSPOSTSIG.F.
References getomega().
1.8.5 
