Runtime Error: unitsystem cgs unrecognized
Hey all! When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true. Best wishes, Patrick I add the field with: def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}") and then call: rad_profile.add_fields("magB") Traceback: File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
Hi Patrick, How old is your file, by chance? FLASH3? FLASH2? Also, if you don't mind, can you paste the script somewhere too? Thanks, John Z On Nov 15, 2012, at 3:41 PM, Patrick Rieser wrote:
Hey all!
When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true.
Best wishes, Patrick
I add the field with:
def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}")
and then call:
rad_profile.add_fields("magB")
Traceback:
File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi John, The files are from FLASH3. Ah and by the way, I am using the yt-dev version, updated some hours ago. Here is the main script: http://paste.yt-project.org/show/2877/ and here the field is defined: http://paste.yt-project.org/show/2878/ I did split my program in way too many modules, so if I should upload all scripts somewhere, just tell me! Best wishes, Patrick Am 15.11.2012 21:43, schrieb John ZuHone:
Hi Patrick,
How old is your file, by chance? FLASH3? FLASH2?
Also, if you don't mind, can you paste the script somewhere too?
Thanks,
John Z
On Nov 15, 2012, at 3:41 PM, Patrick Rieser wrote:
Hey all!
When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true.
Best wishes, Patrick
I add the field with:
def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}")
and then call:
rad_profile.add_fields("magB")
Traceback:
File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Patrick, Ok, so I have reproduced your error. It's very odd, actually. I also do FLASH simulations with magnetic fields, and I have the unitsystem = 'none' in my sims, and yt seems to catch this ok, but it's not catching the 'cgs' in yours. I can verify that I get this error whether I try to use your field or the magnetic fields already defined in yt. I'll see what I can do about this. Best, John Z On Nov 15, 2012, at 3:52 PM, Patrick Rieser wrote:
Hi John,
The files are from FLASH3. Ah and by the way, I am using the yt-dev version, updated some hours ago.
Here is the main script: http://paste.yt-project.org/show/2877/
and here the field is defined: http://paste.yt-project.org/show/2878/
I did split my program in way too many modules, so if I should upload all scripts somewhere, just tell me!
Best wishes, Patrick
Am 15.11.2012 21:43, schrieb John ZuHone:
Hi Patrick,
How old is your file, by chance? FLASH3? FLASH2?
Also, if you don't mind, can you paste the script somewhere too?
Thanks,
John Z
On Nov 15, 2012, at 3:41 PM, Patrick Rieser wrote:
Hey all!
When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true.
Best wishes, Patrick
I add the field with:
def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}")
and then call:
rad_profile.add_fields("magB")
Traceback:
File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
John and Patrick, If you put a breakpoint at line 267 of fields.py in the FLASH frontend, can you see what's going wrong there? There should be an explicit catch for cgs units and I don't see why it's breaking for you. Nathan On Nov 15, 2012, at 6:56 PM, John ZuHone <jzuhone@gmail.com> wrote:
Hi Patrick,
Ok, so I have reproduced your error. It's very odd, actually. I also do FLASH simulations with magnetic fields, and I have the unitsystem = 'none' in my sims, and yt seems to catch this ok, but it's not catching the 'cgs' in yours. I can verify that I get this error whether I try to use your field or the magnetic fields already defined in yt.
I'll see what I can do about this.
Best,
John Z
On Nov 15, 2012, at 3:52 PM, Patrick Rieser wrote:
Hi John,
The files are from FLASH3. Ah and by the way, I am using the yt-dev version, updated some hours ago.
Here is the main script: http://paste.yt-project.org/show/2877/
and here the field is defined: http://paste.yt-project.org/show/2878/
I did split my program in way too many modules, so if I should upload all scripts somewhere, just tell me!
Best wishes, Patrick
Am 15.11.2012 21:43, schrieb John ZuHone:
Hi Patrick,
How old is your file, by chance? FLASH3? FLASH2?
Also, if you don't mind, can you paste the script somewhere too?
Thanks,
John Z
On Nov 15, 2012, at 3:41 PM, Patrick Rieser wrote:
Hey all!
When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true.
Best wishes, Patrick
I add the field with:
def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}")
and then call:
rad_profile.add_fields("magB")
Traceback:
File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Patrick and Nathan, I found it. The 'if' construct that Nathan is referring to was split up into different 'if's, and the only one that had the 'else' was the one for the 'none' case, which is why it wasn't breaking for me. I changed these to 'elif's. If you update the development version to the latest revision it should be fixed. Best, John On Nov 15, 2012, at 10:09 PM, Nathan Goldbaum wrote:
John and Patrick,
If you put a breakpoint at line 267 of fields.py in the FLASH frontend, can you see what's going wrong there? There should be an explicit catch for cgs units and I don't see why it's breaking for you.
Nathan
On Nov 15, 2012, at 6:56 PM, John ZuHone <jzuhone@gmail.com> wrote:
Hi Patrick,
Ok, so I have reproduced your error. It's very odd, actually. I also do FLASH simulations with magnetic fields, and I have the unitsystem = 'none' in my sims, and yt seems to catch this ok, but it's not catching the 'cgs' in yours. I can verify that I get this error whether I try to use your field or the magnetic fields already defined in yt.
I'll see what I can do about this.
Best,
John Z
On Nov 15, 2012, at 3:52 PM, Patrick Rieser wrote:
Hi John,
The files are from FLASH3. Ah and by the way, I am using the yt-dev version, updated some hours ago.
Here is the main script: http://paste.yt-project.org/show/2877/
and here the field is defined: http://paste.yt-project.org/show/2878/
I did split my program in way too many modules, so if I should upload all scripts somewhere, just tell me!
Best wishes, Patrick
Am 15.11.2012 21:43, schrieb John ZuHone:
Hi Patrick,
How old is your file, by chance? FLASH3? FLASH2?
Also, if you don't mind, can you paste the script somewhere too?
Thanks,
John Z
On Nov 15, 2012, at 3:41 PM, Patrick Rieser wrote:
Hey all!
When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true.
Best wishes, Patrick
I add the field with:
def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}")
and then call:
rad_profile.add_fields("magB")
Traceback:
File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Oops, sorry for the mistake! Thanks for catching it John. Nathan On Nov 15, 2012, at 7:16 PM, John ZuHone <jzuhone@gmail.com> wrote:
Patrick and Nathan,
I found it. The 'if' construct that Nathan is referring to was split up into different 'if's, and the only one that had the 'else' was the one for the 'none' case, which is why it wasn't breaking for me. I changed these to 'elif's.
If you update the development version to the latest revision it should be fixed.
Best,
John
On Nov 15, 2012, at 10:09 PM, Nathan Goldbaum wrote:
John and Patrick,
If you put a breakpoint at line 267 of fields.py in the FLASH frontend, can you see what's going wrong there? There should be an explicit catch for cgs units and I don't see why it's breaking for you.
Nathan
On Nov 15, 2012, at 6:56 PM, John ZuHone <jzuhone@gmail.com> wrote:
Hi Patrick,
Ok, so I have reproduced your error. It's very odd, actually. I also do FLASH simulations with magnetic fields, and I have the unitsystem = 'none' in my sims, and yt seems to catch this ok, but it's not catching the 'cgs' in yours. I can verify that I get this error whether I try to use your field or the magnetic fields already defined in yt.
I'll see what I can do about this.
Best,
John Z
On Nov 15, 2012, at 3:52 PM, Patrick Rieser wrote:
Hi John,
The files are from FLASH3. Ah and by the way, I am using the yt-dev version, updated some hours ago.
Here is the main script: http://paste.yt-project.org/show/2877/
and here the field is defined: http://paste.yt-project.org/show/2878/
I did split my program in way too many modules, so if I should upload all scripts somewhere, just tell me!
Best wishes, Patrick
Am 15.11.2012 21:43, schrieb John ZuHone:
Hi Patrick,
How old is your file, by chance? FLASH3? FLASH2?
Also, if you don't mind, can you paste the script somewhere too?
Thanks,
John Z
On Nov 15, 2012, at 3:41 PM, Patrick Rieser wrote:
Hey all!
When I try to call my own added field magB I get a runtime error saying that the cgs unitsystem is unrecognized. When I insert pf["unitsystem"] == "cgs" directly before the call where the traceback starts, it returns true.
Best wishes, Patrick
I add the field with:
def _magB(field,data): return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) add_field("magB", function=_magB, units=r"\rm{Gauss}")
and then call:
rad_profile.add_fields("magB")
Traceback:
File "rad_profile_magb.py", line 100, in <module> rad_profile.add_fields("magB") File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 170, in add_fields self._unlazy_add_fields(fields, weight, accumulation) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 151, in _unlazy_add_fields accumulation, self._args, check_cut = False) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 55, in save_state tr = func(*args, **kwargs) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 264, in _bin_field source_data = self._get_field(source, field, check_cut) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/profiles.py", line 204, in _get_field data.append(source[field][pointI].ravel().astype('float64')) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Auswertung/yt-bit/Plot/common_fields.py", line 7, in _magB return (na.sqrt(data["Bx"]**2 + data["By"]**2 + data["Bz"]**2)) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 328, in __getitem__ self.get_data(key) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2593, in get_data if self._generate_field(field): File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 358, in _generate_field self[field] = self.pf.field_info[field](self) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 378, in __call__ dd = self._function(self, data) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 357, in _Bx factor = GetMagRescalingFactor(data.pf) File "/home/patrick/Programme/yt-x86_64/src/yt-hg/yt/frontends/flash/fields.py", line 353, in GetMagRescalingFactor "value %s is unrecognized" % pf['unitsystem']) RuntimeError: Runtime parameter unitsystem withvalue cgs is unrecognized
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (3)
-
John ZuHone -
Nathan Goldbaum -
Patrick Rieser