== Prototype ==
What running ParNCL will look like.
Now:
prompt> ncl average.ncl
Future:
prompt> mpiexec -np 4 parncl average.ncl
We agreed that the following simple NCL script will be the first thing we try to do with ParCAL
{{{
diri = "/fusion/group/climate/jacob/"
fili = "b40.1850.track1.1deg.006.cam1.T.0100-0110.nc" ; one file with 132 months of data.
f = addfile(diri+fili, "r") ; open the file for reading
T = f->T ; read Temperature field T with dimensions
;[time|132]x[lev|26]x[lat|192]x[lon|288]
Tavg=dim_avg_n(T,0) ;average T over the 0th dimension which is time.
printVarSummary(Tavg)
}}}
The above with graphics output
{{{
diri = "/fusion/group/climate/jacob/"
fili = "b40.1850.track1.1deg.006.cam1.T.0100-0110.nc" ; one file with 132 months of data.
f = addfile(diri+fili, "r") ; open the file for reading
T = f->T ; read Temperature field T with dimensions
;[time|132]x[lev|26]x[lat|192]x[lon|288]
Tavg=dim_avg_n(T,0) ;average T over the 0th dimension which is time.
wks = gsn_open_wks("ps","parvis_1") ; open a PS file
plot = gsn_csm_contour_map(wks,Tavg(0,:,:),False)
}}}