Opening multiple files and NCL wrf functions

Dear all,

I'm trying to open multiple wrf output files and concatenate the as one, following the example named "Spanning Multiple Files":

===============================
; WRF ARW open multiple files
; get list of all files and open as "one big file"
 all_files = systemfunc ("ls /Users/lcana/Dust/wrfout_d01_2004-03*")
 a = addfiles (all_files, "r")
; choose how files are combined and read in variable across files
 ListSetType (a, "cat") ; concatenate or "merge" (default)
===============================

Later, I call different wrf_user functions as the following:

===============================
times = wrf_user_list_times(a) ; get times in the file
ntimes = dimsizes(times) ; number of times in the file

hgt = wrf_user_getvar(a,"ter",it) ; get terrain height
   hgt@description = "Terrain height"
   hgt@units = "m"
====================================================

and so on. When I execute the script, I can see the following output:

Variable: times
Type: string
Total Size: 88 bytes
            11 values
Number of Dimensions: 1
Dimensions and sizes: [11]
Coordinates:
Number Of Attributes: 2
  description : times in file
  _FillValue : missing
(0) 2004-03-02_12:00:00
(1) 2004-03-02_18:00:00
(2) 2004-03-03_00:00:00
(3) 2004-03-03_06:00:00
(4) 2004-03-03_12:00:00
(5) 2004-03-03_18:00:00
(6) 2004-03-04_00:00:00
(7) 2004-03-04_06:00:00
(8) 2004-03-04_12:00:00
(9) 2004-03-04_18:00:00
(10) 2004-03-05_00:00:00

It looks OK, so I assume all things are right. But then there is a problem at this line:

cont_hgt=wrf_contour(a,wks,hgt,opts_h) ; create terrain plot for overlay, line 118

fatal:Argument type mismatch on argument (0) of (wrf_contour) can not coerce
fatal:Execute: Error occurred at or near line 118 in file Viento_en_niveles_z_brisas.ncl

What is wrong here? wrf_contour function is different?

Thanks in advance,

Luis Cana
ULPGC

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 29 08:43:46 2012