Re: Superscripts / fcodes do not work

The first thing to check is your .hluresfile. Maybe you should try copying
the original one from your local computer to your home directory on the new
machine.

saji

On Tue, Nov 30, 2010 at 4:51 PM, Helen Parish <hparish@ess.ucla.edu> wrote:

> Does anyone know why the NCL statements that should create superscripts are
> not working in my script ?. I have used these statements before on my local
> computer where they seemed to work fine. However, I am running on another
> computer (a large multi-user machine), and they no longer seem to produce
> superscripts, but just print the codes exactly the way I typed them in (this
> machine is using NCL version 5.2.1).
>
> I had to replace "gsn_csm.ncl" with my own version, since I use a
> non-standard set of pressure levels. However, I suspect that this is not
> where the codes for doing superscripts are set. I think they may be set
> within "gsn_code.ncl", for which I use the standard version.
>
> Does anyone know what I need to do to get superscripts to print properly ?.
>
> I enclose the script I used and the plot which it produced.
>
> Thanks,
> Helen.
>
> 1) The plot
>
>
> 2) The script I used :
>
> ;***********************
> ; fourtest.ncl
> ;***********************
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> ;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "/p/lscratchd/parish3/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ;***********************
>
> begin
>
> diri = "./"
> fili = systemfunc("cd "+diri+" ; ls hera*.cam2.h0*nc")
> print(fili)
> nfili = dimsizes(fili)
>
> lev_p = (/ 91500., 90000., 87500., 80000., 67500., 50000., 30000.,
> 20000., 10000., 5000., 3000., 1500., 750., 390., 200., 100., 50.,25., 14.,
> 7., 3.5, 1.9, 0.95, 0.4, 0.1, 0.03 /)
>
> do nf = 0, nfili-1
>
> fi = addfile (diri+fili(nf), "r")
>
> nfa = 100.0
> nfb = nfa+nf
>
> diro = "./"
> filo = "uhelen"+nfb+".nc"
> print(filo)
> system ("/bin/rm -f "+diro+filo) ; remove any pre-exist file
> fo = addfile (diro+filo, "c")
>
> hyam = fi->hyam ; read hybrid info
> hybm = fi->hybm
> hyai = fi->hyai ; read hybrid info
> hybi = fi->hybi
> PS = fi->PS
> P0mb = 0.01*fi->P0
> U = fi->U
> time = fi->time
> lat = fi->lat
> lon = fi->lon
>
> lev_p!0 = "lev_p" ; variable and dimension name the
> same
> lev_p&lev_p = lev_p ; create coordinate variable
> lev_p@long_name = "pressure" ; attach some attributes
> lev_p@units = "hPa"
> lev_p@positive = "down"
>
> Up = vinth2p (U, hyam, hybm, lev_p ,PS, 1, P0mb, 2, False)
> copy_VarAtts(U, Up)
> fo->U = Up ; write to netCDF file
> fo->time = time ; write to netCDF file
> fo->lat = lat ; write to netCDF file
> fo->lon = lon ; write to netCDF file
>
> end do
>
> diri = "./"
> fili2 = systemfunc("cd "+diri+" ; ls uhelen*nc")
> print(fili2)
> nfili = dimsizes(fili2)
>
> f2 = addfiles (fili2,"r")
> U2 = addfiles_GetVar(f2,fili2,"U")
> lat2 = addfiles_GetVar(f2,fili2,"lat")
> lon2 = addfiles_GetVar(f2,fili2,"lon")
> lev_p2 = addfiles_GetVar(f2,fili2,"lev_p")
> time2 = addfiles_GetVar(f2,fili2,"time")
>
> printVarSummary( U2 )
> printVarSummary( lat2 )
> printVarSummary( lon2 )
> printVarSummary( time2 )
> printVarSummary( lev_p2 )
>
> dimt2 = dimsizes( U2 )
> ntim2 = dimt2(0)
> klvl2 = dimt2(1)
> nlat2 = dimt2(2)
> mlon2 = dimt2(3)
>
> print(nlat2)
>
> uzon = dim_avg_Wrap(U2) ; uzon(time,lev,lat)
>
> printVarSummary(uzon)
>
> ;***********************
> ; Create Plot
> ;***********************
> d = 1
> sm = 3
> pct = 0.1
>
> nfrq= ntim2/2
> LAT = 65.0 ; latitude 65 degrees N
>
> ; SPC(lev,freq)
> SPC = new( (/klvl2,nfrq/), typeof(U2),"No_FillValue")
>
> do kl=0,klvl2-1
> dof = specx_anal(uzon(lev_p|kl, {lat|LAT}, time|:),d,sm,pct)
> SPC(kl,:) = dof@spcx
> end do
>
> SPC!0 = "lev_p"
> SPC&lev_p = uzon&lev_p
>
> period = 1./dof@frq
> period!0 = "Period" ;
> SPC!1 = "period"
> SPC&period = period
>
> printVarSummary(SPC)
> printMinMax(SPC, True)
>
> wks = gsn_open_wks ("pdf", "spec.alt.65lat.daily.164.7" ) ; open
> workstation
> gsn_define_colormap(wks,"gui_default") ; choose colormap
>
>
> res = True
> res@cnFillOn = True
> res@lbLabelAutoStride = True
> res@gsnMaximize = True ; if [ps, eps, pdf] make large
> res@gsnSpreadColors = True ; span color map
>
> res@gsnCenterString = "Variance (m~S~2~N~ s~S~-2~N~)"
> res@tiXAxisString = "Period (Earth days)"
> res@tiYAxisString = "Pressure (mbar)"
>
> res@pmTitleZone = 3
> res@pmLabelBarOrthogonalPosF = 0.05
>
> res@lbTitleOn = True
> res@lbTitleDirection = "Across" ; title direction
> res@lbTitlePosition = "Right"
> res@lbTitleString = "(m2/s2)"
> res@lbTitleFontHeightF = 0.02
> res@lbTitleOffsetF = 0.02
>
>
> res@cnFillOn = True
> res@cnLineLabelsOn = False
> res@lbOrientation = "Horizontal"
> res@lbLabelAutoStride = True
> res@trYReverse = True
> res@gsnXAxisIrregular2Linear = True
> res@gsnYAxisIrregular2Linear = True
>
> plot = gsn_csm_pres_hgt(wks,SPC(lev_p | 0:20, period | :),res)
>
> end
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Nov 30 01:42:27 2010