[ncl-talk] Why the source grid be recognized as curvilinear grid?
林祥
xianglin72 at icloud.com
Sun Dec 18 10:14:18 MST 2016
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 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161219/84f941b5/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: esmf_regrid_6_sic.ncl
Type: application/octet-stream
Size: 7735 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161219/84f941b5/attachment.obj
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161219/84f941b5/attachment-0001.html
More information about the ncl-talk
mailing list