[ferret_users] Nested loops and variable filename output

To: <ferret_users@xxxxxxxx>

Subject: [ferret_users] Nested loops and variable filename output

From: <Benjamin.Ng@xxxxxxxx>

Date: Wed, 22 Nov 2017 01:31:10 +0000

Accept-language: en-GB, en-AU, en-US

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=W0qLB3sl; arc=pass (i=2 dkim=pass dkdomain=csiro.au); spf=fail (google.com: domain of prvs=492ea0156=benjamin.ng@xxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=prvs=492ea0156=benjamin.ng@xxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=csiro.au

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=W0qLB3sl; arc=pass (i=2 dkim=pass dkdomain=csiro.au); spf=fail (google.com: domain of prvs=492ea0156=benjamin.ng@xxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=prvs=492ea0156=benjamin.ng@xxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=csiro.au

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=W0qLB3sl; arc=pass (i=1 spf=pass spfdomain=csiro.au dkim=pass dkdomain=csiro.au dmarc=pass fromdomain=csiro.au); spf=fail (google.com: domain of prvs=492ea0156=benjamin.ng@xxxxxxxx does not designate 65.127.216.230 as permitted sender) smtp.mailfrom=prvs=492ea0156=Benjamin.Ng@xxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@csiro.au header.s=dkim header.b=W0qLB3sl; spf=pass (google.com: domain of prvs=492ea0156=benjamin.ng@xxxxxxxx designates 2405:b000:e00:257::7:41 as permitted sender) smtp.mailfrom=prvs=492ea0156=Benjamin.Ng@xxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=csiro.au

Arc-message-signature: i=4; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=list-archive:list-help:list-post:list-id:mailing-list:precedence :mime-version:content-language:accept-language:date:thread-index :thread-topic:subject:to:from:arc-authentication-results :arc-message-signature:arc-authentication-results :arc-message-signature:message-id:arc-authentication-results :arc-message-signature:dkim-signature:arc-authentication-results; bh=zVjqPYOHpqpuIoBfR1rbSnwGsErTnwwbbDYFxqjGTAg=; b=VFXee5lKZV1b64daOMn3jWxIwOJht9YiMVHnP7yHJL7AFyZfBEBXtojpV2a7qyEaoW GK++dLsXxRnMWWXxy4nb6uBl76lCnFsJ0YQsIhk/HptdIkykTi0AqGzWWCp4Ke3uyoKC VsokC/X+GH7sBZWaNiTy1p96Rudg2UMXNnsq7wRmE/4Z5YRFZzIIGSnq6JGvfzBoKsIF HnN/I8uKmnxByI0pV/T7vdxrVnYAcZqpHolU/JqF0O+BL2kDkS757LeXp3fjo8W9/bzz ejdKj4IHKJ5MlDuUPPg1j/Cpjb/mqSrRoXrj595AbJtuAjbTQAnQ3BqQ+Z4Aqb44mVf9 Juqg==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:date:thread-index :thread-topic:subject:to:from:dkim-signature :arc-authentication-results:arc-message-signature :arc-authentication-results:arc-message-signature:message-id :arc-authentication-results; bh=zVjqPYOHpqpuIoBfR1rbSnwGsErTnwwbbDYFxqjGTAg=; b=FS1b0NVuPJQsNThyZQVbQFmUGE+tvrgGBFoDaqqDWXG3WyiyHJ7QwWYU31+noTca2y G3xToKZ+VxrUQLCDdQuHYT/fi8T3U2cK1hW75xi3PS/7FrIeX5+cpPeQy66u2oAwLiOX fXkwhJgUQSkM7RmEQQ15kEkJB3AzqiyXOfdBPypbwXNaqehwdIJeLvavV6W1yKOx6LCr +bfx7JOzbqAXtL5e98XPY5aY8OtaMYqa9tLRbkhwsMec1eo52yL7KbLUotB+u83Sbydm L4eM48fx7NtQpLH1vLII7/8lu7pNTv5NYjrafKZUO0TckwcQuP4t2v5i3BhELvLVnNcL p+9w==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:message-id:date :thread-index:thread-topic:subject:to:from:dkim-signature :arc-authentication-results:arc-message-signature :arc-authentication-results; bh=zVjqPYOHpqpuIoBfR1rbSnwGsErTnwwbbDYFxqjGTAg=; b=rKD5EKLJQVD0McHO52M6wLLqIQXWoK9ovboZXD+tVcJrOQMbY3H4VSaYSFaVXvCzyH rBGz8uL5bfqIREkXzDk67eCK64WnUTX6z424/mC3TJqA0z4ovQt1qxAqWDAxm3orXhVN S7PERIeGYgPYFB7PeFjsiG3QNYFrIPreGazyFfSz9A/WxGT3sQ6h7DUrzglc2THTWWh3 i0NJRM5mvhT4thcwXCKnf+a9OFuxB07L6Z9TEJBQ2HA8agKPgRijktlfKcinlUZvE1J8 nVoPr+gaVFT+RmfUMJSqEEOw+4elUEFIbN/q/Ecb9jHPoR7C1Z7fKVSjO3HuZnfYQoSp K7cw==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:message-id:date :thread-index:thread-topic:subject:to:from:dkim-signature :arc-authentication-results; bh=zVjqPYOHpqpuIoBfR1rbSnwGsErTnwwbbDYFxqjGTAg=; b=MJgiU8vVPqBTl5uWvEKmrnSmV9gPHnRmKSdgmTBw8eakhPHKK3pKSsKms3jxctzJ5T 0Bj+vW1i1Qw/eZJrXbeJUMok2lf8qdW60vUHVgH5Qhb+ZsX/+82t9CJf0m2BjKqoB935 kGS9iL6vKdZynFcbBPUYboxoCNpRGOX/jTpMDPCY0SwagIbNFnB6Kueyph7x4GPT5qn4 ugsWmijK4hJIi7xPNWsvzggExqISk+uIVAsOM7le7JWjRUR+w9ohF5neKFUhnX3ApNET ABFW1isYzozQQkydMWawNEDpXqTCEfi7U6RXxSPB5WQnnUR30a5//LRkDOMFdO8ZEIuE OmTA==

Arc-seal: i=4; a=rsa-sha256; t=1511314307; cv=pass; d=google.com; s=arc-20160816; b=jPcruMooMk5f5zIz4B0SicCe2hoIPNcmxEibenLvXAroXCyfTKWwuq+ZvvGc3FL1PG BwGnk0XMdZAF7Czu7PoUA9iay86/zx3gz8cYkXz3E7cdGfIfFmFInonzE3PMMgFPoQBf wYEz2RbIS5aXFmmoiCcRYnVZdYWBReqqUuu3APz0NSWMoA5RFefolq7gUMQ8op/nh05+ RQ22w/Z+Sl4qZU6VsZZq7G2kZ1ZK7OJxdMXAjj1hM+pgD+NK9wyJ6h50bH/ZArpGzd4N H2jAgJuF7VVbBOy8I1U42IdHUiDf969/DZrJn77C9GcfGYo1M5//eA2R/kl6G70Ql5wR R2fA==

Arc-seal: i=3; a=rsa-sha256; t=1511314306; cv=pass; d=google.com; s=arc-20160816; b=ZvxVACsE/ezlIioiKRi6xS2thuW+9aV/r6VProbg1AMdgji0KNC+Rc4N9dsLz+oNTq MbKbWwaA0MVOsm11rVQmlPXJRsDocwbmU7sP0GA9tOYHHEUVoI3Y711Ytf4A0oGTFAao hCrHrLPYYBAnn7T1/J3s+Q82gWC/TdafNKfrgbxs40WGEv2tZSU3SFBOLV09ElX9Dv8L eaq6DGQwTAbIHcmOu4ezc/fOQvMvQQlIt2F3FbYeSPQWpjHB6IcDf6A0qCM6aG/aF839 0MY8l2BFkdhfR3GS68wAW7CALADjUmzPMs9xLBEfnspnbiL77fQ+vgPS24FiTMlmbaOK lVNA==

Arc-seal: i=2; a=rsa-sha256; t=1511314304; cv=pass; d=google.com; s=arc-20160816; b=IdqCsLnghhwEBiRlICNIuJ3I/Xr/262pXfEALEuAABZL+ePpdKEaxrymfjphjg5+44 KOohBimuDzWrGNufgavrVOXgTBIij01Ze2KKlNpYW1Ci9/6G9W/UjRZOwXSIBcsSWTtP 1E/cb6y7h2jXHUsVwez9MF4V4EkRPW2KTh7NcbjVJL92s1CM+pApL8KS7W8jk66kO0JY jMO18rex3fIjXKou5Sj+YQnrtH0bXhf19vUUswmut0FLIfWAFWv4Tr/0/8G7KdKXRAWQ EH9ITuLgcHHhJyZb+lZEutqZA43nvePNB3DIzJkpaL5myj6c4FFuvT6k9IuxQew8L+Kw UJFg==

Arc-seal: i=1; a=rsa-sha256; t=1511314299; cv=none; d=google.com; s=arc-20160816; b=D5b+hrjIsDcCEHTzudLJzgdpytRda89axUkw1rgc85M/I50/ulzU6hhIS2Nxwk5Xos xjrO56oYHJTbE4ltHEq4893wf5D62zQmjygAjRYzT6SFMb31l2j4ZrWLfdFoYX4ttTRQ WInrkyH6dYKh7H7TJGcUidWVwiCac3Pqp2L1Gz9Y1U9+KkbP4aWE8xzao3I2fJ2LS82w qqj/VDnUIwA9edalSuFo+OUTiYufaaemEyBY1W06KDxXlaY1w5RRu9an1F9FUgWiTTWH laAvEHjeDTHWb9ZaUI9TzxjNGwhWY4LAgGxmnwCyzEQUb4B+ZNHYvF59t6tAzyvfqNzl 8z2A==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=message-id:from:to:subject:thread-topic:thread-index:date :accept-language:content-language:mime-version:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=zVjqPYOHpqpuIoBfR1rbSnwGsErTnwwbbDYFxqjGTAg=; b=abgscLikcLZwHXgQYwTUYboKu9HXZ17p7jWA0dV2cirQnLA6DspmV9ry9hxYKqori/ 4l/8e2/mE2KMadXy5EFcTmOGmk1SyRerhCWTi1CBkDckKKf9SUupSfslU8WoPjzP+You iqTts8izJQ38foMjGMfCAETVNWcVTrV+2ZJu8DVp0ROtYeagq6LyunUP9EVd0rDq9oGU KLw3cvjgrjVG3iFKv2pY1TKgXnCS72q9232EKVO+29ubVPscRhYqdOtLBePJ44eMqaZD fVdb4wrHqozJ/Yj2h5NanIOuk8OV/QMMke02FyPr7Pg4dpTk0HbolTRyE5coe4cNzjlr BU2w==

List-archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/>

List-help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@noaa.gov>

List-id: <ferret_users.noaa.gov>

List-post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@noaa.gov>

Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx

Sender: owner-ferret_users@xxxxxxxx

Thread-index: AdNjMAxz+JcnC1u9Rv6p/OXAKLDwLQ==

Thread-topic: Nested loops and variable filename output

Dear Ferret users,

I’m trying to loop over some ocean data to calculate the thermocline and then output the results whilst retaining the original filename formatting. I’m using the following code:

define axis/f=1:103:1 ffile

define grid/f=ffile gfile

file/grid=gfile/format=DELIMITED/var="filenames" filenames.txt

file/grid=gfile/format=DELIMITED/var="output" output.txt

repeat/range=1:103:1/name=nf ( ;\

define symbol ff=`nf` ;\

set data `filenames[n=($ff),d=1]` ;\

let dtdz = thetao[z=@ddc] ;\

let dtdz_min = dtdz[z=@min] ;\

let zero_at_min = dtdz - dtdz_min ;\

let zt = zero_at_min[z=@loc:0] ;\

save/i=1/j=1/clob/ilimits=1:360/jlimits=1:181/ncformat=netcdf4/file=`output[n=($ff),d=2]` zt ;\

repeat/j=1:181 ( save/ncformat=netcdf4/app/file=`output[n=($ff),d=2]` zt ;\ ) ;\

cancel data 3 ;\

)

When running this script in Ferret 6.95 I get a signal = 11 crash

yes? go calculate_zt.jnl

define axis/f=1:103:1 ffile

define grid/f=ffile gfile

file/grid=gfile/format=DELIMITED/var="filenames" filenames.txt

file/grid=gfile/format=DELIMITED/var="output" output.txt

repeat/range=1:103:1/name=nf ( ;\

define symbol ff=`nf` ;\

set data `filenames[n=($ff),d=1]` ;\

let dtdz = thetao[z=@ddc] ;\

let dtdz_min = dtdz[z=@min] ;\

let zero_at_min = dtdz - dtdz_min ;\

let zt = zero_at_min[z=@loc:0] ;\

save/i=1/j=1/clob/ilimits=1:360/jlimits=1:181/ncformat=netcdf4/file=`output[n=($ff),d=2]` zt ;\

repeat/j=1:181 ( save/ncformat=netcdf4/app/file=`output[n=($ff),d=2]` zt ;\ ) ;\

cancel data 3 ;\

)

!-> repeat/range=1:103:1/name=nf ( ;define symbol ff=`nf` ;set data `filenames[n=($ff),d=1]` ;let dtdz = thetao[z=@ddc] ;let dtdz_min = dtdz[z=@min] ;let zero_at_min = dtdz - dtdz_min ;let zt = zero_at_min[z=@loc:0] ;save/i=1/j=1/clob/ilimits=1:360/jlimits=1:181/ncformat=netcdf4/file=`output[n=($ff),d=2]` zt ;repeat/j=1:181 ( save/ncformat=netcdf4/app/file=`output[n=($ff),d=2]` zt ;\ ) ;cancel data 3 ;)

!-> REPEAT: NF:1

!-> define symbol ff=1

!-> set data thetao_Omon_FIO-ESM_historical_r1i1p1_185001-185512_r360x181.nc

!-> LIST/FORMAT=CDF/i=1/j=1/clob/ilimits=1:360/jlimits=1:181/ncformat=netcdf4/file=zt_Omon_FIO-ESM_historical_r1i1p1_185001-185512_r360x181.nc zt

LISTing to file zt_Omon_FIO-ESM_historical_r1i1p1_185001-185512_r360x181.nc

!-> repeat/j=1:181 ( save/ncformat=netcdf4/app/file=`output[n=($ff),d=2]` zt ;\ )

!-> REPEAT: J=1

!-> LIST/FORMAT=CDF/ncformat=netcdf4/app/file=zt_Omon_FIO-ESM_historical_r1i1p1_185001-185512_r360x181.nc zt

LISTing to file zt_Omon_FIO-ESM_historical_r1i1p1_185001-185512_r360x181.nc

!-> REPEAT: J=2

**ERROR Ferret crash; signal = 11

I can avoid this crash by changing the last couple of lines to use a static filename and the loop counter but this is not ideal. E.g.:

save/i=1/j=1/clob/ilimits=1:360/jlimits=1:181/ncformat=netcdf4/file=zt_00($ff).nc zt ;\

repeat/j=1:181 ( save/ncformat=netcdf4/app/file=zt_00($ff).nc zt ) ;\

Cheers,

Ben