code: def set_base_initial_condition(model, monomer, value):
      
      # Build up monomer pattern dict
      sites_dict = {}
      for site in monomer.sites:
          if site in monomer.site_states:
              if site == 'loc' and 'cytoplasm' in monomer.site_states['loc']:
                  sites_dict['loc'] = 'cytoplasm'
              else:
                  sites_dict[site] = monomer.site_states[site][0]
          else:
              sites_dict[site] = None
      mp = monomer(**sites_dict)
      pname = monomer.name + '_0'
      try:
          p = model.parameters[pname]
          p.value = value
      except KeyError:
          p = Parameter(pname, value)
          model.add_component(p)
          model.initial(mp, p)
docs: Set an initial condition for a monomer in its 'default' state.