Hi all,
My code was working under version 4.2.0.a033, but it stopped working
after NCL on the server was upgraded to 4.2.0.a034. The code is below,
trying to plot a map of Neoproterozoic (~600Ma) time. The orography file
was created by myself from a database somewhere also from NCAR. The
problem will appear when the code is running to the line
"mpDataSetName" : "./neo_2x2"
then it will stuck there, without giving any error message. I also tried
another orography file, the same thing happened again.
I have no idea how to fix this because there is no error message. Can
any one help me on this. Many of my scripts are dependent on this function.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
in = addfile("neo_orography_2x2.nc","r")
oro=in->oro(0,:,:)
latmp=in->lat
lonmp=in->lon
paleo_outline(oro,latmp,lonmp,1.,"neo_2x2")
wks = gsn_open_wks("ps","oro_test"); open a ps file
gsn_define_colormap(wks,"gui_default")
vpXF = 0.2 ; make plot bigger
vpYF = 0.8
vpWidthF = 0.6
vpHeightF = 0.6
mapangle = 30
bp = create "baseplot" mapPlotClass wks
"mpProjection" : "Mollweide"
"mpPerimOn" : True
"mpLimitMode" : "MaximalArea"
"mpDataSetName" : "./neo_2x2"
"mpGeophysicalLineThicknessF" : 2.0
end create
draw(bp)
frame(wks)
end
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Apr 01 2007 - 14:43:25 MDT
This archive was generated by hypermail 2.2.0 : Mon Apr 02 2007 - 08:01:14 MDT