[ncl-talk] Cross sections other than west to east...

Thu Feb 20 10:53:43 MST 2020

Hi Paul,

The line numbers in any error messages should tell exactly where the
problem occurs -- the *why* might be elsewhere earlier in the code. In any
case, the error message in your first note indicated there was a
size/dimensions mismatch. Perhaps some strategically placed
print/printVarSummary statements will reveal why.

All that aside, I don't know what might have changed between NCL versions
that would break your code. What version were you using?

Rick

On Thu, Feb 20, 2020 at 10:29 AM Suffern Paul <paul.suffern at ntsb.gov> wrote:

> Hi Rick,
>
>
>
> Thanks for the advice… I switched it to the “wrf_user_ll_xy” instead but
> unfortunately that didn’t fix the issue… I’ve narrowed it down to this part
> of the script I think:
>
>
>
>     ;*******
>
>     ;*******  Loop over cross sections
>
>     ;*******
>
>     do ip = 0, 0           ; we are just going to do one crossection for
> now
>
>
>
>       ; Interpolate data vertically (in z)
>
>       px = p
>
>       do ll = 0,nlevs-1
>
>          w_plane (ll,:) = dim_avg_n( wrf_interp_3d_z( w ,px, pnew(ll)),
> idim )
>
>          th_plane(ll,:) = dim_avg_n( wrf_interp_3d_z( T ,px, pnew(ll)),
> idim )
>
>          th_plane(ll,:) = dim_avg_n( wrf_interp_3d_z( th,px, pnew(ll)),
> idim )
>
>       end do
>
>
>
> For whatever reason it is having in the “wrf_interp_3d_z” part and the
> dimensions… so I don’t know if I made an error in the “w_plane” etc… part
> of it, or something changed in the “wrf_interp_3d_z” in version 6.6 of NCL.
>
>
>
> Thanks again,
>
>
>
> Paul
>
>
>
> *From:* Rick Brownrigg <brownrig at ucar.edu>
> *Sent:* Thursday, February 20, 2020 11:37 AM
> *To:* Suffern Paul <paul.suffern at ntsb.gov>
> *Cc:* ncl-talk at ucar.edu
> *Subject:* Re: [ncl-talk] Cross sections other than west to east...
>
>
>
> [CAUTION] This email originated from outside of the organization. Do not
> click any links or open attachments unless you recognize the sender and
> know the content is safe.
>
> Hi Paul,
>
>
>
> Unfortunately, line 219, referenced in the error message, does not line up
> with anything in the script you sent that could possibly generate that
> error, and note that the filename referenced in the error is different that
> what you attached.
>
>
>
> "Guess"ing as to where the problem might be, the only real place where
> lat/lon get referenced is line 115, in a call to wrf_user_ll_ij.  If that's
> the case, then I noticed that the docs say that function is deprecated as
> of version 6.6.0, and to use wrf_user_ll_xy instead.  I don't know what the
> difference between them is, but perhaps that might resolve your issue
>
>
>
>     http://ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_ll_to_ij.shtml
>
>
>
> Hope that helps...
>
> Rick
>
>
>
> On Thu, Feb 20, 2020 at 9:08 AM Suffern Paul via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
> Good morning,
>
>
>
> I recently upgraded to NCL 6.6.2 (yes it takes me a while to upgrade to
> most current routine/versions). For the past few years I’ve used the
> attached script to make cross sections and it works out superbly well until
> this latest upgrade. I’ve noticed I can make west to east cross section no
> issues when adjusting the:
>
>
>
>   lats  =  (/60.5042, 60.5042/)
>
>   lons  =  (/-150.500, -149.833/)
>
>
>
> “lons” portion of the file. When I used the “lats” portion of the file
> still within the domain the file crashes giving this Error:
>
>
>
> fatal:Dimension size mismatch on subscript #1, left-hand and right-hand
> side dimensions do not match
>
> fatal:["Execute.c":8637]:Execute: Error occurred at or near line 219 in
> file XsectionPresRiLatLon_abt.ncl
>
>
>
> Again, I can change the “lons” file to whatever I want within the domain
> with no issues… Before this upgrade to NCL 6.6.2 I could adjust lats and
> lons to any point in the domain I wanted with no issues.
>
>
>
> Thanks for any ideas!
>
>
>
> Paul
>
>
>
> *Paul Suffern*
>
> Senior Meteorologist, National Transportation Safety Board
>
> Operational Factors Division (AS-30)
>
> Office: 202-314-6338
>
> Cell: 202-384-2162
>
> Email: *paul.suffern at ntsb.gov <paul.suffern at ntsb.gov>*
>
>
>
> CONFIDENTIALITY NOTICE - THIS E-MAIL TRANSMISSION MAY CONTAIN INFORMATION
> THAT IS PRIVILEGED, CONFIDENTIAL, PROPRIETARY, SUBJECT TO COPYRIGHT, AND/OR
> EXEMPT FROM DISCLOSURE UNDER APPLICABLE LAW. IT IS FOR THE USE OF INTENDED
> RECIPIENTS ONLY. If you are not an intended recipient of this message,
> please notify the original sender immediately by forwarding what you
> received and then delete all copies of the correspondence and attachments
> from your computer system. Any use, distribution, or disclosure of this
> message by unintended recipients is not authorized and may be unlawful.
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> CONFIDENTIALITY NOTICE - THIS E-MAIL TRANSMISSION MAY CONTAIN INFORMATION
> THAT IS PRIVILEGED, CONFIDENTIAL, PROPRIETARY, SUBJECT TO COPYRIGHT, AND/OR
> EXEMPT FROM DISCLOSURE UNDER APPLICABLE LAW. IT IS FOR THE USE OF INTENDED
> RECIPIENTS ONLY. If you are not an intended recipient of this message,
> please notify the original sender immediately by forwarding what you
> received and then delete all copies of the correspondence and attachments
> from your computer system. Any use, distribution, or disclosure of this
> message by unintended recipients is not authorized and may be unlawful.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200220/b2cd08a0/attachment.html>