WRF Time Dimension when displaying Mean (err: 3subscr vs 2expect)

From: Jack Ritchie <jritchie_at_nyahnyahspammersnyahnyah>
Date: Wed Sep 18 2013 - 11:21:40 MDT

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