Hi, Pedro! I wold suggest to replace if(mod(istep,iout1d) == 0) by if(iout1d > 0.and.mod(istep,iout1d) == 0), and also for iout0d, 2d, 3d, save. When no output is needed, we can just set iout1d=0 in input.nml. This should be more user-freindly than setting very large numbers for no output.
Hi, Pedro! I wold suggest to replace if(mod(istep,iout1d) == 0) by if(iout1d > 0.and.mod(istep,iout1d) == 0), and also for iout0d, 2d, 3d, save. When no output is needed, we can just set iout1d=0 in input.nml. This should be more user-freindly than setting very large numbers for no output.