92c92,126
< !nilo1  deleted..below
---
> 
>   ! parameters
>   !INTEGER, PARAMETER :: nband=8
> 
>   ! hk is the fractional extra-terrestrial solar flux in each
>   ! of the 8 bands. the sum of hk is 0.47074. (table 3)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(nband) :: hk3 = (/ &
>   !     0.00057_r8, 0.00367_r8, 0.00083_r8, 0.00417_r8, 0.00600_r8, &
>   !     0.00556_r8, 0.05913_r8, 0.39081_r8 /)
> 
>   ! zk is the ozone absorption coefficient. unit: /(cm-atm)stp (table 3)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(nband) :: zk3 = (/ &
>   !     30.47_r8, 187.2_r8, 301.9_r8, 42.83_r8, 7.09_r8, 1.25_r8, 0.0345_r8, 0.0572_r8 /)
> 
>   ! wk is the water vapor absorption coefficient. unit: cm**2/g (table 3)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(nband) :: wk3 = (/ 0.0_r8,0.0_r8,0.0_r8,0.0_r8,0.0_r8,0.0_r8,&
>   !                                                       0.0_r8,0.00075_r8/)
> 
>   ! ry is the extinction coefficient for rayleigh scattering. unit: /mb. (table 3)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(nband) :: ry3 = (/ &
>   !     0.00604_r8, 0.00170_r8, 0.00222_r8, 0.00132_r8, 0.00107_r8, 0.00091_r8, &
>   !     0.00055_r8, 0.00012_r8 /)
> 
>   ! coefficients for computing the extinction coefficients of ice,
>   ! water, and rain particles, independent of spectral band. (table 4)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(2) :: aib3 = (/  3.33e-4_r8, 2.52_r8 /)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(2) :: awb3 = (/ -6.59e-3_r8, 1.65_r8 /)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(2) :: arb3 = (/  3.07e-3_r8, 0.00_r8 /)
> 
>   ! coefficients for computing the asymmetry factor of ice, water,
>   ! and rain particles, independent of spectral band. (table 6)
> 
>   !REAL(KIND=r8), PARAMETER, DIMENSION(3) :: aig3 = (/ 0.74625_r8, 0.0010541_r8, -0.00000264_r8 /)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(3) :: awg3 = (/ 0.82562_r8, 0.0052900_r8, -0.00014866_r8 /)
>   !REAL(KIND=r8), PARAMETER, DIMENSION(3) :: arg3 = (/ 0.883_r8,   0.0_r8,      0.0_r8/)
1185d1218
< !nilo2
1187,1197c1220,1243
<      REAL(KIND=r8), PARAMETER, DIMENSION(11) :: tau_c = (/ &
<             2.432_r8,2.1_r8,2.1_r8,1.901_r8,1.818_r8,1.76_r8,1.562_r8, &
<             1.0_r8,0.568_r8, 0.281_r8,0.129_r8 /)
< 
<         REAL(KIND=r8), PARAMETER, DIMENSION(11) :: ssa_c = (/ &
<           0.653_r8,0.78_r8,0.78_r8,0.858_r8,0.886_r8,0.892_r8,0.903_r8, &
<         0.891_r8,0.836_r8,0.765_r8,0.701_r8 /)
< 
<        REAL(KIND=r8), PARAMETER, DIMENSION(11) :: asym_c = (/ &
<              0.726_r8,0.686_r8,0.686_r8,0.666_r8,0.658_r8,0.656_r8,0.65_r8, &
<             0.637_r8,0.632_r8,0.66_r8,0.776_r8 /)
---
>     !tar    REAL(KIND=r8), PARAMETER, DIMENSION(11) :: tau_c = (/ &
>     !tar         2.432_r8,2.1_r8,2.1_r8,1.901_r8,1.818_r8,1.76_r8,1.562_r8, &
>     !tar         1.0_r8,0.568_r8, 0.281_r8,0.129_r8 /)
> 
>     REAL(KIND=r8), PARAMETER, DIMENSION(8) :: tau_c = (/ &
>          2.09880_r8,1.76676_r8,1.08667_r8,0.87224_r8,0.56770_r8, &
>          0.25307_r8, 0.28098_r8,0.12890_r8 /)
> 
>     !tar     REAL(KIND=r8), PARAMETER, DIMENSION(11) :: ssa_c = (/ &
>     !tar          0.653_r8,0.78_r8,0.78_r8,0.858_r8,0.886_r8,0.892_r8,0.903_r8, &
>     !tar          0.891_r8,0.836_r8,0.765_r8,0.701_r8 /)
> 
>     REAL(KIND=r8), PARAMETER, DIMENSION(8) :: ssa_c = (/ &
>          0.77936_r8,0.89032_r8,0.89711_r8,0.86887_r8,0.83504_r8, &
>          0.75306_r8,0.76509_r8,0.69953_r8 /)
> 
>     !tar    REAL(KIND=r8), PARAMETER, DIMENSION(11) :: asym_c = (/ &
>     !tar         0.726_r8,0.686_r8,0.686_r8,0.666_r8,0.658_r8,0.656_r8,0.65_r8, &
>     !tar         0.637_r8,0.632_r8,0.66_r8,0.776_r8 /)
> 
>     REAL(KIND=r8), PARAMETER, DIMENSION(8) :: asym_c = (/ &
>          0.68684_r8,0.65619_r8,0.63822_r8,0.63599_r8,0.63209_r8, &
>          0.68143_r8,0.65978_r8,0.77775_r8 /)
> 
1201,1211c1247,1271
<          REAL(KIND=r8), PARAMETER, DIMENSION(11) :: tau_m = (/ &
<             1.33_r8,1.241_r8,1.241_r8,1.192_r8,1.173_r8,1.161_r8,1.117_r8, &
<             1.00_r8,0.906_r8,0.799_r8,0.603_r8 /)
< 
<          REAL(KIND=r8), PARAMETER, DIMENSION(11) :: ssa_m = (/ &
<               0.84_r8,0.921_r8,0.921_r8,0.962_r8,0.977_r8,0.98_r8,0.987_r8, &
<               0.989_r8,0.987_r8,0.986_r8,0.861_r8 /)
< 
<         REAL(KIND=r8), PARAMETER, DIMENSION(11) :: asym_m = (/ &
<              0.774_r8,0.757_r8,0.757_r8,0.746_r8,0.745_r8,0.744_r8,0.743_r8, &
<              0.743_r8,0.756_r8,0.772_r8,0.802_r8 /)
---
>     !tar     REAL(KIND=r8), PARAMETER, DIMENSION(11) :: tau_m = (/ &
>     !tar        1.33_r8,1.241_r8,1.241_r8,1.192_r8,1.173_r8,1.161_r8,1.117_r8, &
>     !tar        1.00_r8,0.906_r8,0.799_r8,0.603_r8 /)
> 
>     REAL(KIND=r8), PARAMETER, DIMENSION(8) :: tau_m = (/ &
>          1.23760_r8,1.16595_r8,1.01222_r8,0.97019_r8,0.90512_r8, &
>          0.76353_r8,0.80001_r8,0.60126_r8 /)
> 
> 
>     !tar     REAL(KIND=r8), PARAMETER, DIMENSION(11) :: ssa_m = (/ &
>     !tar          0.84_r8,0.921_r8,0.921_r8,0.962_r8,0.977_r8,0.98_r8,0.987_r8, &
>     !tar          0.989_r8,0.987_r8,0.986_r8,0.861_r8 /)
> 
>     REAL(KIND=r8), PARAMETER, DIMENSION(8) :: ssa_m = (/ &
>          0.92032_r8,0.97951_r8,0.99000_r8,0.98782_r8,0.98620_r8, &
>          0.96099_r8,0.98319_r8,0.86221_r8/)
> 
> 
>     !tar    REAL(KIND=r8), PARAMETER, DIMENSION(11) :: asym_m = (/ &
>     !tar         0.774_r8,0.757_r8,0.757_r8,0.746_r8,0.745_r8,0.744_r8,0.743_r8, &
>     !tar         0.743_r8,0.756_r8,0.772_r8,0.802_r8 /)
> 
>     REAL(KIND=r8), PARAMETER, DIMENSION(8) :: asym_m = (/ &
>          0.75748_r8,0.74395_r8,0.74600_r8,0.74975_r8,0.75541_r8, &
>          0.77743_r8,0.77184_r8,0.80229_r8 /)
1383d1442
< !nilo3----FIXED
2431,2432c2490
<     !nilo4
<      ! parameters
---
>     INTEGER, PARAMETER :: nband2 = 5
2434,2452c2492,2502
<     INTEGER, PARAMETER :: nk=10
<     INTEGER, PARAMETER :: nband2=3
< 
<     ! water vapor absorption coefficient for 10 k-intervals. unit: cm^2/gm (table 2)
<     REAL(KIND=r8), PARAMETER, DIMENSION(nk) :: xk2 = (/    &
<          0.0010_r8, 0.0133_r8, 0.0422_r8, 0.1334_r8, 0.4217_r8, &
<          1.334_r8,  5.623_r8,  31.62_r8,  177.8_r8,  1000.0_r8 /)
< 
<     ! water vapor k-distribution function,
<     ! the sum of hk is 0.52926. unit: fraction (table 2)
<     ! --- new coefficients (tarasova and fomin, 2000)
<     REAL(KIND=r8), PARAMETER, DIMENSION(nband2,nk) :: hk2 = RESHAPE( &
<          SOURCE = (/ &
<          0.19310_r8, 0.06924_r8, 0.00310_r8, 0.05716_r8, 0.01960_r8, 0.00637_r8, &
<          0.02088_r8, 0.00795_r8, 0.00526_r8, 0.02407_r8, 0.01716_r8, 0.00641_r8, &
<          0.01403_r8, 0.01118_r8, 0.00542_r8, 0.00582_r8, 0.01377_r8, 0.00312_r8, &
<          0.00246_r8, 0.02008_r8, 0.00368_r8, 0.00163_r8, 0.00265_r8, 0.00346_r8, &
<          0.00101_r8, 0.00282_r8, 0.00555_r8, 0.00041_r8, 0.00092_r8, 0.00098_r8 /) , &
<          SHAPE = (/nband2 , nk/) )
---
>     INTEGER, PARAMETER, DIMENSION(nband2) :: nk = (/ &
>          1,3,1,4,3 /)
>     !
>     !       
>     REAL(KIND=r8), PARAMETER, DIMENSION(nband2,4) :: hk2 = RESHAPE( &
>          SOURCE = (/ &  
>          0.006887_r8,0.20318_r8,0.01032_r8,0.09006_r8,0.01381_r8, &
>          0.0_r8,0.06211_r8,0.0_r8,0.0187_r8,0.009589_r8, &
>          0.0_r8,0.05358_r8,0.0_r8,0.01906_r8,0.01336_r8, &
>          0.0_r8,0.0_r8,0.0_r8,0.03427_r8,0.0_r8 /) , &
>          SHAPE = (/nband2 , 4/) )
2456,2457c2506,2507
<    REAL(KIND=r8), PARAMETER, DIMENSION(nband2) :: ry2 = (/ 0.0000156_r8, 0.0000018_r8, 0.000000_r8 /)
< 
---
>     REAL(KIND=r8), PARAMETER, DIMENSION(nband2) :: ry2 = (/ & 
>          0.5382_r8,0.5382_r8,0.0_r8,0.0598_r8,0.0_r8 /)  
2460c2510
<     ! ice, water, and rain particles (table 5) 
---
>     ! ice, water, and rain particles 
2463d2512
<   !nilo below from NEW!
2467c2516
<          2.52_r8,   2.52_r8,   2.52_r8, 2.52_r8, 2.52_r8 /) )
---
>          2.52_r8,   2.52_r8,   2.52_r8, 2.52_r8, 2.52_r8 /) )        
2471c2520
<          1.68_r8, 1.72_r8,   1.91_r8,    1.85_r8,   2.16_r8 /) )
---
>          1.68_r8, 1.72_r8,   1.91_r8,    1.85_r8,   2.16_r8 /) ) 
2475,2476c2524,2525
<          0.0_r8    , 0.0_r8    , 0.0_r8    , 0.0_r8  , 0.0_r8  /) )
< 
---
>          0.0_r8    , 0.0_r8    , 0.0_r8    , 0.0_r8  , 0.0_r8  /) )  
>     !
2481c2530
<    REAL(KIND=r8), PARAMETER, DIMENSION(nband2,3) :: aia2  = RESHAPE( &
---
>     REAL(KIND=r8), PARAMETER, DIMENSION(nband2,3) :: aia2  = RESHAPE( &
2483,2485c2532,2534
<          -0.00000260_r8,  0.00215346_r8,  0.08938331_r8, &
<          0.00000746_r8,  0.00073709_r8,  0.00299387_r8, &
<          0.00000000_r8, -0.00000134_r8, -0.00001038_r8 /) )
---
>          -.00000260_r8,-.00000260_r8, .00215346_r8, .00215346_r8, .08938331_r8, &
>          .00000746_r8,  .00000746_r8, .00073709_r8, .00073709_r8,  .00299387_r8, &
>          .00000000_r8, .00000000_r8,-.00000134_r8, -.00000134_r8, -.00001038_r8 /) )
2488,2490c2537,2539
<          0.00000007_r8,-0.00019934_r8, 0.01209318_r8, &
<          0.00000845_r8, 0.00088757_r8, 0.01784739_r8, &
<          -0.00000004_r8,-0.00000650_r8,-0.00036910_r8 /) )
---
>          .00000003_r8, .00000007_r8,.01209318_r8, -.00019934_r8, .01209318_r8, &
>          .00000354_r8,.00000845_r8,.01784739_r8, .00088757_r8, .01784739_r8, &
>          -.000000017_r8, -.00000004_r8,-.00036910_r8, -.00000650_r8, -.00036910_r8 /) ) 
2493,2500c2542,2549
<          0.029_r8,  0.342_r8,    0.466_r8, &
<          0.0000_r8,  0.000_r8,    0.000_r8, &
<          0.0000_r8,  0.000_r8,    0.000_r8 /) )
< 
< 
<     ! coefficients for computing the asymmetry factor of
<     !     ice, water, and rain particles
< 
---
>          .029_r8, .029_r8,      .342_r8,  .342_r8,  .466_r8, &
>          .0000_r8,  .0000_r8,     .000_r8,  .000_r8,  .000_r8, &
>          .0000_r8, .0000_r8,     .000_r8,  .000_r8,  .000_r8 /) ) 
>     !
>     !c-----coefficients for computing the asymmetry factor of 
>     !c     ice, water, and rain particles 
>     !
>     !       
2503,2506c2552,2554
<          0.74935228_r8, 0.76098937_r8, 0.84090400_r8, &
<          0.00119715_r8, 0.00141864_r8, 0.00126222_r8, &
<          -0.00000367_r8,-0.00000396_r8,-0.00000385_r8 /) )
< 
---
>          .74935228_r8, .74935228_r8, .84090400_r8, .76098937_r8, .84090400_r8, &
>          .00119715_r8, .00119715_r8, .00126222_r8 , .00141864_r8, .00126222_r8, &
>          -.00000367_r8,-.00000367_r8, -.00000385_r8 , -.00000396_r8, -.00000385_r8 /) )   
2509,2512c2557,2559
<          0.79375035_r8, 0.74513197_r8, 0.83530748_r8, &
<          0.00832441_r8, 0.01370071_r8, 0.00257181_r8, &
<          -0.00023263_r8,-0.00038203_r8, 0.00005519_r8 /) )
< 
---
>          .79375035_r8, .79375035_r8, .83530748_r8 ,.74513197_r8, .83530748_r8, &
>          .00832441_r8, .00832441_r8,.00257181_r8 ,.01370071_r8, .00257181_r8, &
>          -.00023263_r8, -.00023263_r8,.00005519_r8, -.00038203_r8,.00005519_r8  /) )
2515,2518c2562,2564
<          0.891_r8,  0.948_r8,    0.971_r8, &
<          0.0000_r8,  0.000_r8,    0.000_r8, &
<          0.0000_r8,  0.000_r8,    0.000_r8 /) )
< 
---
>          .891_r8,  .891_r8,      .948_r8,  .948_r8,    .971_r8, &
>          .0000_r8, .0000_r8,     .000_r8,  .000_r8,     .000_r8, &
>          .0000_r8, .0000_r8,     .000_r8,   .000_r8,    .000_r8 /) )
