 
         name = record_name + "_fsa";
         long_name = record.long_name + " (Flux surface average.)";
+        err = nc_def_var( ncid_out, name.data(), NC_DOUBLE, 2, dim_ids1d,
+            &id1d[name]);
+        err = nc_put_att_text( ncid_out, id1d[name], "long_name", long_name.size(),
+            long_name.data());
+        name = record_name + "_std_fsa";
+        long_name = record.long_name + " (Flux surface average standard deviation on outboard midplane.)";
         err = nc_def_var( ncid_out, name.data(), NC_DOUBLE, 2, dim_ids1d,
             &id1d[name]);
         err = nc_put_att_text( ncid_out, id1d[name], "long_name", long_name.size(),

                    }
                    err = nc_put_vara_double( ncid_out, id0d.at(record_name+"_ifs_norm"),
                        start2d_out, count2d, &result );
                }
                else
                {

            long_name.data());

        name = record_name + "_cta2d";
        long_name = record.long_name + " (Convoluted toroidal average on 2d plane.)";
        err = nc_def_var( ncid_out, name.data(), NC_DOUBLE, 3, dim_ids,
            &id2d[name]);
