Re: rect to curv

Dear Mahmoud,

In the future, please do not post your questions multiple times. You posted your three times question over the weekend, and we do not always check our email then.

It sounds like you downloaded a version of NCL that doesn't have "ESMF_RegridWeightGen" bundled with it. As I mentioned below, you need to be running NCL V6.1.0 or
later to use ESMF regridding. All of the V6.1.2 NCL binaries have ESMF_RegridWeightGen bundled with them, with the exception of a CentOS binary that only had gcc V4.1.2 (ESMF needs a later version of gcc to build).

--Mary

On Feb 10, 2013, at 11:14 AM, mahmoud elkarim wrote:

> dear mary and ncl users
>
> when i was trying to use ESMF_regrid_1.ncl to convert from rectilinear to curvilinear http://www.ncl.ucar.edu/Applications/ESMF.shtml i found this error
>
> sh: ESMF_RegridWeightGen: command not found
> fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 2291 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 3166 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 57 in file ESMF_regrid_1.ncl
>
> could you help me
> thank you for your attention
>
> mahmoud gaber
> Nile forecast centre (NFC)
> minister of water resources and irrigations (MWRI)
> cornich el-nil emababa giza- egypt
>
>
>
>
>
> On Thu, Feb 7, 2013 at 6:44 PM, Mary Haley <haley@ucar.edu> wrote:
> Dear Mahmoud,
>
> Please, in the future, include ncl-talk on the response, and not just me. It's important that people see the follow-up questions and answers.
>
> If you want to do regridding with the ESMF routines, you will have to download V6.1.0 or V6.1.2, which was just released today.
>
> If you don't want to upgrade, then you can try the "rgrid2rcm" function:
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/rgrid2rcm.shtml
>
>
> --Mary
>
> On Feb 6, 2013, at 3:21 PM, mahmoud elkarim wrote:
>
>> dear mary
>> thank you so much for helping me but i'm using ncl version 6.00 and i found problem in
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" for that i'm trying to install the latest version that you told about and
>>
>> i wish to install it successfully
>>
>>
>> regards
>> mahmoud
>> On Sat, Feb 2, 2013 at 12:25 AM, Mary Haley <haley@ucar.edu> wrote:
>> Dear Mahmoud Gaber,
>>
>> You are trying to subtract two 2-dimensional arrays which are different sizes. prec1 is 127 x 129, and prec is 111 x 130. NCL cannot subtract two multi-dimensional arrays that are different sizes.
>>
>> I see that you were trying to regrid, so I'll try to help you with that.
>>
>> Your GPCP grid is rectilinear (because it has 1-dimensional coordinate arrays), and your REGCM grid is curvilinear (because it has 2-dimensional coordinate variables on the file).
>>
>> So, using code from the "ESMF_rect_to_curv.ncl" template file in http://www.ncl.ucar.edu/Applications/Templates, you can replace this code:
>>
>> Opt = True
>> Opt@WgtFileName = "30day_SRFmeanDJF.nc" ; default is "weights_file.nc"
>> ;---Destination file options
>> Opt@DstGridType = "0.5x0.5" ; 5x5 degree grid
>> Opt@DstLLCorner = (/ -60.d, 0.d/)
>> Opt@DstURCorner = (/ 60.d, 355.d/)
>>
>> Opt@SrcForceOverwrite = True
>> Opt@DstForceOverwrite = True
>> Opt@WgtForceOverwrite = True
>> ;;Opt@PrintTimings = True
>> ;;Opt@Debug = True
>> ;
>>
>> ; temp_regrid = ESMF_regrid(p,Opt)
>>
>> with this code (UNTESTED):
>>
>> ;---Set up regridding options
>> Opt = True
>>
>> ;---"bilinear" is the default. "patch" and "conserve" are other options.
>> Opt@InterpMethod = "bilinear"
>>
>> Opt@WgtFileName = "gpcp_to_regcm_wgts.nc"
>>
>> Opt@SrcRegional = True
>> Opt@DstRegional = True ;;--Change (maybe)
>>
>> Opt@ForceOverwrite = True
>> Opt@PrintTimings = True
>> Opt@Debug = True
>>
>> prec_regrid = ESMF_regrid(prec,Opt) ; Do the regridding
>>
>> printVarSummary(prec_regrid)
>>
>> Note that you may also want to set these options, which can make the regridding go faster:
>>
>> Opt@DstLLCorner = (/ -60.d, 0.d/) ; lower left corner, (/lat,lon/)
>> Opt@DstURCorner = (/ 60.d, 355.d/) ; upper right corner (/lat,/lon/)
>>
>> Use values that are close to the limits of your REGCM grid.
>>
>> prec_regrid should now be the same size as "prec1".
>>
>> --Mary
>>
>>
>> On Feb 1, 2013, at 2:39 PM, mahmoud elkarim wrote:
>>
>>> dear all
>>> find attached ncl scripits
>>> I'm beginner user in ncl scripts
>>> when i was trying to find precipitation bias which is the difference between regcm output and GPCP observed data i found this error
>>> (fatal:Minus: Dimension size, for dimension number 0, of operands does not match, can't continue)
>>> and i don't know how i can regrid gpcp to regcm
>>>
>>> rgcm file dump is
>>> x = 129 ;
>>> y = 127 ;
>>> m10 = 1 ;
>>> m2 = 1 ;
>>> lev = 2 ;
>>> time = UNLIMITED ; // (1 currently)
>>> nb2 = 2 ;
>>> variables:
>>> float xlon(y, x) ;
>>> xlon:standard_name = "longitude" ;
>>> xlon:long_name = "Longitude at cross points" ;
>>> xlon:units = "degrees_east" ;
>>> xlon:_CoordinateAxisType = "Lon" ;
>>> float xlat(y, x) ;
>>> xlat:standard_name = "latitude" ; smw:cell_methods = "time: point" ;
>>> float tpr (time, y, x) ;
>>> tpr:standard_name = "precipitation_flux" ;
>>> tpr:long_name = "Total precipitation" ;
>>> tpr:units = "kg m-2 day-1" ;
>>> tpr:coordinates = "xlon xlat" ;
>>> tpr:cell_methods = "time: mean (interval: 3 hours)" ;
>>>
>>> xlat:long_name = "Latitude at cross points" ;
>>> xlat:units = "degrees_north" ;
>>> xlat:_CoordinateAxisType = "Lat" ;netcdf gpccDJFmeanbox {
>>>
>>>
>>> and gpgp dump is
>>> dimensions:
>>> lon = 130 ;
>>> lat = 111 ;
>>> time = UNLIMITED ; // (1 currently)
>>> variables:
>>> double lon(lon) ;
>>> lon:standard_name = "longitude" ;
>>> lon:long_name = "longitude" ;
>>> lon:units = "degrees_east" ;
>>> lon:axis = "X" ;
>>> double lat(lat) ;
>>> lat:standard_name = "latitude" ;
>>> lat:long_name = "latitude" ;
>>> lat:units = "degrees_north" ;
>>> lat:axis = "Y" ;
>>> double time(time) ;
>>> time:standard_name = "time" ;
>>> time:units = "days since 1901-01-01 00:00:00" ;
>>> time:calendar = "proleptic_gregorian" ;
>>> float p(time, lat, lon) ;
>>> p:long_name = "full data precipitation version6" ;
>>> p:units = "mm/month" ;
>>> p:code = 16 ;
>>> p:_FillValue = -99999.99f ;
>>>
>>>
>>>
>>> i ask if any one could help me
>>> thanks for attention
>>>
>>>
>>>
>>>
>>> > Mahmoud Gaber
>>> > Nile Forecasting center (NFC)
>>> > minister of water resources and irrigation (MWRI)
>>> <gpcp.ncl>_______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
>
>
>
> <ESMF_regrid_1.ncl>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Feb 11 09:23:07 2013