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)