[ncl-talk] Question regarding the transfer of integer values from shell to NCL

[Prashanth Bhalachandran]

Marston, 
The problem is I have several arrays where time is one of the variables as you will see in the attached code. 

For example, 
 u10  = a->u10 is [76 times] x [801 lats] x [810 lon]
  u   = a->u.  Is [76 times] x [11 levels] [801 lats] x [810 lon]

Using the input argument from the shell script, I am reading these variables at particular times. That is, the text file contains the name of the file and the particular time. 

Therefore, when I use it like this : 
speed(:,:)   = wind_speed(u10(time,0,:,:),v10(time,0,:,:)) * toknots 

Or any such line, there will be an error primarily because the time input from the shell script is a string. As you can see from the printVarSummary of the received time from the shell script. 

Variable: time
Type: string
Total Size: 8 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes:	[1]
Coordinates: 
Number Of Attributes: 1
  _FillValue :	-2147483647
-2147483647



> On Oct 19, 2017, at 10:42 AM, Marston Johnston <shejo284 at gmail.com> wrote:
> 
> It would help if you did a printVarSummary() on the problem array.
>  
>  
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> Marston S. Ward, PhD
> Department of Earth Sciences
> University of Gothenburg, Sweden
> Email: marston.johnston at gu.se <mailto:marston.johnston at gu.se>
> SkypeID: marston.johnston 
> Phone: +46-31-7864901 
> Only the fruitful thing is true!
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>  
>  
> From: ncl-talk <ncl-talk-bounces at ucar.edu <mailto:ncl-talk-bounces at ucar.edu>> on behalf of Prashanth Bhalachandran <prashanth.bhalachandran at gmail.com <mailto:prashanth.bhalachandran at gmail.com>>
> Date: Thursday, 19 October 2017 at 19:31
> To: <ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>>
> Subject: [ncl-talk] Question regarding the transfer of integer values from shell to NCL
>  
> Dear NCL team, 
> Greetings. 
>  
> I have a shell script that does the following : Read from a text file which has a list of dates (yyyyddmm) format and times. For example, the first column in the text represents the date and the second the time. 
> 20121030        13     73.34     51.77    -21.57
> 
> 
> My shell script reads in the first and second column and stores those values and passes off as arguments to my NCL script. My NCL script reads the first argument as the filename and the second as time. I use this time as an index in several of my arrays in my NCL script since my variable dimensions are time x lev x lat x lon. 
> The problem here is that shell reads these values as strings and I am unable to get it as an integer, which is compulsory if I want to use it as an array index. Please see the below error. 
> Can one of you please guide me as to how I pass this value as an integer? 
>  
> Many thanks, 
> Prashanth 
>  
> Error: 
> fatal:Illegal subscript. Subscripts must be integer when not using coordinate indexing
> 
> 
> My code (I’m attaching the entire code in case you want to have a look, but I am pasting the necessary portions here) : 
>  
> #!bin/bash
> file="/scratch/conte/s/sbhalach/DATA/RI_compute/CatRW2030.txt"
> while read f1 f2 f3 f4;
> do
>  
>   export dirname=$(echo $f1)00post
>   echo $dirname
>   cd $dirname
>   export filename=${f1}.nc
>   export time=$f2
>  
> #*************************************************************************
>   cat > varcalc.ncl << EOF 
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;  
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;  
>  
>  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>  load "/scratch/lustreD/s/sbhalach/DATA/func_center.ncl"
>  load "/scratch/lustreD/s/sbhalach/DATA/func_rtheta.ncl"
>  
>  a    = addfile(filename,"r")
>  u10  = a->u10
>  v10  = a->v10
>  slp  = a->slp
>  u    = a->u
>  v    = a->v
> allvars   = new((/40/),float)
> 
> 
> slpsub       = slp(time,0,:,:) 
> ; Note that the time here is taken as string and hence the error. I obviously tried the toint() and strongpoint() functions but then the function only returns missing values. 
> .
> .
> . 
> ;;;;;;;;;;;;;;;;;;;;;; Write to a netcdf file ;;;;;;;;;;;;;;;;;;
>  
>   system("rm -f $f1_$f2_allvars.nc")
>   ncdf = addfile("$f1_$f2_allvars.nc","c")
>   ncdf->allvars = allvars
>  
> EOF
>  
> #*************************************************************************
>  
>  ncl -n 'filename = "$filename"' 'time = "$time"' varcalc.ncl
>  mv $f1_$f2_allvars.nc ../
>  cd ../
>  done <"$file”
> 
> 
> 
> 
> 
> 
> _______________________________________________ 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 <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171019/5f7cb9f0/attachment.html>