def cmd():
sdat = StagyyData(conf.core.path)<EOL>conf.plates.plot = set_of_vars(conf.plates.plot)<EOL>if not conf.plates.vzcheck:<EOL><INDENT>conf.scaling.dimensional = True<EOL>conf.scaling.factors['<STR_LIT>'] = '<STR_LIT:M>'<EOL>main_plates(sdat)<EOL><DEDENT>else:<EOL><INDENT>seuil_memz = <NUM_LIT:0><EOL>nb_plates = []<EOL>time = []<EOL>ch2o = []<EOL>istart, iend = None, None<EOL>for step in sdat.walk.filter(fields=['<STR_LIT:T>']):<EOL><INDENT>if conf.plates.timeprofile:<EOL><INDENT>time.append(step.timeinfo.loc['<STR_LIT:t>'])<EOL>ch2o.append(step.timeinfo.loc[<NUM_LIT:0>])<EOL><DEDENT>istart = step.isnap if istart is None else istart<EOL>iend = step.isnap<EOL>limits, nphi, dvphi, seuil_memz, vphi_surf =detect_plates_vzcheck(step, seuil_memz)<EOL>water_profile = np.mean(step.fields['<STR_LIT>'][<NUM_LIT:0>, :, :, <NUM_LIT:0>], axis=<NUM_LIT:0>)<EOL>limits.sort()<EOL>sizeplates = [limits[<NUM_LIT:0>] + nphi - limits[-<NUM_LIT:1>]]<EOL>for lim in range(<NUM_LIT:1>, len(limits)):<EOL><INDENT>sizeplates.append(limits[lim] - limits[lim - <NUM_LIT:1>])<EOL><DEDENT>lim = len(limits) * [max(dvphi)]<EOL>fig = plt.figure()<EOL>plt.subplot(<NUM_LIT>)<EOL>plt.axis([<NUM_LIT:0>, nphi,<EOL>np.min(vphi_surf) * <NUM_LIT>, np.max(vphi_surf) * <NUM_LIT>])<EOL>plt.plot(vphi_surf)<EOL>plt.subplot(<NUM_LIT>)<EOL>plt.axis(<EOL>[<NUM_LIT:0>, nphi,<EOL>np.min(dvphi) * <NUM_LIT>, np.max(dvphi) * <NUM_LIT>])<EOL>plt.plot(dvphi)<EOL>plt.scatter(limits, lim, color='<STR_LIT>')<EOL>plt.subplot(<NUM_LIT>)<EOL>plt.hist(sizeplates, <NUM_LIT:10>, (<NUM_LIT:0>, nphi / <NUM_LIT:2>))<EOL>plt.subplot(<NUM_LIT>)<EOL>plt.plot(water_profile)<EOL>misc.saveplot(fig, '<STR_LIT>', step.isnap)<EOL>nb_plates.append(len(limits))<EOL><DEDENT>if conf.plates.timeprofile:<EOL><INDENT>for i in range(<NUM_LIT:2>, len(nb_plates) - <NUM_LIT:3>):<EOL><INDENT>nb_plates[i] = (nb_plates[i - <NUM_LIT:2>] + nb_plates[i - <NUM_LIT:1>] +<EOL>nb_plates[i] + nb_plates[i + <NUM_LIT:1>] +<EOL>nb_plates[i + <NUM_LIT:2>]) / <NUM_LIT:5><EOL><DEDENT>figt = plt.figure()<EOL>plt.subplot(<NUM_LIT>)<EOL>plt.axis([time[<NUM_LIT:0>], time[-<NUM_LIT:1>], <NUM_LIT:0>, np.max(nb_plates)])<EOL>plt.plot(time, nb_plates)<EOL>plt.subplot(<NUM_LIT>)<EOL>plt.plot(time, ch2o)<EOL>misc.saveplot(figt, '<STR_LIT>'.format(istart, iend))<EOL><DEDENT><DEDENT>
Implementation of plates subcommand.

    Other Parameters:
        conf.plates
        conf.scaling
        conf.plot
        conf.core