anchor.diff:                             cfc12, cfc22adj
   
   ! ------- Local -------
  -      integer(kind=im) :: jt, jp, igc, ipr, iprsm 
  +      integer(kind=im) :: jt, jp, igc, ipr, iprsm
         real(kind=rb) :: sumk, sumk1, sumk2, sumk3, sumk4, sumk5, sumf1, sumf2
   
   
  
positive.code: ! Input
        integer, intent(in) :: rrtmg_unit
  
  ! Local                                    
        character*80 errmess
        ! logical, external  :: wrf_dm_on_monitor
  
  !     Arrays fracrefao and fracrefbo are the Planck fractions for the lower
  !     and upper atmosphere.
  !     Planck fraction mapping levels: 
  !     Lower: P=1053.63 mb, T= 294.2 K
  !     Upper: P=0.353 mb, T = 262.11 K
  
  !     The array KAO contains absorption coefs at the 16 chosen g-values 
  !     for a range of pressure levels > ~100mb and temperatures.  The first
  !     index in the array, JT, which runs from 1 to 5, corresponds to 
  !     different temperatures.  More specifically, JT = 3 means that the 
  !     data are for the corresponding TREF for this  pressure level, 
  !     JT = 2 refers to the temperatureTREF-15, JT = 1 is for TREF-30, 
  !     JT = 4 is for TREF+15, and JT = 5 is for TREF+30.  The second 
  !     index, JP, runs from 1 to 13 and refers to the corresponding 
  !     pressure level in PREF (e.g. JP = 1 is for a pressure of 1053.63 mb).  
  !     The third index, IG, goes from 1 to 16, and tells us which 
  !     g-interval the absorption coefficients are for.
  
  !     The array KBO contains absorption coefs at the 16 chosen g-values 
  !     for a range of pressure levels < ~100mb and temperatures. The first 
  !     index in the array, JT, which runs from 1 to 5, corresponds to 
  !     different temperatures.  More specifically, JT = 3 means that the 
  !     data are for the reference temperature TREF for this pressure 
  !     level, JT = 2 refers to the temperature TREF-15, JT = 1 is for
  !     TREF-30, JT = 4 is for TREF+15, and JT = 5 is for TREF+30.  
  !     The second index, JP, runs from 13 to 59 and refers to the JPth
  !     reference pressure level (see taumol.f for the value of these
  !     pressure levels in mb).  The third index, IG, goes from 1 to 16,
  !     and tells us which g-interval the absorption coefficients are for.
  
  !     The array KAO_Mxx contains the absorption coefficient for 
  !     a minor species at the 16 chosen g-values for a reference pressure
  !     level below 100~ mb.   The first index refers to temperature 
  !     in 7.2 degree increments.  For instance, JT = 1 refers to a 
  !     temperature of 188.0, JT = 2 refers to 195.2, etc. The second index 
  !     runs over the g-channel (1 to 16).
  
  !     The array KBO_Mxx contains the absorption coefficient for 
  !     a minor species at the 16 chosen g-values for a reference pressure
  !     level above 100~ mb.   The first index refers to temperature 
  !     in 7.2 degree increments.  For instance, JT = 1 refers to a 
  !     temperature of 188.0, JT = 2 refers to 195.2, etc. The second index 
  !     runs over the g-channel (1 to 16).
  
  !     The array FORREFO contains the coefficient of the water vapor
  !     foreign-continuum (including the energy term).  The first 
  !     index refers to reference temperature (296,260,224,260) and 
  !     pressure (970,475,219,3 mbar) levels.  The second index 
  !     runs over the g-channel (1 to 16).
  
  !     The array SELFREFO contains the coefficient of the water vapor
  
negative.code:   ! --------------------
    function initialize_scalar(seed) result(twister)
      integer(kind=im),       intent(in   ) :: seed
      type(randomNumberSequence)                :: twister
  
      integer(kind=im) :: i
      ! See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. In the previous versions,
      !   MSBs of the seed affect only MSBs of the array state[].
      !   2002/01/09 modified by Makoto Matsumoto
  
      twister%state(0) = iand(seed, -1_im)
      do i = 1,  blockSize - 1 ! ubound(twister%state)
         twister%state(i) = 1812433253_im * ieor(twister%state(i-1), &
  