[gradsusr] [EXTERNAL] Station time series ok on Open GrADS 2.0.1 but not on Open GrADS 2.2.1

Thu Aug 20 10:50:48 EDT 2020

Thank you for the reply, Jennifer.

Here are the first rows of ctl file:
more RETEA_PIEMO_20200810.ctl
dset    ^RETEA_PIEMO_20200810.dat
dtype   station
options sequential little_endian
stnmap  ^RETEA_PIEMO_20200810.map

ga-> q file
File 1 : Dati RETE REGIONALE
  Descriptor: RETEA_PIEMO_20200810.ctl
  Binary: RETEA_PIEMO_20200810.dat
  Type = Station Data

ga-> set lat 45
LAT set to 45 45
ga-> set lon 8
LON set to 8 8
ga-> q dims
Default file number is: 1
X is fixed     Lon = 8  X = 8
Y is fixed     Lat = 45  Y = 45
Z is fixed     Lev = 500  Z = 500
T is fixed     Time = 00Z10AUG2020  T = 1
E is fixed     Ens = 1  E = 1
ga-> d t2m(stid=90139)
Data Request Error: Invalid grid coordinates
  Varying dimension 1 decreases: 45 to 45
  Error ocurred getting variable 't2m'
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = t2m(stid=90139)

Daniele

----- Messaggio originale -----
Da: "Adams, Jennifer M. (GSFC-610.2)[ADNET SYSTEMS INC]" <jennifer.m.adams at nasa.gov>
A: "GrADS Users Forum" <gradsusr at gradsusr.org>
Inviato: Giovedì, 20 agosto 2020 15:47:44
Oggetto: Re: [gradsusr] [EXTERNAL] Station time series ok on Open GrADS 2.0.1 but not on Open GrADS 2.2.1

Hi, Daniele -- 

I am able to read your station data with my build of GraDS 2.2.1. Here is the code that is issuing that error message: 

/* Round varying dimensions 'outwards' to integral grid units. */ 

for (i=0;i<5;i++) { 

if (i==pst->idim || i==pst->jdim) { 

dmin[i] = floor(dmin[i]+0.0001); 

dmax[i] = ceil(dmax[i]-0.0001); 

if (dmax[i]<=dmin[i]) { 

gaprnt (0,"Data Request Error: Invalid grid coordinates\n"); 

snprintf(pout,1255," Varying dimension %i decreases: %g to %g\n",i,dmin[i],dmax[i]); 

gaprnt (0,pout); 

snprintf(pout,1255," Error ocurred getting variable '%s'\n",vnam); 

gaprnt (0,pout); 

return (NULL); 

} 

} 

} 

This should only appear for gridded variables, but you are using station data. What is the output from ‘q file’? Does it say “Type = Station Data” ? Is there some weird white space in your ‘dtype station’ entry in the descriptor file? 

Also, what is the output from ‘q dims’ ? Does it show X and Y as fixed? 

--Jennifer 

-- 

Jennifer Miletta Adams 

ADNET Systems, Inc. 

NASA/GSFC, Code 610.2 

Building 32, Room S159 

(301) 614-6070 

On 8/20/20, 5:09 AM, "gradsusr on behalf of Daniele Gandini" <gradsusr-bounces at gradsusr.org on behalf of daniele.gandini at arpa.piemonte.it> wrote: 

Dear GrADS Users, 

sorry for the long mail but the question is important for us. 

We have been using Open GrADS 2.0.1 on a Red Hat Enterprise Linux Server release 6.5 

since many years with excellent results. 

Now we are going to update both Linux server with Fedora 31 operating system 

and Open GrADS version 2.2.1. 

With Open GrADS 2.2.1 we aren't able to plot station time series 

of the attached files written in Fortran 90 language. They are located in Piedmont, north-western Italy. 

============================================================================================================== 

grads -cl "open RETEA_PIEMO_20200810.ctl" 

Welcome to the OpenGrADS Bundle Distribution 

-------------------------------------------- 

For additional information enter "grads --manual". 

Starting "/usr/local/bin/Linux/Versions/2.2.1.oga.1/x86_64/grads -cl 'open RETEA_PIEMO_20200810.ctl' " ... 

ga-> set t 1 last 

Time values set: 2020:8:10:0 2020:8:10:23 

ga-> set lat 45 

LAT set to 45 45 

ga-> set lon 8 

LON set to 8 8 

ga-> d t2m(stid=90139) 

Data Request Error: Invalid grid coordinates 

Varying dimension 1 decreases: 45 to 45 

Error ocurred getting variable 't2m' 

Error ocurred at column 1 

DISPLAY error: Invalid expression 

Expression = t2m(stid=90139) 

===================================================================================================================== 

Instead with Open GrADS 2.0.1: 

grads -cl "open RETEA_PIEMO_20200810.ctl" 

Welcome to the OpenGrADS Bundle Distribution 

-------------------------------------------- 

For additional information enter "grads -h". 

Starting "/opt/grads-2.0.1.oga.1/Contents/Linux/Versions/2.0.1.oga.1/x86_64/grads -cl 'open RETEA_PIEMO_20200810.ctl' " ... 

ga-> set t 1 last 

Time values set: 2020:8:10:0 2020:8:10:23 

ga-> set lat 45 

LAT set to 45 45 

ga-> set lon 8 

LON set to 8 8 

ga-> d t2m(stid=90139) 

ga-> 

all is ok. 

Do you know which can be the problem? 

In the following rows you can see how we installed Open GrADS 2.2.1 

** Install & Configure OpenGrads *** 

[root at xmeteo6 ~]# dnf install bc gcc-gfortran gdal gdal-python gdal-devel ImageMagick perl-Env 

[root at xmeteo6 ~]# wget https://urldefense.proofpoint.com/v2/url?u=https-3A__sourceforge.net_projects_opengrads_files_grads2_2.2.1.oga.1_Linux-2520-252864-2520Bits-2529_opengrads-2D2.2.1.oga.1-2Dbundle-2Dx86-5F64-2Dpc-2Dlinux-2Dgnu-2Dglibc-5F2.17.tar.gz&d=DwIFaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=M2mFZsxo0ZoztoK2Lyoypu1kKkBVKyEoM9CXGWO42zk&m=QH8IXCI1bvP5NhkYHNRGHX3oSBywJY33joPLKgsPT_w&s=R6i5oyzqUAy827sOdyGGE7_0ODCVRb8s-qQRgq7i8io&e= -P /opt 

[root at xmeteo6 ~]# cd /opt 

[root at xmeteo6 opt]# tar xvzf opengrads-2.2.1.oga.1-bundle-x86_64-pc-linux-gnu-glibc_2.17.tar.gz 

[root at xmeteo6 opt]# cp -rf opengrads-2.2.1.oga.1/Contents/* /usr/local/bin/ 

[root at xmeteo6 opt]# rm –rvf opengrads-2.2.1.oga.1/ 

[root at xmeteo6 opt]# rm -rvf opengrads-2.2.1.oga.1-bundle-x86_64-pc-linux-gnu-glibc_2.17.tar.gz 

[root at xmeteo6 opt]# grads -p ==> to test the installation 

Test Installation 

root at xmeteo6 opt]# which perl; echo $? 

/usr/bin/perl 

0 

[root at xmeteo6 opt]# which gcc; echo $? 

/usr/bin/gcc 

0 

[root at xmeteo6 opt]# which gfortran; echo $? 

/usr/bin/gfortran 

0 

[root at xmeteo6 opt]# which bc; echo $? 

/usr/bin/bc 

0 

[root at xmeteo6 opt]# which curl; echo $? 

/usr/bin/curl 

0 

[root at xmeteo6 opt]# which wget; echo $? 

/usr/bin/wget 

0 

[root at xmeteo6 opt]# which convert; echo $? 

/usr/bin/convert 

0 

[root at xmeteo6 opt]# which gdal_rasterize; echo $? 

/usr/bin/gdal_rasterize 

0 

[root at xmeteo6 opt]# which gedit; echo $? 

/usr/bin/gedit 

0 

[root at xmeteo6 opt]# which grads; echo $? 

/usr/local/bin/grads 

0 

For each of these command lines, check status number display on the terminal: 

if you see 0, then the package is installed; but if you have 1, the package is not installed. 

Thank you in advance. 

Daniele Gandini 

Servizio Meteorologico ARPA Piemonte 

Via Pio VII 9 - 10135 Turin (Italy) 

https://urldefense.proofpoint.com/v2/url?u=http-3A__www.arpa.piemonte.it_&d=DwIFaQ&c=ApwzowJNAKKw3xye91w7BE1XMRKi2LN9kiMk5Csz9Zk&r=M2mFZsxo0ZoztoK2Lyoypu1kKkBVKyEoM9CXGWO42zk&m=QH8IXCI1bvP5NhkYHNRGHX3oSBywJY33joPLKgsPT_w&s=ZNLHyGD51KpDvhGoJg0ZbKndndhwLtssQ4SXL7l2gr4&e= 

Daniele.Gandini at arpa.piemonte.it 

-- 
Questo messaggio è stato analizzato con Libraesva ESG ed è risultato non infetto. 

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
 https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&e=457fc64d&h=99e2e735&f=y&p=n