[Thread Prev][Thread Next][Index]

[ferret_users] .jnl file



Respected all
                  I wrote a .jnl file to extract some of the global attribute parameters(alt,day,month,etc ).Now I want to extract all the parameters for particular range of one parameters.The .jnl file is as follows:
**********************************************88
def symb start=1
def symb end=37
def axis/t=1:37:1/unit=hours tax
repeat/range=($start):($end)/name=yy (;\
use mday.`yy`.nc;\
let lat1=..botlat;\
let lat2=..edmaxlat;\
let lon1=..botlon;\
let lon2=..edmaxlon;\
let alt1=..edmaxalt;\
let month=..month;\
let edmax=..edmax;\
let lct=..botlct;\
let day=..day;\
let lat11=lat1+0*t[gt=tax];\
let lat12=lat2+0*t[gt=tax];\
let lon11=lon1+0*t[gt=tax];\
let lon12=lon2+0*t[gt=tax];\
let alt=alt1+0*t[gt=tax];\
let month1=month+0*t[gt=tax];\
let edmax1=edmax+0*t[gt=tax];\
let lct1=lct+0*t[gt=tax];\
let day1=day+0*t[gt=tax];\
let lat=max(lat11,lat12);\
let ed135=if elec_dens[x=70:700:0.001] lt 0 then 0 else elec_dens[x=70:700:0.001]*10^6;\
let ed_tec1=ed135[x=70:200@din]+0*t[gt=tax];\
let ed_tec2=ed135[x=200:350@din]+0*t[gt=tax];\
if `lat le 30` then ;\
list/file=mday9_ma21.nc/format=cdf/clob/l=`yy`/llimits=1:37 lat11,lat12,lon11,lon12,alt,month1,day1,lct1,edmax1,ed_tec1,ed_tec2;\
else ;\
list/app/file=mday9_ma22.nc/format=cdf/clob/l=`yy` lat11,lat12,lon11,lon12,alt,month1,day1,lct1,edmax1,ed_tec1,ed_tec2;\
endif ;\
can da/all ;\
)
************************************************8
 the error is as follows
!-> repeat/range=1:37/name=yy (;use mday`yy`.nc;let lat1=..botlat;let lat2=..edmaxlat;let lon1=..botlon;let lon2=..edmaxlon;let alt1=..edmaxalt;let month=..month;let edmax=..edmax;let lct=..botlct;let day=..day;let lat11=lat1+0*t[gt=tax];let lat12=lat2+0*t[gt=tax];let lon11=lon1+0*t[gt=tax];let lon12=lon2+0*t[gt=tax];let alt=alt1+0*t[gt=tax];let month1=month+0*t[gt=tax];let edmax1=edmax+0*t[gt=tax];let lct1=lct+0*t[gt=tax];let day1=day+0*t[gt=tax];let lat=max(lat11,lat12);let ed135=if elec_dens[x=70:700:0.001] lt 0 then 0 else elec_dens[x=70:700:0.001]*10^6;let ed_tec1=ed135[x=70:200@din]+0*t[gt=tax];let ed_tec2=ed135[x=200:350@din]+0*t[gt=tax];if `lat le 30` then ;list/file=mday9_ma21.nc/format=cdf/clob/l=`yy`/llimits=1:37 lat11,lat12,lon11,lon12,alt,month1,day1,lct1,edmax1,ed_tec1,ed_tec2;else ;list/app/file=mday9_ma22.nc/format=cdf/clob/l=`yy` lat11,lat12,lon11,lon12,alt,month1,day1,lct1,edmax1,ed_tec1,ed_tec2;endif ;can da/all ;)
!-> REPEAT: YY:1
 !-> SET DAT/FORM=CDF mday1.nc
 **ERROR: invalid command: grave accent doesnt evaluate to scalar
          lat le 30
if `lat le 30` then
Command file, command group, or REPEAT execution aborted
********************************
Ples. anyone help me....
 thanks  in advance
Gopal

[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement