area_hi2lores, missing output

From: Erik Noble <enoble_at_nyahnyahspammersnyahnyah>
Date: Fri, 12 Jun 2009 23:58:16 -0400

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

Picture_25.png
Received on Fri Jun 12 2009 - 21:58:16 MDT

This archive was generated by hypermail 2.2.0 : Mon Jun 15 2009 - 08:53:05 MDT