su_uvrad.F90

Go to the documentation of this file.

SUBROUTINE su_uvrad ( KUV )

!**** *SU_UVRAD*   - INITIALIZE COMMON YOE_UVRAD

!     PURPOSE.
!     --------
!           INITIALIZE YOE_UVRAD, THE COMMON THAT CONTAINS COEFFICIENTS
!           NEEDED TO RUN THE ULTRAVIOLET RADIATION PROCESSOR

!**   INTERFACE.
!     ----------
!        *CALL* *SU_UVRAD* FROM *SUECRAD*


!        IMPLICIT ARGUMENTS :
!        --------------------
!        COMMON YOE_UVRAD


!     METHOD.
!     -------
!        SEE DOCUMENTATION

!     EXTERNALS.
!     ----------

!     REFERENCE.
!     ----------
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS

!     AUTHOR.
!     -------
!        JEAN-JACQUES MORCRETTE *ECMWF*

!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 2005-10-04

!     ------------------------------------------------------------------

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook

USE yomlun   , ONLY : nulout
USE yoe_uvrad, ONLY : ipuv, juvlam, &
              & rk250 , rtuv1, rtuv2, rsuvb0, rayuvb, &
              & rtauva, rpiuva, rcguva, rcieas, ruvlam, rfcaer, rfcozo, &
              & rasa  , rasb, rasc, rasd, rase, rasf, rfa0, rfa1, &
              & rfb0  , rfb1, rfb2, rfb3, rfc0, rfc1, rfc2, rfc3, &
              & rfd0  , rfd1, rfd2, rfd3, jcop, rxpo, rxpl 

IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KUV 


INTEGER(KIND=JPIM) :: J, JAE
!INTEGER(KIND=JPIM) :: IACOP(24), IBCOP(120), ICCOP(600), IDCOP(3000)
INTEGER(KIND=JPIM) :: IACOP(24), IBCOP(120), ICCOP(600)

REAL(KIND=JPRB) :: ZBEG, ZEXP, ZHRES, ZRES, ZUVLAM
REAL(KIND=JPRB) :: ZAK250(24)  ,ZATUV1(24)  ,ZATUV2(24)  ,ZASUVB(24)  ,ZARAYUV(24)
REAL(KIND=JPRB) :: ZBK250(120) ,ZBTUV1(120) ,ZBTUV2(120) ,ZBSUVB(120) ,ZBRAYUV(120)
REAL(KIND=JPRB) :: ZCK250(600) ,ZCTUV1(600) ,ZCTUV2(600) ,ZCSUVB(600) ,ZCRAYUV(600)
!REAL(KIND=JPRB) :: ZDK250(3000),ZDTUV1(3000),ZDTUV2(3000),ZDSUVB(3000),ZDRAYUV(3000)

REAL(KIND=JPRB) :: ZHOOK_HANDLE
!     ----------------------------------------------------------------
!#include "abor1.intfb.h"

!     ----------------------------------------------------------------
IF (lhook) CALL dr_hook('SU_UVRAD',0,zhook_handle)

! multiplying factors for sensitivity studies

rfcaer=1.0_jprb
rfcozo=1.0_jprb

!-- explicit computations from 280. to 400. nm by (400-280)/(KUV-1) nm

!-- correspondence table with intervals in Slingo (1989) and Fu (1996)

iacop = (/ 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4 /) 
ibcop = (/ 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,&
         & 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,&
         & 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,&
         & 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 /)
iccop = (/ 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,&
         & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,&
         & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,&
         & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,&
         & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,&
         & 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,&
         & 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,&
         & 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,&
         & 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,&
         & 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,&
         & 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,&
         & 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,&
         & 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,&
         & 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,&
         & 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,&
         & 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,&
         & 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,&
         & 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,&
         & 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,&
         & 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 /)
!IDCOP = (/ 

!          O3 absorption parameters from integrating Molina & Molina, 1986

! RK250: ozone cross-section ( x E-20 cm2 molecule-1 )

! revised 20051031 (code "read_integrate_kurucz_rayleigh_O3")
! check norms: 0.1262176396E+05  0.1293728681E+05  0.1350083624E+05
zak250 = (/ &                                                            
 & 0.3053413e+03_jprb, 0.1737485e+03_jprb, 0.1045573e+03_jprb, 0.5265653e+02_jprb, 0.2620346e+02_jprb, &
 & 0.1350141e+02_jprb, 0.6771311e+01_jprb, 0.3435206e+01_jprb, 0.1629208e+01_jprb, 0.7546285e+00_jprb, &
 & 0.3856211e+00_jprb, 0.1501495e+00_jprb, 0.7608691e-01_jprb, 0.2646207e-01_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)
zatuv1 = (/ &                                                              
 & 0.1955606e-03_jprb, 0.4455332e-03_jprb, 0.6856191e-03_jprb, 0.1009841e-02_jprb, 0.1181320e-02_jprb, &
 & 0.1593266e-02_jprb, 0.1902945e-02_jprb, 0.2377848e-02_jprb, 0.3512989e-02_jprb, 0.4511920e-02_jprb, &
 & 0.4672713e-02_jprb, 0.6355066e-02_jprb, 0.7100463e-02_jprb, 0.1112103e-01_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)
zatuv2 = (/ &                                                              
 &-0.3391206e-06_jprb, 0.1269089e-05_jprb, 0.3833753e-05_jprb, 0.5239359e-05_jprb, 0.1149881e-04_jprb, &
 & 0.1214144e-04_jprb, 0.1456489e-04_jprb, 0.1640689e-04_jprb, 0.8903441e-05_jprb, 0.9585654e-05_jprb, &
 & 0.1444280e-04_jprb, 0.2252139e-04_jprb, 0.7899461e-05_jprb, 0.1340369e-04_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)

! check norms: 0.1262176396E+05  0.1293728681E+05  0.1350083624E+05
zbk250 = (/ &                                                            
 & 0.3746235e+03_jprb, 0.3374782e+03_jprb, 0.3109353e+03_jprb, 0.2882899e+03_jprb, 0.2577510e+03_jprb, &
 & 0.2259643e+03_jprb, 0.2085487e+03_jprb, 0.1854031e+03_jprb, 0.1612119e+03_jprb, 0.1453486e+03_jprb, &
 & 0.1310256e+03_jprb, 0.1133554e+03_jprb, 0.1020065e+03_jprb, 0.8859325e+02_jprb, 0.7909053e+02_jprb, &
 & 0.6807025e+02_jprb, 0.6024819e+02_jprb, 0.5107318e+02_jprb, 0.4535628e+02_jprb, 0.3949330e+02_jprb, &
 & 0.3422062e+02_jprb, 0.2996915e+02_jprb, 0.2577570e+02_jprb, 0.2279167e+02_jprb, 0.2034549e+02_jprb, &
 & 0.1721408e+02_jprb, 0.1564205e+02_jprb, 0.1300914e+02_jprb, 0.1160289e+02_jprb, 0.9688561e+01_jprb, &
 & 0.8512331e+01_jprb, 0.7888720e+01_jprb, 0.6400194e+01_jprb, 0.5981788e+01_jprb, 0.4811029e+01_jprb, &
 & 0.4636916e+01_jprb, 0.3695886e+01_jprb, 0.3476093e+01_jprb, 0.2543994e+01_jprb, 0.2691098e+01_jprb, &
 & 0.2131005e+01_jprb, 0.1743577e+01_jprb, 0.1933369e+01_jprb, 0.1154597e+01_jprb, 0.1144737e+01_jprb, &
 & 0.1227887e+01_jprb, 0.6772851e+00_jprb, 0.8242575e+00_jprb, 0.7724384e+00_jprb, 0.3580751e+00_jprb, &
 & 0.5599410e+00_jprb, 0.5106917e+00_jprb, 0.2370067e+00_jprb, 0.3617942e+00_jprb, 0.2703357e+00_jprb, &
 & 0.1238728e+00_jprb, 0.2085226e+00_jprb, 0.2350139e+00_jprb, 0.9726696e-01_jprb, 0.1077734e+00_jprb, &
 & 0.1266662e+00_jprb, 0.5216565e-01_jprb, 0.4033553e-01_jprb, 0.7404691e-01_jprb, 0.8175474e-01_jprb, &
 & 0.3928745e-01_jprb, 0.2581425e-01_jprb, 0.3220713e-01_jprb, 0.1852549e-01_jprb, 0.1710368e-01_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)
zbtuv1 = (/ &                                                              
 & 0.1304379e-03_jprb, 0.2381045e-03_jprb, 0.1312193e-03_jprb, 0.1382507e-03_jprb, 0.4377790e-03_jprb, &
 & 0.2744315e-03_jprb, 0.3904921e-03_jprb, 0.4057293e-03_jprb, 0.3944389e-03_jprb, 0.6213526e-03_jprb, &
 & 0.5230166e-03_jprb, 0.7115036e-03_jprb, 0.8224138e-03_jprb, 0.7839015e-03_jprb, 0.6757493e-03_jprb, &
 & 0.1013929e-02_jprb, 0.9388676e-03_jprb, 0.1139504e-02_jprb, 0.9152941e-03_jprb, 0.1010882e-02_jprb, &
 & 0.1092914e-02_jprb, 0.1116954e-02_jprb, 0.1157353e-02_jprb, 0.1312172e-02_jprb, 0.1258207e-02_jprb, &
 & 0.1565364e-02_jprb, 0.1304358e-02_jprb, 0.1921441e-02_jprb, 0.1597455e-02_jprb, 0.1585083e-02_jprb, &
 & 0.1804721e-02_jprb, 0.1717419e-02_jprb, 0.2201039e-02_jprb, 0.1650979e-02_jprb, 0.2436943e-02_jprb, &
 & 0.1758300e-02_jprb, 0.2866617e-02_jprb, 0.2292916e-02_jprb, 0.3437942e-02_jprb, 0.2268845e-02_jprb, &
 & 0.3270049e-02_jprb, 0.3506425e-02_jprb, 0.2534346e-02_jprb, 0.5591768e-02_jprb, 0.3701606e-02_jprb, &
 & 0.3035392e-02_jprb, 0.6216419e-02_jprb, 0.3456239e-02_jprb, 0.4134974e-02_jprb, 0.8414933e-02_jprb, &
 & 0.3396263e-02_jprb, 0.4222453e-02_jprb, 0.8249625e-02_jprb, 0.3627255e-02_jprb, 0.6302904e-02_jprb, &
 & 0.1011540e-01_jprb, 0.4498586e-02_jprb, 0.4007648e-02_jprb, 0.8463577e-02_jprb, 0.7236953e-02_jprb, &
 & 0.5225470e-02_jprb, 0.1114658e-01_jprb, 0.1362563e-01_jprb, 0.6584077e-02_jprb, 0.4911811e-02_jprb, &
 & 0.1365494e-01_jprb, 0.9170420e-02_jprb, 0.7847505e-02_jprb, 0.8155373e-02_jprb, 0.2917714e-01_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)
zbtuv2 = (/ &                                                              
 &-0.2310595e-05_jprb,-0.1588555e-05_jprb,-0.3803504e-06_jprb, 0.2203029e-05_jprb,-0.2347744e-05_jprb, &
 & 0.3014092e-05_jprb,-0.2353008e-05_jprb, 0.3071326e-05_jprb, 0.2498274e-05_jprb, 0.6234744e-06_jprb, &
 & 0.3393139e-05_jprb, 0.5542024e-05_jprb, 0.9758300e-06_jprb, 0.5219182e-05_jprb, 0.3721842e-05_jprb, &
 & 0.2097421e-05_jprb, 0.4340026e-05_jprb, 0.5201248e-05_jprb, 0.8369003e-05_jprb, 0.8713780e-05_jprb, &
 & 0.8890532e-05_jprb, 0.1285612e-04_jprb, 0.1431880e-04_jprb, 0.1206355e-04_jprb, 0.8875655e-05_jprb, &
 & 0.1178909e-04_jprb, 0.1027358e-04_jprb, 0.6029656e-05_jprb, 0.1471264e-04_jprb, 0.2339166e-04_jprb, &
 & 0.1449467e-04_jprb, 0.1392035e-04_jprb, 0.9490171e-05_jprb, 0.1963981e-04_jprb, 0.1545245e-04_jprb, &
 & 0.2217010e-04_jprb, 0.9402666e-05_jprb, 0.1555344e-04_jprb, 0.1678093e-04_jprb, 0.1418057e-04_jprb, &
 & 0.1018282e-04_jprb, 0.1016892e-04_jprb, 0.6648122e-05_jprb,-0.1761027e-05_jprb, 0.1343048e-04_jprb, &
 & 0.5945775e-05_jprb, 0.5228630e-05_jprb, 0.1189961e-04_jprb, 0.1062544e-04_jprb, 0.6783912e-05_jprb, &
 & 0.1611166e-04_jprb, 0.9911433e-05_jprb, 0.4732061e-05_jprb, 0.2326484e-04_jprb, 0.6846333e-05_jprb, &
 & 0.5566848e-05_jprb, 0.2998335e-04_jprb, 0.4369986e-05_jprb, 0.2288655e-04_jprb, 0.4555816e-04_jprb, &
 & 0.2229065e-04_jprb,-0.5176777e-05_jprb, 0.9089374e-05_jprb, 0.2211234e-04_jprb,-0.5195374e-04_jprb, &
 &-0.5204193e-04_jprb, 0.1285348e-03_jprb,-0.3511670e-05_jprb, 0.6086764e-04_jprb,-0.2923827e-03_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)

! check norms: 0.1262176396E+05  0.1293728681E+05  0.1350083624E+05
zck250(  1:300) = (/ &
 & 0.3912691e+03_jprb, 0.3839502e+03_jprb, 0.3743189e+03_jprb, 0.3712760e+03_jprb, 0.3682782e+03_jprb, &
 & 0.3607995e+03_jprb, 0.3484512e+03_jprb, 0.3362245e+03_jprb, 0.3306447e+03_jprb, 0.3257288e+03_jprb, &
 & 0.3205932e+03_jprb, 0.3156011e+03_jprb, 0.3106752e+03_jprb, 0.3061641e+03_jprb, 0.3025993e+03_jprb, &
 & 0.2990364e+03_jprb, 0.2958928e+03_jprb, 0.2922895e+03_jprb, 0.2815917e+03_jprb, 0.2710969e+03_jprb, &
 & 0.2643130e+03_jprb, 0.2602894e+03_jprb, 0.2555913e+03_jprb, 0.2508463e+03_jprb, 0.2458966e+03_jprb, &
 & 0.2400022e+03_jprb, 0.2328428e+03_jprb, 0.2266644e+03_jprb, 0.2237314e+03_jprb, 0.2211230e+03_jprb, &
 & 0.2176822e+03_jprb, 0.2129666e+03_jprb, 0.2082285e+03_jprb, 0.2050648e+03_jprb, 0.2015646e+03_jprb, &
 & 0.1964826e+03_jprb, 0.1907146e+03_jprb, 0.1843511e+03_jprb, 0.1790203e+03_jprb, 0.1740547e+03_jprb, &
 & 0.1699693e+03_jprb, 0.1658762e+03_jprb, 0.1631024e+03_jprb, 0.1594428e+03_jprb, 0.1561083e+03_jprb, &
 & 0.1522514e+03_jprb, 0.1489648e+03_jprb, 0.1456738e+03_jprb, 0.1427900e+03_jprb, 0.1400463e+03_jprb, &
 & 0.1368885e+03_jprb, 0.1338666e+03_jprb, 0.1308455e+03_jprb, 0.1279194e+03_jprb, 0.1252291e+03_jprb, &
 & 0.1221635e+03_jprb, 0.1173480e+03_jprb, 0.1134466e+03_jprb, 0.1104490e+03_jprb, 0.1075005e+03_jprb, &
 & 0.1052409e+03_jprb, 0.1034875e+03_jprb, 0.1014312e+03_jprb, 0.9917348e+02_jprb, 0.9690936e+02_jprb, &
 & 0.9363521e+02_jprb, 0.9064002e+02_jprb, 0.8742842e+02_jprb, 0.8607183e+02_jprb, 0.8419236e+02_jprb, &
 & 0.8235930e+02_jprb, 0.8037249e+02_jprb, 0.7882393e+02_jprb, 0.7689006e+02_jprb, 0.7481278e+02_jprb, &
 & 0.7260259e+02_jprb, 0.7000890e+02_jprb, 0.6741465e+02_jprb, 0.6606238e+02_jprb, 0.6431214e+02_jprb, &
 & 0.6274430e+02_jprb, 0.6096603e+02_jprb, 0.5956343e+02_jprb, 0.5820831e+02_jprb, 0.5665523e+02_jprb, &
 & 0.5500309e+02_jprb, 0.5351154e+02_jprb, 0.5118981e+02_jprb, 0.5002563e+02_jprb, 0.4881731e+02_jprb, &
 & 0.4748677e+02_jprb, 0.4631187e+02_jprb, 0.4524318e+02_jprb, 0.4422579e+02_jprb, 0.4326350e+02_jprb, &
 & 0.4201139e+02_jprb, 0.4046981e+02_jprb, 0.3921751e+02_jprb, 0.3804230e+02_jprb, 0.3728167e+02_jprb, &
 & 0.3646048e+02_jprb, 0.3568511e+02_jprb, 0.3469347e+02_jprb, 0.3368040e+02_jprb, 0.3223647e+02_jprb, &
 & 0.3120151e+02_jprb, 0.3046987e+02_jprb, 0.2973342e+02_jprb, 0.2916021e+02_jprb, 0.2869675e+02_jprb, &
 & 0.2819547e+02_jprb, 0.2704570e+02_jprb, 0.2633426e+02_jprb, 0.2524253e+02_jprb, 0.2440977e+02_jprb, &
 & 0.2368355e+02_jprb, 0.2307297e+02_jprb, 0.2260014e+02_jprb, 0.2230349e+02_jprb, 0.2196699e+02_jprb, &
 & 0.2159911e+02_jprb, 0.2100203e+02_jprb, 0.2042105e+02_jprb, 0.1959785e+02_jprb, 0.1850330e+02_jprb, &
 & 0.1799117e+02_jprb, 0.1739673e+02_jprb, 0.1696323e+02_jprb, 0.1671992e+02_jprb, 0.1652047e+02_jprb, &
 & 0.1630084e+02_jprb, 0.1608451e+02_jprb, 0.1578843e+02_jprb, 0.1524567e+02_jprb, 0.1440076e+02_jprb, &
 & 0.1386541e+02_jprb, 0.1329988e+02_jprb, 0.1272576e+02_jprb, 0.1254554e+02_jprb, 0.1245282e+02_jprb, &
 & 0.1218871e+02_jprb, 0.1178283e+02_jprb, 0.1145067e+02_jprb, 0.1133450e+02_jprb, 0.1119256e+02_jprb, &
 & 0.1078407e+02_jprb, 0.1020251e+02_jprb, 0.9528469e+01_jprb, 0.9329364e+01_jprb, 0.9126274e+01_jprb, &
 & 0.8931937e+01_jprb, 0.8718665e+01_jprb, 0.8575605e+01_jprb, 0.8437309e+01_jprb, 0.8310787e+01_jprb, &
 & 0.8215147e+01_jprb, 0.8127766e+01_jprb, 0.8001044e+01_jprb, 0.7630502e+01_jprb, 0.7218388e+01_jprb, &
 & 0.6842024e+01_jprb, 0.6619379e+01_jprb, 0.6379956e+01_jprb, 0.6211814e+01_jprb, 0.6073732e+01_jprb, &
 & 0.6047354e+01_jprb, 0.6167841e+01_jprb, 0.6213631e+01_jprb, 0.5978124e+01_jprb, 0.5663016e+01_jprb, &
 & 0.5312896e+01_jprb, 0.4966990e+01_jprb, 0.4632250e+01_jprb, 0.4511959e+01_jprb, 0.4437874e+01_jprb, &
 & 0.4521541e+01_jprb, 0.4730288e+01_jprb, 0.4850666e+01_jprb, 0.4597253e+01_jprb, 0.4254414e+01_jprb, &
 & 0.3974671e+01_jprb, 0.3772243e+01_jprb, 0.3622352e+01_jprb, 0.3596173e+01_jprb, 0.3597771e+01_jprb, &
 & 0.3587278e+01_jprb, 0.3569527e+01_jprb, 0.3534352e+01_jprb, 0.3400931e+01_jprb, 0.3269153e+01_jprb, &
 & 0.3038523e+01_jprb, 0.2772412e+01_jprb, 0.2532108e+01_jprb, 0.2400635e+01_jprb, 0.2269651e+01_jprb, &
 & 0.2318103e+01_jprb, 0.2555933e+01_jprb, 0.2783417e+01_jprb, 0.2859555e+01_jprb, 0.2929581e+01_jprb, &
 & 0.2773406e+01_jprb, 0.2471774e+01_jprb, 0.2117337e+01_jprb, 0.1887952e+01_jprb, 0.1665876e+01_jprb, &
 & 0.1545364e+01_jprb, 0.1530765e+01_jprb, 0.1569366e+01_jprb, 0.1765859e+01_jprb, 0.2076064e+01_jprb, &
 & 0.2118826e+01_jprb, 0.2013361e+01_jprb, 0.1904165e+01_jprb, 0.1817854e+01_jprb, 0.1717719e+01_jprb, &
 & 0.1528353e+01_jprb, 0.1301476e+01_jprb, 0.1075097e+01_jprb, 0.9913148e+00_jprb, 0.9390067e+00_jprb, &
 & 0.9410631e+00_jprb, 0.9908390e+00_jprb, 0.1077318e+01_jprb, 0.1253445e+01_jprb, 0.1470490e+01_jprb, &
 & 0.1500942e+01_jprb, 0.1410377e+01_jprb, 0.1291770e+01_jprb, 0.1121325e+01_jprb, 0.9340122e+00_jprb, &
 & 0.7967983e+00_jprb, 0.6878130e+00_jprb, 0.6101190e+00_jprb, 0.6281103e+00_jprb, 0.6616779e+00_jprb, &
 & 0.6857364e+00_jprb, 0.7056597e+00_jprb, 0.7477958e+00_jprb, 0.9194159e+00_jprb, 0.1109930e+01_jprb, &
 & 0.1092774e+01_jprb, 0.9152249e+00_jprb, 0.7291639e+00_jprb, 0.6087120e+00_jprb, 0.4870125e+00_jprb, &
 & 0.4089467e+00_jprb, 0.3685752e+00_jprb, 0.3364085e+00_jprb, 0.3355197e+00_jprb, 0.3394912e+00_jprb, &
 & 0.4048120e+00_jprb, 0.4981886e+00_jprb, 0.6072237e+00_jprb, 0.6562844e+00_jprb, 0.7088529e+00_jprb, &
 & 0.6887763e+00_jprb, 0.6012469e+00_jprb, 0.5178619e+00_jprb, 0.4263444e+00_jprb, 0.3389390e+00_jprb, &
 & 0.2766929e+00_jprb, 0.2419911e+00_jprb, 0.2081429e+00_jprb, 0.2205261e+00_jprb, 0.2407883e+00_jprb, &
 & 0.2829715e+00_jprb, 0.3323930e+00_jprb, 0.3840338e+00_jprb, 0.4009112e+00_jprb, 0.4141605e+00_jprb, &
 & 0.3841738e+00_jprb, 0.3191549e+00_jprb, 0.2538516e+00_jprb, 0.2141375e+00_jprb, 0.1749193e+00_jprb, &
 & 0.1442347e+00_jprb, 0.1308030e+00_jprb, 0.1175839e+00_jprb, 0.1152005e+00_jprb, 0.1144295e+00_jprb, &
 & 0.1332690e+00_jprb, 0.1750619e+00_jprb, 0.2049352e+00_jprb, 0.2269636e+00_jprb, 0.2466895e+00_jprb, &
 & 0.2630055e+00_jprb, 0.2741632e+00_jprb, 0.2762103e+00_jprb, 0.2284683e+00_jprb, 0.1768542e+00_jprb, &
 & 0.1404621e+00_jprb, 0.1127350e+00_jprb, 0.8942173e-01_jprb, 0.8041714e-01_jprb, 0.7391840e-01_jprb, &
 & 0.7558672e-01_jprb, 0.8552747e-01_jprb, 0.9933794e-01_jprb, 0.1231135e+00_jprb, 0.1472332e+00_jprb /)
zck250(301:600) = (/ &
 & 0.1526061e+00_jprb, 0.1421443e+00_jprb, 0.1296576e+00_jprb, 0.1091694e+00_jprb, 0.8614825e-01_jprb, &
 & 0.6917399e-01_jprb, 0.5783544e-01_jprb, 0.4674525e-01_jprb, 0.4370108e-01_jprb, 0.4142633e-01_jprb, &
 & 0.3910948e-01_jprb, 0.3637619e-01_jprb, 0.3495317e-01_jprb, 0.4133272e-01_jprb, 0.4971278e-01_jprb, &
 & 0.5383842e-01_jprb, 0.5446874e-01_jprb, 0.6006002e-01_jprb, 0.8558061e-01_jprb, 0.1183642e+00_jprb, &
 & 0.1234464e+00_jprb, 0.1012597e+00_jprb, 0.8177673e-01_jprb, 0.7039373e-01_jprb, 0.5842135e-01_jprb, &
 & 0.4994100e-01_jprb, 0.4395068e-01_jprb, 0.3678789e-01_jprb, 0.3423710e-01_jprb, 0.3205213e-01_jprb, &
 & 0.2882322e-01_jprb, 0.2490090e-01_jprb, 0.2239739e-01_jprb, 0.2488575e-01_jprb, 0.2809792e-01_jprb, &
 & 0.3163878e-01_jprb, 0.3529353e-01_jprb, 0.3743464e-01_jprb, 0.3154792e-01_jprb, 0.2476330e-01_jprb, &
 & 0.2125607e-01_jprb, 0.2063931e-01_jprb, 0.1970776e-01_jprb, 0.1711069e-01_jprb, 0.1448247e-01_jprb, &
 & 0.1419092e-01_jprb, 0.1697305e-01_jprb, 0.1922265e-01_jprb, 0.1809877e-01_jprb, 0.1633014e-01_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)
zctuv1(  1:300) = (/ &
 &-0.7300886e-04_jprb, 0.2165305e-04_jprb, 0.1518885e-03_jprb, 0.1801729e-03_jprb, 0.2046781e-03_jprb, &
 & 0.2321395e-03_jprb, 0.2666132e-03_jprb, 0.2905016e-03_jprb, 0.2372348e-03_jprb, 0.1789425e-03_jprb, &
 & 0.1385595e-03_jprb, 0.1176325e-03_jprb, 0.1025867e-03_jprb, 0.1330112e-03_jprb, 0.1624851e-03_jprb, &
 & 0.1525305e-03_jprb, 0.8864058e-04_jprb, 0.4904664e-04_jprb, 0.1451164e-03_jprb, 0.2617269e-03_jprb, &
 & 0.3550630e-03_jprb, 0.4328231e-03_jprb, 0.5061485e-03_jprb, 0.5026669e-03_jprb, 0.4904191e-03_jprb, &
 & 0.4514644e-03_jprb, 0.3765797e-03_jprb, 0.3052066e-03_jprb, 0.2502422e-03_jprb, 0.1951159e-03_jprb, &
 & 0.2224293e-03_jprb, 0.3463365e-03_jprb, 0.4626359e-03_jprb, 0.4548639e-03_jprb, 0.4312663e-03_jprb, &
 & 0.3981121e-03_jprb, 0.3601595e-03_jprb, 0.3408163e-03_jprb, 0.4272911e-03_jprb, 0.5288183e-03_jprb, &
 & 0.5513737e-03_jprb, 0.4471895e-03_jprb, 0.3864898e-03_jprb, 0.3626459e-03_jprb, 0.3489821e-03_jprb, &
 & 0.4398179e-03_jprb, 0.5966221e-03_jprb, 0.7405845e-03_jprb, 0.6992127e-03_jprb, 0.6250672e-03_jprb, &
 & 0.5602797e-03_jprb, 0.5163022e-03_jprb, 0.4840466e-03_jprb, 0.5107349e-03_jprb, 0.5462615e-03_jprb, &
 & 0.5881982e-03_jprb, 0.6573852e-03_jprb, 0.7159264e-03_jprb, 0.7562820e-03_jprb, 0.7965609e-03_jprb, &
 & 0.8163875e-03_jprb, 0.8190289e-03_jprb, 0.8226782e-03_jprb, 0.8284381e-03_jprb, 0.8345755e-03_jprb, &
 & 0.8056005e-03_jprb, 0.7495604e-03_jprb, 0.7066258e-03_jprb, 0.7724642e-03_jprb, 0.9057206e-03_jprb, &
 & 0.8502832e-03_jprb, 0.6413947e-03_jprb, 0.5219530e-03_jprb, 0.6165869e-03_jprb, 0.7758490e-03_jprb, &
 & 0.8803623e-03_jprb, 0.9650903e-03_jprb, 0.1049895e-02_jprb, 0.1079331e-02_jprb, 0.1114359e-02_jprb, &
 & 0.1041519e-02_jprb, 0.8810123e-03_jprb, 0.7704380e-03_jprb, 0.8720518e-03_jprb, 0.1053372e-02_jprb, &
 & 0.1148837e-02_jprb, 0.1189115e-02_jprb, 0.1224080e-02_jprb, 0.1139373e-02_jprb, 0.1039175e-02_jprb, &
 & 0.9717263e-03_jprb, 0.9659483e-03_jprb, 0.9509932e-03_jprb, 0.8842079e-03_jprb, 0.8138584e-03_jprb, &
 & 0.8657698e-03_jprb, 0.1011626e-02_jprb, 0.1128555e-02_jprb, 0.1102516e-02_jprb, 0.1043855e-02_jprb, &
 & 0.1035198e-02_jprb, 0.1093726e-02_jprb, 0.1145715e-02_jprb, 0.1101793e-02_jprb, 0.1023007e-02_jprb, &
 & 0.1049129e-02_jprb, 0.1125353e-02_jprb, 0.1193913e-02_jprb, 0.1135826e-02_jprb, 0.1054100e-02_jprb, &
 & 0.1004690e-02_jprb, 0.9727073e-03_jprb, 0.9889501e-03_jprb, 0.1206755e-02_jprb, 0.1413361e-02_jprb, &
 & 0.1483673e-02_jprb, 0.1389585e-02_jprb, 0.1303692e-02_jprb, 0.1202618e-02_jprb, 0.1067236e-02_jprb, &
 & 0.1081125e-02_jprb, 0.1193359e-02_jprb, 0.1297538e-02_jprb, 0.1362035e-02_jprb, 0.1438721e-02_jprb, &
 & 0.1498571e-02_jprb, 0.1605460e-02_jprb, 0.1667381e-02_jprb, 0.1561088e-02_jprb, 0.1404859e-02_jprb, &
 & 0.1300112e-02_jprb, 0.1232071e-02_jprb, 0.1186359e-02_jprb, 0.1317165e-02_jprb, 0.1565706e-02_jprb, &
 & 0.1757936e-02_jprb, 0.1998126e-02_jprb, 0.2219445e-02_jprb, 0.1934279e-02_jprb, 0.1663420e-02_jprb, &
 & 0.1494387e-02_jprb, 0.1478608e-02_jprb, 0.1506739e-02_jprb, 0.1668795e-02_jprb, 0.1919090e-02_jprb, &
 & 0.1843774e-02_jprb, 0.1546209e-02_jprb, 0.1246735e-02_jprb, 0.1458146e-02_jprb, 0.1869647e-02_jprb, &
 & 0.1982769e-02_jprb, 0.1804524e-02_jprb, 0.1709744e-02_jprb, 0.1773561e-02_jprb, 0.1859586e-02_jprb, &
 & 0.1826165e-02_jprb, 0.1656768e-02_jprb, 0.1557886e-02_jprb, 0.1705018e-02_jprb, 0.1915113e-02_jprb, &
 & 0.2075051e-02_jprb, 0.2126461e-02_jprb, 0.2189845e-02_jprb, 0.2264029e-02_jprb, 0.2330736e-02_jprb, &
 & 0.2142399e-02_jprb, 0.1641383e-02_jprb, 0.1314965e-02_jprb, 0.1430123e-02_jprb, 0.1671123e-02_jprb, &
 & 0.2000943e-02_jprb, 0.2398306e-02_jprb, 0.2805397e-02_jprb, 0.2760983e-02_jprb, 0.2603031e-02_jprb, &
 & 0.2254001e-02_jprb, 0.1707438e-02_jprb, 0.1358394e-02_jprb, 0.1617800e-02_jprb, 0.2065673e-02_jprb, &
 & 0.2466817e-02_jprb, 0.2769094e-02_jprb, 0.3017141e-02_jprb, 0.3025354e-02_jprb, 0.2981591e-02_jprb, &
 & 0.2718787e-02_jprb, 0.2352677e-02_jprb, 0.2036865e-02_jprb, 0.2123214e-02_jprb, 0.2251885e-02_jprb, &
 & 0.2545599e-02_jprb, 0.2963245e-02_jprb, 0.3421649e-02_jprb, 0.3761743e-02_jprb, 0.4168039e-02_jprb, &
 & 0.3792586e-02_jprb, 0.2691818e-02_jprb, 0.1859473e-02_jprb, 0.1698663e-02_jprb, 0.1596011e-02_jprb, &
 & 0.1866468e-02_jprb, 0.2472971e-02_jprb, 0.3382006e-02_jprb, 0.4041891e-02_jprb, 0.4838723e-02_jprb, &
 & 0.5173073e-02_jprb, 0.4874604e-02_jprb, 0.4306427e-02_jprb, 0.3231780e-02_jprb, 0.2000433e-02_jprb, &
 & 0.1861439e-02_jprb, 0.2225228e-02_jprb, 0.2649365e-02_jprb, 0.3034028e-02_jprb, 0.3543490e-02_jprb, &
 & 0.4232267e-02_jprb, 0.5139813e-02_jprb, 0.6346525e-02_jprb, 0.6436912e-02_jprb, 0.6250926e-02_jprb, &
 & 0.5587344e-02_jprb, 0.4808163e-02_jprb, 0.3922214e-02_jprb, 0.2997367e-02_jprb, 0.2173786e-02_jprb, &
 & 0.2030102e-02_jprb, 0.2242679e-02_jprb, 0.2624961e-02_jprb, 0.3565763e-02_jprb, 0.5075541e-02_jprb, &
 & 0.6281868e-02_jprb, 0.6886248e-02_jprb, 0.7266843e-02_jprb, 0.6076686e-02_jprb, 0.4795772e-02_jprb, &
 & 0.4356297e-02_jprb, 0.4466772e-02_jprb, 0.4258595e-02_jprb, 0.2897404e-02_jprb, 0.1843119e-02_jprb, &
 & 0.2009932e-02_jprb, 0.3199593e-02_jprb, 0.4951469e-02_jprb, 0.6041670e-02_jprb, 0.7461846e-02_jprb, &
 & 0.8546823e-02_jprb, 0.8710152e-02_jprb, 0.8794515e-02_jprb, 0.8287750e-02_jprb, 0.7636299e-02_jprb, &
 & 0.5824848e-02_jprb, 0.4316517e-02_jprb, 0.3105479e-02_jprb, 0.2465729e-02_jprb, 0.1831760e-02_jprb, &
 & 0.2086637e-02_jprb, 0.3358716e-02_jprb, 0.4856066e-02_jprb, 0.5922231e-02_jprb, 0.7184053e-02_jprb, &
 & 0.8372539e-02_jprb, 0.9001604e-02_jprb, 0.9555919e-02_jprb, 0.8104750e-02_jprb, 0.6525334e-02_jprb, &
 & 0.4792899e-02_jprb, 0.3547894e-02_jprb, 0.2745707e-02_jprb, 0.3341994e-02_jprb, 0.4115321e-02_jprb, &
 & 0.4877749e-02_jprb, 0.5827051e-02_jprb, 0.7167797e-02_jprb, 0.7492950e-02_jprb, 0.7709746e-02_jprb, &
 & 0.8481136e-02_jprb, 0.9578651e-02_jprb, 0.1096287e-01_jprb, 0.1099204e-01_jprb, 0.1081739e-01_jprb, &
 & 0.8763371e-02_jprb, 0.6211870e-02_jprb, 0.4962033e-02_jprb, 0.3850856e-02_jprb, 0.2970878e-02_jprb, &
 & 0.2839165e-02_jprb, 0.3280359e-02_jprb, 0.3735243e-02_jprb, 0.4371049e-02_jprb, 0.5292864e-02_jprb, &
 & 0.6320171e-02_jprb, 0.7608543e-02_jprb, 0.9269996e-02_jprb, 0.9811061e-02_jprb, 0.1011115e-01_jprb, &
 & 0.1063044e-01_jprb, 0.1127774e-01_jprb, 0.1068954e-01_jprb, 0.6261593e-02_jprb, 0.3448325e-02_jprb /)
zctuv1(301:600) = (/ &
 & 0.3132930e-02_jprb, 0.4512714e-02_jprb, 0.6235344e-02_jprb, 0.7328861e-02_jprb, 0.9068431e-02_jprb, &
 & 0.1026526e-01_jprb, 0.1043364e-01_jprb, 0.1075774e-01_jprb, 0.1175188e-01_jprb, 0.1299149e-01_jprb, &
 & 0.1349927e-01_jprb, 0.1335356e-01_jprb, 0.1329088e-01_jprb, 0.1371228e-01_jprb, 0.1416486e-01_jprb, &
 & 0.1314228e-01_jprb, 0.1087734e-01_jprb, 0.8150859e-02_jprb, 0.4523423e-02_jprb, 0.2050808e-02_jprb, &
 & 0.1681349e-02_jprb, 0.2845095e-02_jprb, 0.4677435e-02_jprb, 0.6798010e-02_jprb, 0.1090589e-01_jprb, &
 & 0.1337907e-01_jprb, 0.1275446e-01_jprb, 0.1242592e-01_jprb, 0.1394778e-01_jprb, 0.1606850e-01_jprb, &
 & 0.1529291e-01_jprb, 0.1153001e-01_jprb, 0.8916092e-02_jprb, 0.6729881e-02_jprb, 0.4765625e-02_jprb, &
 & 0.5167450e-02_jprb, 0.7806076e-02_jprb, 0.1046871e-01_jprb, 0.9795097e-02_jprb, 0.8494707e-02_jprb, &
 & 0.7421230e-02_jprb, 0.6606390e-02_jprb, 0.6140483e-02_jprb, 0.8835159e-02_jprb, 0.1348886e-01_jprb, &
 & 0.2018839e-01_jprb, 0.2876065e-01_jprb, 0.4566265e-01_jprb, 0.3430709e-01_jprb, 0.2582433e-01_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)
zctuv2(  1:300) = (/ &
 & 0.5449728e-05_jprb, 0.1073871e-05_jprb,-0.4617351e-05_jprb,-0.4662597e-05_jprb,-0.4245903e-05_jprb, &
 &-0.3501344e-05_jprb,-0.2327032e-05_jprb,-0.1160485e-05_jprb,-0.9958612e-06_jprb,-0.9097969e-06_jprb, &
 &-0.6708045e-06_jprb,-0.3015507e-06_jprb,-0.6392786e-08_jprb,-0.3021240e-06_jprb,-0.6042436e-06_jprb, &
 & 0.1495235e-06_jprb, 0.2286339e-05_jprb, 0.3876945e-05_jprb, 0.3135600e-05_jprb, 0.1962560e-05_jprb, &
 & 0.2706349e-06_jprb,-0.2130587e-05_jprb,-0.4417644e-05_jprb,-0.4455218e-05_jprb,-0.4247367e-05_jprb, &
 &-0.2599687e-05_jprb, 0.1002400e-05_jprb, 0.4013848e-05_jprb, 0.4290779e-05_jprb, 0.4239120e-05_jprb, &
 & 0.2425512e-05_jprb,-0.1615081e-05_jprb,-0.5240423e-05_jprb,-0.4114007e-05_jprb,-0.2223535e-05_jprb, &
 & 0.7356294e-06_jprb, 0.4315885e-05_jprb, 0.7074211e-05_jprb, 0.3583437e-05_jprb,-0.7265943e-06_jprb, &
 &-0.2152072e-05_jprb, 0.8515158e-06_jprb, 0.2627439e-05_jprb, 0.3447826e-05_jprb, 0.3982830e-05_jprb, &
 & 0.3435070e-05_jprb, 0.2089333e-05_jprb, 0.6670214e-06_jprb,-0.5065600e-06_jprb,-0.1649474e-05_jprb, &
 &-0.5022627e-06_jprb, 0.2711703e-05_jprb, 0.5396061e-05_jprb, 0.5160576e-05_jprb, 0.4445397e-05_jprb, &
 & 0.4654078e-05_jprb, 0.6341880e-05_jprb, 0.7375423e-05_jprb, 0.5832266e-05_jprb, 0.3559810e-05_jprb, &
 & 0.1810887e-05_jprb, 0.4885454e-06_jprb,-0.4923044e-06_jprb, 0.8285694e-06_jprb, 0.2424262e-05_jprb, &
 & 0.4601969e-05_jprb, 0.6532422e-05_jprb, 0.7979189e-05_jprb, 0.5542080e-05_jprb, 0.6684266e-06_jprb, &
 & 0.5030713e-06_jprb, 0.3990714e-05_jprb, 0.6064054e-05_jprb, 0.5041105e-05_jprb, 0.3097796e-05_jprb, &
 & 0.2510567e-05_jprb, 0.2758335e-05_jprb, 0.2725572e-05_jprb, 0.1902094e-05_jprb, 0.5449911e-06_jprb, &
 & 0.1825130e-05_jprb, 0.5139096e-05_jprb, 0.7489894e-05_jprb, 0.6139778e-05_jprb, 0.3478090e-05_jprb, &
 & 0.2995152e-05_jprb, 0.3705460e-05_jprb, 0.4955684e-05_jprb, 0.5840223e-05_jprb, 0.6808994e-05_jprb, &
 & 0.7637068e-05_jprb, 0.8072566e-05_jprb, 0.8472739e-05_jprb, 0.8785633e-05_jprb, 0.9082861e-05_jprb, &
 & 0.8783771e-05_jprb, 0.8007679e-05_jprb, 0.7444976e-05_jprb, 0.8667521e-05_jprb, 0.1001028e-04_jprb, &
 & 0.9730282e-05_jprb, 0.7358927e-05_jprb, 0.5605692e-05_jprb, 0.9000273e-05_jprb, 0.1483381e-04_jprb, &
 & 0.1546519e-04_jprb, 0.1342223e-04_jprb, 0.1142913e-04_jprb, 0.1142773e-04_jprb, 0.1189336e-04_jprb, &
 & 0.1269055e-04_jprb, 0.1521436e-04_jprb, 0.1622249e-04_jprb, 0.1450292e-04_jprb, 0.1258747e-04_jprb, &
 & 0.1118601e-04_jprb, 0.1046087e-04_jprb, 0.1042888e-04_jprb, 0.1318807e-04_jprb, 0.1757355e-04_jprb, &
 & 0.1569589e-04_jprb, 0.8876419e-05_jprb, 0.3477785e-05_jprb, 0.6092179e-05_jprb, 0.1169360e-04_jprb, &
 & 0.1280658e-04_jprb, 0.1156909e-04_jprb, 0.1077928e-04_jprb, 0.1150577e-04_jprb, 0.1265658e-04_jprb, &
 & 0.1216864e-04_jprb, 0.1071351e-04_jprb, 0.9272525e-05_jprb, 0.9172681e-05_jprb, 0.9266210e-05_jprb, &
 & 0.8021223e-05_jprb, 0.5337233e-05_jprb, 0.2790878e-05_jprb, 0.5715482e-05_jprb, 0.8468879e-05_jprb, &
 & 0.1234670e-04_jprb, 0.1693622e-04_jprb, 0.1945495e-04_jprb, 0.1515071e-04_jprb, 0.8224375e-05_jprb, &
 & 0.1201112e-04_jprb, 0.2331155e-04_jprb, 0.3526882e-04_jprb, 0.2870291e-04_jprb, 0.1563162e-04_jprb, &
 & 0.1140909e-04_jprb, 0.1557432e-04_jprb, 0.1762633e-04_jprb, 0.1500678e-04_jprb, 0.1178058e-04_jprb, &
 & 0.1138556e-04_jprb, 0.1392922e-04_jprb, 0.1566068e-04_jprb, 0.1500365e-04_jprb, 0.1378532e-04_jprb, &
 & 0.1177548e-04_jprb, 0.9777851e-05_jprb, 0.7775072e-05_jprb, 0.8716728e-05_jprb, 0.9806249e-05_jprb, &
 & 0.1316521e-04_jprb, 0.1904226e-04_jprb, 0.2291855e-04_jprb, 0.2240074e-04_jprb, 0.2072333e-04_jprb, &
 & 0.1811899e-04_jprb, 0.1469667e-04_jprb, 0.1133881e-04_jprb, 0.1322723e-04_jprb, 0.1647023e-04_jprb, &
 & 0.1916838e-04_jprb, 0.2205816e-04_jprb, 0.2387456e-04_jprb, 0.2335926e-04_jprb, 0.2209182e-04_jprb, &
 & 0.1890222e-04_jprb, 0.1381362e-04_jprb, 0.9271030e-05_jprb, 0.5131299e-05_jprb, 0.1653106e-05_jprb, &
 & 0.4911985e-05_jprb, 0.1272395e-04_jprb, 0.1980745e-04_jprb, 0.2038804e-04_jprb, 0.2024396e-04_jprb, &
 & 0.2059572e-04_jprb, 0.2124187e-04_jprb, 0.2100664e-04_jprb, 0.1594810e-04_jprb, 0.8547316e-05_jprb, &
 & 0.9426703e-05_jprb, 0.1696874e-04_jprb, 0.2096693e-04_jprb, 0.1660292e-04_jprb, 0.1010086e-04_jprb, &
 & 0.8478976e-05_jprb, 0.9815730e-05_jprb, 0.1096572e-04_jprb, 0.1022864e-04_jprb, 0.8534034e-05_jprb, &
 & 0.7089896e-05_jprb, 0.6877433e-05_jprb, 0.7622835e-05_jprb, 0.1016432e-04_jprb, 0.1179338e-04_jprb, &
 & 0.1105352e-04_jprb, 0.9418647e-05_jprb, 0.7047521e-05_jprb, 0.2836491e-05_jprb,-0.3273142e-05_jprb, &
 &-0.5378550e-05_jprb,-0.4812351e-05_jprb,-0.4587549e-05_jprb, 0.2554223e-06_jprb, 0.6261604e-05_jprb, &
 & 0.1238181e-04_jprb, 0.1535396e-04_jprb, 0.1632936e-04_jprb, 0.1245009e-04_jprb, 0.7953743e-05_jprb, &
 & 0.7041203e-05_jprb, 0.8162968e-05_jprb, 0.8718669e-05_jprb, 0.4926565e-05_jprb,-0.3380653e-05_jprb, &
 &-0.8394464e-05_jprb,-0.4490310e-05_jprb, 0.1038037e-05_jprb, 0.1348030e-04_jprb, 0.2415614e-04_jprb, &
 & 0.2358109e-04_jprb, 0.1487318e-04_jprb, 0.8895374e-05_jprb, 0.7548401e-05_jprb, 0.6078768e-05_jprb, &
 & 0.6853832e-05_jprb, 0.9282053e-05_jprb, 0.1039655e-04_jprb, 0.1021246e-04_jprb, 0.8594036e-05_jprb, &
 & 0.6369831e-05_jprb, 0.6854452e-05_jprb, 0.7271157e-05_jprb, 0.6951013e-05_jprb, 0.6097368e-05_jprb, &
 & 0.1010198e-04_jprb, 0.1231559e-04_jprb, 0.1322275e-04_jprb, 0.1709857e-04_jprb, 0.2124202e-04_jprb, &
 & 0.1979320e-04_jprb, 0.1126303e-04_jprb, 0.4723688e-06_jprb,-0.7984352e-06_jprb,-0.7715917e-07_jprb, &
 &-0.1499937e-05_jprb,-0.3995374e-05_jprb,-0.4728345e-05_jprb, 0.8784824e-05_jprb, 0.2100014e-04_jprb, &
 & 0.2723090e-04_jprb, 0.2794500e-04_jprb, 0.2646355e-04_jprb, 0.2066187e-04_jprb, 0.1373170e-04_jprb, &
 & 0.8347183e-05_jprb, 0.3109514e-05_jprb,-0.3910176e-05_jprb, 0.5399752e-05_jprb, 0.2236864e-04_jprb, &
 & 0.2737938e-04_jprb, 0.1373889e-04_jprb,-0.4865972e-05_jprb,-0.7329967e-05_jprb,-0.7027517e-05_jprb, &
 & 0.5835257e-05_jprb, 0.1728207e-04_jprb, 0.2252474e-04_jprb, 0.3336825e-04_jprb, 0.4236085e-04_jprb, &
 & 0.3230816e-04_jprb, 0.8722150e-05_jprb,-0.9830307e-05_jprb,-0.4446554e-05_jprb, 0.7582591e-05_jprb, &
 & 0.1569976e-04_jprb, 0.1360494e-04_jprb, 0.1223443e-04_jprb, 0.2941493e-04_jprb, 0.5059241e-04_jprb, &
 & 0.4548364e-04_jprb, 0.1569561e-04_jprb, 0.3601336e-05_jprb, 0.4300248e-04_jprb, 0.6058728e-04_jprb /)
zctuv2(301:600) = (/ &
 & 0.5357265e-04_jprb, 0.3029424e-04_jprb, 0.1513925e-05_jprb,-0.1018678e-04_jprb,-0.2847171e-04_jprb, &
 &-0.3043390e-04_jprb,-0.9641391e-05_jprb, 0.1707224e-04_jprb, 0.1344189e-04_jprb, 0.3707683e-05_jprb, &
 & 0.7901584e-05_jprb, 0.2753039e-04_jprb, 0.4065997e-04_jprb, 0.8414642e-05_jprb,-0.2501079e-04_jprb, &
 &-0.1891941e-04_jprb, 0.1456770e-04_jprb, 0.3651696e-04_jprb, 0.2586492e-04_jprb, 0.1050168e-04_jprb, &
 & 0.5334877e-06_jprb,-0.1962072e-04_jprb,-0.5089322e-04_jprb,-0.8441173e-04_jprb,-0.1526737e-03_jprb, &
 &-0.1566662e-03_jprb,-0.7777984e-04_jprb, 0.1420665e-04_jprb, 0.2271288e-05_jprb,-0.2627546e-04_jprb, &
 & 0.1128679e-04_jprb, 0.1177677e-03_jprb, 0.1880122e-03_jprb, 0.1648969e-03_jprb, 0.1299399e-03_jprb, &
 & 0.6704369e-04_jprb,-0.2064946e-04_jprb,-0.9502823e-04_jprb,-0.4572404e-04_jprb, 0.3763052e-04_jprb, &
 & 0.8176474e-04_jprb, 0.7521586e-04_jprb, 0.6640679e-04_jprb, 0.4663960e-04_jprb, 0.1536782e-05_jprb, &
 &-0.1030770e-03_jprb,-0.2760625e-03_jprb,-0.6169242e-03_jprb,-0.3990865e-03_jprb,-0.2450004e-03_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, &
 & 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb, 0.0000000e+00_jprb /)

!ZDK250 = (/ &
!ZDTUV1 = (/ &
!ZDTUV2 = (/ &


                                    
! revised 20051031
! TSI UV: 0.1057502799E+03 integrated over the 280-400 nm interval
zasuvb = (/ &                                                            
 & 0.1251796e+01_jprb, 0.1800686e+01_jprb, 0.3127783e+01_jprb, 0.3057467e+01_jprb, 0.3225982e+01_jprb, &
 & 0.3601007e+01_jprb, 0.3833149e+01_jprb, 0.3772373e+01_jprb, 0.3740471e+01_jprb, 0.5108631e+01_jprb, &
 & 0.4926397e+01_jprb, 0.4104165e+01_jprb, 0.4356608e+01_jprb, 0.5318757e+01_jprb, 0.5635120e+01_jprb, &
 & 0.4680448e+01_jprb, 0.5237057e+01_jprb, 0.6325728e+01_jprb, 0.5219252e+01_jprb, 0.5970084e+01_jprb, &
 & 0.4724881e+01_jprb, 0.4951663e+01_jprb, 0.5338753e+01_jprb, 0.6442023e+01_jprb /)

! TSI UV: 0.1057502799E+03 integrated over the 280-400 nm interval
zbsuvb = (/ &                                                            
 & 0.6848311e-01_jprb, 0.2498671e+00_jprb, 0.3945846e+00_jprb, 0.3491516e+00_jprb, 0.1897091e+00_jprb, &
 & 0.8228300e-01_jprb, 0.2997684e+00_jprb, 0.4743167e+00_jprb, 0.4136641e+00_jprb, 0.5306536e+00_jprb, &
 & 0.7163181e+00_jprb, 0.6737111e+00_jprb, 0.5750798e+00_jprb, 0.6513190e+00_jprb, 0.5113551e+00_jprb, &
 & 0.6462190e+00_jprb, 0.5669341e+00_jprb, 0.6326718e+00_jprb, 0.4581728e+00_jprb, 0.7534690e+00_jprb, &
 & 0.5087756e+00_jprb, 0.7125608e+00_jprb, 0.5832244e+00_jprb, 0.7415306e+00_jprb, 0.6798909e+00_jprb, &
 & 0.7188407e+00_jprb, 0.6775800e+00_jprb, 0.8080093e+00_jprb, 0.8375965e+00_jprb, 0.5589809e+00_jprb, &
 & 0.7169390e+00_jprb, 0.8965086e+00_jprb, 0.7290298e+00_jprb, 0.8052218e+00_jprb, 0.6854493e+00_jprb, &
 & 0.7742989e+00_jprb, 0.6145189e+00_jprb, 0.9867508e+00_jprb, 0.6230114e+00_jprb, 0.7737931e+00_jprb, &
 & 0.8929048e+00_jprb, 0.6839315e+00_jprb, 0.6537453e+00_jprb, 0.6464708e+00_jprb, 0.8634182e+00_jprb, &
 & 0.9398837e+00_jprb, 0.1117516e+01_jprb, 0.1016506e+01_jprb, 0.9115749e+00_jprb, 0.1123150e+01_jprb, &
 & 0.9910813e+00_jprb, 0.9517728e+00_jprb, 0.1039072e+01_jprb, 0.9489145e+00_jprb, 0.9955564e+00_jprb, &
 & 0.9388580e+00_jprb, 0.6909630e+00_jprb, 0.7680316e+00_jprb, 0.8123949e+00_jprb, 0.8939177e+00_jprb, &
 & 0.9223812e+00_jprb, 0.8122940e+00_jprb, 0.8338263e+00_jprb, 0.9790789e+00_jprb, 0.8090277e+00_jprb, &
 & 0.1097064e+01_jprb, 0.1042524e+01_jprb, 0.9627741e+00_jprb, 0.1277082e+01_jprb, 0.9393126e+00_jprb, &
 & 0.1239610e+01_jprb, 0.1066807e+01_jprb, 0.9717703e+00_jprb, 0.1129375e+01_jprb, 0.1227557e+01_jprb, &
 & 0.1157243e+01_jprb, 0.9689299e+00_jprb, 0.8114618e+00_jprb, 0.6022537e+00_jprb, 0.1140560e+01_jprb, &
 & 0.1067406e+01_jprb, 0.1023137e+01_jprb, 0.1176517e+01_jprb, 0.9235301e+00_jprb, 0.1046466e+01_jprb, &
 & 0.1342558e+01_jprb, 0.1303960e+01_jprb, 0.1248713e+01_jprb, 0.1120003e+01_jprb, 0.1310493e+01_jprb, &
 & 0.1126627e+01_jprb, 0.1383549e+01_jprb, 0.9327253e+00_jprb, 0.8592345e+00_jprb, 0.9171154e+00_jprb, &
 & 0.1090722e+01_jprb, 0.1083875e+01_jprb, 0.1342697e+01_jprb, 0.1462411e+01_jprb, 0.9903797e+00_jprb, &
 & 0.1264847e+01_jprb, 0.1143896e+01_jprb, 0.6993406e+00_jprb, 0.6519147e+00_jprb, 0.9648829e+00_jprb, &
 & 0.9177119e+00_jprb, 0.9760993e+00_jprb, 0.8564172e+00_jprb, 0.8915388e+00_jprb, 0.1309896e+01_jprb, &
 & 0.1252062e+01_jprb, 0.1518355e+01_jprb, 0.1037630e+01_jprb, 0.3918341e+00_jprb, 0.1138872e+01_jprb, &
 & 0.1384648e+01_jprb, 0.7875828e+00_jprb, 0.9430014e+00_jprb, 0.1651466e+01_jprb, 0.1675324e+01_jprb /)

! TSI UV: 0.1057502799E+03 integrated over the 280-400 nm interval
zcsuvb(  1:300) = (/ &
 & 0.1158150e-01_jprb, 0.4307253e-02_jprb, 0.8679530e-02_jprb, 0.1612620e-01_jprb, 0.2778863e-01_jprb, &
 & 0.3575870e-01_jprb, 0.3593124e-01_jprb, 0.5114494e-01_jprb, 0.6681306e-01_jprb, 0.6021919e-01_jprb, &
 & 0.8154579e-01_jprb, 0.6622351e-01_jprb, 0.8759850e-01_jprb, 0.7109950e-01_jprb, 0.8811727e-01_jprb, &
 & 0.7734752e-01_jprb, 0.6886194e-01_jprb, 0.7562198e-01_jprb, 0.5066989e-01_jprb, 0.7665027e-01_jprb, &
 & 0.5342708e-01_jprb, 0.5731519e-01_jprb, 0.3119192e-01_jprb, 0.2780394e-01_jprb, 0.1997097e-01_jprb, &
 & 0.1199952e-01_jprb, 0.3781212e-02_jprb, 0.1343422e-01_jprb, 0.2034198e-01_jprb, 0.3272607e-01_jprb, &
 & 0.4846868e-01_jprb, 0.5698028e-01_jprb, 0.7398895e-01_jprb, 0.4847769e-01_jprb, 0.7185277e-01_jprb, &
 & 0.9515611e-01_jprb, 0.1128814e+00_jprb, 0.9564205e-01_jprb, 0.7717521e-01_jprb, 0.9346189e-01_jprb, &
 & 0.4597404e-01_jprb, 0.3428929e-01_jprb, 0.1077923e+00_jprb, 0.1154901e+00_jprb, 0.1101183e+00_jprb, &
 & 0.9279595e-01_jprb, 0.9740494e-01_jprb, 0.7512583e-01_jprb, 0.1420005e+00_jprb, 0.1233264e+00_jprb, &
 & 0.1506065e+00_jprb, 0.1314347e+00_jprb, 0.1466181e+00_jprb, 0.1625762e+00_jprb, 0.1250826e+00_jprb, &
 & 0.1025410e+00_jprb, 0.1180653e+00_jprb, 0.1372793e+00_jprb, 0.1567150e+00_jprb, 0.1591106e+00_jprb, &
 & 0.1509978e+00_jprb, 0.1523881e+00_jprb, 0.1050467e+00_jprb, 0.8648758e-01_jprb, 0.8015965e-01_jprb, &
 & 0.1484698e+00_jprb, 0.1179774e+00_jprb, 0.1492430e+00_jprb, 0.1174614e+00_jprb, 0.1181674e+00_jprb, &
 & 0.1251424e+00_jprb, 0.1001071e+00_jprb, 0.9280812e-01_jprb, 0.1512694e+00_jprb, 0.4202820e-01_jprb, &
 & 0.1189204e+00_jprb, 0.1582238e+00_jprb, 0.9405126e-01_jprb, 0.1425538e+00_jprb, 0.1324698e+00_jprb, &
 & 0.1597341e+00_jprb, 0.1710007e+00_jprb, 0.8407424e-01_jprb, 0.5323930e-01_jprb, 0.9888584e-01_jprb, &
 & 0.6824519e-01_jprb, 0.1120127e+00_jprb, 0.1361283e+00_jprb, 0.1293499e+00_jprb, 0.1869357e+00_jprb, &
 & 0.1177258e+00_jprb, 0.9133044e-01_jprb, 0.4441319e-01_jprb, 0.9906802e-01_jprb, 0.1056354e+00_jprb, &
 & 0.1822143e+00_jprb, 0.1713972e+00_jprb, 0.6409213e-01_jprb, 0.1762541e+00_jprb, 0.1595112e+00_jprb, &
 & 0.4608663e-01_jprb, 0.6802742e-01_jprb, 0.1393778e+00_jprb, 0.1646635e+00_jprb, 0.9062016e-01_jprb, &
 & 0.1500694e+00_jprb, 0.1700569e+00_jprb, 0.1527697e+00_jprb, 0.1529162e+00_jprb, 0.8674857e-01_jprb, &
 & 0.4436884e-01_jprb, 0.1046963e+00_jprb, 0.1463353e+00_jprb, 0.8535615e-01_jprb, 0.2024678e+00_jprb, &
 & 0.1631171e+00_jprb, 0.1501661e+00_jprb, 0.1783627e+00_jprb, 0.1560721e+00_jprb, 0.9381256e-01_jprb, &
 & 0.1548443e+00_jprb, 0.1065238e+00_jprb, 0.1752939e+00_jprb, 0.1560214e+00_jprb, 0.8720748e-01_jprb, &
 & 0.1526506e+00_jprb, 0.2124056e+00_jprb, 0.1212607e+00_jprb, 0.1718601e+00_jprb, 0.6066382e-01_jprb, &
 & 0.1302736e+00_jprb, 0.1977617e+00_jprb, 0.1057830e+00_jprb, 0.1346538e+00_jprb, 0.1091080e+00_jprb, &
 & 0.1827814e+00_jprb, 0.1490962e+00_jprb, 0.1923628e+00_jprb, 0.1350661e+00_jprb, 0.1487027e+00_jprb, &
 & 0.1854675e+00_jprb, 0.1422262e+00_jprb, 0.1845453e+00_jprb, 0.1926982e+00_jprb, 0.1326593e+00_jprb, &
 & 0.9339384e-01_jprb, 0.7605023e-01_jprb, 0.1398333e+00_jprb, 0.1053312e+00_jprb, 0.1443723e+00_jprb, &
 & 0.4935243e-01_jprb, 0.1231221e+00_jprb, 0.1592260e+00_jprb, 0.1695701e+00_jprb, 0.2156683e+00_jprb, &
 & 0.1872401e+00_jprb, 0.2222135e+00_jprb, 0.1774713e+00_jprb, 0.1779839e+00_jprb, 0.1315997e+00_jprb, &
 & 0.1245990e+00_jprb, 0.1306088e+00_jprb, 0.1570354e+00_jprb, 0.1659742e+00_jprb, 0.1508124e+00_jprb, &
 & 0.1550767e+00_jprb, 0.1582761e+00_jprb, 0.1203409e+00_jprb, 0.1615870e+00_jprb, 0.2099411e+00_jprb, &
 & 0.1686777e+00_jprb, 0.1619672e+00_jprb, 0.5918952e-01_jprb, 0.1481797e+00_jprb, 0.1474352e+00_jprb, &
 & 0.1851919e+00_jprb, 0.1679120e+00_jprb, 0.1998436e+00_jprb, 0.1389546e+00_jprb, 0.8239684e-01_jprb, &
 & 0.1047535e+00_jprb, 0.9949701e-01_jprb, 0.1487465e+00_jprb, 0.1276272e+00_jprb, 0.1338947e+00_jprb, &
 & 0.1867179e+00_jprb, 0.2201701e+00_jprb, 0.1779457e+00_jprb, 0.2388429e+00_jprb, 0.1630742e+00_jprb, &
 & 0.1068201e+00_jprb, 0.8450688e-01_jprb, 0.1127901e+00_jprb, 0.1274841e+00_jprb, 0.1914103e+00_jprb, &
 & 0.2002285e+00_jprb, 0.9097867e-01_jprb, 0.1278886e+00_jprb, 0.1344432e+00_jprb, 0.2202541e+00_jprb, &
 & 0.1482008e+00_jprb, 0.1469951e+00_jprb, 0.2015719e+00_jprb, 0.1878379e+00_jprb, 0.2082990e+00_jprb, &
 & 0.1017584e+00_jprb, 0.1451019e+00_jprb, 0.7703697e-01_jprb, 0.1780647e+00_jprb, 0.1819695e+00_jprb, &
 & 0.1600806e+00_jprb, 0.1267505e+00_jprb, 0.1588537e+00_jprb, 0.9667448e-01_jprb, 0.1113860e+00_jprb, &
 & 0.9133881e-01_jprb, 0.1729763e+00_jprb, 0.1181825e+00_jprb, 0.1292425e+00_jprb, 0.1347307e+00_jprb, &
 & 0.1512866e+00_jprb, 0.1808238e+00_jprb, 0.1699806e+00_jprb, 0.2193174e+00_jprb, 0.1420098e+00_jprb, &
 & 0.1576642e+00_jprb, 0.1784694e+00_jprb, 0.1526431e+00_jprb, 0.2519739e+00_jprb, 0.1991331e+00_jprb, &
 & 0.2155590e+00_jprb, 0.2423780e+00_jprb, 0.2025711e+00_jprb, 0.2256757e+00_jprb, 0.2313318e+00_jprb, &
 & 0.1832516e+00_jprb, 0.2300989e+00_jprb, 0.2196651e+00_jprb, 0.2103950e+00_jprb, 0.1730957e+00_jprb, &
 & 0.1979747e+00_jprb, 0.1740626e+00_jprb, 0.1905000e+00_jprb, 0.1610193e+00_jprb, 0.1880183e+00_jprb, &
 & 0.2203149e+00_jprb, 0.2412104e+00_jprb, 0.2545475e+00_jprb, 0.1662324e+00_jprb, 0.2408452e+00_jprb, &
 & 0.2558326e+00_jprb, 0.1969197e+00_jprb, 0.2039628e+00_jprb, 0.1443392e+00_jprb, 0.1900270e+00_jprb, &
 & 0.2009350e+00_jprb, 0.1847773e+00_jprb, 0.1558351e+00_jprb, 0.1924818e+00_jprb, 0.2177436e+00_jprb, &
 & 0.2086823e+00_jprb, 0.1804881e+00_jprb, 0.2275202e+00_jprb, 0.2078572e+00_jprb, 0.2145247e+00_jprb, &
 & 0.1762646e+00_jprb, 0.2206831e+00_jprb, 0.2270376e+00_jprb, 0.1269940e+00_jprb, 0.1979352e+00_jprb, &
 & 0.2014534e+00_jprb, 0.2058836e+00_jprb, 0.2000303e+00_jprb, 0.1874826e+00_jprb, 0.2007064e+00_jprb, &
 & 0.1420998e+00_jprb, 0.2326092e+00_jprb, 0.2102334e+00_jprb, 0.2060209e+00_jprb, 0.1478948e+00_jprb, &
 & 0.6166791e-01_jprb, 0.9824482e-01_jprb, 0.2134419e+00_jprb, 0.1737989e+00_jprb, 0.1438094e+00_jprb, &
 & 0.9751125e-01_jprb, 0.9984105e-01_jprb, 0.1418594e+00_jprb, 0.2413629e+00_jprb, 0.1874570e+00_jprb, &
 & 0.1152467e+00_jprb, 0.1874475e+00_jprb, 0.1836579e+00_jprb, 0.1809469e+00_jprb, 0.1450959e+00_jprb, &
 & 0.2264130e+00_jprb, 0.1111494e+00_jprb, 0.1591455e+00_jprb, 0.1895232e+00_jprb, 0.2076865e+00_jprb /)
zcsuvb(301:600) = (/ &
 & 0.2530582e+00_jprb, 0.1691832e+00_jprb, 0.1673257e+00_jprb, 0.1651305e+00_jprb, 0.1676837e+00_jprb, &
 & 0.1801740e+00_jprb, 0.1678922e+00_jprb, 0.1142739e+00_jprb, 0.1588322e+00_jprb, 0.1911218e+00_jprb, &
 & 0.2137500e+00_jprb, 0.1589394e+00_jprb, 0.1462175e+00_jprb, 0.1521134e+00_jprb, 0.1628060e+00_jprb, &
 & 0.2311198e+00_jprb, 0.1948069e+00_jprb, 0.1652501e+00_jprb, 0.1950638e+00_jprb, 0.1928383e+00_jprb, &
 & 0.1185256e+00_jprb, 0.1315500e+00_jprb, 0.1481399e+00_jprb, 0.1733126e+00_jprb, 0.2374996e+00_jprb, &
 & 0.2374471e+00_jprb, 0.1786689e+00_jprb, 0.2087820e+00_jprb, 0.2633031e+00_jprb, 0.2088632e+00_jprb, &
 & 0.2220497e+00_jprb, 0.1871010e+00_jprb, 0.2358347e+00_jprb, 0.1865678e+00_jprb, 0.2109712e+00_jprb, &
 & 0.2466503e+00_jprb, 0.2067040e+00_jprb, 0.2054577e+00_jprb, 0.1018214e+00_jprb, 0.2021408e+00_jprb, &
 & 0.2510176e+00_jprb, 0.2619543e+00_jprb, 0.2186540e+00_jprb, 0.2692426e+00_jprb, 0.2762131e+00_jprb, &
 & 0.1443544e+00_jprb, 0.2331391e+00_jprb, 0.2068319e+00_jprb, 0.1624159e+00_jprb, 0.1925712e+00_jprb, &
 & 0.2315444e+00_jprb, 0.2675816e+00_jprb, 0.2617563e+00_jprb, 0.2397498e+00_jprb, 0.2389780e+00_jprb, &
 & 0.2017490e+00_jprb, 0.2164805e+00_jprb, 0.1560795e+00_jprb, 0.2334687e+00_jprb, 0.2590291e+00_jprb, &
 & 0.2519819e+00_jprb, 0.2069263e+00_jprb, 0.1131312e+00_jprb, 0.1501324e+00_jprb, 0.2495983e+00_jprb, &
 & 0.2326745e+00_jprb, 0.2538886e+00_jprb, 0.2257110e+00_jprb, 0.2008390e+00_jprb, 0.2162623e+00_jprb, &
 & 0.2815691e+00_jprb, 0.2071696e+00_jprb, 0.2931835e+00_jprb, 0.2567119e+00_jprb, 0.1889232e+00_jprb, &
 & 0.2792606e+00_jprb, 0.2335529e+00_jprb, 0.1652571e+00_jprb, 0.2136573e+00_jprb, 0.2655147e+00_jprb, &
 & 0.2021998e+00_jprb, 0.2492024e+00_jprb, 0.2237493e+00_jprb, 0.9668919e-01_jprb, 0.1970893e+00_jprb, &
 & 0.9243014e-01_jprb, 0.1704992e+00_jprb, 0.1867496e+00_jprb, 0.1740344e+00_jprb, 0.1877484e+00_jprb, &
 & 0.1485550e+00_jprb, 0.9281972e-01_jprb, 0.1346586e+00_jprb, 0.8005841e-01_jprb, 0.1461620e+00_jprb, &
 & 0.1703585e+00_jprb, 0.2748065e+00_jprb, 0.2135569e+00_jprb, 0.2206926e+00_jprb, 0.2611456e+00_jprb, &
 & 0.3002764e+00_jprb, 0.2001613e+00_jprb, 0.1852186e+00_jprb, 0.1820937e+00_jprb, 0.1996564e+00_jprb, &
 & 0.1438404e+00_jprb, 0.2192468e+00_jprb, 0.2426639e+00_jprb, 0.2698561e+00_jprb, 0.1475303e+00_jprb, &
 & 0.1369558e+00_jprb, 0.2283586e+00_jprb, 0.1959572e+00_jprb, 0.3457803e+00_jprb, 0.2694656e+00_jprb, &
 & 0.1617160e+00_jprb, 0.1227123e+00_jprb, 0.2226935e+00_jprb, 0.1845919e+00_jprb, 0.2318165e+00_jprb, &
 & 0.2967969e+00_jprb, 0.2018478e+00_jprb, 0.1869166e+00_jprb, 0.2005376e+00_jprb, 0.1603670e+00_jprb, &
 & 0.2054483e+00_jprb, 0.2567736e+00_jprb, 0.2826306e+00_jprb, 0.2930521e+00_jprb, 0.3046532e+00_jprb, &
 & 0.2922131e+00_jprb, 0.2425442e+00_jprb, 0.2202092e+00_jprb, 0.3031651e+00_jprb, 0.2458283e+00_jprb, &
 & 0.2005188e+00_jprb, 0.3134617e+00_jprb, 0.2558356e+00_jprb, 0.2571902e+00_jprb, 0.2217072e+00_jprb, &
 & 0.2084442e+00_jprb, 0.2559362e+00_jprb, 0.2557498e+00_jprb, 0.2190597e+00_jprb, 0.1808136e+00_jprb, &
 & 0.2219874e+00_jprb, 0.3175022e+00_jprb, 0.2378181e+00_jprb, 0.2659648e+00_jprb, 0.2672206e+00_jprb, &
 & 0.2852412e+00_jprb, 0.2738525e+00_jprb, 0.2187990e+00_jprb, 0.1703537e+00_jprb, 0.1783808e+00_jprb, &
 & 0.1845599e+00_jprb, 0.2636097e+00_jprb, 0.3560195e+00_jprb, 0.3115435e+00_jprb, 0.2678166e+00_jprb, &
 & 0.1147860e+00_jprb, 0.1107158e+00_jprb, 0.2578273e+00_jprb, 0.2698551e+00_jprb, 0.1795411e+00_jprb, &
 & 0.2748137e+00_jprb, 0.2236250e+00_jprb, 0.9111431e-01_jprb, 0.1095607e+00_jprb, 0.1601207e+00_jprb, &
 & 0.2477373e+00_jprb, 0.2842807e+00_jprb, 0.1589658e+00_jprb, 0.1025690e+00_jprb, 0.1235626e+00_jprb, &
 & 0.1114009e+00_jprb, 0.2359608e+00_jprb, 0.2825647e+00_jprb, 0.3238992e+00_jprb, 0.1368962e+00_jprb, &
 & 0.2139942e+00_jprb, 0.2490571e+00_jprb, 0.1949766e+00_jprb, 0.2244310e+00_jprb, 0.2014157e+00_jprb, &
 & 0.2007936e+00_jprb, 0.2693866e+00_jprb, 0.2962085e+00_jprb, 0.2980880e+00_jprb, 0.2782205e+00_jprb, &
 & 0.3581445e+00_jprb, 0.3251818e+00_jprb, 0.3322366e+00_jprb, 0.2409279e+00_jprb, 0.2059198e+00_jprb, &
 & 0.1942687e+00_jprb, 0.2585852e+00_jprb, 0.1726274e+00_jprb, 0.2225615e+00_jprb, 0.1423369e+00_jprb, &
 & 0.2277843e+00_jprb, 0.2619196e+00_jprb, 0.3190190e+00_jprb, 0.2141752e+00_jprb, 0.2419485e+00_jprb, &
 & 0.3032540e+00_jprb, 0.2357463e+00_jprb, 0.1822032e+00_jprb, 0.1373469e+00_jprb, 0.2853460e+00_jprb, &
 & 0.1031687e+00_jprb, 0.1750166e+00_jprb, 0.1704368e+00_jprb, 0.1019661e+00_jprb, 0.1487524e+00_jprb, &
 & 0.1650886e+00_jprb, 0.1164109e+00_jprb, 0.1054042e+00_jprb, 0.1594288e+00_jprb, 0.1055822e+00_jprb, &
 & 0.1176607e+00_jprb, 0.1849293e+00_jprb, 0.2224796e+00_jprb, 0.1783373e+00_jprb, 0.2614759e+00_jprb, &
 & 0.1500050e+00_jprb, 0.2488365e+00_jprb, 0.1675993e+00_jprb, 0.1911893e+00_jprb, 0.1600818e+00_jprb, &
 & 0.9722187e-01_jprb, 0.1919937e+00_jprb, 0.2495219e+00_jprb, 0.2255473e+00_jprb, 0.2118145e+00_jprb, &
 & 0.1620823e+00_jprb, 0.1034425e+00_jprb, 0.2278996e+00_jprb, 0.2416400e+00_jprb, 0.1213528e+00_jprb, &
 & 0.2125401e+00_jprb, 0.1247030e+00_jprb, 0.2620897e+00_jprb, 0.1179802e+00_jprb, 0.1742258e+00_jprb, &
 & 0.1820025e+00_jprb, 0.3155541e+00_jprb, 0.2765884e+00_jprb, 0.2444874e+00_jprb, 0.2912637e+00_jprb, &
 & 0.2198656e+00_jprb, 0.2843479e+00_jprb, 0.2241306e+00_jprb, 0.2015474e+00_jprb, 0.3221704e+00_jprb, &
 & 0.2833514e+00_jprb, 0.3601597e+00_jprb, 0.2911269e+00_jprb, 0.2972085e+00_jprb, 0.2865082e+00_jprb, &
 & 0.2065020e+00_jprb, 0.2025947e+00_jprb, 0.2778561e+00_jprb, 0.2096391e+00_jprb, 0.1410383e+00_jprb, &
 & 0.9857938e-01_jprb, 0.6269173e-01_jprb, 0.2530281e-01_jprb, 0.7657863e-01_jprb, 0.1286815e+00_jprb, &
 & 0.2086445e+00_jprb, 0.2141411e+00_jprb, 0.1877717e+00_jprb, 0.2910647e+00_jprb, 0.2372502e+00_jprb, &
 & 0.2744186e+00_jprb, 0.2735835e+00_jprb, 0.2960536e+00_jprb, 0.2737275e+00_jprb, 0.2668652e+00_jprb, &
 & 0.3166273e+00_jprb, 0.1598142e+00_jprb, 0.1693170e+00_jprb, 0.9992833e-01_jprb, 0.4189600e-01_jprb, &
 & 0.7145026e-01_jprb, 0.1279196e+00_jprb, 0.1842234e+00_jprb, 0.2720304e+00_jprb, 0.2873777e+00_jprb, &
 & 0.3626651e+00_jprb, 0.2933673e+00_jprb, 0.3022911e+00_jprb, 0.3212442e+00_jprb, 0.3718986e+00_jprb, &
 & 0.3289628e+00_jprb, 0.3477682e+00_jprb, 0.3787437e+00_jprb, 0.3112315e+00_jprb, 0.3086174e+00_jprb /)

!ZDSUVB = (/ &
                                               
!          Rayleigh cross-sections from Bodhaine et al., 1999
! NB: values in cm2, so tau = RAYUVB * Pressure * Avogadro / ( m_a * grav )
!     Pressure: N cm-1 s-2
!     Avogadro: Avog number (mol-1)
!     m_a     : dry air mean molecular weight ( g mol-1 )
!     grav    : gravity ( cm s-2) 

! revised 20051031
!Rayl: 0.3293000611E-25 sun-weighted over the 280-400 nm interval
zarayuv = (/ &                                                           
 & 0.7306844e-25_jprb, 0.6736580e-25_jprb, 0.6322517e-25_jprb, 0.5858758e-25_jprb, 0.5445114e-25_jprb, &
 & 0.5087498e-25_jprb, 0.4744089e-25_jprb, 0.4430921e-25_jprb, 0.4145744e-25_jprb, 0.3880453e-25_jprb, &
 & 0.3640141e-25_jprb, 0.3414811e-25_jprb, 0.3210803e-25_jprb, 0.3019948e-25_jprb, 0.2842767e-25_jprb, &
 & 0.2681506e-25_jprb, 0.2527365e-25_jprb, 0.2386328e-25_jprb, 0.2259065e-25_jprb, 0.2130931e-25_jprb, &
 & 0.2022461e-25_jprb, 0.1907690e-25_jprb, 0.1816129e-25_jprb, 0.1714296e-25_jprb /)

!Rayl: 0.3293000611E-25 sun-weighted over the 280-400 nm interval
zbrayuv = (/ &                                                           
 & 0.7553307e-25_jprb, 0.7443457e-25_jprb, 0.7335825e-25_jprb, 0.7226265e-25_jprb, 0.7125960e-25_jprb, &
 & 0.6991365e-25_jprb, 0.6898332e-25_jprb, 0.6798154e-25_jprb, 0.6684466e-25_jprb, 0.6591285e-25_jprb, &
 & 0.6496839e-25_jprb, 0.6396853e-25_jprb, 0.6313266e-25_jprb, 0.6216223e-25_jprb, 0.6126177e-25_jprb, &
 & 0.6032645e-25_jprb, 0.5953125e-25_jprb, 0.5853245e-25_jprb, 0.5775356e-25_jprb, 0.5693963e-25_jprb, &
 & 0.5606095e-25_jprb, 0.5534949e-25_jprb, 0.5444878e-25_jprb, 0.5378845e-25_jprb, 0.5302978e-25_jprb, &
 & 0.5231385e-25_jprb, 0.5155581e-25_jprb, 0.5083396e-25_jprb, 0.5011746e-25_jprb, 0.4939369e-25_jprb, &
 & 0.4866159e-25_jprb, 0.4809942e-25_jprb, 0.4740003e-25_jprb, 0.4675913e-25_jprb, 0.4614712e-25_jprb, &
 & 0.4555732e-25_jprb, 0.4489270e-25_jprb, 0.4430247e-25_jprb, 0.4366891e-25_jprb, 0.4312100e-25_jprb, &
 & 0.4253825e-25_jprb, 0.4196292e-25_jprb, 0.4146292e-25_jprb, 0.4089423e-25_jprb, 0.4035688e-25_jprb, &
 & 0.3982210e-25_jprb, 0.3931992e-25_jprb, 0.3881491e-25_jprb, 0.3831836e-25_jprb, 0.3782540e-25_jprb, &
 & 0.3735359e-25_jprb, 0.3685556e-25_jprb, 0.3638909e-25_jprb, 0.3593377e-25_jprb, 0.3547793e-25_jprb, &
 & 0.3502970e-25_jprb, 0.3455926e-25_jprb, 0.3412171e-25_jprb, 0.3372662e-25_jprb, 0.3331015e-25_jprb, &
 & 0.3291422e-25_jprb, 0.3249210e-25_jprb, 0.3210885e-25_jprb, 0.3170751e-25_jprb, 0.3128713e-25_jprb, &
 & 0.3093283e-25_jprb, 0.3055725e-25_jprb, 0.3020082e-25_jprb, 0.2981950e-25_jprb, 0.2946112e-25_jprb, &
 & 0.2911369e-25_jprb, 0.2875635e-25_jprb, 0.2842758e-25_jprb, 0.2809044e-25_jprb, 0.2775960e-25_jprb, &
 & 0.2742853e-25_jprb, 0.2711246e-25_jprb, 0.2676895e-25_jprb, 0.2647165e-25_jprb, 0.2615410e-25_jprb, &
 & 0.2587241e-25_jprb, 0.2555613e-25_jprb, 0.2524211e-25_jprb, 0.2495965e-25_jprb, 0.2469929e-25_jprb, &
 & 0.2439028e-25_jprb, 0.2412799e-25_jprb, 0.2385009e-25_jprb, 0.2358339e-25_jprb, 0.2331176e-25_jprb, &
 & 0.2306411e-25_jprb, 0.2278344e-25_jprb, 0.2251960e-25_jprb, 0.2230511e-25_jprb, 0.2205796e-25_jprb, &
 & 0.2178428e-25_jprb, 0.2154956e-25_jprb, 0.2130606e-25_jprb, 0.2108947e-25_jprb, 0.2085232e-25_jprb, &
 & 0.2061758e-25_jprb, 0.2039777e-25_jprb, 0.2017122e-25_jprb, 0.1995736e-25_jprb, 0.1972342e-25_jprb, &
 & 0.1952774e-25_jprb, 0.1930078e-25_jprb, 0.1910372e-25_jprb, 0.1890537e-25_jprb, 0.1869344e-25_jprb, &
 & 0.1849739e-25_jprb, 0.1830610e-25_jprb, 0.1811516e-25_jprb, 0.1791227e-25_jprb, 0.1772642e-25_jprb, &
 & 0.1754530e-25_jprb, 0.1739054e-25_jprb, 0.1715780e-25_jprb, 0.1699967e-25_jprb, 0.1682695e-25_jprb /)

!Rayl: 0.3293000611E-25 sun-weighted over the 280-400 nm interval
zcrayuv(  1:300) = (/ &
 & 0.7616638e-25_jprb, 0.7596217e-25_jprb, 0.7566579e-25_jprb, 0.7545073e-25_jprb, 0.7520894e-25_jprb, &
 & 0.7498396e-25_jprb, 0.7474992e-25_jprb, 0.7449660e-25_jprb, 0.7426870e-25_jprb, 0.7405151e-25_jprb, &
 & 0.7382129e-25_jprb, 0.7359474e-25_jprb, 0.7336698e-25_jprb, 0.7312532e-25_jprb, 0.7293129e-25_jprb, &
 & 0.7270750e-25_jprb, 0.7246929e-25_jprb, 0.7225889e-25_jprb, 0.7201327e-25_jprb, 0.7179669e-25_jprb, &
 & 0.7159132e-25_jprb, 0.7137091e-25_jprb, 0.7112686e-25_jprb, 0.7092680e-25_jprb, 0.7072331e-25_jprb, &
 & 0.7050776e-25_jprb, 0.7027431e-25_jprb, 0.7005233e-25_jprb, 0.6984647e-25_jprb, 0.6963897e-25_jprb, &
 & 0.6943879e-25_jprb, 0.6922317e-25_jprb, 0.6899903e-25_jprb, 0.6880721e-25_jprb, 0.6858853e-25_jprb, &
 & 0.6837090e-25_jprb, 0.6817655e-25_jprb, 0.6795759e-25_jprb, 0.6775398e-25_jprb, 0.6756201e-25_jprb, &
 & 0.6737114e-25_jprb, 0.6711798e-25_jprb, 0.6694924e-25_jprb, 0.6673852e-25_jprb, 0.6654869e-25_jprb, &
 & 0.6633041e-25_jprb, 0.6614539e-25_jprb, 0.6595648e-25_jprb, 0.6575899e-25_jprb, 0.6556557e-25_jprb, &
 & 0.6535321e-25_jprb, 0.6515833e-25_jprb, 0.6496196e-25_jprb, 0.6476328e-25_jprb, 0.6457961e-25_jprb, &
 & 0.6441038e-25_jprb, 0.6419320e-25_jprb, 0.6400914e-25_jprb, 0.6382210e-25_jprb, 0.6362626e-25_jprb, &
 & 0.6343807e-25_jprb, 0.6325313e-25_jprb, 0.6305017e-25_jprb, 0.6287837e-25_jprb, 0.6271078e-25_jprb, &
 & 0.6250013e-25_jprb, 0.6232995e-25_jprb, 0.6213345e-25_jprb, 0.6199295e-25_jprb, 0.6177485e-25_jprb, &
 & 0.6158079e-25_jprb, 0.6138453e-25_jprb, 0.6123282e-25_jprb, 0.6104834e-25_jprb, 0.6085159e-25_jprb, &
 & 0.6068199e-25_jprb, 0.6050634e-25_jprb, 0.6031437e-25_jprb, 0.6016957e-25_jprb, 0.5996980e-25_jprb, &
 & 0.5980282e-25_jprb, 0.5962179e-25_jprb, 0.5947632e-25_jprb, 0.5930616e-25_jprb, 0.5910388e-25_jprb, &
 & 0.5893103e-25_jprb, 0.5879395e-25_jprb, 0.5857085e-25_jprb, 0.5842694e-25_jprb, 0.5827528e-25_jprb, &
 & 0.5809260e-25_jprb, 0.5791326e-25_jprb, 0.5774750e-25_jprb, 0.5757388e-25_jprb, 0.5740871e-25_jprb, &
 & 0.5725365e-25_jprb, 0.5709284e-25_jprb, 0.5695706e-25_jprb, 0.5675827e-25_jprb, 0.5660968e-25_jprb, &
 & 0.5644708e-25_jprb, 0.5629105e-25_jprb, 0.5610676e-25_jprb, 0.5597107e-25_jprb, 0.5578470e-25_jprb, &
 & 0.5561693e-25_jprb, 0.5547664e-25_jprb, 0.5532939e-25_jprb, 0.5515983e-25_jprb, 0.5500726e-25_jprb, &
 & 0.5488026e-25_jprb, 0.5466552e-25_jprb, 0.5453631e-25_jprb, 0.5435469e-25_jprb, 0.5421854e-25_jprb, &
 & 0.5406856e-25_jprb, 0.5390386e-25_jprb, 0.5376519e-25_jprb, 0.5362250e-25_jprb, 0.5343698e-25_jprb, &
 & 0.5331479e-25_jprb, 0.5315746e-25_jprb, 0.5301434e-25_jprb, 0.5286928e-25_jprb, 0.5268597e-25_jprb, &
 & 0.5257406e-25_jprb, 0.5240588e-25_jprb, 0.5227215e-25_jprb, 0.5212194e-25_jprb, 0.5196392e-25_jprb, &
 & 0.5181440e-25_jprb, 0.5168025e-25_jprb, 0.5153224e-25_jprb, 0.5140691e-25_jprb, 0.5122811e-25_jprb, &
 & 0.5110644e-25_jprb, 0.5097039e-25_jprb, 0.5081770e-25_jprb, 0.5066004e-25_jprb, 0.5054125e-25_jprb, &
 & 0.5039219e-25_jprb, 0.5023992e-25_jprb, 0.5009734e-25_jprb, 0.4998271e-25_jprb, 0.4982583e-25_jprb, &
 & 0.4969272e-25_jprb, 0.4957173e-25_jprb, 0.4941637e-25_jprb, 0.4930264e-25_jprb, 0.4915094e-25_jprb, &
 & 0.4900671e-25_jprb, 0.4884935e-25_jprb, 0.4873834e-25_jprb, 0.4859938e-25_jprb, 0.4846769e-25_jprb, &
 & 0.4834488e-25_jprb, 0.4820086e-25_jprb, 0.4807528e-25_jprb, 0.4794680e-25_jprb, 0.4781788e-25_jprb, &
 & 0.4766772e-25_jprb, 0.4754919e-25_jprb, 0.4741631e-25_jprb, 0.4728048e-25_jprb, 0.4716431e-25_jprb, &
 & 0.4703830e-25_jprb, 0.4689923e-25_jprb, 0.4677535e-25_jprb, 0.4665251e-25_jprb, 0.4652005e-25_jprb, &
 & 0.4638987e-25_jprb, 0.4627368e-25_jprb, 0.4614312e-25_jprb, 0.4600859e-25_jprb, 0.4587120e-25_jprb, &
 & 0.4576507e-25_jprb, 0.4563778e-25_jprb, 0.4551894e-25_jprb, 0.4540353e-25_jprb, 0.4527886e-25_jprb, &
 & 0.4514944e-25_jprb, 0.4502190e-25_jprb, 0.4491759e-25_jprb, 0.4478702e-25_jprb, 0.4466893e-25_jprb, &
 & 0.4453446e-25_jprb, 0.4442118e-25_jprb, 0.4429979e-25_jprb, 0.4417654e-25_jprb, 0.4406395e-25_jprb, &
 & 0.4393794e-25_jprb, 0.4382254e-25_jprb, 0.4371032e-25_jprb, 0.4359890e-25_jprb, 0.4347318e-25_jprb, &
 & 0.4337158e-25_jprb, 0.4325351e-25_jprb, 0.4312833e-25_jprb, 0.4302253e-25_jprb, 0.4289430e-25_jprb, &
 & 0.4278239e-25_jprb, 0.4267910e-25_jprb, 0.4255298e-25_jprb, 0.4244454e-25_jprb, 0.4233537e-25_jprb, &
 & 0.4221707e-25_jprb, 0.4210167e-25_jprb, 0.4198555e-25_jprb, 0.4189473e-25_jprb, 0.4176732e-25_jprb, &
 & 0.4165337e-25_jprb, 0.4155153e-25_jprb, 0.4144353e-25_jprb, 0.4134188e-25_jprb, 0.4122111e-25_jprb, &
 & 0.4110894e-25_jprb, 0.4101121e-25_jprb, 0.4090171e-25_jprb, 0.4079204e-25_jprb, 0.4068998e-25_jprb, &
 & 0.4056811e-25_jprb, 0.4046853e-25_jprb, 0.4035664e-25_jprb, 0.4025730e-25_jprb, 0.4014375e-25_jprb, &
 & 0.4004650e-25_jprb, 0.3994781e-25_jprb, 0.3983609e-25_jprb, 0.3973505e-25_jprb, 0.3963121e-25_jprb, &
 & 0.3953201e-25_jprb, 0.3941969e-25_jprb, 0.3932154e-25_jprb, 0.3922008e-25_jprb, 0.3911375e-25_jprb, &
 & 0.3901542e-25_jprb, 0.3890660e-25_jprb, 0.3881347e-25_jprb, 0.3870995e-25_jprb, 0.3861014e-25_jprb, &
 & 0.3850985e-25_jprb, 0.3841576e-25_jprb, 0.3831229e-25_jprb, 0.3821930e-25_jprb, 0.3811757e-25_jprb, &
 & 0.3802202e-25_jprb, 0.3791729e-25_jprb, 0.3782096e-25_jprb, 0.3772830e-25_jprb, 0.3762521e-25_jprb, &
 & 0.3752196e-25_jprb, 0.3743804e-25_jprb, 0.3733339e-25_jprb, 0.3724655e-25_jprb, 0.3714239e-25_jprb, &
 & 0.3704893e-25_jprb, 0.3695292e-25_jprb, 0.3686195e-25_jprb, 0.3676398e-25_jprb, 0.3667089e-25_jprb, &
 & 0.3657894e-25_jprb, 0.3649280e-25_jprb, 0.3638438e-25_jprb, 0.3630178e-25_jprb, 0.3620675e-25_jprb, &
 & 0.3610680e-25_jprb, 0.3602397e-25_jprb, 0.3592707e-25_jprb, 0.3583906e-25_jprb, 0.3574758e-25_jprb, &
 & 0.3565372e-25_jprb, 0.3556517e-25_jprb, 0.3547120e-25_jprb, 0.3538869e-25_jprb, 0.3530207e-25_jprb, &
 & 0.3519870e-25_jprb, 0.3511801e-25_jprb, 0.3502476e-25_jprb, 0.3493818e-25_jprb, 0.3486291e-25_jprb, &
 & 0.3476539e-25_jprb, 0.3466662e-25_jprb, 0.3458988e-25_jprb, 0.3450374e-25_jprb, 0.3441917e-25_jprb, &
 & 0.3432756e-25_jprb, 0.3424206e-25_jprb, 0.3415836e-25_jprb, 0.3407089e-25_jprb, 0.3398823e-25_jprb, &
 & 0.3391132e-25_jprb, 0.3381723e-25_jprb, 0.3372699e-25_jprb, 0.3364352e-25_jprb, 0.3356600e-25_jprb, &
 & 0.3348421e-25_jprb, 0.3340288e-25_jprb, 0.3330966e-25_jprb, 0.3322651e-25_jprb, 0.3314748e-25_jprb /)
zcrayuv(301:600) = (/ &
 & 0.3305984e-25_jprb, 0.3298072e-25_jprb, 0.3289859e-25_jprb, 0.3282017e-25_jprb, 0.3273561e-25_jprb, &
 & 0.3265128e-25_jprb, 0.3257629e-25_jprb, 0.3249522e-25_jprb, 0.3241397e-25_jprb, 0.3233116e-25_jprb, &
 & 0.3225721e-25_jprb, 0.3217845e-25_jprb, 0.3209827e-25_jprb, 0.3202007e-25_jprb, 0.3193858e-25_jprb, &
 & 0.3186312e-25_jprb, 0.3177970e-25_jprb, 0.3169526e-25_jprb, 0.3162512e-25_jprb, 0.3154191e-25_jprb, &
 & 0.3147949e-25_jprb, 0.3139129e-25_jprb, 0.3130702e-25_jprb, 0.3123660e-25_jprb, 0.3115790e-25_jprb, &
 & 0.3108486e-25_jprb, 0.3101754e-25_jprb, 0.3092851e-25_jprb, 0.3085713e-25_jprb, 0.3078730e-25_jprb, &
 & 0.3070664e-25_jprb, 0.3062619e-25_jprb, 0.3055719e-25_jprb, 0.3048145e-25_jprb, 0.3040600e-25_jprb, &
 & 0.3033070e-25_jprb, 0.3025883e-25_jprb, 0.3018769e-25_jprb, 0.3011228e-25_jprb, 0.3004097e-25_jprb, &
 & 0.2996710e-25_jprb, 0.2989633e-25_jprb, 0.2982268e-25_jprb, 0.2974796e-25_jprb, 0.2967974e-25_jprb, &
 & 0.2961358e-25_jprb, 0.2953785e-25_jprb, 0.2945691e-25_jprb, 0.2939451e-25_jprb, 0.2931466e-25_jprb, &
 & 0.2925504e-25_jprb, 0.2918363e-25_jprb, 0.2911492e-25_jprb, 0.2904339e-25_jprb, 0.2896760e-25_jprb, &
 & 0.2890574e-25_jprb, 0.2883687e-25_jprb, 0.2876630e-25_jprb, 0.2869037e-25_jprb, 0.2862616e-25_jprb, &
 & 0.2855907e-25_jprb, 0.2848499e-25_jprb, 0.2842363e-25_jprb, 0.2836590e-25_jprb, 0.2828614e-25_jprb, &
 & 0.2822013e-25_jprb, 0.2815331e-25_jprb, 0.2808206e-25_jprb, 0.2801999e-25_jprb, 0.2795129e-25_jprb, &
 & 0.2788733e-25_jprb, 0.2781790e-25_jprb, 0.2775226e-25_jprb, 0.2768396e-25_jprb, 0.2761947e-25_jprb, &
 & 0.2755280e-25_jprb, 0.2748942e-25_jprb, 0.2742622e-25_jprb, 0.2736305e-25_jprb, 0.2729839e-25_jprb, &
 & 0.2723205e-25_jprb, 0.2716291e-25_jprb, 0.2710531e-25_jprb, 0.2703178e-25_jprb, 0.2697369e-25_jprb, &
 & 0.2690815e-25_jprb, 0.2684563e-25_jprb, 0.2678317e-25_jprb, 0.2671934e-25_jprb, 0.2666264e-25_jprb, &
 & 0.2660071e-25_jprb, 0.2652226e-25_jprb, 0.2647622e-25_jprb, 0.2640349e-25_jprb, 0.2634147e-25_jprb, &
 & 0.2628049e-25_jprb, 0.2622270e-25_jprb, 0.2615982e-25_jprb, 0.2610121e-25_jprb, 0.2603948e-25_jprb, &
 & 0.2597863e-25_jprb, 0.2592153e-25_jprb, 0.2585788e-25_jprb, 0.2579637e-25_jprb, 0.2574628e-25_jprb, &
 & 0.2567483e-25_jprb, 0.2561915e-25_jprb, 0.2555600e-25_jprb, 0.2549877e-25_jprb, 0.2545187e-25_jprb, &
 & 0.2537348e-25_jprb, 0.2531826e-25_jprb, 0.2526576e-25_jprb, 0.2520277e-25_jprb, 0.2514406e-25_jprb, &
 & 0.2508984e-25_jprb, 0.2502242e-25_jprb, 0.2497412e-25_jprb, 0.2491734e-25_jprb, 0.2485537e-25_jprb, &
 & 0.2479957e-25_jprb, 0.2474433e-25_jprb, 0.2468558e-25_jprb, 0.2462352e-25_jprb, 0.2456774e-25_jprb, &
 & 0.2451092e-25_jprb, 0.2445473e-25_jprb, 0.2440110e-25_jprb, 0.2434430e-25_jprb, 0.2428878e-25_jprb, &
 & 0.2423375e-25_jprb, 0.2418300e-25_jprb, 0.2412748e-25_jprb, 0.2407345e-25_jprb, 0.2401573e-25_jprb, &
 & 0.2395626e-25_jprb, 0.2390419e-25_jprb, 0.2384824e-25_jprb, 0.2379713e-25_jprb, 0.2374113e-25_jprb, &
 & 0.2368813e-25_jprb, 0.2363609e-25_jprb, 0.2357891e-25_jprb, 0.2352318e-25_jprb, 0.2346731e-25_jprb, &
 & 0.2342027e-25_jprb, 0.2336675e-25_jprb, 0.2330989e-25_jprb, 0.2326043e-25_jprb, 0.2320905e-25_jprb, &
 & 0.2315493e-25_jprb, 0.2310179e-25_jprb, 0.2305187e-25_jprb, 0.2299280e-25_jprb, 0.2294418e-25_jprb, &
 & 0.2288918e-25_jprb, 0.2284270e-25_jprb, 0.2279157e-25_jprb, 0.2274044e-25_jprb, 0.2269143e-25_jprb, &
 & 0.2263537e-25_jprb, 0.2258695e-25_jprb, 0.2253686e-25_jprb, 0.2248540e-25_jprb, 0.2243068e-25_jprb, &
 & 0.2238817e-25_jprb, 0.2233570e-25_jprb, 0.2228893e-25_jprb, 0.2223422e-25_jprb, 0.2217755e-25_jprb, &
 & 0.2213169e-25_jprb, 0.2208572e-25_jprb, 0.2203020e-25_jprb, 0.2198316e-25_jprb, 0.2194407e-25_jprb, &
 & 0.2188110e-25_jprb, 0.2183839e-25_jprb, 0.2178888e-25_jprb, 0.2174415e-25_jprb, 0.2169765e-25_jprb, &
 & 0.2164350e-25_jprb, 0.2159686e-25_jprb, 0.2154432e-25_jprb, 0.2150311e-25_jprb, 0.2144809e-25_jprb, &
 & 0.2141063e-25_jprb, 0.2135743e-25_jprb, 0.2131350e-25_jprb, 0.2126439e-25_jprb, 0.2121757e-25_jprb, &
 & 0.2116993e-25_jprb, 0.2112348e-25_jprb, 0.2107558e-25_jprb, 0.2103303e-25_jprb, 0.2098430e-25_jprb, &
 & 0.2093762e-25_jprb, 0.2089122e-25_jprb, 0.2084616e-25_jprb, 0.2079919e-25_jprb, 0.2075581e-25_jprb, &
 & 0.2070484e-25_jprb, 0.2066188e-25_jprb, 0.2061627e-25_jprb, 0.2057537e-25_jprb, 0.2052658e-25_jprb, &
 & 0.2048236e-25_jprb, 0.2043740e-25_jprb, 0.2039009e-25_jprb, 0.2034790e-25_jprb, 0.2030402e-25_jprb, &
 & 0.2026373e-25_jprb, 0.2021292e-25_jprb, 0.2017653e-25_jprb, 0.2012671e-25_jprb, 0.2008244e-25_jprb, &
 & 0.2003679e-25_jprb, 0.1999900e-25_jprb, 0.1995156e-25_jprb, 0.1990708e-25_jprb, 0.1986896e-25_jprb, &
 & 0.1982515e-25_jprb, 0.1977577e-25_jprb, 0.1973389e-25_jprb, 0.1969472e-25_jprb, 0.1965130e-25_jprb, &
 & 0.1961169e-25_jprb, 0.1956593e-25_jprb, 0.1952552e-25_jprb, 0.1948516e-25_jprb, 0.1944288e-25_jprb, &
 & 0.1939545e-25_jprb, 0.1935074e-25_jprb, 0.1931325e-25_jprb, 0.1927054e-25_jprb, 0.1922953e-25_jprb, &
 & 0.1919001e-25_jprb, 0.1914619e-25_jprb, 0.1910454e-25_jprb, 0.1906483e-25_jprb, 0.1902818e-25_jprb, &
 & 0.1898249e-25_jprb, 0.1894532e-25_jprb, 0.1889855e-25_jprb, 0.1886213e-25_jprb, 0.1882224e-25_jprb, &
 & 0.1877868e-25_jprb, 0.1874022e-25_jprb, 0.1869706e-25_jprb, 0.1865394e-25_jprb, 0.1861920e-25_jprb, &
 & 0.1857844e-25_jprb, 0.1854381e-25_jprb, 0.1849851e-25_jprb, 0.1846284e-25_jprb, 0.1842192e-25_jprb, &
 & 0.1838483e-25_jprb, 0.1834272e-25_jprb, 0.1830332e-25_jprb, 0.1826590e-25_jprb, 0.1822672e-25_jprb, &
 & 0.1818896e-25_jprb, 0.1814961e-25_jprb, 0.1810901e-25_jprb, 0.1807128e-25_jprb, 0.1803498e-25_jprb, &
 & 0.1799625e-25_jprb, 0.1795890e-25_jprb, 0.1791572e-25_jprb, 0.1788046e-25_jprb, 0.1784347e-25_jprb, &
 & 0.1780655e-25_jprb, 0.1776813e-25_jprb, 0.1773257e-25_jprb, 0.1769180e-25_jprb, 0.1765592e-25_jprb, &
 & 0.1761816e-25_jprb, 0.1758429e-25_jprb, 0.1754198e-25_jprb, 0.1750965e-25_jprb, 0.1747065e-25_jprb, &
 & 0.1743662e-25_jprb, 0.1739582e-25_jprb, 0.1736024e-25_jprb, 0.1732771e-25_jprb, 0.1729447e-25_jprb, &
 & 0.1724799e-25_jprb, 0.1721488e-25_jprb, 0.1718100e-25_jprb, 0.1714394e-25_jprb, 0.1710822e-25_jprb, &
 & 0.1707207e-25_jprb, 0.1703568e-25_jprb, 0.1700035e-25_jprb, 0.1696602e-25_jprb, 0.1692917e-25_jprb, &
 & 0.1689509e-25_jprb, 0.1686000e-25_jprb, 0.1682677e-25_jprb, 0.1679007e-25_jprb, 0.1675450e-25_jprb /)

!ZDRAYUV = (/


!          Liquid water cloud optical parameters from Slingo, 1989, JAS 46, 120
rasa = (/ 3.094e-02_jprb, 2.944e-02_jprb, 3.308e-02_jprb, 2.801e-02_jprb /)
rasb = (/ 1.252_jprb    , 1.270_jprb    , 1.246_jprb    , 1.293_jprb     /)
rasc = (/ 7.90e-07_jprb ,-6.50e-07_jprb ,-3.00e-07_jprb , 1.00e-06_jprb  /)
rasd = (/ 3.69e-07_jprb , 4.33e-07_jprb , 2.36e-07_jprb , 0.00e+00_jprb  /)
rase = (/ 0.844_jprb    , 0.841_jprb    , 0.839_jprb    , 0.836_jprb     /)
rasf = (/ 1.558e-03_jprb, 1.680e-03_jprb, 1.946e-03_jprb, 2.153e-03_jprb /)

!          Ice water cloud optical parameters from Fu, 1996, J. Clim., 9, 2070 &
!          extinction coefficient
rfa0 = (/-0.236447e-03_jprb,-0.266955e-03_jprb,-0.293599e-03_jprb,-0.258858e-03_jprb /)
rfa1 = (/ 0.253817e+01_jprb, 0.254179e+01_jprb, 0.254540e+01_jprb, 0.253815e+01_jprb /)
!                  single scattering albedo
rfb0 = (/-0.269916e-06_jprb,-0.100570e-06_jprb, 0.182299e-06_jprb,-0.198529e-06_jprb /)
rfb1 = (/ 0.212909e-06_jprb, 0.160441e-06_jprb, 0.112454e-06_jprb, 0.939480e-07_jprb /)
rfb2 = (/-0.265397e-09_jprb,-0.205663e-09_jprb,-0.108031e-09_jprb,-0.254540e-09_jprb /)
rfb3 = (/ 0.112983e-11_jprb, 0.888595e-12_jprb, 0.499801e-12_jprb, 0.110876e-11_jprb /)
!          assymetry factor
rfc0 = (/ 0.733260e+00_jprb, 0.737809e+00_jprb, 0.740926e+00_jprb, 0.743546e+00_jprb /)
rfc1 = (/ 0.918317e-03_jprb, 0.897515e-03_jprb, 0.890216e-03_jprb, 0.908674e-03_jprb /)
rfc2 = (/-0.422974e-06_jprb,-0.217099e-06_jprb,-0.187919e-06_jprb,-0.465326e-06_jprb /)
rfc3 = (/-0.107976e-07_jprb,-0.116090e-07_jprb,-0.116760e-07_jprb,-0.105786e-07_jprb /)
!          forward Delta function
rfd0 = (/ 0.119114e+00_jprb, 0.118072e+00_jprb, 0.117402e+00_jprb, 0.116918e+00_jprb /)
rfd1 = (/ 0.271746e-03_jprb, 0.271772e-03_jprb, 0.266303e-03_jprb, 0.259549e-03_jprb /)
rfd2 = (/ 0.167483e-05_jprb, 0.169253e-05_jprb, 0.177039e-05_jprb, 0.186749e-05_jprb /)
rfd3 = (/-0.114603e-07_jprb,-0.115557e-07_jprb,-0.118611e-07_jprb,-0.122603e-07_jprb /)

!          Aerosols (from GADS)
! continental average 80%
! maritime clean      95%
! desert              50% 
! urban               80% 
rtauva(1,:) = (/ & 
 & 2.2150_jprb , 1.0985_jprb , 1.1208_jprb , 2.2535_jprb , 1.03858_jprb , 1.12044_jprb /)  
rtauva(2,:) = (/ & 
 & 1.9353_jprb , 1.0754_jprb , 1.0955_jprb , 1.9756_jprb , 1.03858_jprb , 1.12044_jprb /)  
rtauva(3,:) = (/ & 
 & 1.6822_jprb , 1.0506_jprb , 1.0716_jprb , 1.7074_jprb , 1.03858_jprb , 1.12044_jprb /)  
rtauva(4,:) = (/ & 
 & 1.4664_jprb , 1.0347_jprb , 1.0499_jprb , 1.4806_jprb , 1.03858_jprb , 1.12044_jprb /)

rpiuva(1, :)= (/ &
 &  .8736_jprb , 0.9868_jprb , 0.6758_jprb , 0.7806_jprb, .9401905_jprb , 0.9989_jprb /)  
rpiuva(2, :)= (/ &
 &  .9253_jprb , 0.9969_jprb , 0.7080_jprb , 0.8226_jprb, .9401905_jprb , 0.9990_jprb /)  
rpiuva(3, :)= (/ &
 &  .9322_jprb , 0.9983_jprb , 0.7512_jprb , 0.8278_jprb, .9401905_jprb , 0.9990_jprb /)  
rpiuva(4, :)= (/ &
 &  .9315_jprb , 0.9985_jprb , 0.7849_jprb , 0.8259_jprb, .9401905_jprb , 0.9991_jprb /)  

rcguva(1, :)= (/ &
 & 0.7445_jprb , 0.8093_jprb , 0.8242_jprb , 0.7285_jprb, .7008249_jprb , 0.7567_jprb /)  
rcguva(2, :)= (/ &
 & 0.7317_jprb , 0.8061_jprb , 0.8068_jprb , 0.7137_jprb, .7008249_jprb , 0.7675_jprb /)  
rcguva(3, :)= (/ &
 & 0.7262_jprb , 0.8046_jprb , 0.7861_jprb , 0.7086_jprb, .7008249_jprb , 0.7743_jprb /)  
rcguva(4, :)= (/ &
 & 0.7209_jprb , 0.7993_jprb , 0.7700_jprb , 0.7040_jprb, .7008249_jprb , 0.7782_jprb /)  

!- Definition of the CIE erythemasl action spectrum

rxpo = (/ 0.000_jprb, -0.094_jprb, -0.015_jprb /)
rxpl = (/ 298.0_jprb,  328.0_jprb,  139.0_jprb /)

rk250(:) =0._jprb
rtuv1(:) =0._jprb
rtuv2(:) =0._jprb
rsuvb0(:)=0._jprb
rayuvb(:)=0._jprb
jcop(:)  =0
juvlam(:)=0

zres=120._jprb/kuv
zbeg=280._jprb-zres
zhres=zres*0.5_jprb

DO j=1,kuv
  ruvlam(j)=280._jprb+(j-1)*zres
  ipuv(j)  =int((ruvlam(j)-zbeg)/zres+0.0001_jprb)
  zuvlam=ruvlam(j)+zhres
  IF( zuvlam <= rxpl(1) ) THEN
    rcieas(j) = 10._jprb**rxpo(1)
  ELSEIF (zuvlam <= rxpl(2) ) THEN
    zexp= rxpo(2)*(zuvlam-rxpl(1))
    rcieas(j) = 10._jprb**zexp
  ELSE
    zexp= rxpo(3)*(zuvlam-rxpl(3))
    rcieas(j) = 10._jprb**zexp
  ENDIF
  IF (ruvlam(j) == 305._jprb .OR. ruvlam(j) == 310._jprb .OR. &
    & ruvlam(j) == 320._jprb .OR. ruvlam(j) == 380._jprb) juvlam(j) = 1

  WRITE(nulout,'("Eryth.Action Spectrum:",I3,F8.2,E12.5,F10.2,I3)') &
    & ipuv(j),zuvlam,rcieas(j),ruvlam(j),juvlam(j)
ENDDO

IF (kuv == 24) THEN
  DO j=1,kuv
    jcop(j)=iacop(j)
    rk250(j)=zak250(j)
    rtuv1(j)=zatuv1(j)
    rtuv2(j)=zatuv2(j)
    rsuvb0(j)=zasuvb(j)
    rayuvb(j)=zarayuv(j)
    juvlam(j)=1
  ENDDO
ELSEIF (kuv == 120) THEN
  DO j=1,kuv
    jcop(j)=ibcop(j)
    rk250(j)=zbk250(j)
    rtuv1(j)=zbtuv1(j)
    rtuv2(j)=zbtuv2(j)
    rsuvb0(j)=zbsuvb(j)
    rayuvb(j)=zbrayuv(j)
  ENDDO
ELSEIF (kuv == 600) THEN
  DO j=1,kuv
    jcop(j)=iccop(j)
    rk250(j)=zck250(j)
    rtuv1(j)=zctuv1(j)
    rtuv2(j)=zctuv2(j)
    rsuvb0(j)=zcsuvb(j)
    rayuvb(j)=zcrayuv(j)
  ENDDO
ELSE
  stop
ENDIF
 
   
WRITE(nulout,'("SU_UVRAD O3:",25I3)') (jcop(j),j=1,kuv)
DO j=1,kuv
  WRITE(nulout,'("SU_UVRAD O3:",I3,5E12.5)') j,rk250(j),rtuv1(j),rtuv2(j),rsuvb0(j),rayuvb(j)
ENDDO
DO j=1,4
  WRITE(nulout,'("CLOUDSLINGO:",I3,6E12.5)') j,rasa(j),rasb(j),rasc(j),rasd(j),rase(j),rasf(j)
ENDDO
DO j=1,4
  WRITE(nulout,'("ICECLDFUA B:",I3,6E12.5)') j,rfa0(j),rfa1(j),rfb0(j),rfb1(j),rfb2(j),rfb3(j)
ENDDO
DO j=1,4
  WRITE(nulout,'("ICECLDFUC D:",I3,8E12.5)') j,rfc0(j),rfc1(j),rfc2(j),rfc3(j),&
     &                                        rfd0(j),rfd1(j),rfd2(j),rfd3(j)
ENDDO
DO j=1,4
  DO jae=1,6
    WRITE(nulout,'("AEROSOLS   :",2I3,3E12.5)') j,jae,rtauva(j,jae),rpiuva(j,jae),rcguva(j,jae)
  ENDDO
ENDDO
  
!     ----------------------------------------------------------------

IF (lhook) CALL dr_hook('SU_UVRAD',1,zhook_handle)
END SUBROUTINE su_uvrad