 
         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
                {

                    dg::blas1::pointwiseDot( transferH2dX, volX2d, transferH2dX); //multiply by sqrt(g)
                    poloidal_average( transferH2dX, t1d, false); //average over eta
                    dg::blas1::scal( t1d, 4*M_PI*M_PI*f0); //
                    dg::blas1::copy( 0., fsa1d); //get rid of previous nan in fsa1d (nasty bug)
                    if( record_name[0] != 'j')
                        dg::blas1::pointwiseDivide( t1d, dvdpsip, fsa1d );
                    else
                        dg::blas1::copy( t1d, fsa1d);
                    //3. Interpolate fsa on 2d plane : <f>
                    dg::blas2::gemv(fsa2rzmatrix, fsa1d, transferH2d); //fsa on RZ grid
                }
