wiki:ParNCLDeveloperInfo/TestingParNCL

Version 1 (modified by jayesh, 7 years ago) (diff)

Some info on testing ParNCL

Testing ParNCL

You can test ParNCL by using a simple script provided below. The data files used in the script are available with the ParGAL source (<PARGAL_SRC_DIR>/GridFiles/CAM/EulSpec/datasets/eul26x48x96.t[0-2].nc)

begin
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

end

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