* In ./src/core_init_atmosphere/mpas_init_atm_thompson_aerosols, get the dimension of

  index_nifa and index_nwfa from the pool lbc_state instead of state. Corrected the
  call to subroutine init_atm_thompson_aerosols_lbc accordingly.

* In ./src/core_init_atmosphere/mpas_init_atm_thompson_aerosols.F, replaced the index
  name of nifa and nwfa from index_nifa (index_nwfa) to index_lbc_nifa (index_lbc_nwfa)
  from the pool lbc_state.

Author: ldfowler58

src/core_init_atmosphere/mpas_init_atm_cases.F

@@ -331,7 +331,7 @@ subroutine init_atm_setup_case(domain, stream_manager)
                                       diag, lbc_state, block_ptr % dimensions, block_ptr % configs)
 
                call mpas_get_time(start_time, dateTimeString=timeStart)
-               call init_atm_thompson_aerosols_lbc(timeString, timeStart, block_ptr, mesh, diag, state, lbc_state)
+               call init_atm_thompson_aerosols_lbc(timeString, timeStart, block_ptr, mesh, diag, lbc_state)
 
                block_ptr => block_ptr % next
             end do

src/core_init_atmosphere/mpas_init_atm_thompson_aerosols.F

@@ -719,13 +719,12 @@ subroutine init_hinterp_gocart(block,mesh)
  end subroutine init_hinterp_gocart
 
 !=================================================================================================================
- subroutine init_atm_thompson_aerosols_lbc(timestamp,timestart,block,mesh,diag,state,lbc_state)
+ subroutine init_atm_thompson_aerosols_lbc(timestamp,timestart,block,mesh,diag,lbc_state)
 !=================================================================================================================
 
 !input arguments:
  character(len=StrKIND),intent(in):: timestart,timestamp
  type(mpas_pool_type),intent(in):: diag
- type(mpas_pool_type),intent(in):: state
 
 !inout arguments:
  type(block_type),intent(inout),target:: block
@@ -771,8 +770,8 @@ subroutine init_atm_thompson_aerosols_lbc(timestamp,timestart,block,mesh,diag,st
  call mpas_pool_get_dimension(mesh,'nVertLevels'  ,nVertLevels  )
  call mpas_pool_get_dimension(mesh,'nMonths'      ,nMonths      )
 
- call mpas_pool_get_dimension(state,'index_nifa',index_nifa)
- call mpas_pool_get_dimension(state,'index_nwfa',index_nwfa)
+ call mpas_pool_get_dimension(lbc_state,'index_lbc_nifa',index_nifa)
+ call mpas_pool_get_dimension(lbc_state,'index_lbc_nwfa',index_nwfa)
 
  call mpas_pool_get_array(diag,'pressure_base',pressure)