[ncl-talk] Can ESMF_regrid be used for output Gaussian grid?

Rashed Mahmood rashidcomsis at gmail.com
Tue Feb 28 18:26:29 MST 2017


Hi everyone,

My data is 1*1 degree resolution which is Not global (latitute span:
45N-89N) and I want to convert it into a Gaussian grid as below. I thought
ESMF_regrid would be easy But I am not sure what to specify for attribute
"Opt at DstGridType       = "??????" for a Gaussian destination grid?

So I tried using area_conserve_remap as in following simple code:

begin
   ifil   = addfile("test_file.nc","r")
   cod    = ifil->cod
   printVarSummary(cod)

   NLATT  = 64
   MLONT  =  2*NLATT
   LATT   = latGau(NLATT, "lat", "latitude" , "degrees_north")
   LONT   =  fspan(-180,180,128)
   TLAT   = LATT({45:90})

   opt       = True
   opt at NLATi = 178
   opt at NALTo = NLATT
   gcod      = area_conserve_remap_Wrap (cod&lon, cod&lat, cod ,LONT, TLAT,
opt)

end


However, I am getting error as below, which I understand is due not global
lats but not sure how to fix it? The sample data file attached.

Copyright (C) 1995-2015 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.3.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.

Variable: cod
Type: float
Total Size: 64980 bytes
            16245 values
Number of Dimensions: 2
Dimensions and sizes:   [lat | 45] x [lon | 361]
Coordinates:
            lat: [  45..  89]
            lon: [-180.. 180]
Number Of Attributes: 3
  _FillValue :  -999
  remap :       remapped via ESMF_regrid_with_weights: Bilinear
  missing_value :       -999
 Error in "cremapbin":
 Could not map global lat array into grid array
Abort (core dumped)

any thoughts?



Thank you,
Rashed
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170228/db947d39/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_file.nc
Type: application/x-netcdf
Size: 68672 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170228/db947d39/attachment.nc 


More information about the ncl-talk mailing list