Re: Regrid error : area_conserve_remap

Hi. Ok. I missed the part on the example page that says this function
assumes the input is global. It is not.
To make this work, the page says I need to declare the "dimension of the
*global* latitude set of which the grid is a subset."
My model data is on a 20km by 20km regional grid. Do I need to make this
declaration in "degrees" only.

opt = True
opt@NLATi = ???

On Sat, Sep 24, 2011 at 7:20 PM, Dennis Shea <shea@ucar.edu> wrote:

> area_conserve_remap requires that *both* the input and output grid be
> global. The input grid
>
> Variable: U
> [snip]
>
> lat: [-18.40244..35.33525]
> lon: [-35.74922..35.74922]
>
> is not global
>
>
> On 9/24/11 2:30 PM, Erik Noble wrote:
>
>> Hi. I am trying to remap a hi resolution file to a 2.5 degree global
>> grid, similar to NCEP Reanalysis II, so that I can compare data. I'm
>> following the examples using the area_conserve_remap function.
>> http://www.ncl.ucar.edu/**Document/Functions/**Contributed/area_conserve_
>> **remap_Wrap.shtml<http://www.ncl.ucar.edu/Document/Functions/Contributed/area_conserve_remap_Wrap.shtml>
>> I keep getting this error:
>>
>> Error in "cremapbin":
>> Could not map global lat array into grid array
>> Segmentation fault
>>
>> I am not sure where I am creating this error. Could I please have some
>> help with this?
>> Sincerely,
>> Erik
>>
>>
>> [eunoble@borg01t008] Experiment_01 $ ncl test_regrid.ncl
>> Copyright (C) 1995-2011 - All Rights Reserved
>> University Corporation for Atmospheric Research
>> NCAR Command Language Version 6.0.0
>> The use of this software is governed by a License Agreement.
>> See http://www.ncl.ucar.edu/ for more details.
>> (0) ****************************************************************
>> ***
>> (0) These are the WRF Model output files you are extracting from.
>> (0) number of FILES = 97
>> (0) ****************************************************************
>> ***
>> (0) Extracting WRF Model Output Variables
>>
>> Variable: U
>> Type: float
>> Total Size: 6844488 bytes
>> 1711122 values
>> Number of Dimensions: 4
>> Dimensions and sizes: [time | 1] x [lev | 14] x [lat | 311] x [lon |
>> 393]
>> Coordinates:
>> time: [0..0]
>> lev: [1000..50]
>> lat: [-18.40244..35.33525]
>> lon: [-35.74922..35.74922]
>> Number Of Attributes: 6
>> _FillValue : 9.96921e+36
>> description : y-wind component
>> units : m s-1
>> coordinates : lat lon
>> stagger :
>> FieldType : 104
>> Error in "cremapbin":
>> Could not map global lat array into grid array
>> Segmentation fault
>>
>>
>> ;;;;;;;;;;;;;;;; NCL SCRIPT
>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;**;;;;;;;;;;;;;;;;;;;
>>
>> ;*******************************************************************
>> ; Load Libaries
>> ;*******************************************************************
>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/gsn_csm.ncl"
>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/contributed.**ncl"
>>
>> begin
>> ;*******************************************************************
>> ; Change these parts
>> ;*******************************************************************
>> DATADir = "./" ;input directory
>> Files = "wrfout_d01*_wind_diagnostics.**nc<http://wind_diagnostics.nc>
>> <http://wind_diagnostics.nc>" ;input files
>>
>> ncdiro = "./" ;output directory
>> NCFILE1 = "wind_diagnostics" ;output file
>> ;*******************************************************************
>> ; Read in WRF Files
>> ;*******************************************************************
>> ;FILES = systemfunc (" ls -1 " + DATADir + Files) ; file paths
>> FILES = systemfunc (" ls -1 " + Files) ; file paths
>> numFILES = dimsizes(FILES)
>>
>> print("*********************************************************
>> **********")
>> print("These are the WRF Model output files you are extracting from.")
>> print("number of FILES = " + numFILES)
>> ;print(FILES)
>>
>> print("*********************************************************
>> **********")
>> print ("Extracting WRF Model Output Variables")
>> ;****************************************
>> ; DO LOOP: Extraction
>> ;****************************************
>> ; Use wrf_user_getvar function to extract each variable.
>> ; "-1" means all times, which doesn't matter here
>> ; b/c there is only 1 time in each file
>> ;**************************************
>> deltaTime = ispan(0,numFILES-1,1)
>> do ifil = 0,0,1 ;0,numFILES-1
>> a = addfile(FILES(ifil)+".nc","r") ;read individual wrfout file
>> U = a->U
>> V = a->V
>> printVarSummary(U)
>>
>> ;**************************************************
>> ; interpolate to new grid (NCEP R2 2.5 GRID --> 73,144)
>> ;*************************************************
>> loni = U&lon
>> lati = U&lat
>> newlat = fspan(-90.,90,73)
>> newlon = fspan(0.,359.,144)
>> newlat@units = "degrees_north"
>> newlon@units = "degrees_east"
>>
>>
>> Uo = area_conserve_remap_Wrap (loni,lati,U, newlon,newlat, False)
>> printVarSummary(Uo)
>>
>> end do
>> end
>>
>>
>>
>>
>>
>> ______________________________**_________________________
>> Erik Noble
>> Center for Climate Systems Research, Columbia University
>> NASA Goddard Institute for Space Studies
>> 2880 Broadway, New York, NY-10025, USA
>>
>> Tel: +1 (212) 678-5512; Fax: +1 (212) 678-5552
>> E-mail: enoble@giss.nasa.gov <mailto:enoble@giss.nasa.gov>,
>> en2192@columbia.edu <mailto:en2192@columbia.edu>
>>
>> WEB: http://www.giss.nasa.gov/**staff/enoble.html<http://www.giss.nasa.gov/staff/enoble.html>
>> GISS: http://www.giss.nasa.gov
>> CCSR: http://www.ccsr.columbia.edu
>> ______________________________**_________________________
>>
>>
>>
>> ______________________________**_________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Sep 24 18:27:47 2011