<div dir="ltr">Hello:<div><br></div><div>Have you defined var3 and z prior? If so, you can check the dimensions with printVarSummary(var_name).</div><div><br></div><div>I don't see why you need the "0" as the first index, since it is just a placeholder.</div><div><br></div><div>One thing to try would be just to set var3 = var3d(:,:,:,:) just to see that you are accessing var3d correctly.</div><div><br></div><div>Barry</div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Oct 11, 2018 at 6:41 PM Schuster Alina <<a href="mailto:alina-galchenko@yandex.ru">alina-galchenko@yandex.ru</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear colleagues!<br>
<br>
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<br>
<br>
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.<br>
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 246 in file /home/alina/NCL64/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl<br>
<br>
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 29 in file TKE.ncl<br>
<br>
I don't know how to fix it.<br>
<br>
Line 246 contains:<br>
if ( nd .eq. 4 ) then<br>
var3(0,:,:,:,:) = var3d(:,:,:,:)<br>
z(0,:,:,:,:) = z_in(:,:,:,:)<br>
end if<br>
<br>
Don't know what's wrong. It can be the problem is with the variable TKE and its 4 dimensions:<br>
Variable: tke<br>
Type: float<br>
Total Size: 57600 bytes<br>
14400 values<br>
Number of Dimensions: 4<br>
Dimensions and sizes: [Time | 60] x [bottom_top_stag | 60] x [south_north | 2] x [west_east | 2]<br>
Coordinates:<br>
Time: [ 0.. 59]<br>
Number Of Attributes: 6<br>
FieldType : 104<br>
MemoryOrder : XYZ<br>
description : TKE from PBL<br>
units : J<br>
stagger : Z<br>
coordinates : XLONG XLAT XTIME<br>
<br>
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).<br>
Thanks in advance!<br>
<br>
Best Regards,<br>
Alina Schuster<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Associate Scientist, Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>