Method do_intitialization

SRC/domain/pattern/drm/H5DRMLoadPattern.cpp:320–722

void H5DRMLoadPattern::do_intitialization()
{
    Domain *theDomain = this->getDomain();

    if (theDomain == 0)
    {
        return;
    }

    if (MPI_local_rank == 0)
    {
        H5DRMout << "initializing - filename : " << dataset_fname << "\n";
    }

    //===========================================================================
    // Open the specified file for read only
    //===========================================================================
    hid_t file_access_plist   = H5Pcreate(H5P_FILE_ACCESS);


    ih5_fname = H5Fopen(dataset_fname.c_str(), H5F_ACC_RDONLY, file_access_plist);

    if (ih5_fname < 0)
    {
        H5DRMerror << "Error opening HDF5 file: " << dataset_fname << endl;
        return;
    }


    //===========================================================================
    // MAP DRM data to local domain
    //===========================================================================

    H5DRMout << "Initialize" << endln;

    ID internal;
    Matrix xyz;
    Vector drmbox_x0;
    read_double_dataset_into_vector(ih5_fname, "DRM_Metadata/drmbox_x0", drmbox_x0);
    read_double_dataset_into_matrix(ih5_fname, "DRM_Data/xyz", xyz);
    read_int_dataset_into_id(ih5_fname, "DRM_Data/internal", internal);
    read_int_dataset_into_id(ih5_fname, "DRM_Data/data_location", station_id2data_pos);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/drmbox_xmax", drmbox_xmax);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/drmbox_xmin", drmbox_xmin);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/drmbox_ymax", drmbox_ymax);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/drmbox_ymin", drmbox_ymin);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/drmbox_zmax", drmbox_zmax);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/drmbox_zmin", drmbox_zmin);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/dt", dt);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/tstart", tstart);
    read_scalar_double_dataset_into_double(ih5_fname, "DRM_Metadata/tend", tend);

    N_timesteps = int((tend - tstart) / dt);

    last_integration_time = theDomain->getCurrentTime();

    int NDRM_points = xyz.noRows();