Enclosed you can find MATLAB codes and results for this study: https://doi.org/10.1038/s41598-020-80279-0 Two challenges to optimizing transcranial direct current stimulation (tDCS) are selecting between, often similar, electrode montages and accounting for inter-individual differences in response. These two factors are related by how tDCS montage determines current flow through the brain considered across or within individuals. MRI-based computational head models (CHMs) predict how brain anatomy determines electric field (EF) patterns for a given tDCS montage. Because conventional tDCS produces diffuse brain current flow, stimulation outcomes may be understood as modulation of global networks. Therefore, we developed a network-led, rather than region-led, approach. We specifically considered two common “frontal” tDCS montages that nominally target the dorsolateral prefrontal cortex; asymmetric “unilateral” (anode/cathode: F4/Fp1) and symmetric “bilateral” (F4/F3) electrode montages. We showed that cathode location significantly affects EFs in the limbic network. Furthermore, using a finer parcellation of large-scale networks, we found significant differences in some of the main nodes within a network, even if there is no difference at the network level.