[Thread Prev][Thread Next][Index]

symbols in subscript unrecognized?



Hi ferret users

Here is another thing i don't quite understand. I have a script which
gets some parameters (i.e. names of variables) and defines some symbols
for those. It then calls a subscript and hands some parameters to this
one. Now look at the following script(s) and the error it produces:

   go calc_correlation_map ($var_reg) PC_reg
    !-> go calc_correlation_map rvar PC_reg
   !
---------------------------------------------------------------------------
   ! Description: calculates correlation map based of ERA data & indexes

   ! arguments  : $1 : variable name            [MSL]
   !            : $2 : correlation type         [point, NAO, SOI,
NAO_JONES,
   !                                             NAO_NCEP, 
NAO_NCEP_DJFM
   !                                             PminEpoint,SminEpoint
   !                                             PminEwGL, SminEwGL,
   !                                             SH_AOI, PC_reg         
]
   !            : $3 : longitude
   !            : $4 : latitude
   !            : $5 : lag                      [0.]
   !            : $6 : file name
   !
   !
--------------------------------------------------------------------------------
 
   define symbol var_cc = $1%MSL%          ! Processed Variable to
correlate to
    !-> define symbol var_cc = rvar                ! Processed Variable
to correlate to
   define symbol type = $2%point%          ! Type of regression
    !-> define symbol type = PC_reg                ! Type of regression
   define symbol lon = $3%20w%             ! longitude for point
regression
    !-> define symbol lon = 20w            ! longitude for point
regression
   define symbol lat = $4%60N%             ! latitude for point
regression
    !-> define symbol lat = 60N            ! latitude for point
regression
   define symbol lag = $5%0%               ! lag
    !-> define symbol lag = 0              ! lag
   define symbol file_name = $6%file%      ! file (regridding for
NAO_NCEP only)
    !-> define symbol file_name = file     ! file (regridding for
NAO_NCEP only)
    
 
   ! -------------------------------------------------------------------
   ! calculate normalized variable
     go calc_normal_var ($var_cc)
    !-> go calc_normal_var rvar
   !
---------------------------------------------------------------------------
   ! Description: calculate normalized variable var
   !
   ! Comment    : this script uses a number of working variable
($var)_dummy
   !              1) Use an expilicit region for the calculation.
   !              2) File must be set befor running the script
   !
   ! Output     : Var_norm       ->  normalized variable var
   !            : Var_std        ->  standard deviation of variable var
   !
   ! author     : ChA  08/97 [christof@climate.unibe.ch]
   ! adapted by : DM 09/1999 [mueller@climate.unibe.ch]
   !
   ! Usage      : GO calc_normal_var varibale_name
   !
   ! arguments  : $1 : variable name            [MSL]
   !
   !
---------------------------------------------------------------------------
   !cancel mode verify
   ! - define symbol var
   define symbol var = $1%MSL%
     !-> define symbol var = rvar
    
   show grid ($var)
    !-> show grid rvar
    **ERROR: invalid command: grid for rvar is not available.
             (LOAD this variable and try SHOW GRID again.)
   show grid rvar
   Command file, command group, or REPEAT execution aborted
   yes? show grid ($var)
    !-> show grid rvar
       GRID GIC1
    name       axis              # pts   start                end
    LONGITUDE LONGITUDE          144mr   0E                   2.5W
    LATITUDE  LATITUDE            29 r   20N                  90N
    normal    Z
    TIME1     T (MON)            472 r   23485                23956
   yes? 
   

So it seems, that the command show grid ($var) does work when back in
ferret typed in manually, but doesn't in the subscript.
I invoke the show grid command, because the following command in my
script, 
   define symbol min_L = `($var),RETURN=LSTART`
causes ferret tho hang for several minutes and then to terminate itself
without any further error massage.

The variable rvar is defined as follows: 
   define symbol fn=pres.mon.mean.jan1958_apr1997.des.20_90
   define symbol klimavar=pres_des
   ...
   use "/usr/users/mueller/dataset/mout/($fn).data.cdf"
   let rvar=($klimavar)[d=2]


which seems to be correct, because the following datasets are open:

   yes? show data
        currently SET data sets:
       1>
/usr/users/mueller/dataset/Wallace/AO_index_jan1958_apr1997.cdf
    name     title                             I         J        
K         L
    AOI      AO Index                         ...       ...      
...       1:472
    
       2>
/usr/users/mueller/dataset/mout/pres.mon.mean.jan1958_apr1997.des.20_90.data.cdf
    name     title                             I         J        
K         L
    PRES_DES SLP_DESEASONED                   1:144     1:29     
...       1:472
    TREND_STEIG                               1:144     1:29     
...       ...
    TREND_ABSN                                1:144     1:29     
...       ...
    PRES_DES_MEAN                             1:144     1:29     
...       ...
    
       3>
/usr/users/mueller/dataset/ncep/pres.mon.mean.jan1958_apr1997.cdf
    name     title                             I         J        
K         L
    PRES     Monthly mean surface pressure    1:144     1:73     
...       1:472
    
       4>
/usr/users/mueller/dataset/mout/weight_my/pres.mon.mean.jan1958_apr1997.des.20_90.DeTrend.DeMean.cdf       
(default)
    name     title                             I         J        
K         L
    EXPVAR   Explained Variance               ...       ...      
1:472     ...
    DEXPVAR  Error  Explained Variance        ...       ...      
1:472     ...
    EOF                                       1:144     1:29     
1:472     ...
    PC                                        ...       ...      
1:472     1:472
    PRES_DES_CALC                                1:144     1:29     
...       1:472
 
   yes? 


Is there anybody who can tell me what's happening (or why it isn't)? And
even more: Does anybody know a workaround for this?

Thank you in advance!

Daniel


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement