<div dir="ltr"><div><div><div><div><div><div><div>[1] General background<br></div><div><br>NCL is a strongly typed language. The default rules will not allow <br><br></div> a = (/1 4, 5/) ; a(3)<br></div> b = (/2,5,6,9,10/) ; b(5)<br><br></div> a = b <br><br></div>left hand side is size 3; right hand size is size 5. This violate an NCL 'strongly typed' rule and would yield<br><br>Fatal:Dimension sizes of left hand side do not match right hand side<br><br>-----<br><br><br></div>FYI: NCL (6.1.1 onward) syntax contains the := reassignment operator which over-rides the default rule.<br><br></div> a := b ; this would be allowed via reassignment<br><br></div>---- 6.1.1 documentation; reassignment allows different sizes & types to be reassigned.<br><div><p>
<b>New reassignment operator</b></p><p>
There's a new operator, ":=", that allows you to reassign variables
without having to delete them first.<br>
For example:
</p><p>
</p><pre> x = fspan(-10.5,20.8,0.1) ; array of floats
; x = "This is a string" ; this won't work
x <b>:=</b> "This is a string" ; single string
</pre>
<p>
This operator works for attributes as well, but not for coordinate arrays.
</p><p>
<a name="m_892631072938679669_NewFunction6.1.1"></a>
</p><p>
++++++++++++++++++++++++++++++<wbr>+++++++++</p><p><br></p><p>[2]<b> </b>o, how does [1] affect your code?<br></p><div><br></div><div><div>fout->U250 = (/u250/)</div><div>fout->U500 = (/u500/)</div><div>fout->V250 = (/v250/)</div><div>fout->V500 = (/v500/)<br><br></div><div>It is my speculation, that the predefined 'U250' dimension sizes for the file referenced by 'fout' <br><br>dimSizes = (/-1,ny,nx,ny+1,nx+1,19/)<br><br>are of different sizes than the (/u250/) variable dimension sizes.<br><br></div><div>do a <br><br>printVarSummary(u250)<br><br></div><div>or <br><br></div><div>dim250 = dimsizes(u250)<br></div><div>print(dim250)<br><br></div><div>To check the dimension sizes.<br><br><br></div><div>Possibly, even though the script failed, a netCDF file was created. <br>Do a 'ncdump -h' on the file. Look at the dimension sizes.<br><br></div><div>NOTE: the := will *not* work for this issue.<br><br></div><div>Good luck<br></div><div><br><br></div></div><div><br></div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 5, 2017 at 1:31 PM, Sébastien Marinier <span dir="ltr"><<a href="mailto:marinier.sebastien@gmail.com" target="_blank">marinier.sebastien@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear NCL users,<div><br></div><div>I've been using NCL for post-processing some WRF output files: the script "crops" the domain and interpolates some variables at some pressure levels. When I run my script, I get this error:</div><div><br></div><div><div>fatal:Dimension sizes of left hand side do not match right hand side</div><div>fatal:["Execute.c":8575]:<wbr>Execute: Error occurred at or near line 140 in file test_crop_3d.ncl</div><div><br></div><div>fatal:Dimension sizes of left hand side do not match right hand side</div><div>fatal:["Execute.c":8575]:<wbr>Execute: Error occurred at or near line 141 in file test_crop_3d.ncl</div><div><br></div><div>fatal:Dimension sizes of left hand side do not match right hand side</div><div>fatal:["Execute.c":8575]:<wbr>Execute: Error occurred at or near line 142 in file test_crop_3d.ncl</div><div><br></div><div>fatal:Dimension sizes of left hand side do not match right hand side</div><div>fatal:["Execute.c":8575]:<wbr>Execute: Error occurred at or near line 143 in file test_crop_3d.ncl</div></div><div><br></div><div>I looked at this kind of error in the archives, but didn't find any solution. Lines 140 to 143 are the following lines:</div><div><br></div><div><div>fout->U250 = (/u250/)</div><div>fout->U500 = (/u500/)</div><div>fout->V250 = (/v250/)</div><div>fout->V500 = (/v500/)</div></div><div><br></div><div>Here's the full script:</div><div><br></div><div><div>;-----------------------------<wbr>------------------------------<wbr>-</div><div>;Getting coordinates and indices to crop the domain</div><div>;-----------------------------<wbr>------------------------------<wbr>-</div><div>fname_coords = "./<a href="http://wrfout_constants.nc" target="_blank">wrfout_constants.nc</a>"</div><div>f_coords = addfile(fname_coords,"r")</div><div><br></div><div>lon = f_coords->XLONG(0,:,:)</div><div>lat = f_coords->XLAT(0,:,:)</div><div><br></div><div>latN = 60.0</div><div>latS = 40.0</div><div>lonW = -140.0</div><div>lonE = -85.0</div><div><br></div><div>hpa = 0.01</div><div><br></div><div>ij = region_ind(lat,lon,latS,latN,<wbr>lonW,lonE)</div><div><br></div><div>iStrt = ij(0)</div><div>iLast = ij(1)</div><div>jStrt = ij(2)</div><div>jLast = ij(3)</div><div><br></div><div>new_lat = lat(iStrt:iLast,jStrt:jLast)</div><div>new_lon = lon(iStrt:iLast,jStrt:jLast)</div><div><br></div><div>;-----------------------------<wbr>------------------------------<wbr>-</div><div>;Opening files</div><div>;-----------------------------<wbr>------------------------------<wbr>-</div><div>f = addfile("wrf3d_d01_2013-06-01_<wbr>00:00:00","r")</div><div><br></div><div>;-----------------------------<wbr>------------------------------<wbr>--</div><div>;Create new netCDF file</div><div>;-----------------------------<wbr>------------------------------<wbr>--</div><div>dims = dimsizes(new_lat)</div><div>ny = dims(0)</div><div>nx = dims(1)</div><div><br></div><div>i=0</div><div><br></div><div>time = f->Times(i,:)</div><div>tk = f->TK(i,:,iStrt:iLast,jStrt:<wbr>jLast)</div><div>p = f->P(i,:,iStrt:iLast,jStrt:<wbr>jLast)</div><div>z_tmp = f->Z(i,:,iStrt:iLast,jStrt:<wbr>jLast)</div><div>qvapor = f->QVAPOR(i,:,iStrt:iLast,<wbr>jStrt:jLast)</div><div>u_tmp = f->U(i,:,iStrt:iLast,jStrt:<wbr>jLast)</div><div>v_tmp = f->V(i,:,iStrt:iLast,jStrt:<wbr>jLast)</div><div><br></div><div>;Unstagger some variables</div><div>z = wrf_user_unstagger(z_tmp,z_<wbr>tmp@stagger)</div><div>;u = wrf_user_unstagger(u_tmp,u_<wbr>tmp@stagger)</div><div>;v = wrf_user_unstagger(v_tmp,v_<wbr>tmp@stagger)</div><div>u = u_tmp</div><div>v = v_tmp</div><div><br></div><div>;Compute RH</div><div>rh = wrf_rh(qvapor,p,tk)</div><div><br></div><div>;Interpolation at pressure levels</div><div>;250 mb</div><div>u250 = wrf_user_intrp3d(u,p*hpa,"h",<wbr>250,0.,False)</div><div>v250 = wrf_user_intrp3d(v,p*hpa,"h",<wbr>250,0.,False)</div><div>z250 = wrf_user_intrp3d(z,p*hpa,"h",<wbr>250,0.,False)</div><div>t250 = wrf_user_intrp3d(tk,p*hpa,"h",<wbr>250,0.,False)</div><div>rh250 = wrf_user_intrp3d(rh,p*hpa,"h",<wbr>250,0.,False)</div><div>;500 mb</div><div>u500 = wrf_user_intrp3d(u,p*hpa,"h",<wbr>500,0.,False)</div><div>v500 = wrf_user_intrp3d(v,p*hpa,"h",<wbr>500,0.,False)</div><div>z500 = wrf_user_intrp3d(z,p*hpa,"h",<wbr>500,0.,False)</div><div>t500 = wrf_user_intrp3d(tk,p*hpa,"h",<wbr>500,0.,False)</div><div>rh500 = wrf_user_intrp3d(rh,p*hpa,"h",<wbr>500,0.,False)</div><div>;sfc</div><div>psfc = p(0,:,:)</div><div>slp = wrf_slp(z,tk,p,qvapor)</div><div>rh0 = rh(0,:,:)</div><div><br></div><div>;Output Files</div><div>diro = "./"</div><div>filo = "wrf_lvls_"+time+".nc"</div><div>system("/bin/rm -f "+diro+filo) ; remove if exists</div><div>fout = addfile(diro+filo,"c") ; open output file</div><div><br></div><div>setfileoption(fout,"<wbr>DefineMode",True)</div><div><br></div><div>;Global attributes</div><div>fAtt = True</div><div>fAtt@title = "Cropped and interpolated (sfc, 500 mb and 250 mb) data files from Liu et al. (2016) simulations"</div><div>fAtt@Conventions = "None"</div><div>fAtt@creation_date = systemfunc("date")</div><div>fileattdef(fout,fAtt)</div><div><br></div><div>;Dimensions of coordinates</div><div>dimNames = (/"Time","south_north","west_<wbr>east","south_north_stag","<wbr>west_east_stag","DateStrLen"/)</div><div>dimSizes = (/-1,ny,nx,ny+1,nx+1,19/)</div><div>dimUnlim = (/True,False,False,False,<wbr>False,False/)</div><div>filedimdef(fout,dimNames,<wbr>dimSizes,dimUnlim)</div><div><br></div><div>;... of variables</div><div>filevardef(fout,"Times",<wbr>typeof(time),getvardims(time))</div><div>filevardef(fout,"XLAT",typeof(<wbr>new_lat),getvardims(new_lat))</div><div>filevardef(fout,"XLONG",<wbr>typeof(new_lon),getvardims(<wbr>new_lon))</div><div>filevardef(fout,"U250",typeof(<wbr>u250),getvardims(u250))</div><div>filevardef(fout,"U500",typeof(<wbr>u500),getvardims(u500))</div><div>filevardef(fout,"V250",typeof(<wbr>v250),getvardims(v250))</div><div>filevardef(fout,"V500",typeof(<wbr>v500),getvardims(v500))</div><div>filevardef(fout,"Z250",typeof(<wbr>z250),getvardims(z250))</div><div>filevardef(fout,"Z500",typeof(<wbr>z500),getvardims(z500))</div><div>filevardef(fout,"T250",typeof(<wbr>t250),getvardims(t250))</div><div>filevardef(fout,"T500",typeof(<wbr>t500),getvardims(t500))</div><div>filevardef(fout,"RH250",<wbr>typeof(rh250),getvardims(<wbr>rh250))</div><div>filevardef(fout,"RH500",<wbr>typeof(rh500),getvardims(<wbr>rh500))</div><div>filevardef(fout,"RH0",typeof(<wbr>rh0),getvardims(rh0))</div><div>filevardef(fout,"PSFC",typeof(<wbr>psfc),getvardims(psfc))</div><div>filevardef(fout,"SLP",typeof(<wbr>slp),getvardims(slp))</div><div><br></div><div>;Copy attributes</div><div>;filevarattdef(fout,"Times",<wbr>time)</div><div>filevarattdef(fout,"XLAT",new_<wbr>lat)</div><div>filevarattdef(fout,"XLONG",<wbr>new_lon)</div><div>filevarattdef(fout,"U250",<wbr>u250)</div><div>filevarattdef(fout,"U500",<wbr>u500)</div><div>filevarattdef(fout,"V250",<wbr>v250)</div><div>filevarattdef(fout,"V500",<wbr>v500)</div><div>filevarattdef(fout,"Z250",<wbr>z250)</div><div>filevarattdef(fout,"Z500",<wbr>z500)</div><div>filevarattdef(fout,"T250",<wbr>t250)</div><div>filevarattdef(fout,"T500",<wbr>t500)</div><div>filevarattdef(fout,"RH250",<wbr>rh250)</div><div>filevarattdef(fout,"RH500",<wbr>rh500)</div><div>filevarattdef(fout,"RH0",rh0)</div><div>filevarattdef(fout,"PSFC",<wbr>psfc)</div><div>filevarattdef(fout,"SLP",slp)</div><div><br></div><div>setfileoption(fout,"<wbr>DefineMode",False)</div><div><br></div><div>;Output variables</div><div>fout->Times = (/time/)</div><div>fout->XLAT = (/new_lat/)</div><div>fout->XLONG = (/new_lon/)</div><div>fout->U250 = (/u250/)</div><div>fout->U500 = (/u500/)</div><div>fout->V250 = (/v250/)</div><div>fout->V500 = (/v500/)</div><div>fout->Z250 = (/z250/)</div><div>fout->Z500 = (/z500/)</div><div>fout->T250 = (/t250/)</div><div>fout->T500 = (/t500/)</div><div>fout->RH250 = (/rh250/)</div><div>fout->RH500 = (/rh500/)</div><div>fout->RH0 = (/rh0/)</div><div>fout->PSFC = (/psfc/)</div><div>fout->SLP = (/slp/)</div></div><div>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>-------------</div><div><br></div><div>I think there is a problem with the staggered coordinates, but I can't put the finger on it...</div><div><br></div><div>Thanks for the help,</div><div>Sebastien.</div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>