[MITgcm-devel] exf_interp_read

[Martin Losch]

Hi there,

something has happened to exf_interp_read

After updating the code I get this:

mpiexec /opt/pce/bin/xd1launcher ../build/mitgcmuv
EXF_INTERP_READ: exf_interp_bufferSize too small
EXF_INTERP_READ: exf_interp_bufferSize too small
EXF_INTERP_READ: exf_interp_bufferSize too small
EXF_INTERP_READ: exf_interp_bufferSize too small
STOP

in EXF_PARAMS.h I see
C  To read input data without dynamical allocation  
(EXF_INTERP_USE_DYNALLOC undef),
C  buffer size currently set to 20000 (allows to read-in a 2x2 global  
data set)
       INTEGER    exf_interp_bufferSize
       PARAMETER( exf_interp_bufferSize = 10000 )

wouldn't it make sense to set the default really to 20000 (as  
promised in the comment) which is larger than 192x94=18048 (the  
"standard" size of ncep/core fields that we use for the cs510 for  
example)?

Plus with the fields that I downloaded from Dimitiris, I have a  
runoff file which is 360x180 = 64800.
I am not happy about having a local copy of EXF_PARAMS.h in every  
single experiment that uses EXF_USE_INTERPOLATION. How do you feel  
about this? Couldn't we have a better default (more like 64800)?


Martin