Zonal regridding seg faults unless regridTool='regrid2'

[katemarvel]

I have zonally averaged data I need to re-grid to a common grid g.

print data.getAxisList()

[ id: lat
Designated a latitude axis.
units: degrees_north
Length: 145
First: -90.0
Last: 90.0
Other axis attributes:
long_name: latitude
standard_name: latitude
axis: Y
Python id: 0x35911d0
,
id: axis_2
Designated a longitude axis.
units: degrees_east
Length: 1
First: 180.0
Last: 180.0
Other axis attributes:
axis: X
modulo: 360.0
topology: circular
Python id: 0x3587e90
]
The following code leads to a seg fault:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl)

But this code works fine:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='regrid2')

Finally, trying

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='esmf')

gives the following error:
NameError: ESMC_FieldRegridStore() failed with rc = 6. Please check the log files (named "PET*.ESMF_LogFile").

[doutriaux1]

Kate,

I don't think ESMF can handle zonal grids. So you're right falling back to regrid2 is right.

Shouldn't seg fault though, it should automatically fallback on regrid2, could you please be so kind to submit a bug on github? Thx!

C.

From: katemarvel <notifications@github.zerozr99.workers.devmailto:notifications@github.com>
Reply-To: UV-CDAT/uvcdat <reply@reply.github.commailto:reply@reply.github.com>
Date: Friday, February 8, 2013 9:14 AM
To: UV-CDAT/uvcdat <uvcdat@noreply.github.commailto:uvcdat@noreply.github.com>
Subject: [uvcdat] Zonal regridding seg faults unless regridTool='regrid2' (#118)

I have zonally averaged data I need to re-grid to a common grid g.

print data.getAxisList()

[ id: lat
Designated a latitude axis.
units: degrees_north
Length: 145
First: -90.0
Last: 90.0
Other axis attributes:
long_name: latitude
standard_name: latitude
axis: Y
Python id: 0x35911d0
,
id: axis_2
Designated a longitude axis.
units: degrees_east
Length: 1
First: 180.0
Last: 180.0
Other axis attributes:
axis: X
modulo: 360.0
topology: circular
Python id: 0x3587e90
]
The following code leads to a seg fault:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl)

But this code works fine:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='regrid2')

Finally, trying

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='esmf')

gives the following error:
NameError: ESMC_FieldRegridStore() failed with rc = 6. Please check the log files (named "PET*.ESMF_LogFile").

—
Reply to this email directly or view it on GitHubhttps://github.com//issues/118.

[katemarvel]

I thought that was submitting a bug on github? I did what Paul told me to do.

Is there any documentation at all for the new regridder? I have no idea what ESGF is.

Kate Marvel
Program for Climate Model Diagnosis
and Intercomparison
Lawrence Livermore National Laboratory
P.O. Box 808, L-103
Livermore, CA 94551 USA
Phone: (925) 422-2697

From: Charles <notifications@github.zerozr99.workers.devmailto:notifications@github.com>
Reply-To: UV-CDAT/uvcdat <reply@reply.github.commailto:reply@reply.github.com>
Date: Friday, February 8, 2013 9:17 AM
To: UV-CDAT/uvcdat <uvcdat@noreply.github.commailto:uvcdat@noreply.github.com>
Cc: "Marvel, Katherine D." <marvel1@llnl.govmailto:marvel1@llnl.gov>
Subject: Re: [uvcdat] Zonal regridding seg faults unless regridTool='regrid2' (#118)

Kate,

I don't think ESMF can handle zonal grids. So you're right falling back to regrid2 is right.

Shouldn't seg fault though, it should automatically fallback on regrid2, could you please be so kind to submit a bug on github? Thx!

C.

From: katemarvel <notifications@github.zerozr99.workers.devmailto:notifications@github.zerozr99.workers.devmailto:notifications@github.com>
Reply-To: UV-CDAT/uvcdat <reply@reply.github.commailto:reply@reply.github.commailto:reply@reply.github.com>
Date: Friday, February 8, 2013 9:14 AM
To: UV-CDAT/uvcdat <uvcdat@noreply.github.commailto:uvcdat@noreply.github.commailto:uvcdat@noreply.github.com>
Subject: [uvcdat] Zonal regridding seg faults unless regridTool='regrid2' (#118)

I have zonally averaged data I need to re-grid to a common grid g.

print data.getAxisList()

[ id: lat
Designated a latitude axis.
units: degrees_north
Length: 145
First: -90.0
Last: 90.0
Other axis attributes:
long_name: latitude
standard_name: latitude
axis: Y
Python id: 0x35911d0
,
id: axis_2
Designated a longitude axis.
units: degrees_east
Length: 1
First: 180.0
Last: 180.0
Other axis attributes:
axis: X
modulo: 360.0
topology: circular
Python id: 0x3587e90
]
The following code leads to a seg fault:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl)

But this code works fine:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='regrid2')

Finally, trying

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='esmf')

gives the following error:
NameError: ESMC_FieldRegridStore() failed with rc = 6. Please check the log files (named "PET*.ESMF_LogFile").

—
Reply to this email directly or view it on GitHubhttps://github.com//issues/118.

—
Reply to this email directly or view it on GitHubhttps://github.com//issues/118#issuecomment-13300747.

[doutriaux1]

Oh yes you're right…doh…

From: katemarvel <notifications@github.zerozr99.workers.devmailto:notifications@github.com>
Reply-To: UV-CDAT/uvcdat <reply@reply.github.commailto:reply@reply.github.com>
Date: Friday, February 8, 2013 9:19 AM
To: UV-CDAT/uvcdat <uvcdat@noreply.github.commailto:uvcdat@noreply.github.com>
Cc: Charles Doutriaux <doutriaux1@llnl.govmailto:doutriaux1@llnl.gov>
Subject: Re: [uvcdat] Zonal regridding seg faults unless regridTool='regrid2' (#118)

I thought that was submitting a bug on github? I did what Paul told me to do.

Is there any documentation at all for the new regridder? I have no idea what ESGF is.

Kate Marvel
Program for Climate Model Diagnosis
and Intercomparison
Lawrence Livermore National Laboratory
P.O. Box 808, L-103
Livermore, CA 94551 USA
Phone: (925) 422-2697

From: Charles <notifications@github.zerozr99.workers.devmailto:notifications@github.zerozr99.workers.devmailto:notifications@github.com>
Reply-To: UV-CDAT/uvcdat <reply@reply.github.commailto:reply@reply.github.commailto:reply@reply.github.com>
Date: Friday, February 8, 2013 9:17 AM
To: UV-CDAT/uvcdat <uvcdat@noreply.github.commailto:uvcdat@noreply.github.commailto:uvcdat@noreply.github.com>
Cc: "Marvel, Katherine D." <marvel1@llnl.govmailto:marvel1@llnl.govmailto:marvel1@llnl.gov>
Subject: Re: [uvcdat] Zonal regridding seg faults unless regridTool='regrid2' (#118)

Kate,

I don't think ESMF can handle zonal grids. So you're right falling back to regrid2 is right.

Shouldn't seg fault though, it should automatically fallback on regrid2, could you please be so kind to submit a bug on github? Thx!

C.

From: katemarvel <notifications@github.zerozr99.workers.devmailto:notifications@github.zerozr99.workers.devmailto:notifications@github.commailto:notifications@github.com>
Reply-To: UV-CDAT/uvcdat <reply@reply.github.commailto:reply@reply.github.commailto:reply@reply.github.commailto:reply@reply.github.com>
Date: Friday, February 8, 2013 9:14 AM
To: UV-CDAT/uvcdat <uvcdat@noreply.github.commailto:uvcdat@noreply.github.commailto:uvcdat@noreply.github.commailto:uvcdat@noreply.github.com>
Subject: [uvcdat] Zonal regridding seg faults unless regridTool='regrid2' (#118)

I have zonally averaged data I need to re-grid to a common grid g.

print data.getAxisList()

[ id: lat
Designated a latitude axis.
units: degrees_north
Length: 145
First: -90.0
Last: 90.0
Other axis attributes:
long_name: latitude
standard_name: latitude
axis: Y
Python id: 0x35911d0
,
id: axis_2
Designated a longitude axis.
units: degrees_east
Length: 1
First: 180.0
Last: 180.0
Other axis attributes:
axis: X
modulo: 360.0
topology: circular
Python id: 0x3587e90
]
The following code leads to a seg fault:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl)

But this code works fine:

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='regrid2')

Finally, trying

gl = cdms.createZonalGrid(g)
regridded = data.regrid(gl,regridTool='esmf')

gives the following error:
NameError: ESMC_FieldRegridStore() failed with rc = 6. Please check the log files (named "PET*.ESMF_LogFile").

—
Reply to this email directly or view it on GitHubhttps://github.com//issues/118.

—
Reply to this email directly or view it on GitHubhttps://github.com//issues/118#issuecomment-13300747.

—
Reply to this email directly or view it on GitHubhttps://github.com//issues/118#issuecomment-13300867..

[doutriaux1]

@pletzer do you know what's going on here? moved for 2.2 release

[doutriaux1]

#1060 fixes this