<div dir="ltr">Kyle's right about reading error messages, specifically this one: <div><b><span style="font-size:13px">fatal:Undefined identifier: (cv_inv_calendar) is undefined, can't continue</span><br style="font-size:13px"><span style="font-size:13px">fatal:["Execute.c":8578]:</span><u style="font-size:13px"></u><span style="font-size:13px">Execu</span><span style="font-size:13px">te: Error occurred at or near line 103 in file IR_download_and_crop.ncl</span></b><br style="font-size:13px"><br style="font-size:13px"></div><div>It's <b>cd</b>_inv_calendar()</div><div><br></div><div><br></div><div>However there is still a problem. The below works on 6.0 and 6.2.0 but I get the same error as you in 6.2.1 on the lines using ut_inv_calendar. </div><div><br></div><div>This seems to be a backward compatible issue with the deprecated function and using cd_inv_calendar should solve the problem. </div><div><br></div><div><br></div><div>Alan. </div><div><br></div><div><br></div><div><br></div><div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div><br></div><div>begin </div><div>;# Set intital string variables</div><div>time_array=new((/240/),"double")</div><div>time_array@units = "hours ref 1-1-1 00:00:0.0"</div><div>;######## Begin the time loop ########</div><div>beg=15</div><div>day_ind = 0</div><div>hr_ind = 0</div><div>;########## SET THE TIMES ######## </div><div><br></div><div><div> time_array(48*day_ind+hr_ind*2) = cd_inv_calendar(2006,7,beg+day_ind+1,hr_ind+1,0,0,time_array@units, 0)</div><div> time_array(48*day_ind+hr_ind*2+1) = cd_inv_calendar(2006,7,beg+day_ind+1,hr_ind+1,30,0,time_array@units, 0)</div><div><br></div><div>print( time_array(48*day_ind+hr_ind*2) )</div><div>print( time_array(48*day_ind+hr_ind*2+1) ) </div></div><div><br></div><div><br></div><div> time_array(48*day_ind+hr_ind*2) = ut_inv_calendar(2006,7,beg+day_ind+1,hr_ind+1,0,0,time_array@units, 0)</div><div> time_array(48*day_ind+hr_ind*2+1) = ut_inv_calendar(2006,7,beg+day_ind+1,hr_ind+1,30,0,time_array@units, 0)</div><div><br></div><div>print( time_array(48*day_ind+hr_ind*2) )</div><div>print( time_array(48*day_ind+hr_ind*2+1) ) </div><div><br></div><div>end</div><br style="font-size:13px"><br style="font-size:13px"></div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 27, 2015 at 8:19 AM, Kyle Griffin <span dir="ltr"><<a href="mailto:ksgriffin2@wisc.edu" target="_blank">ksgriffin2@wisc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p dir="ltr">Conrado,</p>
<p dir="ltr">Please read your error messages more carefully! They tell you exactly what is going wrong. In this case, it tells me that you did not spell the "cd" part of the function name correctly...fix that so the name is as Dennis suggested and let us know how that goes.<br></p>
<p dir="ltr">Kyle</p>
<div class="gmail_quote"><div><div class="h5">On Feb 27, 2015 6:58 AM, <<a href="mailto:jorge.conrado@cptec.inpe.br" target="_blank">jorge.conrado@cptec.inpe.br</a>> wrote:<br type="attribution"></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><br>
<br>
<br>
Hi,<br>
<br>
I don't have experience with the NCL and I downloaded a script (IR_download_and_crop.ncl) in the Web to download images from the ftp address:<br>
<br>
<br>
<a href="ftp://disc2.nascom.nasa.gov/data/s4pa/TRMM_ANCILLARY/MERG/" target="_blank">ftp://disc2.nascom.nasa.gov/<u></u>data/s4pa/TRMM_ANCILLARY/MERG/</a><br>
<br>
<br>
This the script:<br>
<br>
<br>
<br>
<br>
;########## IR_download_and_crop.ncl #############<br>
;# This program downloads data from the NASA Merged IR database and crops to<br>
;# the specified region and time period. A NETCDF file of the data is made<br>
;# and the orignial fortran binary data is deleted.<br>
<br>
;# 1. Set inital variables<br>
;# 2. Read in 5 days of data over W. Africa write to NETCDF at 4km, 30 min<br>
<br>
load "$NCARG_ROOT/lib/ncarg/<u></u>nclscripts/csm/gsn_code.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/<u></u>nclscripts/csm/gsn_csm.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/<u></u>nclscripts/csm/skewt_func.ncl"<br>
load "$NCARG_ROOT/lib/ncarg/<u></u>nclscripts/csm/contributed.<u></u>ncl"<br>
<br>
begin<br>
<br>
;######## SET INITIAL VARS ########<br>
<br>
;# Reference information for the Merged IR Dataset<br>
;# For 2006...<br>
;# Jul begins on day 182 and ends on 212<br>
;# Aug begins on day 213 and ends on 243<br>
;# Sep begins on day 244 and ends on 273<br>
<br>
;# Set intital string variables<br>
hr=(/"00","01","02","03","04",<u></u>"05","06","07","08","09","10",<u></u>"11","12","13","14","15","16",<u></u>"17","18","19","20","21","22",<u></u>"23"/)<br>
<br>
day_31=(/"01","02","03","04","<u></u>05","06","07","08","09","10","<u></u>11","12","13","14","15","16","<u></u>17","18","19","20","21","22","<u></u>23","24","25","26","27","28","<u></u>29","30","31"/)<br>
<br>
day_30=(/"01","02","03","04","<u></u>05","06","07","08","09","10","<u></u>11","12","13","14","15","16","<u></u>17","18","19","20","21","22","<u></u>23","24","25","26","27","28","<u></u>29","30"/)<br>
<br>
time_array=new((/240/),"<u></u>double")<br>
<br>
time_array@units = "hours ref 1-1-1 00:00:0.0"<br>
<br>
;# For the July brightness temp. array a byte array is storred to reduce memory requirements and storage space.<br>
;# Convert to another data type and add 75.<br>
<br>
tb_array=new((/240,825,1925/),<u></u>"float")<br>
<br>
;######## Begin the time loop ########<br>
<br>
beg=15<br>
<br>
do day_ind=0,4<br>
do hr_ind=0,23<br>
<br>
file_name_download = "<a href="ftp://disc2.nascom.nasa.gov/data/s4pa/TRMM_ANCILLARY/MERG/2006/" target="_blank">ftp://disc2.nascom.nasa.gov/<u></u>data/s4pa/TRMM_ANCILLARY/MERG/<u></u>2006/</a>"+(beg+day_ind+182)+"/<u></u>merg_200607"+day_31(beg+day_<u></u>ind)+hr(hr_ind)+"_4km-pixel.Z"<br>
system("wget "+file_name_download)<br>
<br>
;####### READ THE FORTRAN BINARY ##############<br>
<br>
file_name_read_zip = "merg_200607"+day_31(beg+day_<u></u>ind)+hr(hr_ind)+"_4km-pixel.Z"<br>
file_name_read = "merg_200607"+day_31(beg+day_<u></u>ind)+hr(hr_ind)+"_4km-pixel"<br>
<br>
print(file_name_read)<br>
<br>
system("gunzip "+file_name_read_zip)<br>
<br>
;# Set Fortran read options<br>
<br>
setfileoption("bin","<u></u>ReadByteOrder","LittleEndian")<br>
ir_temp = byte2flt(fbindirread (file_name_read ,0, (/2,3298,9896/),"byte"))+75.<br>
<br>
system("rm -f "+file_name_read)<br>
<br>
;######## CROP THE IR GRID FOR WEST AFR #######<br>
;# LARGE: 319.998E(8795)-59.97576E(1648)<u></u>,40.0013N(549)-19.9772S(2197)<br>
;# LIMIT: 334.9874E(9207)-44.98637E(<u></u>1236),25.00667N(961)-4.<u></u>982578S(1785)<br>
<br>
lon=fspan(0.0182,0.0182+0.<u></u>036378335*9896,9896)<br>
lat=fspan(59.982,59.982-0.<u></u>036383683*3298,3298)<br>
<br>
lon@units = "degrees-east"<br>
lat@units = "degrees-north"<br>
<br>
ir_temp!2 = "lon"<br>
ir_temp!1 = "lat"<br>
ir_temp!0 = "time"<br>
<br>
ir_temp_re_0=new((/825,1925/),<u></u>"float")<br>
ir_temp_re_1=new((/825,1925/),<u></u>"float")<br>
<br>
ir_temp_re_0(:,:687)=ir_temp(<u></u>lat|961:1785,lon|9207:9894,<u></u>time|0)<br>
ir_temp_re_0(:,688:)=ir_temp(<u></u>lat|961:1785,lon|:1236,time|0)<br>
<br>
ir_temp_re_1(:,:687)=ir_temp(<u></u>lat|961:1785,lon|9207:9894,<u></u>time|1)<br>
ir_temp_re_1(:,688:)=ir_temp(<u></u>lat|961:1785,lon|:1236,time|1)<br>
<br>
lon_out = new((/1925/),"float")<br>
<br>
lon_out(:687)=lon(9207:9894)<br>
lon_out(688:)=lon(0:1236)<br>
<br>
lat_out = lat(961:1785)<br>
<br>
;########## PLACE IN ARRAY #######<br>
<br>
tb_array(48*day_ind+hr_ind*2,:<u></u>,:)=ir_temp_re_0 ;# First half hour<br>
tb_array(48*day_ind+hr_ind*2+<u></u>1,:,:)=ir_temp_re_1 ;# Second half hour<br>
<br>
;########## SET THE TIMES ########<br>
<br>
time_array(48*day_ind+hr_ind*<u></u>2) = ut_inv_calendar(2006,7,beg+<u></u>day_ind+1,hr_ind+1,0,0,time_<u></u>array@units, 0)<br>
time_array(48*day_ind+hr_ind*<u></u>2+1) = ut_inv_calendar(2006,7,beg+<u></u>day_ind+1,hr_ind+1,30,0,time_<u></u>array@units, 0)<br>
<br>
end do<br>
end do<br>
<br>
;######## CREATE THE NETCDF FILE ########<br>
<br>
ntim = 240<br>
nlat = 825<br>
nlon = 1925<br>
<br>
diro = "./" ; Output directory<br>
filo = "<a href="http://ir_test_jul_16_to_jul_20.nc" target="_blank">ir_test_jul_16_to_jul_20.nc</a>" ; Output file<br>
system("/bin/rm -f " + diro + filo) ; remove if exists<br>
fout = addfile (diro + filo, "c") ; open output file<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; explicitly declare file definition mode. Improve efficiency.<br>
;=============================<u></u>==============================<u></u>========<br>
<br>
setfileoption(fout,"<u></u>DefineMode",True)<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; create global attributes of the file<br>
;=============================<u></u>==============================<u></u>========<br>
<br>
fAtt = True ; assign file attributes<br>
fAtt@title = "NCL Efficient Approach to netCDF Creation"<br>
fAtt@source_file = "<a href="http://original-file.nc" target="_blank">original-file.nc</a>"<br>
fAtt@Conventions = "None"<br>
fAtt@creation_date = systemfunc ("date")<br>
fileattdef( fout, fAtt ) ; copy file attributes<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; predefine the coordinate variables and their dimensionality<br>
; Note: to get an UNLIMITED record dimension, we set the dimensionality<br>
; to -1 (or the actual size) and set the dimension name to True.<br>
;=============================<u></u>==============================<u></u>========<br>
<br>
dimNames = (/"time", "lat", "lon"/)<br>
dimSizes = (/ -1 , nlat, nlon/)<br>
dimUnlim = (/ True , False, False/)<br>
filedimdef(fout,dimNames,<u></u>dimSizes,dimUnlim)<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; predefine the the dimensionality of the variables to be written out<br>
;=============================<u></u>==============================<u></u>========<br>
; Here we are using NCL functions to facilitate defining<br>
; each variable's dimension name(s) and type.<br>
; The following could be replaced with explicit, user defined dimension<br>
; names different from those associated with the variable in memory.<br>
; Say, PS(time,lat,lon) in the NCL script. They could be redefined for the file via:<br>
; filevardef(fout, "PS" ,typeof(PS) ,(/"TIME","latitude","<u></u>longitude"/))<br>
;=============================<u></u>==============================<u></u>========<br>
<br>
filevardef(fout, "time" ,(/"double"/),(/"time"/))<br>
filevardef(fout, "lat" ,(/"float"/),(/"lat"/))<br>
filevardef(fout, "lon" ,(/"float"/),(/"lon"/))<br>
filevardef(fout, "TB" ,(/"float"/),(/"time","lat","<u></u>lon"/))<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; Copy attributes associated with each variable to the file<br>
; All attributes associated with each variable will be copied.<br>
;=============================<u></u>==============================<u></u>=========<br>
<br>
filevarattdef(fout,"time" ,time_array) ; copy time attributes<br>
filevarattdef(fout,"lat" ,lat_out) ; copy lat attributes<br>
filevarattdef(fout,"lon" ,lon_out) ; copy lon attributes<br>
filevarattdef(fout,"TB" ,tb_array) ; copy PS attributes<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; explicitly exit file definition mode. **NOT REQUIRED**<br>
;=============================<u></u>==============================<u></u>========<br>
<br>
setfileoption(fout,"<u></u>DefineMode",False)<br>
<br>
;=============================<u></u>==============================<u></u>========<br>
; output only the data values since the dimensionality and such have<br>
; been predefined. The "(/", "/)" syntax tells NCL to only output the<br>
; data values to the predefined locations on the file.<br>
;=============================<u></u>==============================<u></u>=========<br>
<br>
fout->time = (/time_array/)<br>
fout->lat = (/lat_out/)<br>
fout->lon = (/lon_out/)<br>
fout->TB = (/tb_array/)<br>
<br>
end<br>
<br>
<br>
<br>
I run it and the script download the images, but this the line<br>
<br>
<br>
time_array(48*day_ind+hr_ind*<u></u>2) = ut_inv_calendar(2006,7,beg+<u></u>day_ind+1,hr_ind+1,0,0,time_<u></u>array@units, 0)<br>
time_array(48*day_ind+hr_ind*<u></u>2+1) =<br>
<br>
<br>
I had the message:<br>
<br>
<br>
fatal:ut_inv_calendar: Invalid specification string<br>
fatal:["Execute.c":8578]:<u></u>Execute: Error occurred at or near line 103 in file IR_download_and_crop.ncl<br>
<br>
<br>
<br>
I already change the ut_inv_calendar to cd_inv_calendar but I had the same error.<br>
<br>
<br>
fatal:Undefined identifier: (cv_inv_calendar) is undefined, can't continue<br>
fatal:["Execute.c":8578]:<u></u>Execute: Error occurred at or near line 103 in file IR_download_and_crop.ncl<br>
<br>
<br>
<br>
<br>
fatal:ut_inv_calendar: Invalid specification string<br>
fatal:["Execute.c":8578]:<u></u>Execute: Error occurred at or near line 103 in file IR_download_and_crop.ncl<br>
<br>
<br>
I don't know how can I solve this error.<br>
<br>
<br>
Thank for all that help me with the question I asked yesterday.<br>
<br>
<br>
Conrado<br>
<br>
<br>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>