!-------------------------------------- LICENCE BEGIN ------------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
! version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
!without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
!See the above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
***S/P - PREINTR - DETERMINES THE INTERPRETION POINTS OF H2O AND CO2
*
#include "phy_macros_f.h"
subroutine preintr (inpr, dir, qq, rhc, il1, il2, ilg, lay) 1
*
#include "impnone.cdk"
*
integer ilg, lay, il1, il2, jendr, k, i, j, l, lp1
real rrmco2
real dir(ilg,lay), qq(ilg,lay), rhc(ilg,lay), standr(5)
integer inpr(ilg,lay)
*
*Authors
* J. Li, M. Lazare, CCCMA, rt code for gcm4
* (Ref: J. Li, H. W. Barker, 2005:
* JAS Vol. 62, no. 2, pp. 286\226309)
* P. Vaillancourt, D. Talbot, RPN/CMC;
* adapted for CMC/RPN physics (May 2006)
*
*Revisions
* 001
*
*Object
* This subroutine determines the interpretion points for the ratio
* of h2o and co2.
*
*Arguments
*
* inpr number of the ratio level for the standard 5 ratios
* dir interpolation factor for mass ratio of h2o / co2
* between two neighboring standard input ratios
* q: water vapor mass mixing ratio
* rhc the ratio of the h2o mass mixing to co2 mass mixing
*----------------------------------------------------------------------
*
*Implicites
#include "tracegases.cdk"
**
data standr / .06, .24, 1., 4., 16. /
c
rrmco2 = 1.0 / (rmco2 + 1.e-10)
jendr = 5
do 400 k = 1, lay
do 400 i = il1, il2
inpr(i,k) = 0
rhc(i,k) = qq(i,k) * rrmco2
c
do 200 j = 1, jendr
if (rhc(i,k) .gt. standr(j)) then
inpr(i,k) = inpr(i,k) + 1
endif
200 continue
c
l = inpr(i,k)
lp1 = l + 1
if (l .ge. 1 .and. l .lt. 5) then
dir(i,k) = (rhc(i,k) - standr(l)) /
1 (standr(lp1) - standr(l))
else
c
c----------------------------------------------------------------------
c dir is not used with values of {0,5} in tlinehc, but we
c initialize here to avoid problems with nan when used
c in multitasking mode.
c----------------------------------------------------------------------
c
dir(i,k) = 0.0
endif
c
400 continue
c
return
end