[MITgcm-support] exf /obcs not opening forcing files
Martin Losch
mlosch at awi-bremerhaven.de
Mon Sep 25 02:50:09 EDT 2006
Hi Jill,
can you post all your "data.*" files and the content of you "code"
directory?
there is some "historical" connection between exf and ecco, but you
can move all of the exf-relevant stuff of ECCO_CPPOPTIONS.h to a
EXF_OPTIONS.h file
Martin
On Sep 25, 2006, at 7:03 AM, jschwarz at awi-bremerhaven.de wrote:
> Hi all,
>
> Another question from the beginner's corner: With exf and obcs
> packages listed in /code/package.config, and useOBCS=TRUE in
> data.pkg, and make reporting no errors; exf and obcs packages both
> look fine.. I am unable to convice exf to even attempt to open my
> open boundary initialisation files or the daily forcing data
> files. In the runtime output blurb, data.obcs is read, the open
> boundaries are reported, velocities at the boundary are initialised
> to junk (e+32) and t, s to the reference values. As data.exf is
> processed, the blurb actually reports that all the input filenames
> are ' '. This is despite setting ALLOW_OBCS_PRESCRIBE and careful
> checking that i'm reading and writing at precision 32bits, and
> specifying input files (which do exist) as in /verification and
> Matt's examples.
>
> It seems that exf only attempts to initialise from external forcing
> files if it finds #define ALLOW_ATM_TEMP/WIND etc; but i only found
> references to these settings in the ecco options file - surely ecco
> isn't required for exf?
>
> And also; when i switch on seaice, i get the following compile error:
> **************************
> In file budget.f:1396
>
> ATEMP(I,J,bi,bj)=MAX(273.16d0+MIN_ATEMP,ATEMP(I,J,bi,bj))
> 1
> Error: Statement function at (1) is recursive
> In file budget.f:1397
>
> LWDOWN(I,J,bi,bj)=MAX(MIN_LWDOWN,LWDOWN(I,J,bi,bj))
> 1
> Error: Statement function at (1) is recursive
> In file budget.f:1448
>
> & +D1*UG(I,J)*ATEMP(I,J,bi,bj)+D1I*UG(I,J)*AQH(I,J,bi,bj)
> 1
> Error: Function 'aqh' at (1) has no IMPLICIT type
> In file budget.f:1454
>
> & +D1*UG(I,J)*ATEMP(I,J,bi,bj)+D1I*UG(I,J)*AQH(I,J,bi,bj)
> 1
> Error: Function 'aqh' at (1) has no IMPLICIT type
> In file budget.f:1517
>
> QSWI(I,J,bi,bj)=-(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj)
> 1
> Error: Function 'swdown' at (1) has no IMPLICIT type
> make: *** [budget.o] Error 1
> ******************************************************
> .. any ideas why? (This is using checkpoint58)
>
> At the risk of being verbose.. i'm listing the initialisation
> section of the output text (a 3step trial with no seaice, using all
> dummy input files - no time variation, just constant fields, not
> that it matters because the files never get opened..) at the end
> of this.
> Any help appreciated!
> cheers,
> jill.
More information about the MITgcm-support
mailing list