[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