Testing ParNCL

You can test ParNCL by using a simple script provided below. The data files used in the script are available here

fili = systemfunc("ls eul26x48x96.t[0-2].nc")  ; 3 files

setfileoption("nc","SuppressClose",False)   ; alter default mode

f = addfiles(fili, "r") ; add the files

t = f[:]->T ; read temp

t_avg = dim_avg_n(t, 0) ; calc dim_avg

print(t_avg) ; print result


For correctness you can compare the result manually with output from the serial version of NCL.