[ncl-talk] Why cannot run my script for a long time

Guido Cioni guidocioni at gmail.com
Mon Sep 11 08:30:30 MDT 2017


Be sure to include ncl-talk when you're replying. 
The script is not enough. Which kind of error are you getting on the machine? 

Also the setup that you're describing is not clear to me.

> On 11. Sep 2017, at 16:28, aida farughy <farughy at gmail.com> wrote:
> 
> Yes of course
> It seems to be a problem with memory machin but I run this script on server.
> Please help me.
> Thanks
>  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>  l =addfile("/home/aida/Desktop/Newfolder/dec10-14.nc <http://dec10-14.nc/>","r")
>  z =l->msl(:,::-1,:)
>  loni =l->longitude
>  lati =l->latitude(::-1)
>  time =l->time
>  z :=(z*z at scale_factor+z@add_offset)/100
>  ntime =dimsizes(time)
>  nlat =dimsizes(lati)
>  nlon =dimsizes(loni)
>  
> ;;;;;;;;;;;;;;;
>     wks =gsn_open_wks("png","/home/aida/Desktop/dec2014.png")
>     res =True
>     res at gsnMaximize =True
>     res at gsnFrame =False
>     res at gsnDraw =False
>     res at tiMainString ="dec"
>     res at mpLimitMode ="LatLon"
>     res at mpMinLonF =10.0
>     res at mpMaxLonF =90.0
>     res at mpMinLatF =10.0
>     res at mpMaxLatF =70.0
>     res at mpOutlineOn =True
>     res at mpInlandWaterFillColor ="PaleTurquoise3"
>     res at mpOceanFillColor ="PaleTurquoise3"
>     map =gsn_csm_map(wks,res)
>     res_poly =True
>     res_poly at gsLineThicknessF =4.
>     dum =new((/ntime/),graphic)
>     dume =new((/ntime/),graphic)
>     t1 =new((/750/),integer)
>     xp =new((/750/),float)
>     yp =new((/750/),float)
>    ; xi =new((/750,750,750/),float)
> ;;;;;;;;;;
>  p =0
>  do gg =0,ntime-1
>       do j =26,nlat-7
>       do i =3,nlon-7
>      ; print(i+" "+j)
>     ...
>     ...
>     delete([/temp1D,temp,hkh,ilon,ilat,inds,m,minm,m1,h1,hk,z1z,mz1z/]) 
>  end do
> end do    
>   end do
> ;;;;;;;;;;;;;;;;;;;;;
> 
> ic =0
> do j =1,p
> lat1 =yp(j)
> lon1 =xp(j)
> t11 =t1(j)
> if (.not.ismissing(lat1)).and.(.not.ismissing(lon1)) then
>  k =j+1
> do while (k.le.p)
> ;print(lat1)  ;;;for check
>  res_poly at gsLineColor ="Black"
>  dum(ic) =gsn_add_polyline(wks,map,(/xp(j),xp(k)/),(/yp(j),yp(k)/),res_poly)
>  ic :=ic+1
>  j =j+1
>  k =k+1
> ...
> ...
>  delete([/gx,dx,lat2,lat1,lon2,lon1/])
> else
>  delete([/lat2,lon2/])
>  k =k+1
> end if
> end do
> ;print(j)
> end do
> ;;;;;;plot;;;;;;
> ;mkres =True
> ;mkres at gsMarkerIndex =16
> ;mkres at gsMarkerSizeF =0.01
> ;mkres at gsMarkerColor ="black"
> ;pl =gsn_add_polymarker(wks,map,xp,yp,mkres)
>  ;res at gsnAddCyclic =False
>  ;res3 =res
>  ;res3 at cnFillOn =True
>  ;res3 at cnFillPalette ="BlWhRe"
>  ;con =gsn_csm_contour(wks,l->msl(400,::-1,:),res)
>  ;zon =gsn_csm_contour(wks,z(89,:,:),res3)
>  ;overlay(map,con)
>  ;overlay(map,zon)
>  draw(map)
>  frame(wks)
>  delete([/p,xp,yp,ic/])
> 
> 
> 
> 
> On Mon, Sep 11, 2017 at 6:25 PM, Guido Cioni <guidocioni at gmail.com <mailto:guidocioni at gmail.com>> wrote:
> Hi,
> we need more details.
> 
> It seems to be a problem related to the machine and not NCL itself. The program will exit only upon errors or memory overload, not time limits...
> Also, please provide the script that you're using.
> Cheers 
>  
>> On 11. Sep 2017, at 15:51, aida farughy <farughy at gmail.com <mailto:farughy at gmail.com>> wrote:
>> 
>> 
>> Hi
>> When I running NCL for a long time, I get time out error and my plot doesn't go to the end of time.
>> How can I increase the array lenght to overcome this problem?
>> Best Regards
>> Aida
>> Sent from my iPad
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
> 
> 
> Guido Cioni
> http://guidocioni.altervista <http://guidocioni.altervista/>.org
> 
> 


Guido Cioni
http://guidocioni.altervista <http://guidocioni.altervista/>.org

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170911/2b46afd9/attachment.html>


More information about the ncl-talk mailing list