[ncl-talk] Sub: linin2

[dale zuri]

Hi Rick,

begin
;;;;;;;;;;;;;;;;;ERA;;;;;;;;;;;;;;;;;;;;;;;;;;;
;  iz = addfile("psa_mask_cfsv2.nc","r")
  iz = addfile("tmp2m.2012050218.time.nc","r")
  print(iz)
  printVarSummary(iz)

P1=iz->TMP_P0_L103_GGA0(:,15:40,180:255)
  printVarSummary(P1)
  printMinMax(P1,0)
  ds = addfile("nlwrs.sfc.gauss.1948.nc","r")
  pp=ds->nlwrs(:,:,:)
  p = pp(:,{76:51},{168:239})

  printVarSummary(p)
;****************************************************
 psa_gauss = linint2_Wrap(P1&lon_0,P1&lat_0,P1, False, p&lon,p&lat, 0)

  printVarSummary(psa_gauss)
  printMinMax(psa_gauss,0)
return



I changed the line. Then also the same issue of getting a missing number as
output.
P1=iz->TMP_P0_L103_GGA0(:,15:40,180:255)
 p = pp(:,{76:51},{168:239})

Variable: P1
Type: float
Total Size: 9555936 bytes
            2388984 values
Number of Dimensions: 3
Dimensions and sizes:    [forecast_time0 | 1209] x [lat_0 | 26] x [lon_0 |
76]
Coordinates:
            forecast_time0: [6..7254]
            lat_0: [75.11737..51.49582]
            lon_0: [168.7498..239.0622]

Variable: p
Type: float
Total Size: 723216 bytes
            180804 values
Number of Dimensions: 3
Dimensions and sizes:    [time | 366] x [lat | 13] x [lon | 38]
Coordinates:
            time: [1297320..1306080]
            lat: [75.2351..52.3799]
            lon: [168.75..238.125]


On Mon, Nov 13, 2017 at 8:31 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:

> Hi Dale,
>
> It would really help if you could show the results of those various
> printVarSummary()'s And include one for pp/p -- are you certain your
> coordinate subscripting is correct in this line:
>
>  p = pp(:,{51:75},{168:239})
>
> Rick
>
>
> On Mon, Nov 13, 2017 at 9:33 PM, dale zuri <dalezuri at gmail.com> wrote:
>
>> Hi all,
>> Please help me understand this issue.
>> I didn't get any error. But the interpolated output are missing numbers.
>>  iz = addfile("tmp2m.2012050218.time.nc","r")
>>   print(iz)
>>   printVarSummary(iz)
>>
>> P1=iz->TMP_P0_L103_GGA0(:,75:51,168:239)
>>   printVarSummary(P1)
>>   printMinMax(P1,0)
>>   ds = addfile("nlwrs.sfc.gauss.1948.nc","r")
>>   pp=ds->nlwrs(:,::-1,:)
>>   p = pp(:,{51:75},{168:239})
>>
>>   printVarSummary(p)
>> ;****************************************************
>>  psa_gauss = linint2_Wrap(P1&lon_0,P1&lat_0,P1, False, p&lon,p&lat, 0)
>>   printVarSummary(psa_gauss)
>>   printMinMax(psa_gauss,0)
>>
>>
>> ;;;;
>> (0)    Temperature: min=1e+20   max=1e+20
>>
>> Thanks
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171113/e02e027a/attachment.html>