<div dir="ltr"><div dir="ltr">Hi Jennifer,<br></div><div dir="ltr"><br></div><div>Thanks for confirming the issue.</div><div><br></div><div>Will the update require a new build of GrADS or will it be possible to overwrite one of the  libraries / scripts in an existing install of GrADS and Open GrADS (I use both)?</div><div><br></div><div>As an interim I've written a function using Zeller's congruence to get the day of the week when ??? is returned by q time.</div><div><br></div><div>I'm still testing and it is quite clunky so I won't post here. However, if anyone gets stuck I'm happy to try and help.</div><div dir="ltr"><br></div><div>Thanks</div><div><br></div><div>Brian</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 21 Sep 2020 at 14:34, <<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Send gradsusr mailing list submissions to<br>
        <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:gradsusr-owner@gradsusr.org" target="_blank">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. display values in a script (Frank Colby)<br>
   2. Re: display values in a script (James T. Potemra)<br>
   3. Re: [EXTERNAL]  Q time in 2021 returns ??? for the day<br>
      (Adams, Jennifer M. (GSFC-610.2)[ADNET SYSTEMS INC])<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sun, 20 Sep 2020 13:10:03 -0400<br>
From: Frank Colby <<a href="mailto:Frank_Colby@uml.edu" target="_blank">Frank_Colby@uml.edu</a>><br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
Subject: [gradsusr] display values in a script<br>
Message-ID: <<a href="mailto:aaa52185-47c4-95f3-f918-e1d0b9e64c64@uml.edu" target="_blank">aaa52185-47c4-95f3-f918-e1d0b9e64c64@uml.edu</a>><br>
Content-Type: text/plain; charset=utf-8; format=flowed<br>
<br>
Dear GradUsers,<br>
<br>
I can open a file, set latitude and longitude and time, and display a <br>
variable.? The value of the variable displays in the command window.<br>
<br>
If I try this in a script, it displays nothing.<br>
<br>
Does anyone know why this is?? The script is below.<br>
<br>
Thanks,<br>
Frank Colby<br>
UMass Lowell<br>
<br>
*********************Script<br>
'set mpdset hires'<br>
filen = 20100219<br>
 ? latt = 42.361<br>
 ? lont = -071.010<br>
* All times in one file 'set t ' t1<br>
say 'ERA5-'filen'-<a href="http://s2.nc" rel="noreferrer" target="_blank">s2.nc</a>'<br>
'sdfopen ERA5-'filen'-<a href="http://s2.nc" rel="noreferrer" target="_blank">s2.nc</a>'<br>
*<br>
'set lat 'latt<br>
'set lon 'lont<br>
*<br>
t1=1<br>
t2=24<br>
while (t1<=t2)<br>
'set t ' t1<br>
'd 'sp<br>
*<br>
say t1<br>
say '<enter>=next, q=quit, b=back, t<enter>n<enter>=n-th step'<br>
say 'i to create an image'<br>
pull ans<br>
*<br>
if (ans ="q"); break; endif;<br>
if (ans = "b")<br>
 ? t1 = t1 -2<br>
endif<br>
if (ans = "t"); say 'Time step =';pull num;<br>
 ? if (num > t2); say 'Timestep too high!'; num=t2-1; endif;<br>
 ? t1 = num - 1<br>
endif<br>
*<br>
t1 = t1 + 1<br>
'c'<br>
endwhile<br>
'reinit'<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sun, 20 Sep 2020 16:06:40 -1000<br>
From: "James T. Potemra" <<a href="mailto:jimp@hawaii.edu" target="_blank">jimp@hawaii.edu</a>><br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>>, Frank Colby<br>
        <<a href="mailto:Frank_Colby@uml.edu" target="_blank">Frank_Colby@uml.edu</a>><br>
Subject: Re: [gradsusr] display values in a script<br>
Message-ID: <<a href="mailto:66df65db-03f0-f173-fb80-04efd0161cf7@hawaii.edu" target="_blank">66df65db-03f0-f173-fb80-04efd0161cf7@hawaii.edu</a>><br>
Content-Type: text/plain; charset=utf-8; format=flowed<br>
<br>
Hi Frank,<br>
<br>
In a script you'll have to explicitly print the output, so right after<br>
<br>
'd 'sp<br>
<br>
add<br>
<br>
say result<br>
<br>
Jim<br>
<br>
On 9/20/20 7:10 AM, Frank Colby wrote:<br>
> Dear GradUsers,<br>
><br>
> I can open a file, set latitude and longitude and time, and display a <br>
> variable.? The value of the variable displays in the command window.<br>
><br>
> If I try this in a script, it displays nothing.<br>
><br>
> Does anyone know why this is?? The script is below.<br>
><br>
> Thanks,<br>
> Frank Colby<br>
> UMass Lowell<br>
><br>
> *********************Script<br>
> 'set mpdset hires'<br>
> filen = 20100219<br>
> ? latt = 42.361<br>
> ? lont = -071.010<br>
> * All times in one file 'set t ' t1<br>
> say 'ERA5-'filen'-<a href="http://s2.nc" rel="noreferrer" target="_blank">s2.nc</a>'<br>
> 'sdfopen ERA5-'filen'-<a href="http://s2.nc" rel="noreferrer" target="_blank">s2.nc</a>'<br>
> *<br>
> 'set lat 'latt<br>
> 'set lon 'lont<br>
> *<br>
> t1=1<br>
> t2=24<br>
> while (t1<=t2)<br>
> 'set t ' t1<br>
> 'd 'sp<br>
> *<br>
> say t1<br>
> say '<enter>=next, q=quit, b=back, t<enter>n<enter>=n-th step'<br>
> say 'i to create an image'<br>
> pull ans<br>
> *<br>
> if (ans ="q"); break; endif;<br>
> if (ans = "b")<br>
> ? t1 = t1 -2<br>
> endif<br>
> if (ans = "t"); say 'Time step =';pull num;<br>
> ? if (num > t2); say 'Timestep too high!'; num=t2-1; endif;<br>
> ? t1 = num - 1<br>
> endif<br>
> *<br>
> t1 = t1 + 1<br>
> 'c'<br>
> endwhile<br>
> 'reinit'<br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 21 Sep 2020 13:29:51 +0000<br>
From: "Adams, Jennifer M. (GSFC-610.2)[ADNET SYSTEMS INC]"<br>
        <<a href="mailto:jennifer.m.adams@nasa.gov" target="_blank">jennifer.m.adams@nasa.gov</a>><br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
Subject: Re: [gradsusr] [EXTERNAL]  Q time in 2021 returns ??? for the<br>
        day<br>
Message-ID: <<a href="mailto:AEA487EE-C1BE-4F31-9A61-FEF262A799E2@nasa.gov" target="_blank">AEA487EE-C1BE-4F31-9A61-FEF262A799E2@nasa.gov</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi, Brian ?<br>
Good catch! That routine to calculate the day of the week has hard-coded years in it:<br>
<br>
/* Return day of week for date/time  0=sunday, 6=saturday */<br>
gaint dayweek (struct dt *dtime) {<br>
struct dt anch;<br>
gaint i,j;<br>
  if (dtime->yr<1950 || dtime->yr>2020) return(7);<br>
  anch.yr = 1950;<br>
  anch.mo = 1;<br>
  anch.dy = 1;<br>
  <a href="http://anch.hr" rel="noreferrer" target="_blank">anch.hr</a> = 0;<br>
  <a href="http://anch.mn" rel="noreferrer" target="_blank">anch.mn</a> = 0;<br>
  i = timdif(&anch,dtime,0);<br>
  i = i/1440;<br>
  j = i/7;<br>
  i = i - j*7;<br>
  return(i);<br>
}<br>
<br>
When year is outside this arbitrary window, it returns ?7?, which is the array index for this set of strings for the day of the week:<br>
<br>
char *dweek[8] = {"Sun","Mon","Tue","Wed","Thu","Fri","Sat","???"};<br>
<br>
I will have to update this in the source code.<br>
--Jennifer<br>
<br>
<br>
<br>
--<br>
Jennifer Miletta Adams<br>
ADNET Systems, Inc.<br>
NASA/GSFC, Code 610.2<br>
Building 32, Room S159<br>
(301) 614-6070<br>
<br>
<br>
From: gradsusr <<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>> on behalf of Brian Gaze <<a href="mailto:brianw.gaze@googlemail.com" target="_blank">brianw.gaze@googlemail.com</a>><br>
Reply-To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
Date: Sunday, September 20, 2020 at 7:12 AM<br>
To: "<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>" <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
Subject: [EXTERNAL] [gradsusr] Q time in 2021 returns ??? for the day<br>
<br>
Hi,<br>
<br>
When doing a q time on a ctl file which contains the following line a string containing the DAY (e.g. SUN) is returned as expected:<br>
<br>
tdef 1 linear 18Z20dec2020 1mo<br>
<br>
However, when doing a q time on ctl files which are for a date in 2021 (as below) the DAY returned is ???<br>
<br>
tdef 1 linear 06Z04jan2021 1mo<br>
<br>
The month, year and time are returned correctly in both cases. It is only a problem with the DAY part of the string.<br>
<br>
I have experienced this behaviour with GrADS 2.2.0 on Ubuntu and Open GrADS v2.1.a2.oga.1 on Windows Server 2016.<br>
<br>
Does anyone know the reason for this?<br>
<br>
PS: If necessary I can post both ctl files in full.<br>
<br>
Thanks<br>
<br>
Brian<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://gradsusr.org/pipermail/gradsusr/attachments/20200921/125bf083/attachment.html" rel="noreferrer" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20200921/125bf083/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
------------------------------<br>
<br>
End of gradsusr Digest, Vol 127, Issue 15<br>
*****************************************<br>
</blockquote></div></div>