[ncl-talk] sres at mpCenterLatF won't accept decimal values?

[Marco Miani]

Dear community,


hope this e-mail finds you all well.


I am working on animating NCL-produced png figures, using satellite view, showing how my tracers travel around the world. To do it, I am observing on board of my satellite and, inside a while loop (looping over every time-step), I update center latitude and longitude (mpCenterLatF) by increasing the value (either linearly, or with a sine function...). Once all the frames (png) are output, I use a python script to stitch them together and get my nice animation in *.avi.

Now, the problem arises when adding a given value to lat, lon to mpCenterLatF inside the loop. It seems that if the added value is an integer (no decimals), all works well. IF, otherwise, I add a float (with decimals) NCL refuses to go on, stops, and complains. Here below, a snippet of my code. This is working in Python (PyNGL), but would crash in "pure" NCL.


The error message I got is: fatal:["NclAtt.c":262]:Attribute assignment type mismatch.

Any clues on how to fix it?

Thanks!

OS: Linux
Version: 6.6.2

Marco



sres at mpProjection               = "Satellite" ; choose map projection

sres at mpLimitMode                 = "angles"
sres at mpLeftAngleF                = 20.
sres at mpRightAngleF               = 20.
sres at mpTopAngleF                 = 20.
sres at mpBottomAngleF              = 20.

sres at pmTickMarkDisplayMode     = "Always"     ; turn on automatic tickmarks
sres at mpGridMaskMode            = "MaskLand"   ; Mask grid over land.
sres at mpAreaMaskingOn           = 1
sres at mpMaskAreaSpecifiers      = mask_specs

sres at mpCenterLonF               = 51       ; choose center lon
sres at mpCenterLatF               = 25       ; choose center lat

sres at mpSatelliteAngle1F         = 1.4*55.2957795130823*asin(1./1.4)/8.
sres at mpSatelliteAngle2F         = 90.
sres at mpSatelliteDistF           = 2.35     ; choose satellite view

i = 0
Nsteps= 60
   do while(i.le.Nsteps)
      ..do some magic
      sres at mpCenterLonF               = 51+i      ; Works only if "i" is integer, NOT float...
      ;sres at mpCenterLatF               = 25.2      ; This won't work! <-------------------------------------- this is the trouble maker!
      sres at mpCenterLatF               = 25        ; This does work
      ...do more magic
      print("")
      print(sprinti("Drawing frame %0.1i ",i))
      i=i+1
   end do
print("All done. Reading for merging...")        ;Next step is to call python script, to generate an animation...


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200602/09dfd3b7/attachment.html>