Re: [ferret_users] How to shade/plot by using scaling equation for data in TDS?

Dear Ferret-users,

I am very thankfull if someone could correct my script, since the result that I've got is always "No valid data" inside the box of the shade window.

The script is followings:

use "http://localhost:8080/thredds/dodsC/A20103212010328.L3m_8D_SST4_4"

! Latitude/Longitude Step of data: 0.041666668

! example region of interest with fine resolution

define axis/x=10e:50e:0.004535/unit=longitude xfine

define axis/y=10s:50s:0.004535/unit=longitude yfine

define grid/x=xfine/y=yfine gfine

    !    apply the Scaling Equation: (Slope*l3m_data) + Intercept = Parameter value
    !    Slope: 7.17185E-4
    !    Intercept: -2.0

LET/units="`l3m_data,RETURN=units`"/TITLE="sst" sst = ((7.17185E-4)*l3m_data[g=gfine])-2.0

SAVE/FILE="/home/dessy/sstx.nc" sst

use "/home/dessy/sstx.nc"

shade sst

Thank you very much for any help and suggestion.

Best wishes,

dessy

From: dessy berlianty <dessyberlianty@xxxxxxxxx>
To: "ferret_users@xxxxxxxx" <ferret_users@xxxxxxxx>
Sent: Tuesday, May 15, 2012 3:16 PM
Subject: [ferret_users] How to shade/plot by using scaling equation for data in TDS?

Dear Ferret-users,

I would like to ask a little bit silly question, I guess.

After successfull in access and read data from TDS, now I have problem like followings when I tried to shade/plot/contour by applied the scaling question as mentioned in global attributes:

yes? set memory/size=300

yes? go "/home/dessy/Documents/JNLFerret/readsmisst.jnl"
!

!    read data SMI/HDF from TDS
!
use "http://localhost:8080/thredds/dodsC/A20103212010328.L3m_8D_SST4_4"
show data/f
     currently SET data sets:
    1> http://localhost:8080/thredds/dodsC/A20103212010328.L3m_8D_SST4_4 (default)

name     title                             I         J         K         L
L3M_DATA                                  1:8640    1:4320    ...       ...
               on grid GBA1 with -1.E+34 & NaN for missing data
             X=0.5:8640.5 Y=0.5:4320.5
L3M_QUAL                                  1:8640    1:4320    ...       ...
               on grid GBA2 with -1.E+34 & NaN for missing data
             X=0.5:8640.5 Y=0.5:4320.5

let sst=(7.17185E-4*L3M_DATA)-2.0
shade/level=(25,31,1) sst
**ERROR Ferret crash; signal = 11
%

Notes from global attributes of OPeNDAP Dataset Access Form:
Scaling Equation: (Slope*l3m_data) + Intercept = Parameter value
Slope: 7.17185E-4
Intercept: -2.0

Is there any clue of command so I could make the figure from this type of data in Ferret directly?

Any help highly appreciated.

Thank you very much.

Best wishes,

dessy