Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Activity
    Guillaume
    @PlanetGus
    Hello ?

    Quick question, I got this line with gmt5:

    psbasemap -R -J $FRAME -B+t"My Title" -O >> $filename

    But when I run the script, I got this error : psbasemap: Offending option +tMy Title ... What am I doing wrong ?

    Oh, and $FRAME="-Ba20f10/a10f10/WeSn"
    Guillaume
    @PlanetGus

    hum...
    I somehow solved the problem :
    FRAME="-Bxa20f10+l\"Longitude \" -Bya10f10+l\"Latitude\""

    And inside the loop
    psbasemap -R -J $FRAME -BWeSn+t"$day" -K -O >> $filename

    But the "longitude" and "latitude" doesn't appear (even if I change +l\"text\" by +l'text' or just +ltext) and there's no error hint

    Leonardo Uieda
    @leouieda
    Hi @PlanetGus thanks for posting! Sorry, I wasn't really checking this chat room yet. We're still trying to figure out the best way to handle user questions.
    This might be because you're setting the frame twice. Have you tried putting $FRAME after the -BWeSn..?
    Guillaume
    @PlanetGus
    Yes ... It doesn't change anything

    Before loop :
    FRAME='-Bxa30f10+lLongitude -Bya10f10+lLatitude -BWeSn'

    Inside loop :
    psbasemap -R -J $FRAME+t"$day"

    Result : I got the map, the frame and the title, but not the "Longitude" / "Latitude" label next to the axis

    Leonardo Uieda
    @leouieda
    @PlanetGus could you post your script here so that I can try to run it?
    Guillaume
    @PlanetGus
    Guillaume
    @PlanetGus

    Sure, It might be a little messy ...
    This one works. It create (main loop) 3 "mean state" map and (sub loop) several daily map of the same quantity.

    My first concern was : Am I able to use a mean value when lowering the resolution for arrow instead of a node value. => It has been dealt with grdfilter

    My second concern was : Can I do the same with the time serie. Apparently not, so bluntly generated a new file from matlab

    My third concern was : Can I add the date to the title. I worked it out with julian day reading and increment through a tier file.

    My last concern (cosmetic one) is : Why the labels don't work. It is not really useful, but as I'm just learning GMT I was curious.

    Eric Fielding
    @EJFielding
    I notice that gmt.soest.hawaii.edu is offline at the moment. Is there any word about how long it might be down?
    Leonardo Uieda
    @leouieda
    Hi @EJFielding it should be fixed now
    We had an outage on the weekend and no one was here to reboot
    Eric Fielding
    @EJFielding
    Thanks. Stuff happens.
    Guillaume
    @PlanetGus

    Hi !
    I try to use psxy to plot a box and it has a strange behaviour....
    The command is this :
    (projection is -Jmxx)

    gmt psxy -R -J -W2.0,black,-. -K -O << EOF >> $filename_track
    145 18
    145 -18
    250 -18
    250 18
    145 18
    EOF

    So it should create a rectangle right ?
    But the result is a strange curved rectangle ... Instead of having a straight line between {145,18} and {250,18} it plots an arch extended north linking the two coordinates.
    Same with {145,-18} and {250,-18} but the arch extends southward.

    Any idea ?

    Leonardo Uieda
    @leouieda
    @PlanetGus can you share an example output?
    Guillaume
    @PlanetGus
    cyclone_19971207.pdf
    Sure ... here it is .
    Guillaume
    @PlanetGus
    Ok ... c'est réglé, il me manquait l'option -A
    Guillaume
    @PlanetGus

    Hi ... I got this problem if anybody is willing to have a look at it :
    http://gmt.soest.hawaii.edu/boards/1/topics/8287

    Thanks,

    Leonardo Uieda
    @leouieda
    Hi @PlanetGus, we're all at UCSD for workshops and summit right now so that might be the reason for the delay
    sorry
    Guillaume
    @PlanetGus
    Oh ok... Thanks @leouieda ! Have a good workshop :)
    SquirrelKnight
    @SquirrelKnight
    Hi. I was wondering if anyone could provide advice for using PSCOUPE. I am currently attempting to plot symbols instead of focal mechanisms, such that the size of the symbols are dependent on the magnitude. Unfortunately, I can only get a uniform size with the -Fsc flag. I have tried with the -Sa command and without, but the results end up being the same.
    Here is an example from my script: awk '$34 == "N" {print $12,$11,$13,$22,$23,$24,$9,0,0,$1}' $focal_file | gmt pscoupe -Gblue $region $projection -Aa$min_lon/$min_lat/$max_lon/$max_lat/90/10/$min_depth/$max_depth -Fsc -Sa0.6 -O -K >> $grd_map # For HASH Normal events
    Please note that is from the command without the -Fsc flag.