<div dir="ltr">Hi Toni,<br><br> yeah, that worked. Thank You very much. But one more error still exists. Could you please check that also.<br><br><br>fatal:Dimension (1) of (vort_hov) is not named and therefore<br>doesn't have an associated coordinate variable<br><br>fatal:["Execute.c":8637]:Execute: Error occurred at or near line 18 in<br>file hov_sv2cv_2015_test.ncl<br><br>1 f = addfile("curv_vort_lat_lon_925hPa_20150101_20151231.nc","r")<br><br>2 time=f->time<br><br>3 lon=f->longitude<br><br>4 lat=f->latitude<br><br>5 u = f->cur_vort(:,:,:)<br><br>6 printVarSummary(u)<br><br>7 vort_hov = new((/dimsizes(time),dimsizes(lon)/),float) ; this is<br>the 2-dimensional variable you use for the Hovmoeller diagram<br><br>8 printVarSummary(vot_hov)<br><br>9 do time_step = 0,dimsizes(time)-1<br><br>10 do lon_step = 0,dimsizes(lon)-1<br><br>11 vort_timeseries = avg(u(time_step,:,lon_step)) ; a loop to<br>average the data by latitude for each time step and longitudinal step<br>12 end do<br>13 end do<br><br>14 res=True<br><br>15 res@cnFillOn = True ; turn on color fill<br><br>16 res@cnFillPalette = "BlWhRe"<br><br>17 res@gsnLeftString="cv"<br><br>18 <font color="#ff0000">plot = gsn_csm_hov(wks, vort_timeseries(:,{60:120}), res )</font><br></div><br><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jan 30, 2019 at 2:05 AM Toni Klemm <<a href="mailto:toni-klemm@tamu.edu">toni-klemm@tamu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">Hi Saran,<div><br></div><div>I think the errors all occurred because of the first one. Change </div><div><br></div><div><blockquote type="cite">9 do time_step = 0,dimsizes(time)</blockquote>and </div><div><blockquote type="cite">10 do lon_step = 0,dimsizes(lon)</blockquote><div><br></div>to</div><div><br></div><div>do time_step = 0,dimsizes(time)<b><u>-1</u></b></div><div>and<br><div><div>
<div dir="auto" style="text-align:start;text-indent:0px;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="text-align:start;text-indent:0px;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="text-align:start;text-indent:0px;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="text-align:start;text-indent:0px;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="text-align:start;text-indent:0px;word-wrap:break-word;line-break:after-white-space"><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px">do lon_step = 0,dimsizes(lon)<b><u>-1</u></b></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px"><font face="Calibri"><br></font></div><div style="text-align:start;text-indent:0px"><font face="Calibri">This should let you run the do-loops without crashing. </font></div><div style="text-align:start;text-indent:0px"><font face="Calibri"><br></font></div><div style="text-align:start;text-indent:0px"><font face="Calibri">To explain the “-1”, the position of the first item in a list in NCL is <i>zero</i>, not one. So the position of the last item in a list is not the length of the list (“</font>dimsizes()”), but “dimsizes()-1”. Without “-1” you were trying to call on a list item that wasn’t there, which is what “subscript <i>out of range</i>” means.</div><div style="text-align:start;text-indent:0px"><font face="Calibri"><br></font></div><div style="text-align:start;text-indent:0px"><font face="Calibri">Best,</font></div><div style="text-align:start;text-indent:0px"><font face="Calibri">Toni</font></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px"><font face="Calibri"><br class="gmail-m_1105898657825178436Apple-interchange-newline"><br></font></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px"><font face="Calibri"><b>Toni Klemm, Ph.D.<br></b>Postdoctoral Research Associate</font></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px">Department of Ecosystem Science and Management</div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px">College of Agriculture and Life Sciences</div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px">Texas A&M University, College Station, TX</div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px"><font face="Calibri">Contributor to the <a href="http://www.eccforum.org" target="_blank">Early Career Climate Forum</a><br><a href="http://www.toni-klemm.de" target="_blank">www.toni-klemm.de</a> | <a href="http://twitter.com/toniklemm" target="_blank">@toniklemm</a></font></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px"><font face="Calibri"><br></font></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-variant-caps:normal;text-align:start;text-indent:0px"><font face="Calibri"><br></font></div></div></div><br class="gmail-m_1105898657825178436Apple-interchange-newline"></div><br class="gmail-m_1105898657825178436Apple-interchange-newline"></div><br class="gmail-m_1105898657825178436Apple-interchange-newline"></div><br class="gmail-m_1105898657825178436Apple-interchange-newline"><br class="gmail-m_1105898657825178436Apple-interchange-newline">
</div>
<div><br><blockquote type="cite"><div>On Jan 29, 2019, at 2:09 PM, Saran Rajendran <<a href="mailto:happysaran.raj@gmail.com" target="_blank">happysaran.raj@gmail.com</a>> wrote:</div><br class="gmail-m_1105898657825178436Apple-interchange-newline"><div><div>Hi Toni<br><br> Thank You for that valuable information. I thoroughly read<br>your code and I understood my mistake. I modified my code as per your<br>correction, but I'm getting some errors again.<br><br>These are the errors:<br><br>fatal:Subscript out of range, error in subscript #2<br><br>fatal:An error occurred reading u<br><br>fatal:["Execute.c":8637]:Execute: Error occurred at or near line 11 in<br>file hov_sv2cv_2015_test.ncl<br><br><br>fatal:Dimension (1) of (vort_timeseries) is not named and therefore<br>doesn't have an associated coordinate variable<br><br>fatal:["Execute.c":8637]:Execute: Error occurred at or near line 18 in<br>file hov_sv2cv_2015_test.ncl<br><br><br><br>And the modified code is<br><br>1 f = addfile("curv_vort_lat_lon_925hPa_20150101_20151231.nc","r")<br><br>2 time=f->time<br><br>3 lon=f->longitude<br><br>4 lat=f->latitude<br><br>5 u = f->cur_vort(:,:,:)<br><br>6 printVarSummary(u)<br><br>7 vort_hov = new((/dimsizes(time),dimsizes(lon)/),float) ; this is<br>the 2-dimensional variable you use for the Hovmoeller diagram<br><br>8 printVarSummary(vot_hov)<br><br>9 do time_step = 0,dimsizes(time)<br><br>10 do lon_step = 0,dimsizes(lon)<br><br>11 vort_timeseries = avg(u(time_step,:,lon_step)) ; a loop to<br>average the data by latitude for each time step and longitudinal step<br>12 end do<br>13 end do<br><br>14 res=True<br><br>15 res@cnFillOn = True ; turn on color fill<br><br>16 res@cnFillPalette = "BlWhRe"<br><br>17 res@gsnLeftString="cv"<br><br>18 plot = gsn_csm_hov(wks, vort_timeseries(:,{60:120}), res )<br><br>On Tue, Jan 29, 2019 at 7:55 PM Toni Klemm <<a href="mailto:toni-klemm@tamu.edu" target="_blank">toni-klemm@tamu.edu</a>> wrote:<br><blockquote type="cite"><br>Hi Saran,<br><br>That’s for two reasons:<br><br>(1) You are only using the first time dimension of your data, but all latitude and all longitude dimensions, so by your own setting you're plotting latitude/longitude, not time/longitude:<br><br>u = f->cur_vort(0,:,:)<br><br><br>(2) You need to prepare the data in the correct way first, so that you can plot it. Your code below only does the plotting, but not the data manipulation so that you can actually plot what you need. Average the latitude dimension (y axis) for each longitudinal grid cell “column” and each time step, and then plot the data.<br>This might work (though, it might not. I didn’t test it. Please read through it and try to understand it. Don’t just plug it into your script. Also, there might be an easier way to do this):<br><br>u = f->cur_vort(:,:,:) ; u contains all your (3-dimensional) curvature vorticity data<br><br>vort_timeseries = new((/dimsizes(time),dimsizes(lon)/),float) ; this is the 2-dimensional variable you use for the Hovmoeller diagram<br><br>do time_step = 0,dimsizes(time)<br><br> do lon_step = 0,dimsizes(lon)<br><br> vort_timeseries(time_step,lon_step) = avg(u(time_step,:,lon_step)) ; a loop to average the data by latitude for each time step and longitudinal step<br><br> end do ; lon_step loop<br><br>end do ; time_step loop<br><br><br>vort_timeseries should be the data you need for your Hovmoeller diagram.<br><br><br>Good luck,<br>Toni<br><br><br>Toni Klemm, Ph.D.<br>Postdoctoral Research Associate<br>Department of Ecosystem Science and Management<br>College of Agriculture and Life Sciences<br>Texas A&M University, College Station, TX<br>Contributor to the Early Career Climate Forum<br><a href="http://www.toni-klemm.de" target="_blank">www.toni-klemm.de</a> | @toniklemm<br><br><br><br><br><br><br>On Jan 29, 2019, at 5:16 AM, Saran Rajendran <<a href="mailto:happysaran.raj@gmail.com" target="_blank">happysaran.raj@gmail.com</a>> wrote:<br><br>Dear NCL users,<br><br> When I tried this script to plot hovmueller diagram<br>(time vs longitude) of curvature vorticity, the code runs. But instead<br>of time in y-axis, it shows latitude. I want time in the y-axis.<br><br><br>f = addfile("curv_vort_lat_lon_925hPa_20150101_20151231.nc","r")<br><br>u = f->cur_vort(0,:,:)<br><br> wks = gsn_open_wks("png","hovmueller_cv_2015")<br><br> res=True<br><br> res@cnFillOn = True ; turn on color fill<br><br>res@cnFillPalette = "BlWhRe"<br><br>res@gsnLeftString="cv"<br><br>plot = gsn_csm_hov(wks, u(:,{60:120}), res )<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" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br><br><br></blockquote><br><br>-- <br>Saran Rajendran<br>Junior Research Fellow<br>Centre for Atmospheric Science<br>Indian Institute of Technology Delhi<br></div></div></blockquote></div><br></div></div></div></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Saran Rajendran<br>Junior Research Fellow<br>Centre for Atmospheric Science<br>Indian Institute of Technology Delhi<br></div>