[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