delorie.com/archives/browse.cgi   search  
Mail Archives: cygwin/2015/11/23/15:39:58

X-Recipient: archive-cygwin AT delorie DOT com
DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id
:list-unsubscribe:list-subscribe:list-archive:list-post
:list-help:sender:mime-version:in-reply-to:references:date
:message-id:subject:from:to:content-type; q=dns; s=default; b=nS
qONZE7s+qC6zIfI9bxiCWDybOWMMcDzyD7NTu/6O7wl6HG7BABolCFR3xBx7EJn2
oEGIDNO2BAT6dmDa6YgT9EM5hAjftn/GTbSY/FIxCb1A94ojoSW1GAPKu7LH6l0y
OUtFRnflfdmPCfTtnb5yDbGE04AMT03rPuAUnr/G8=
DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id
:list-unsubscribe:list-subscribe:list-archive:list-post
:list-help:sender:mime-version:in-reply-to:references:date
:message-id:subject:from:to:content-type; s=default; bh=rk5mFwJ8
FZWnKgtxANCFo7e76u4=; b=uZ6nK0Rwig2bY4faynajw/dqCLM+ZaO5wR4CPeAr
/xshA7GQINFnaxotSSQ3Wv2+yjZgI8IgyvIAVwSAMwPQ81Lb9NnbXaUJ7XYQUWVl
GmD+pA6F2Vyk3+go61/jpubKVRSaKfRQS9gNbcC8VfOwcjDwWrzk6nu2p/tROULU
BhI=
Mailing-List: contact cygwin-help AT cygwin DOT com; run by ezmlm
List-Id: <cygwin.cygwin.com>
List-Subscribe: <mailto:cygwin-subscribe AT cygwin DOT com>
List-Archive: <http://sourceware.org/ml/cygwin/>
List-Post: <mailto:cygwin AT cygwin DOT com>
List-Help: <mailto:cygwin-help AT cygwin DOT com>, <http://sourceware.org/ml/#faqs>
Sender: cygwin-owner AT cygwin DOT com
Mail-Followup-To: cygwin AT cygwin DOT com
Delivered-To: mailing list cygwin AT cygwin DOT com
Authentication-Results: sourceware.org; auth=none
X-Virus-Found: No
X-Spam-SWARE-Status: No, score=-1.4 required=5.0 tests=AWL,BAYES_05,FREEMAIL_ENVFROM_END_DIGIT,FREEMAIL_FROM,RCVD_IN_DNSWL_LOW,SPF_PASS autolearn=ham version=3.3.2
X-HELO: mail-wm0-f48.google.com
MIME-Version: 1.0
X-Received: by 10.194.114.105 with SMTP id jf9mr31989580wjb.102.1448311180539; Mon, 23 Nov 2015 12:39:40 -0800 (PST)
In-Reply-To: <CAE3taFC=tjki1FTfbjyt+KnEGsWmrj63FMeBjgRswRwRSKdAWQ@mail.gmail.com>
References: <CAE3taFC=tjki1FTfbjyt+KnEGsWmrj63FMeBjgRswRwRSKdAWQ AT mail DOT gmail DOT com>
Date: Mon, 23 Nov 2015 15:39:40 -0500
Message-ID: <CAD8GWsvXCC0GB8mxozLrT4aWpfoRhe7mEu9whSi2o6kp59q61w@mail.gmail.com>
Subject: Re: Awk not ouputting results via echo
From: Lee <ler762 AT gmail DOT com>
To: cygwin AT cygwin DOT com
X-IsSubscribed: yes

On 11/23/15, Lester Anderson <arctica1963 AT gmail DOT com> wrote:
> Hello,
>
> Having sorted how to access a separate drive location, I have expanded
> the script to include an awk section which processes new coordinate
> limits which should pass to img2grd to extract the data grids.
> However, the awk scripting does not seem to generate any data - is
> there something I am missing? Modified from a published script

As has already been pointed out, the read isn't going to get envars
set the way you want, so call another script with the proper
variables.  And AWK is going to want to read input, so you need to put
everything inside BEGIN{ "your stuff here" ; exit }
And I use -v and single quotes with awk ,so how about something like

awk  -v lon_min="$lon_min" -v lon_max="$lon_max" -v lat_min="$lat_min" \
     -v lat_max="$lat_max" -v R_d="$R_d" -v R_i="$R_i" \
'BEGIN{R_t=6370;
  ... etc ...
      cmd=sprintf("%s %f %f %f %f %f %f %f %f %f %f\n",
                  "echo ", on_e_min,lon_e_max,lat_e_min,lat_e_max,lon_i_min,
                           lon_i_max,lat_i_min,lat_i_max,lat_av*180./pi,lon_av)
      system(cmd)
  exit
}'


.. you need to write your own secondary script & not use "echo"

Regards,
Lee


>
> #!/bin/bash
> ruta_elev="Q:/geophys/Potential-field datasets/Topography/topo_17.1.img"
> ruta_grav="Q:/geophys/Potential-field datasets/Gravity/grav.img.23.1"
> lon_min=-12
> lon_max=0
> lat_min=28
> lat_max=39
> R_d=167
> R_i=20
>
> echo $lon_min $lon_max $lat_min $lat_max $R_d $R_i |
> awk "{R_t=6370;
>       pi=3.14159;
>       lat_av=(($lat_max+$lat_min)/2)*(pi/180.);
>       lon_av=(($lon_min+$lon_max)/2);
>       d_lat_e=($R_d/R_t)*180./pi;
>       d_lon_e=($R_d/(R_t*cos(lat_av)))*180./pi;
>       d_lat_i=($R_i/R_t)*180./pi;
>       d_lon_i=($R_i/(R_t*cos(lat_av)))*180./pi;
>       lon_e_min=$lon_min-d_lon_e;
>       lon_e_max=$lon_max+d_lon_e;
>       lat_e_min=$lat_min-d_lat_e;
>       lat_e_max=$lat_max+d_lat_e;
>       lon_i_min=$lon_min-d_lon_i;
>       lon_i_max=$lon_max+d_lon_i;
>       lat_i_min=$lat_min-d_lat_i;
>       lat_i_max=$lat_max+d_lat_i;
>       print
> lon_e_min,lon_e_max,lat_e_min,lat_e_max,lon_i_min,lon_i_max,lat_i_min,lat_i_max,lat_av*180./pi,lon_av;
> }" | read lon_e_min lon_e_max lat_e_min lat_e_max lon_i_min lon_i_max
> lat_i_min lat_i_max lat_av lon_av
> proj='m'$lon_av'/'$lat_av'/16'
> echo $lon_e_min $lon_e_max $lat_e_min $lat_e_max $lon_i_min $lon_i_max
> $lat_i_min $lat_i_max $lat_av $lon_av
> #
> img2grd "$ruta_elev"  -V -T1
> -R$lon_e_min/$lon_e_max/$lat_e_min/$lat_e_max -S1 -Gelev.grd -I1m -D
> img2grd "$ruta_grav" -T1 -V
> -R$lon_e_min/$lon_e_max/$lat_e_min/$lat_e_max -S0.1 -GFA.grd -I1m -D
> #
>
> Any pointers on this, thanks
>
> Lester

--
Problem reports:       http://cygwin.com/problems.html
FAQ:                   http://cygwin.com/faq/
Documentation:         http://cygwin.com/docs.html
Unsubscribe info:      http://cygwin.com/ml/#unsubscribe-simple

- Raw text -


  webmaster     delorie software   privacy  
  Copyright © 2019   by DJ Delorie     Updated Jul 2019