Re: area_hi2lores, missing output

I figured this out. It was a memory problem.

On Fri, Jun 12, 2009 at 11:58 PM, Erik Noble <enoble_at_giss.nasa.gov> wrote:

> Hi. Thank you for the help earlier with area_hi2lores and memory. I used
> area_hi2lores on a CMORPH precipitation file. After I use the function and
> then output my regridded data, the data does not graphically show up.
> Does anyone else use CMORPH? If so, I know that each record contains a 1440
> x 480 REAL*4 array of data which is oriented from 0.125E EASTward and from
> 59.875N SOUTHward, with a grid increment of 0.25 degrees of latitude and
> longitude. I don't think I am doing the lat and lon coordinates correctly.
> My command-line terminal history is below. Also a picture is attached when
> I try to look at the output.
>
> Thank you for your time.
> Sincerely,
> Erik
>
> noble:CMORPH enoble$ ncl
> Copyright (C) 1995-2009 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.1.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> ncl 0> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> ncl 1> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ncl 4> a = addfile("20060901_CMORPH_3hr-025deg.nc","r")
> ncl 5> CMORPH_rate = a->CMORPH_RATE
> ncl 6> printVarSummary(CMORPH_rate)
>
>
> Variable: CMORPH_rate
> Type: float
> Total Size: 22118400 bytes
> 5529600 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 8] x [lat | 480] x [lon | 1440]
> Coordinates:
> time: [2006090100..2006090121]
> lat: [59.875..-59.875]
> lon: [0.125..359.875]
> Number Of Attributes: 4
> lonFlip : longitude coordinate variable has been reordered via lonFlip
> long_name : CMORPH
> units : mm/hr
> _FillValue : -9999
> ncl 7> CMORPH_rate=CMORPH_rate(:,::-1,:)
> ncl 8> printVarSummary(CMORPH_rate)
>
>
> Variable: CMORPH_rate
> Type: float
> Total Size: 22118400 bytes
> 5529600 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 8] x [lat | 480] x [lon | 1440]
> Coordinates:
> time: [2006090100..2006090121]
> lat: [-59.875..59.875]
> lon: [0.125..359.875]
> Number Of Attributes: 4
> lonFlip : longitude coordinate variable has been reordered via lonFlip
> long_name : CMORPH
> units : mm/hr
> _FillValue : -9999
> ncl 9> ;********************
> ncl 10> ; area_hi2lores
> ncl 11> nlat = 240
> ncl 12> mlon = 480
> ncl 13> LAT = 59.875 - ispan(0,nlat-1,1)*0.50
> ncl 14> LON = 0.125 + ispan(0,mlon-1,1)*0.50
> ncl 15> lat = CMORPH_rate(0,:,0)
> ncl 16> lon = CMORPH_rate(0,0,:)
> ncl 17> po = area_hi2lores_Wrap (lon,lat, CMORPH_rate , True, 1, LON,
> LAT, False) ; "po" being the example from ncl documentation
>
> ncl 18> printVarSummary(po)
>
>
> Variable: po
> Type: float
> Total Size: 3686400 bytes
> 921600 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 8] x [LAT | 240] x [LON | 480]
> Coordinates:
> time: [2006090100..2006090121]
> LAT: [59.875..-59.625]
> LON: [0.125..239.625]
> Number Of Attributes: 4
> lonFlip : longitude coordinate variable has been reordered via lonFlip
> long_name : CMORPH
> units : mm/hr
> _FillValue : -9999
> ncl 19> ncout ="CMORPH_0901_3hr_.50deg.nc"
> ncl 20> system("/bin/rm -f "+ncout)
> ncl 21> b = addfile(ncout,"c")
> ncl 22> b->CMORPH_RATE = po
> ncl 23> exit
> noble:CMORPH enoble$ ncview CMORPH_0901_3hr_.50deg.nc
> Ncview 1.93f David W. Pierce 12 September 2008
> http://meteora.ucsd.edu:80/~pierce/ncview_home_page.html<http://meteora.ucsd.edu:80/%7Epierce/ncview_home_page.html>
> Copyright (C) 1993 through 2008, David W. Pierce
> Ncview comes with ABSOLUTELY NO WARRANTY; for details type `ncview -w'.
> This is free software licensed under the Gnu General Public License version
> 3; type `ncview -c' for redistribution details.
>
> Note: no Ncview app-defaults file found, using internal defaults
> calculating min and maxes for CMORPH_RATE...
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Jun 14 2009 - 11:10:50 MDT