[ferret_users] [FERRET]memory issues in ferret

[zhen.li@xxxxxxxxxxxxx]

Dear Ferret Users:

Lately, I've often run into the memory limit in ferret. I am trying to 
process a multi-year model output (see the script
below). The script fails at T=@VAR transform for both original, fine grid 
(1/4 x 1/6) and regridded coarse (1x1) grid.

The problem occurs after I executes "std_var_245=salt245[T=@VAR]^0.5", 
because std_var_245 still have the same dimension as salt245 with l=120. 
T=@VAR transform doesn't work for both the fine and coarse grid.

On the other hand, T=@VAR transform does work if I just use it on 
"salt[l=1,T=@VAR] " without the calculation of sigma and saltonsigma. I 
believe that after calculation of rho_un and ZAXREPLACE, the memory has 
been used up. So, T=@VAR transform can not be carried out. I don't know how 
to fix this problem.

Any help is greatly appreciated.

Below is my script.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
use ts_9099  [monthly salinity and temperature data for pacific ocean from 
year 1990 to 1999]
 [SALT     Salinity                                1:700     1:390 
1:43      1:120
 TEMP     Potential temperature           1:700     1:390     1:43      1:120]
Note: the data set is too large, so I have to regrid to 1deg, but even with 
1 deg data set,
I still run into the memory problem,

define axis/x=120E:295E:1/units=degrees x1deg
define axis/Y=45S:20N:1/units=degrees y1deg

let salt_1deg=salt[GX=x1deg@ave, GY=y1deg@ave]  >>>dimension [176 66 43 120]
let temp_1deg=salt[GX=x1deg@ave, GY=y1deg@ave]

let sigma=rho_un(salt_1deg,temp_1deg,0)-1000;
define axis/z=23:26:0.5 axden;

let saltonsigma=ZAXREPLACE(salt_1deg, sigma, z[gz=axden]);

let salt245=saltonsigma[z=24.5];
let std_var_245=salt245[T=@VAR]^0.5

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>



--
This message has been scanned for viruses and dangerous
content by GMAO's MailScanner, and is believed to be clean.