[Thread Prev][Thread Next][Index]
[ferret_users] [FERRET]memory issues in ferret
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 command
"std_var_245=salt245[T=@VAR]^0.5", since
"std_var_245" is still a 3D variable with L=120 (time
dimension). T=@VAR transform doesn't seem to work for both fine and
coarse grid.
On the other hand, T=@VAR transform does work if I just use it on
"salt[k=1,T=@VAR] " without 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 MailScanner and is believe to be clean.
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement