From joaho at statoil.com Thu Feb 2 11:36:59 2017 From: joaho at statoil.com (Joakim Hove) Date: Thu, 2 Feb 2017 11:36:59 +0000 Subject: [Opm] Possible test failures Message-ID: Hello; I have now merged these two PR's: https://github.com/OPM/opm-data/pull/179 https://github.com/OPM/opm-output/pull/160 Currently all normal tests pass, but I have not confirmed that the norne / polymer tests pass. If things fail you know where to direct the flak - and I will come in and try to fix the breakage (i.e. further updates of the reference data). Joakim ------------------------------------------------------------------- The information contained in this message may be CONFIDENTIAL and is intended for the addressee only. Any unauthorised use, dissemination of the information or copying of this message is prohibited. If you are not the addressee, please notify the sender immediately by return e-mail and delete this message. Thank you -------------- next part -------------- An HTML attachment was scrubbed... URL: From bobmerrill at mersep.com Thu Feb 9 15:32:29 2017 From: bobmerrill at mersep.com (bobmerrill at mersep.com) Date: Thu, 9 Feb 2017 10:32:29 -0500 (EST) Subject: [Opm] Flow UNRST binary file format Message-ID: <11745.182.69.141.209.1486654349.squirrel@www.mersep.com> Dear OPM Mailing List: Where can I find the binary format for the UNRST file? Resinsight is a good visualisation tool, but if I want to use Flow for numerical well testing, I need to extract the pressures of a block (or blocks) for every timestep and output the same to a well testing application. The “Results” pane in Resinsight provides a plot of the selected cell’s property vs. time, but I need to create a text file of that data with (time, value) pairs. I need the time values to be output with a precision of 1E-5 to 1E-4 days (approximately 1-10 seconds) for near wellbore effects. Once I get the binary format, I’ll write a Perl or Python script to read it and extract the necessary data. I'd be happy to share it, if it works. Many thanks. Regards, Bob Merrill From joaho at statoil.com Fri Feb 10 14:54:05 2017 From: joaho at statoil.com (Joakim Hove) Date: Fri, 10 Feb 2017 14:54:05 +0000 Subject: [Opm] Flow UNRST binary file format In-Reply-To: <11745.182.69.141.209.1486654349.squirrel@www.mersep.com> References: <11745.182.69.141.209.1486654349.squirrel@www.mersep.com> Message-ID: The reading and writing of binary eclipse files is handled by the ert library: https://github.com/Ensembles/ert which flow links to. The data in a binary Eclipse file are organised in "keywords" (mind you - these keywords are very different from the keywords used in the eclipse input file!). In the ert library a keyword is implemented with the C stucture ecl_kw: https://github.com/Ensembles/ert/blob/master/libecl/src/ecl_kw.c - and the full restart file is organized in the ecl_file structure: https://github.com/Ensembles/ert/blob/master/libecl/src/ecl_file.c So when extracting the data you are interested in you have (at least) two alternatives: 1. You can use the RestartIO class from the the opm-output and write a small C++ program - that way you will only be working with C++/Opm datastructures. 2. You can write an ert program using the low level ecl_xxx datastructures; if going this way you can use the python bindings. Personally I would use the ERT Python bindings; then your script could be something like this: #!/usr/bin/env python block_list = [(1,2,3) , (3,4,5)] # NB: These block numbers should start counting at zero. import sys from ert.ecl import EclGrid from ert.ecl import EclRestartFile case = sys.argv[1] grid = EclGrid("%s.EGRID" % case) rst_file = EclRestartFile("%s.UNRST" % case) with open("pressure.txt", "w") as f: time_list = rst_file.timeList( ) for index,time in enumerate( time_list ); view = rst_file.restartView( seqnum_index = index ) kw = view["PRESSURE"][0] f.write("%s" % time_list[2]) for (i,j,k) in block_list: f.write( "%g " % kw[i,j,k]) f.write("\n") Please just e-mail me directly if you have further questions. Joakim ________________________________ From: Opm on behalf of bobmerrill at mersep.com Sent: Thursday, February 9, 2017 4:32:29 PM To: opm at opm-project.org Subject: [Opm] Flow UNRST binary file format Dear OPM Mailing List: Where can I find the binary format for the UNRST file? Resinsight is a good visualisation tool, but if I want to use Flow for numerical well testing, I need to extract the pressures of a block (or blocks) for every timestep and output the same to a well testing application. The “Results” pane in Resinsight provides a plot of the selected cell’s property vs. time, but I need to create a text file of that data with (time, value) pairs. I need the time values to be output with a precision of 1E-5 to 1E-4 days (approximately 1-10 seconds) for near wellbore effects. Once I get the binary format, I’ll write a Perl or Python script to read it and extract the necessary data. I'd be happy to share it, if it works. Many thanks. Regards, Bob Merrill _______________________________________________ Opm mailing list Opm at opm-project.org http://opm-project.org/cgi-bin/mailman/listinfo/opm ------------------------------------------------------------------- The information contained in this message may be CONFIDENTIAL and is intended for the addressee only. Any unauthorised use, dissemination of the information or copying of this message is prohibited. If you are not the addressee, please notify the sender immediately by return e-mail and delete this message. Thank you -------------- next part -------------- An HTML attachment was scrubbed... URL: From joaho at statoil.com Fri Feb 10 15:11:59 2017 From: joaho at statoil.com (Joakim Hove) Date: Fri, 10 Feb 2017 15:11:59 +0000 Subject: [Opm] Flow UNRST binary file format In-Reply-To: References: <11745.182.69.141.209.1486654349.squirrel@www.mersep.com>, Message-ID: * You can use the RestartIO class from the the opm-output and write a small C++ program - that way you will only be working with C++/Opm datastructures. * You can write an ert program using the low level ecl_xxx datastructures; if going this way you can use the python bindings. A third option - would of course to be to modify opm-output, and add the calls to write your data into the RestartIO implemeentation - the RestartIO::save() https://github.com/OPM/opm-output/blob/master/opm/output/eclipse/RestartIO.cpp#L485 could be augmented with something like (untested again): { std::fstream out("pressure.txt" , ios_base::app); const auto& p = cells.data("PRESSURE") out << sim_time << " "; for (const auto& ijk : {{1,2,3} , {2,3,4}}) { size_t active_index = grid.getActiveIndex( ijk[0] , ijk[1] , ijk[2]); out << p[active_index]; } out.close( ); } That way you would get 'pressure.txt' file as part of your simulation - which is obviously neat, but you would have to maintain change to C++ code yourself, as I do not see us merging something special like this. I guess this will be the "best" solution - and the Python solution is the simplest. Your call. Joakim ------------------------------------------------------------------- The information contained in this message may be CONFIDENTIAL and is intended for the addressee only. Any unauthorised use, dissemination of the information or copying of this message is prohibited. If you are not the addressee, please notify the sender immediately by return e-mail and delete this message. Thank you -------------- next part -------------- An HTML attachment was scrubbed... URL: From tobk at dtu.dk Mon Feb 27 15:01:34 2017 From: tobk at dtu.dk (Tobias Kasper Skovborg Ritschel) Date: Mon, 27 Feb 2017 15:01:34 +0000 Subject: [Opm] Installation from binary on Ubuntu 16.04, undefined symbol Message-ID: Hi I have tried to install OPM from binary packages on Ubuntu 16.04 by following the installation instructions for Ubuntu 14.04 on opm-project.org (http://opm-project.org/?page_id=245). When I type "flow" on the command line or try to run the first tutorial (http://opm-project.org/?page_id=197) I get the following error flow: symbol lookup error: /usr/lib/libptscotch-5.1.so: undefined symbol: ompi_mpi_op_bor I have tried (afterwards) to upgrade prerequisite libraries as instructed in Note 3 on the installation instructions page, but the error remains. I would be very happy if you could help me with this. Let me know if I can provide any other useful information. Kind regards Tobias Ritschel -------------- next part -------------- An HTML attachment was scrubbed... URL: