[MITgcm-support] Matlab query ?
Amitabh Mitra
amitabhmitra1 at rediffmail.com
Mon Jun 13 06:23:06 EDT 2005
Dear sir,
I have run the global_ocean.90x40x15 experiment
for one year with a resolution of 1x1 for the region
80S to 80N.(Data file is Attached).I transferred all the output files
to my PC (in D:) which has free space of 1.37 GB in C: drive.
When I run the program mit_loadglobal.m in Matlab
an error message is displayed showing shortage of
memory space as shown below:
-----------------------------------------------------------------------
Warning: buoyancyRelation not found in D:\MITrun\data
assuming buoyancyRelation = OCEANIC
Reading 36 time levels: 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180
??? Out of memory. Type HELP MEMORY for your options.
Error in ==> C:\matlabR12\work\rdmds.m
On line 131 ==> AA(r0(1):rN(1),r0(2):rN(2),r0(3):rN(3),iter)=A;
Error in ==> C:\matlabR12\work\mit_timesteps.m
On line 16 ==> [dummy, timesteps] = rdmds('D:\MITrun\uVeltave',NaN);
Error in ==> C:\matlabR12\work\mit_loadglobal.m
On line 31 ==> mit_timesteps
----------------------------------------------------------------------
Also I am unable to install Matlab in SGI-Origin.
So how could I plot my output?
Thanking You,
With regards,
Amitabh Mitra
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20050613/17a95cca/attachment.htm>
-------------- next part --------------
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
# tRef = 19*20.,
tRef = 24.0, 23.0, 22.0, 21.0, 20.0, 19.0,
18.0, 17.0, 16.0, 15.0, 14.0, 13.0,
12.0, 11.0, 10.0, 9.0, 8.0, 7.0, 6.0,
# sRef = 19*35.,
sRef = 34.65, 34.75, 34.82, 34.87, 34.90, 34.90,
34.86, 34.78, 34.69, 34.60, 34.58, 34.62,
34.68, 34.72, 34.73, 34.74, 34.73, 34.73,
34.72,
# no_slip_sides=.false.,
# no_slip_bottom=.TRUE.,
viscAr=1.E-3,
viscAh=4200.,
diffKhT=0.0,
diffKrT=3.E-5,
diffKhS=0.0,
diffKrS=3.E-5,
rhonil=1035.,
gravity=9.81,
eosType = 'JMD95Z',
ivdc_kappa=100.,
implicitDiffusion=.TRUE.,
useOldFreezing=.TRUE.,
useRealFreshWaterFlux=.TRUE.,
useCDscheme=.TRUE.,
useNHMTerms=.TRUE.,
# turn on looped cells
hFacMin=.05,
hFacMindr=50.,
# set precision of data files
readBinaryPrec=32,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=500,
cg2dTargetResidual=1.E-13,
&
# Time stepping parameters
&PARM03
# nIter0 = 0,
# nTimeSteps = 30,
# 100 years of integration will yield a reasonable flow field
startTime = 0.,
endTime = 31104000.,
deltaTmom = 800.0,
tauCD = 321428.,
deltaTtracer= 800.0,
deltaTClock = 172800.0,
# if you are using a version later than checkpoint45d on the main branch
# you can uncomment the following line and increase the time step
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the
# asynchronous time stepping
deltaTfreesurf = 172800.0,
abEps = 0.1,
pChkptFreq= 311040000.,
dumpFreq= 311040000.,
dumpFreq= 864000.,
taveFreq= 311040000.,
taveFreq= 864000.,
#monitorFreq=31104000.,
monitorFreq=1.,
# 2 months restoring timescale for temperature
tauThetaClimRelax = 5184000.0,
# 6 months restoring timescale for salinity
tauSaltClimRelax = 15552000.0,
periodicExternalForcing=.TRUE.,
externForcingPeriod=2592000.,
externForcingCycle=31104000.,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.FALSE.,
usingSphericalPolarGrid=.TRUE.,
delR= 10., 10., 10., 10., 20., 25., 25.,
25., 25., 50., 50., 50., 100.,
100., 100., 100., 100., 100.,100.,
phiMin=-80.,
dySpacing=1.,
dxSpacing=1.,
thetamin=0.,
&
# Input datasets
&PARM05
bathyFile= 'bathymetry.bin',
hydrogThetaFile='lev_t.bin',
hydrogSaltFile= 'lev_s.bin',
zonalWindFile= 'trenberth_taux1.bin',
meridWindFile= 'trenberth_tauy1.bin',
thetaClimFile= 'sst1.bin',
saltClimFile= 'lev_sss1.bin',
surfQFile= 'ncep_qnet.bin',
# fresh water flux is turned off, uncomment next line to turn on
# (not recommened together with surface salinity restoring)
# EmPmRFile= 'ncep_emp.bin',
&
More information about the MITgcm-support
mailing list