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

Suffern Paul paul.suffern at ntsb.gov
Thu Feb 20 11:03:48 MST 2020 

Hi Rick,

The version is was using was from back in 2014 or so… yes, a very long time ago… I just don’t/didn’t have time to do updates until was forced too this past month, but since everything was still running smoothly didn’t need too. The hiccup in this script is the only ncl script I’m running that has issues from the port over to NCL 6.6, etc…

Thanks again,

Paul


From: Rick Brownrigg <brownrig at ucar.edu>
Sent: Thursday, February 20, 2020 12:54 PM
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,

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<mailto: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<mailto:brownrig at ucar.edu>>
Sent: Thursday, February 20, 2020 11:37 AM
To: Suffern Paul <paul.suffern at ntsb.gov<mailto:paul.suffern at ntsb.gov>>
Cc: ncl-talk at ucar.edu<mailto: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<mailto: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<mailto: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<mailto: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.
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/8b81dece/attachment.html>

More information about the ncl-talk mailing list