Folks,
I am brand-spanking-new to NCL (as well as to WRF
output) having just installed it (after the, unh!,
globus toolkit) over the weekend.
I wrote out slp files after taking data from original
wrf*.nc files (as the "to slp" calculations were very
slow) to separate *nc files, the "ncdump -h" output
from one of which is below.
Upon opening these files from the routine also found
below, I get a 'subscripts used vs expected' message,
Yet the dimensions appear to be consistent and correct
in the printvarsum output.
It may be that the wrf specific routines expect there
to be a time dimension in all cases. How does one
handle this in the case of e.g. means, variance ...?
Please advise. Thanks!! Jack
=============
NCDUMP -h
=============
netcdf slp00Zbar092000 {
dimensions:
south-north = 290 ;
west-east = 175 ;
variables:
float XLAT(south-north, west-east) ;
XLAT:stagger = "" ;
XLAT:units = "degree_north" ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:MemoryOrder = "XY" ;
XLAT:FieldType = 104 ;
float XLONG(south-north, west-east) ;
XLONG:stagger = "" ;
XLONG:units = "degree_east" ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:MemoryOrder = "XY" ;
XLONG:FieldType = 104 ;
float SLP00(south-north, west-east) ;
SLP00:coordinates = "XLONG XLAT" ;
SLP00:stagger = "" ;
SLP00:units = "hPa" ;
SLP00:description = "Sea Level Pressure" ;
SLP00:MemoryOrder = "XY" ;
SLP00:FieldType = 104 ;
// global attributes:
:TITLE = " OUTPUT FROM WRF V3.4.1 MODEL" ;
:START_DATE = "2000-09-01_00:00:00" ;
:SIMULATION_START_DATE = "1999-10-01_00:00:00" ;
:WEST-EAST_GRID_DIMENSION = 176 ;
:SOUTH-NORTH_GRID_DIMENSION = 291 ;
:BOTTOM-TOP_GRID_DIMENSION = 28 ;
:DX = 10000.f ;
:DY = 10000.f ;
:GRIDTYPE = "C" ;
:DIFF_OPT = 1 ;
:KM_OPT = 4 ;
:DAMP_OPT = 3 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 3 ;
:RA_LW_PHYSICS = 1 ;
:RA_SW_PHYSICS = 1 ;
:SF_SFCLAY_PHYSICS = 1 ;
:SF_SURFACE_PHYSICS = 2 ;
:BL_PBL_PHYSICS = 1 ;
:CU_PHYSICS = 1 ;
:SURFACE_INPUT_SOURCE = 1 ;
:SST_UPDATE = 1 ;
:GRID_FDDA = 0 ;
:GFDDA_INTERVAL_M = 0 ;
:GFDDA_END_H = 0 ;
:GRID_SFDDA = 0 ;
:SGFDDA_INTERVAL_M = 0 ;
:SGFDDA_END_H = 0 ;
:HYPSOMETRIC_OPT = 2 ;
:SF_URBAN_PHYSICS = 0 ;
:SHCU_PHYSICS = 0 ;
:MFSHCONV = 0 ;
:FEEDBACK = 0 ;
:SMOOTH_OPTION = 0 ;
:SWRAD_SCAT = 1.f ;
:W_DAMPING = 1 ;
:MOIST_ADV_OPT = 1 ;
:SCALAR_ADV_OPT = 1 ;
:TKE_ADV_OPT = 1 ;
:DIFF_6TH_OPT = 2 ;
:DIFF_6TH_FACTOR = 0.12f ;
:OBS_NUDGE_OPT = 0 ;
:BUCKET_MM = -1.f ;
:BUCKET_J = -1.f ;
:PREC_ACC_DT = 0.f ;
:OMLCALL = 0 ;
:ISFTCFLX = 0 ;
:ISHALLOW = 0 ;
:DFI_OPT = 0 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 175 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 176 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 290 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 291 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 27 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 28 ;
:GRID_ID = 2 ;
:PARENT_ID = 1 ;
:I_PARENT_START = 9 ;
:J_PARENT_START = 26 ;
:PARENT_GRID_RATIO = 5 ;
:DT = 54.f ;
:CEN_LAT = 40.30994f ;
:CEN_LON = -125.1265f ;
:TRUELAT1 = 36.583f ;
:TRUELAT2 = 0.f ;
:MOAD_CEN_LAT = 36.58298f ;
:STAND_LON = -120.646f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 0.f ;
:JULYR = 1999 ;
:JULDAY = 275 ;
:MAP_PROJ = 3 ;
:MMINLU = "USGS" ;
:NUM_LAND_CAT = 24 ;
:ISWATER = 16 ;
:ISLAKE = -1 ;
:ISICE = 24 ;
:ISURBAN = 1 ;
:ISOILWATER = 14 ;
}
..........................................
ncl.xq.fix running on Mac OS X 10.8.4
=============
OUTPUT:
=============
Copyright (C) 1995-2013 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.1.2
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
fatal:Number of subscripts do not match number of dimensions of variable, (3)
subscripts used, (2) subscripts expected
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4019 in file
$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4111 in file
$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 4194 in file
$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 131 in file
zncl_slpanom_wrf_cln.ncl
(0)
(0) --------------------------------END Sum slpbr
Variable: slpbr
Type: float
Total Size: 203000 bytes
50750 values
Number of Dimensions: 2
Dimensions and sizes: [south-north | 290] x [west-east | 175]
Coordinates:
Number Of Attributes: 6
coordinates : XLONG XLAT
stagger :
units : Pa
description : MSLP DAY MEAN SEP 2000
MemoryOrder : XY
FieldType : 104
(0)
(0) MSLP DAY MEAN SEP 2000: min=100970 max=102208
(0) ------
(0)
(0) --------------------------------END Sum XLAT
Variable: XLAT
Type: float
Total Size: 203000 bytes
50750 values
Number of Dimensions: 2
Dimensions and sizes: [south-north | 290] x [west-east | 175]
Coordinates:
Number Of Attributes: 5
FieldType : 104
MemoryOrder : XY
description : LATITUDE, SOUTH IS NEGATIVE
units : degree_north
stagger :
(0)
(0) LATITUDE, SOUTH IS NEGATIVE: min=26.8863 max=51.515
(0) ------
(0)
(0) --------------------------------END Sum XLONG
Variable: XLONG
Type: float
Total Size: 203000 bytes
50750 values
Number of Dimensions: 2
Dimensions and sizes: [south-north | 290] x [west-east | 175]
Coordinates:
Number Of Attributes: 5
FieldType : 104
MemoryOrder : XY
description : LONGITUDE, WEST IS NEGATIVE
units : degree_east
stagger :
(0)
(0) LONGITUDE, WEST IS NEGATIVE: min=-134.872 max=-115.381
(0) ------
=============
CODE
=============
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
f00 = addfile ("slp00Zbar092000.nc","r")
f06 = addfile ("slp06Zbar092000.nc","r")
f12 = addfile ("slp12Zbar092000.nc","r")
f18 = addfile ("slp18Zbar092000.nc","r")
system("/bin/rm -f slpDAYbar092000.nc")
ncbr = addfile ("slpDAYbar092000.nc", "c")
vattbs = getvaratts(f00) ; Any Var, inclu file ptr
vnames = getfilevarnames(f00)
fdmszs = getfiledimsizes(f00 )
fdnams = getvardims(f00 )
; print(" ")
; print("------------------------------------getfiledimsizes:-")
; print(fdmszs+" "+fdnams)
; print("------------------------------------getfiledimsizes:-")
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; We generate plots, but what kind do we prefer?
; type = "x11"
type = "pdf"
; type = "ps"
; type = "ncgm"
; type@wkColorModel = "cmyk" ; make output cmyk
wksbr = gsn_open_wks(type,"plt_slpDAYbar")
gsn_define_colormap(wksbr,"BlAqGrYeOrRe") ; choose colormap
; Set some Basic Plot options
res = True
res@MainTitle = "NCL WRF NWBY"
res@Footer = False
pltres = True
mpres = True
mpres@mpGeophysicalLineColor = "Black"
mpres@mpNationalLineColor = "Black"
mpres@mpUSStateLineColor = "Black"
mpres@mpGridLineColor = "Black"
mpres@mpLimbLineColor = "Black"
mpres@mpPerimLineColor = "Black"
mpres@mpGeophysicalLineThicknessF = 2.0
mpres@mpGridLineThicknessF = 2.0
mpres@mpLimbLineThicknessF = 2.0
mpres@mpNationalLineThicknessF = 2.0
mpres@mpUSStateLineThicknessF = 2.0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
time = -1 ; only likes -1
XLONG = wrf_user_getvar(f00 ,"XLONG",time) ;
XLAT = wrf_user_getvar(f00 , "XLAT",time) ;
slp00 = wrf_user_getvar(f00 ,"SLP00",time) ;
slp06 = wrf_user_getvar(f06 ,"SLP06",time) ;
slp12 = wrf_user_getvar(f12 ,"SLP12",time) ;
slp18 = wrf_user_getvar(f18 ,"SLP18",time) ;
; slpbr := slp00(:,:) * 0
; slp00a := slp00(:,:) * 0
; slp06a := slp06(:,:) * 0
; slp12a := slp12(:,:) * 0
; slp18a := slp18(:,:) * 0
slpbr = slp00(:,:) * 0
slp00a = slp00(:,:) * 0
slp06a = slp06(:,:) * 0
slp12a = slp12(:,:) * 0
slp18a = slp18(:,:) * 0
slpbr = 100 * ( slp00(:,:) + slp06(:,:) + slp12(:,:) + slp18(:,:) ) / 4
slpbr!0 = slp00!0 ;
slpbr!1 = slp00!1
slpbr@FieldType = slp00@FieldType
slpbr@MemoryOrder = slp00@MemoryOrder
slpbr@description = slp00@description
slpbr@units = "Pa"
slpbr@stagger = slp00@stagger
slpbr@coordinates = slp00@coordinates
XLAT!0 = slp00!0 ;
XLAT!1 = slp00!1
XLAT@FieldType = slp00@FieldType
XLAT@MemoryOrder = slp00@MemoryOrder
XLAT@stagger = slp00@stagger
XLONG!0 = slp00!0 ;
XLONG!1 = slp00!1
XLONG@FieldType = slp00@FieldType
XLONG@MemoryOrder = slp00@MemoryOrder
XLONG@stagger = slp00@stagger
optsbr = res
slpbr@description = "MSLP DAY SEP 2000"
slpbr@units = "Pa"
optsbr@cnLineColor = "Black"
optsbr@cnFillOn = True
optsbr@ContourParameters = (/ 100.0 /)
optsbr@cnInfoLabelOrthogonalPosF = 0.07 ; offset second label
information
optsbr@gsnContourLineThicknessesScale = 2.0
contour_slpbrF = wrf_contour(f00,wksbr,slpbr(:,:),optsbr)
optsbr@cnFillOn = False
optsbr@cnLineColor = "grey42"
optsbr@ContourParameters = (/ 100.0 /)
optsbr@cnInfoLabelOrthogonalPosF = 0.07 ; offset second label
information
optsbr@gsnContourLineThicknessesScale = 2.0
contour_slpbrL = wrf_contour(f00,wksbr,slpbr(:,:),optsbr)
plot =
wrf_map_overlays(f00,wksbr,(/contour_slpbrF,contour_slpbrL/),pltres,mpres)
print(" ")
print("--------------------------------END Sum slpbr ")
printVarSummary (slpbr)
printMinMax (slpbr,True)
print("------")
print(" ")
print("--------------------------------END Sum XLAT ")
printVarSummary (XLAT)
printMinMax (XLAT,True)
print("------")
print(" ")
print("--------------------------------END Sum XLONG ")
printVarSummary (XLONG)
printMinMax (XLONG,True)
print("------")
delete(optsbr)
Cheers,
Jack
------------------------------------------------------------------------------
Jack Ritchie / UCSD / SIO / CASPO
------------------------------------------------------------------------------
Mailing Address: Contact:
Scripps Institution of Oceanography Tel.(858)534-7515
University of California, San Diego Fax (858)534-8561
9500 Gilman Drive Nierenberg Hall 426
La Jolla, CA 92093-0224 jritchie@ucsd.edu
------------------------------------------------------------------------------
Express Delivery Address:
8810 Shellback Way
Nierenberg Hall 426
La Jolla, CA 92037
------------------------------------------------------------------------------
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 18 11:21:49 2013
This archive was generated by hypermail 2.1.8 : Tue Oct 01 2013 - 14:41:43 MDT