Hi,
I am setting up a parallel run with GETM 2.2 and it is failing here:
init_time
Time step: 1.50000000000000 seconds
Time format: 2
Start: 2007-12-01 00:00:00
Stop: 2007-12-01 01:00:00
==> 0 day(s) and 3600 seconds ==> 2400
micro time steps
init_domain test
FATAL ERROR: read_par_setup: Global grid sizes do not match
FATAL ERROR: Expected -1 by -1
FATAL ERROR: Read 1225 by 538
Looking at read_par_setup.F90:
if (iextr_read /= iextr .OR. jextr_read /= jextr) then
FATAL 'read_par_setup: Global grid sizes do not match'
FATAL ' Expected ',iextr,' by ',jextr
FATAL ' Read ',iextr_read,' by ',jextr_read
stop
end if
It seems that iextr and jextr are still -1 as initialised and haven’t been update to the domain size. Can you tell me where this should be read from?
The values in par_setup.dat (iextr_read, jextr_read) are the correct ones.
This is the &domain section of the getm.inp:
&domain
vert_cord = 1,
maxdepth = 55.0,
bathymetry = 'topo.nc',
longitude = -2.8,
latitude = 52.0,
f_plane = .true.,
openbdy = .true.,
bdyinfofile = 'bdyinfo.dat',
crit_depth = 0.5,
min_depth = 0.2,
kdum = 25,
ddu = -1.0,
ddl = -1.0,
d_gamma = 40.0,
gamma_surf = .true.,
il = -1,
ih = -1,
jl = -1,
jh = -1,
z0_method = 0,
z0_const = 0.005,
/
Thanks and I apologise for the corporate disclaimer at the end of this email,
Tiago