[ncl-talk] TKE vs time and height, wrf_user_intrp3d doesn't work

[Ehsan Taghizadeh]

 Dear Mary,
I have same problem as Alina. May I answer your question about printVarSummary(z,tke)?Following are the answers:

Variable: tkeType: floatTotal Size: 1008886320 bytes            252221580 valuesNumber of Dimensions: 4Dimensions and sizes:   [Time | 61] x [bottom_top_stag | 31] x [south_north | 342] x [west_east | 390]Coordinates:Number Of Attributes: 6  FieldType :   104  MemoryOrder : XYZ  description : TKE from PBL  units :       m2 s-2  stagger :     Z  coordinates : XLONG XLAT XTIME
Variable: zType: floatTotal Size: 976341600 bytes            244085400 valuesNumber of Dimensions: 4Dimensions and sizes:   [Time | 61] x [bottom_top | 30] x [south_north | 342] x [west_east | 390]Coordinates:Number Of Attributes: 6  coordinates : XLONG XLAT  stagger :  units :       m  description : Height  MemoryOrder : XYZ  FieldType :   104       
It seems the only difference is about (bottom_top) and (bottom_top_stag) dimension names.I think "wrf_user_interp3d" doesn't work for TKE_PBL? I also used this script (http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/LEVELS_INTERP/wrf_Height2.ncl) for "TKE_PBL" variable of WRF, but same error has occured:
"mension size mismatch, dimension (1) of left hand side reference does not have the same size as the right hand side reference after subscripting."
SincerelyEhsan    On Tuesday, October 16, 2018, 9:31:43 PM GMT+3:30, Mary Haley <haley at ucar.edu> wrote:  
 
 Hi Alina,
I'm guessing that the error is in how wrf_user_intrp3d is being called.  If you look at the documentation:
http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_intrp3d.shtml

and read the section on input arguments, it takes that the first two input arrays must be the same size, which in your case is tke and z. Did you do a printVarSummary on z to see what size it is? 
--Mary

On Thu, Oct 11, 2018 at 9:41 AM, Schuster Alina <alina-galchenko at yandex.ru> wrote:


Dear colleagues!

I am using a scm wrf model (an idealised case) and need to plot a variable TKE (turbulent kinetic energy) vs time and height. But I get an error

fatal:Dimension size mismatch, dimension (1) of left hand side reference does not have the same size as the right hand side reference after subscripting.
fatal:["Execute.c":8640]: Execute: Error occurred at or near line 246 in file /home/alina/NCL64/lib/ncarg/ nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8640]: Execute: Error occurred at or near line 29 in file TKE.ncl

I don't know how to fix it.

Line 246 contains:
 if ( nd .eq. 4 ) then
       var3(0,:,:,:,:) = var3d(:,:,:,:)
          z(0,:,:,:,:) = z_in(:,:,:,:)
     end if

Don't know what's wrong. It can be the problem is with the variable TKE and its 4 dimensions:
Variable: tke
Type: float
Total Size: 57600 bytes
            14400 values
Number of Dimensions: 4
Dimensions and sizes: [Time | 60] x [bottom_top_stag | 60] x [south_north | 2] x [west_east | 2]
Coordinates:
            Time: [ 0.. 59]
Number Of Attributes: 6
  FieldType : 104
  MemoryOrder : XYZ
  description : TKE from PBL
  units : J
  stagger : Z
  coordinates : XLONG XLAT XTIME

May be someone could advise me an alternative to the command wrf_user_intrp3d or may be someone has a ready script for TKE vs height and time (gsn_csm_contour).
Thanks in advance!

Best Regards,
Alina Schuster

______________________________ _________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/ mailman/listinfo/ncl-talk



_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181016/23cd9f8d/attachment.html>