[ncl-talk] Why the source grid be recognized as curvilinear grid?
林祥
xianglin72 at icloud.com
Mon Dec 19 01:54:23 MST 2016
i’ve resolved the problem. Thanks!
> 在 2016年12月19日，01:14，林祥 <xianglin72 at icloud.com> 写道：
>
> Hi, all
>
> I’m regridding multiple CMIP5 model output with ESMF_regrid function, among all the data file, lat/lon coordinate array can be 1-dimension or 2 dimension.
>
> Now I met a problem , Please see the below message output from the scripts:
>
>
> (0) 4. /data1/linx/CMIP5_SIC/historical/BNU-ESM/sic_OImon_BNU-ESM_historical_r1i1p1_185001-200512.nc ; ( Loop 4 )
>
>
> Variable: Opt
> Type: logical
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> Number Of Attributes: 9
> Debug : True
> PrintTimings : True
> DstGridLon : <ARRAY of 360 elements>
> DstGridLat : <ARRAY of 180 elements>
> DstGridMask : <ARRAY of 64800 elements>
> SrcGridMask : <ARRAY of 72000 elements>
> ForceOverwrite : True
> DstFileName : destination_grid_file.nc
> SrcFileName : source_grid_file.nc
> (0) True
> (0) Generating interpolation weights from CMIP5 to
> (0) World 1 degree grid using the bilinear method.
> (0) get_src_grid_info: source lat dims = (200,360)
> (0) get_src_grid_info: source lon dims = (200,360)
> (0) get_src_grid_info: source grid type is 'curvilinear'
> (0) curvilinear_to_SCRIP: calculating grid corners...
> (0) curvilinear_to_SCRIP: no lat values are at the poles, so
> (0) calculating grid corners using
> (0) calc_SCRIP_corners_noboundaries...
> (0) calc_SCRIP_corners_noboundaries
> (0) min/max original lat: -81.5003/89.4879
> (0) min/max original lon: 0.0457764/359.949
> (0) calc_SCRIP_corners_noboundaries
> (0) min/max Extlat2d: -82.5003/90.4798
> (0) min/max Extlon2d: -358.806/360.995
> (0) calc_SCRIP_corners_noboundaries
> (0) min/max ExtGridCenter_lat: -82.0003/89.9838
> (0) min/max ExtGridCenter_lon: -179.22/360.397
> (0) =====> CPU Elapsed Time: curvilinear_to_SCRIP: 0.164974 seconds <=====
> (0) get_dst_grid_info: destination lat dims = (180)
> (0) get_dst_grid_info: destination lon dims = (360)
> (0) curvilinear_to_SCRIP: calculating grid corners...
> (0) curvilinear_to_SCRIP: no lat values are at the poles, so
> (0) calculating grid corners using
> (0) calc_SCRIP_corners_noboundaries...
> (0) calc_SCRIP_corners_noboundaries
> (0) min/max original lat: -89.5/89.5
> (0) min/max original lon: -179.5/179.5
> (0) calc_SCRIP_corners_noboundaries
> (0) min/max Extlat2d: -90.5/90.5
> (0) min/max Extlon2d: -180.5/180.5
> (0) calc_SCRIP_corners_noboundaries
> (0) min/max ExtGridCenter_lat: -90/90
> (0) min/max ExtGridCenter_lon: -180/180
> (0) =====> CPU Elapsed Time: rectilinear_to_SCRIP: 0.155979 seconds <=====
> (0) ESMF_regrid_gen_weights: number of processors used: 1
> (0) --------------------------------------------------
> (0) ESMF_regrid_gen_weights: the following command is about to be executed on the system:
> (0) 'ESMF_RegridWeightGen --source source_grid_file.nc --destination destination_grid_file.nc --weight ESMF_regrid_wgt_bilinear.nc --method bilinear -i'
> (0) --------------------------------------------------
> (0) ESMF_regrid_gen_weights: output from 'ESMF_RegridWeightGen':
> (0) Starting weight generation with these inputs:
> (1) Source File: source_grid_file.nc
> (2) Destination File: destination_grid_file.nc
> (3) Weight File: ESMF_regrid_wgt_bilinear.nc
> (4) Source File is in SCRIP format
> (5) Source Grid is a global grid
> (6) Source Grid is a logically rectangular grid
> (7) Destination File is in SCRIP format
> (8) Destination Grid is a global grid
> (9) Destination Grid is a logically rectangular grid
> (10) Regrid Method: bilinear
> (11) Pole option: ALL
> (12) Ignore unmapped destination points
> (13)
> (14) Completed weight generation successfully.
> (15)
> (0) --------------------------------------------------
> (0) ESMF_regrid_gen_weights: 'ESMF_RegridWeightGen' was successful.
> (0) =====> CPU Elapsed Time: ESMF_regrid_gen_weights: 0.0169983 seconds <=====
> (0) ESMF_regrid_with_weights: regridding using interpolation weights ...
> (0) ESMF_regrid_with_weights: warning: destination grid is not
> (0) completely covered by the source grid. This is not an error.
> (0) It just means your destination grid covers a larger area
> (0) than your source grid.
> (0) ESMF_regrid_with_weights: Source Grid:
> (0) rank: 2
> (0) dimensions: 200 360
> (0) original source rank: 3
> (0) latitude min/max: -81.50032043457031/89.48786163330078
> (0) longitude min/max:0.0457763671875/359.9493103027344
> (0) ESMF_regrid_with_weights: Destination Grid:
> (0) dimensions: 180 360
> (0) latitude min/max: -89.5/89.5
> (0) longitude min/max:-179.5/179.5
> (0) ESMF_regrid_with_weights: retrieving interpolation weights ...
> (0) ESMF_regrid_with_weights: calling sparse_matrix_mult to apply weights...
> (0) ESMF_regrid_with_weights: dstData
> (0) Dimensions: 1872 180 360
> (0) minSrcData: 0
> (0) maxSrcData: 99.9994
> (0) minDstData: -2.50832e-12
> (0) maxDstData: 99.9981
> (0) =====> CPU Elapsed Time: ESMF_regrid_with_weights: 122.265 seconds <=====
>
> Variable: thetao_regrid
> Type: float
> Total Size: 485222400 bytes
> 121305600 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1872] x [lat | 180] x [lon | 360]
> Coordinates:
> time: [15.5..56924.5]
> lat: [89.5..-89.5]
> lon: [-179.5..179.5]
> Number Of Attributes: 11
> remap : remapped via ESMF_regrid_with_weights: Bilinear remapping
> standard_name : sea_ice_area_fraction
> long_name : Sea Ice Area Fraction
> comment : fraction of grid cell covered by sea ice.
> units : %
> cell_methods : time: mean
> cell_measures : area: areacello
> history : 2012-04-01T10:38:11Z altered by CMOR: replaced missing value flag (1e+30) with standard missing value (1e+20).
> missing_value : 1e+20
> _FillValue : 1e+20
> associated_files : baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation <http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation> gridspecFile: gridspec_seaIce_fx_BNU-ESM_historical_r0i0p0.nc areacello: areacello_fx_BNU-ESM_historical_r0i0p0.nc
> warning:lbLabelBrOn is not a resource in the given object
> warning:lbLabelBrOn is not a resource in the given object
>
> ; it’s ok for the curvilinear grid to rectilinear grid for Loop 4
>
>
> (0) 5. /data1/linx/CMIP5_SIC/historical/CanCM4/sic_OImon_CanCM4_historical_r1i1p1_196101-200512.nc
>
>
> Variable: Opt
> Type: logical
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> Number Of Attributes: 9
> Debug : True
> PrintTimings : True
> DstGridLon : <ARRAY of 360 elements>
> DstGridLat : <ARRAY of 180 elements>
> DstGridMask : <ARRAY of 64800 elements>
> SrcGridMask : <ARRAY of 8192 elements>
> ForceOverwrite : True
> DstFileName : destination_grid_file.nc
> SrcFileName : source_grid_file.nc
> (0) True
> (0) Generating interpolation weights from CMIP5 to
> (0) World 1 degree grid using the bilinear method. ; ************************************************************
> (0) get_src_grid_info: source lat dims = (64) ; the lat coordinate of source data is 1-dimension
> (0) get_src_grid_info: source lon dims = (128)
> (0) get_src_grid_info: source grid type is ‘curvilinear’ ; but recognized as “curvilinear” and report error as below
> ; *************************************************
>
> fatal:Number of dimensions in parameter (1) of (curvilinear_to_SCRIP) is (1), (2) dimensions were expected
> fatal:["Execute.c":8567]:Execute: Error occurred at or near line 3320 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> fatal:["Execute.c":8567]:Execute: Error occurred at or near line 3427 in file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl
>
> fatal:["Execute.c":8567]:Execute: Error occurred at or near line 147 in file esmf_regrid_6_sic.ncl
>
>
>
> Thank!
>
>
> Lin
>
> <esmf_regrid_6_sic.ncl>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161219/db26e921/attachment.html
More information about the ncl-talk
mailing list