## 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)