Plotting conditional_effects in brms

General
1

I have a model fitted in brms. I am trying to plot the conditional effects,

mod2<-
  bf(
    PCFBPAGI ~ 
      E0 - (((expIPRED^G) * Emax)/((expIPRED^G) +EC50^G)), 
    E0   ~ 1 + (1|i|NID),
    Emax ~ 1 + (1|i|NID),
    EC50 ~ 1 + (1|i|NID),
    G    ~ 1 + (1|i|NID),
    nl=TRUE
  )

fit2<-brm(mod2,
          data = subdat,
          prior = prior_1,
          chains = 3,
          future=T,
          seed=2021,
          iter = 4000,
          warmup = 1000)

conditions <- data.frame(NID = unique(subdat$NID))
rownames(conditions) <- unique(subdat$NID)
me_subdat<-
  conditional_effects(
    fit2, 
    effects = "expIPRED",
    conditions = conditions,
    re_formula = NULL,
    method = "predict"
  )

plot(me_subdat, ncol = 5, points = TRUE)

I have 102 unique values of NID. As you see above, I am using ncol=5, the plot I get is trying to dispaly all panels (all 102 NID values) in one page.

According to hte help for plot in brms ncol is described as:
"Number of plots to display per column for each effect. If NULL (default), ncol is computed internally based on the number of rows of conditions.

But, plot I get is putting all 102 subpanels in 5 columns in one page. How do I define for example 5 columns and 5 rows for the plot to have 25 panels in each page?
Do I have anything wrong in the model defined above that makes the plotting incorrect?
Thanks for help in advance