[MITgcm-support] exf /obcs not opening forcing files

jschwarz at awi-bremerhaven.de jschwarz at awi-bremerhaven.de
Mon Sep 25 01:03:49 EDT 2006


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.
----------------------------------------------------------
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) //                      MITgcm UV
(PID.TID 0000.0001) //                      =========
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // execution environment starting up...
(PID.TID 0000.0001)
(PID.TID 0000.0001) // MITgcmUV version:  checkpoint58
(PID.TID 0000.0001) // Build user:        root
(PID.TID 0000.0001) // Build host:        tasman
(PID.TID 0000.0001) // Build date:        Sat Sep 23 08:11:18 NZST 2006
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Execution Environment parameter file "eedata"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># Example "eedata" file
(PID.TID 0000.0001) ># Lines beginning "#" are comments
(PID.TID 0000.0001) ># nTx - No. threads per process in X
(PID.TID 0000.0001) ># nTy - No. threads per process in Y
(PID.TID 0000.0001) > &EEPARMS
(PID.TID 0000.0001) > nTx=1,
(PID.TID 0000.0001) > nTy=1,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) ># Note: Some systems use & as the
(PID.TID 0000.0001) ># namelist terminator. Other systems
(PID.TID 0000.0001) ># use a / character (as shown here).
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Computational Grid Specification ( see files "SIZE.h" )
(PID.TID 0000.0001) //                                  ( and "eedata"       )
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)      nPx =    1 ; /* No. processes in X */
(PID.TID 0000.0001)      nPy =    1 ; /* No. processes in Y */
(PID.TID 0000.0001)      nSx =    1 ; /* No. tiles in X per process */
(PID.TID 0000.0001)      nSy =    1 ; /* No. tiles in Y per process */
(PID.TID 0000.0001)      sNx =  120 ; /* Tile size in X */
(PID.TID 0000.0001)      sNy =   82 ; /* Tile size in Y */
(PID.TID 0000.0001)      OLx =    3 ; /* Tile overlap distance in X */
(PID.TID 0000.0001)      OLy =    3 ; /* Tile overlap distance in Y */
(PID.TID 0000.0001)      nTx =    1 ; /* No. threads in X per process */
(PID.TID 0000.0001)      nTy =    1 ; /* No. threads in Y per process */
(PID.TID 0000.0001)       Nr =   21 ; /* No. levels in the vertical   */
(PID.TID 0000.0001)       nX =  120 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
(PID.TID 0000.0001)       nY =   82 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */
(PID.TID 0000.0001)   nTiles =    1 ; /* Total no. tiles per process ( = nSx*nSy ) */
(PID.TID 0000.0001)   nProcs =    1 ; /* Total no. processes ( = nPx*nPy ) */
(PID.TID 0000.0001) nThreads =    1 ; /* Total no. threads per process ( = nTx*nTy ) */
(PID.TID 0000.0001) usingMPI =    F ; /* Flag used to control whether MPI is in use */
(PID.TID 0000.0001)                   /*  note: To execute a program with MPI calls */
(PID.TID 0000.0001)                   /*  it must be launched appropriately e.g     */
(PID.TID 0000.0001)                   /*  "mpirun -np 64 ......"                    */
(PID.TID 0000.0001) useCoupler=   F ; /* Flag used to control communications with */
(PID.TID 0000.0001)                   /*  other model components, through a coupler */
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Mapping of tiles to threads
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // -o- Thread   1, tiles (   1:   1,   1:   1)
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile <-> Tile connectvity table
(PID.TID 0000.0001) // ======================================================
(PID.TID 0000.0001) // Tile number: 000001 (process no. = 000001)
(PID.TID 0000.0001) //        WEST: Tile = 000001, Process = 000001, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //        EAST: Tile = 000001, Process = 000001, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //       SOUTH: Tile = 000001, Process = 000001, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001) //       NORTH: Tile = 000001, Process = 000001, Comm = put
(PID.TID 0000.0001) //                bi = 000001, bj = 000001
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model parameter file "data"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) ># | Model parameters |
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># Continuous equation parameters
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#   tRef            - Reference vertical potential temperature          (deg C)
(PID.TID 0000.0001) >#   sRef            - Reference vertical salinity                         (PSU)
(PID.TID 0000.0001) >#   viscAh          - Horizontal eddy viscosity coefficient             (m^2/s)
(PID.TID 0000.0001) >#   viscAz          - Vertical eddy viscosity coefficient               (m^2/s)
(PID.TID 0000.0001) >#   diffKhT         - Horizontal temperature diffusivity                (m^2/s)
(PID.TID 0000.0001) >#   diffKzT         - Vertical temperature diffusivity                  (m^2/s)
(PID.TID 0000.0001) >#   diffKhS         - Horizontal salt diffusivity                       (m^2/s)
(PID.TID 0000.0001) >#   diffKzS         - Vertical salt diffusivity                         (m^2/s)
(PID.TID 0000.0001) >#   f0              - Reference coriolis parameter,
(PID.TID 0000.0001) >#                     south edge of f on beta plane                       (1/s)
(PID.TID 0000.0001) >#   beta            - df/dy                                         (s^-1.m^-1)
(PID.TID 0000.0001) >#   tAlpha          - Linear EOS thermal expansion coefficient           (1/oC)
(PID.TID 0000.0001) >#   sBeta           - Linear EOS haline contraction coefficient         (1/ppt)
(PID.TID 0000.0001) >#   gravity         - Acceleration due to gravity                       (m/s^2)
(PID.TID 0000.0001) >#   gBaro           - Accel. due to gravity used in barotropic equation (m/s^2)
(PID.TID 0000.0001) >#   rigidLid        - Set to true to use rigid lid
(PID.TID 0000.0001) >#   implicitFreeSurface - Set to true to use implicit free surface
(PID.TID 0000.0001) >#   eosType         - Flag for linear or polynomial equation of state
(PID.TID 0000.0001) >#   momAdvection    - On/Off flag for momentum self transport
(PID.TID 0000.0001) >#   momViscosity    - On/Off flag for momentum mixing
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &PARM01
(PID.TID 0000.0001) > tRef= 1.3 , 1.3 , 1.3 , 1.3 , 1.3 ,
(PID.TID 0000.0001) >       1.3 , 1.3 , 1.3 , 1.3 , 1.3 ,
(PID.TID 0000.0001) >       1.3 , 1.3 , 1.3 , 1.3 , 1.3 ,
(PID.TID 0000.0001) >       1.3 , 1.3 , 1.3 , 1.3 , 1.3 ,
(PID.TID 0000.0001) >       1.3 ,
(PID.TID 0000.0001) > sRef= 33.65, 33.75, 33.82, 33.87, 33.90,
(PID.TID 0000.0001) >       33.90, 33.86, 33.78, 33.69, 33.60,
(PID.TID 0000.0001) >       33.58, 33.62, 33.68, 33.72, 33.73,
(PID.TID 0000.0001) >       33.74, 33.73, 33.73, 33.72, 33.72,
(PID.TID 0000.0001) >       33.71,
(PID.TID 0000.0001) > no_slip_sides=.FALSE.,
(PID.TID 0000.0001) > no_slip_bottom=.TRUE.,
(PID.TID 0000.0001) > viscAz=1.93e-5,
(PID.TID 0000.0001) > viscAh=5.E4,
(PID.TID 0000.0001) > diffKhT=0.0,
(PID.TID 0000.0001) > diffKzT=1.46e-5,
(PID.TID 0000.0001) > diffKhS=1.e3,
(PID.TID 0000.0001) > diffKzS=1.46e-5,
(PID.TID 0000.0001) > tAlpha=2.E-4,
(PID.TID 0000.0001) > sBeta =7.4E-4,
(PID.TID 0000.0001) > rigidLid=.FALSE.,
(PID.TID 0000.0001) > implicitFreeSurface=.TRUE.,
(PID.TID 0000.0001) > nonHydrostatic=.TRUE.,
(PID.TID 0000.0001) > eosType='JMD95Z',
(PID.TID 0000.0001) > readBinaryPrec=32,
(PID.TID 0000.0001) > writeBinaryPrec=32,
(PID.TID 0000.0001) > saltStepping=.FALSE.,
(PID.TID 0000.0001) > tempStepping=.FALSE.,
(PID.TID 0000.0001) > momStepping=.TRUE.,
(PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
(PID.TID 0000.0001) > implicitViscosity=.TRUE.,
(PID.TID 0000.0001) > useSingleCpuIO=.TRUE.,
(PID.TID 0000.0001) > allowFreezing=.TRUE.,
(PID.TID 0000.0001) > HeatCapacity_Cp = 3986.D0
(PID.TID 0000.0001) > gravity         = 9.8156,
(PID.TID 0000.0001) > gBaro           = 9.8156,
(PID.TID 0000.0001) > rhoNil          = 1027.D0
(PID.TID 0000.0001) > useCDscheme=.TRUE.,
(PID.TID 0000.0001) > useNHMTerms=.TRUE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Elliptic solver parameters
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#   cg2dMaxIters       - Maximum number of 2d solver iterations
(PID.TID 0000.0001) >#   cg2dTargetResidual - Solver target residual
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &PARM02
(PID.TID 0000.0001) > cg2dMaxIters=500,
(PID.TID 0000.0001) > cg2dTargetResidual=1.D-8
(PID.TID 0000.0001) > cg3dMaxIters=500,
(PID.TID 0000.0001) > cg3dTargetResidual=1.D-8
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Time stepping parameters
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#   startTime         - Integration starting time                (s)
(PID.TID 0000.0001) >#   endTime           - Integration ending time                  (s)
(PID.TID 0000.0001) >#   tauCD             - CD scheme coupling timescale             (s)
(PID.TID 0000.0001) >#   deltaTMom         - Timestep for momemtum equations          (s)
(PID.TID 0000.0001) >#   deltaTtracer      - Tracer timestep                          (s)
(PID.TID 0000.0001) >#   deltaTClock       - Timestep used as model "clock"           (s)
(PID.TID 0000.0001) >#   abEps             - Adams-Bashforth stabilising factor
(PID.TID 0000.0001) >#   pChkPtFreq        - Frequency of permanent check pointing    (s)
(PID.TID 0000.0001) >#   chkPtFreq         - Frequency of rolling check pointing      (s)
(PID.TID 0000.0001) >#   dumpFreq          - Frequency at which model state is stored (s)
(PID.TID 0000.0001) >#   tauThetaClimRelax - Relaxation to climatology time scale     (s)
(PID.TID 0000.0001) >#   tauSaltClimRelax  - Relaxation to climatology time scale     (s)
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &PARM03
(PID.TID 0000.0001) > tauCD=172800.,
(PID.TID 0000.0001) > startTime=0.0,
(PID.TID 0000.0001) > endTime=64800.,
(PID.TID 0000.0001) > deltaTmom=21600.0,
(PID.TID 0000.0001) > deltaTtracer=21600.0,
(PID.TID 0000.0001) > deltaTClock =21600.0,
(PID.TID 0000.0001) > cAdjFreq=0.,
(PID.TID 0000.0001) > abEps=0.1,
(PID.TID 0000.0001) > pChkptFreq=216000.,
(PID.TID 0000.0001) > chkptFreq= 0.,
(PID.TID 0000.0001) > dumpFreq = 21600.,
(PID.TID 0000.0001) > tavefreq = 216000.,
(PID.TID 0000.0001) > tauSaltClimRelax = 172800.,
(PID.TID 0000.0001) > tauThetaClimRelax = 172800.,
(PID.TID 0000.0001) > monitorFreq=1.,
(PID.TID 0000.0001) > forcing_In_AB = .FALSE.
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Gridding parameters
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#   usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
(PID.TID 0000.0001) >#   usingCartesianGrid      - On/Off flag for selecting cartesian coordinates
(PID.TID 0000.0001) >#   delX                    - Zonal grid spacing         (degrees)
(PID.TID 0000.0001) >#   delY                    - Meridional grid spacing    (degrees)
(PID.TID 0000.0001) >#   delZ                    - Vertical grid spacing      (m)
(PID.TID 0000.0001) >#   phiMin                  - Southern boundary latitude (degrees)
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &PARM04
(PID.TID 0000.0001) > usingCartesianGrid=.FALSE.,
(PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
(PID.TID 0000.0001) > delXfile = 'DELX3.ross_1deg_sn_ciso',
(PID.TID 0000.0001) > delYfile = 'DELY3.ross_1deg_sn_ciso',
(PID.TID 0000.0001) > delZ= 5., 5., 20., 20., 25., 25., 50., 50., 100.,
(PID.TID 0000.0001) >       200., 250., 250., 500., 500., 500., 500.,
(PID.TID 0000.0001) >       500., 500., 500., 500., 500.,
(PID.TID 0000.0001) > phimin=-80.,
(PID.TID 0000.0001) > thetamin=140.,
(PID.TID 0000.0001) > rSphere         = 6371.D3
(PID.TID 0000.0001) > hFacMin = 0.3,
(PID.TID 0000.0001) > hFacMinDz = 50.,
(PID.TID 0000.0001) ># hFacInf = 0.05,
(PID.TID 0000.0001) ># hFacSup = 0.95,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Input datasets
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#   bathyFile       - File containing bathymetry
(PID.TID 0000.0001) >#   hydrogThetaFile - File containing initial potential temperature data
(PID.TID 0000.0001) >#   hydrogSaltFile  - File containing initial salinity data
(PID.TID 0000.0001) >#   zonalWindFile   - File containing zonal wind data
(PID.TID 0000.0001) >#   meridWindFile   - File containing meridional wind data
(PID.TID 0000.0001) >#   thetaClimFile   - File containing theta climatology used for relaxation
(PID.TID 0000.0001) >#   saltClimFile    - File containing salt climatology used for relaxation
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &PARM05
(PID.TID 0000.0001) > bathyFile       = 'bathy_etgebda_sn_iso1deg_b.bin',
(PID.TID 0000.0001) ># hydrogThetaFile = 'theta00_g1fd_sn.bin',
(PID.TID 0000.0001) ># hydrogSaltFile  = 'sal00_g1fd_sn.bin',
(PID.TID 0000.0001) ># thetaClimFile   = 'theta00_g1fd_sn.bin',
(PID.TID 0000.0001) ># saltClimFile    = 'sal00_g1fd_sn.bin',
(PID.TID 0000.0001) > hydrogThetaFile='dummy_theta.bin',
(PID.TID 0000.0001) > hydrogSaltFile='dummy_salinity.bin',
(PID.TID 0000.0001) > thetaClimFile='dummy_theta.bin',
(PID.TID 0000.0001) > saltClimFile='dummy_salinity.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001)
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM01
(PID.TID 0000.0001) S/R INI_PARMS ; read PARM01 : OK
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM02
(PID.TID 0000.0001) S/R INI_PARMS ; read PARM02 : OK
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM03
(PID.TID 0000.0001) S/R INI_PARMS ; read PARM03 : OK
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM04
(PID.TID 0000.0001) S/R INI_PARMS ; read PARM04 : OK
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM05
(PID.TID 0000.0001) S/R INI_PARMS ; read PARM05 : OK
(PID.TID 0000.0001)  PACKAGES_BOOT: opening data.pkg
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.pkg
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.pkg"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># Packages
(PID.TID 0000.0001) > &PACKAGES
(PID.TID 0000.0001) >  useGMRedi = .TRUE.,
(PID.TID 0000.0001) >  useKPP    = .TRUE.,
(PID.TID 0000.0001) >  useSEAICE = .FALSE.,
(PID.TID 0000.0001) >  useMNC    = .TRUE.,
(PID.TID 0000.0001) >  useOBCS   = .TRUE.,
(PID.TID 0000.0001) >  useGrdchk = .FALSE.,
(PID.TID 0000.0001) >  useDiagnostics = .FALSE.,
(PID.TID 0000.0001) >#  useEXF    = .TRUE.,
(PID.TID 0000.0001) >#  useCAL    = .TRUE.,
(PID.TID 0000.0001) >#  useCDscheme = .TRUE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001)  PACKAGES_BOOT: finished reading data.pkg
(PID.TID 0000.0001)  MNC_READPARMS: opening file 'data.mnc'
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.mnc
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.mnc"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># MNC output
(PID.TID 0000.0001) > &MNC_01
(PID.TID 0000.0001) > mnc_use_outdir=.TRUE.,
(PID.TID 0000.0001) > mnc_outdir_str='../mnc_out',
(PID.TID 0000.0001) > mnc_outdir_date=.TRUE.,
(PID.TID 0000.0001) > snapshot_mnc=.TRUE.,
(PID.TID 0000.0001) > timeave_mnc=.TRUE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001)  MNC_READPARMS: finished reading data.mnc
(PID.TID 0000.0001)  GM_READPARMS: opening data.gmredi
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.gmredi
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.gmredi"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># GMREDI parameters
(PID.TID 0000.0001) > &GM_PARM01
(PID.TID 0000.0001) >  GM_background_K = 571.0
(PID.TID 0000.0001) >  GM_taper_scheme = 'ldd97'
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001)  GM_READPARMS: finished reading data.gmredi
(PID.TID 0000.0001)  KPP_INIT: opening data.kpp
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.kpp
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.kpp"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ***************
(PID.TID 0000.0001) ># Open boundaries
(PID.TID 0000.0001) ># ***************
(PID.TID 0000.0001) > &KPP_PARM01
(PID.TID 0000.0001) > KPPmixingMaps   = .TRUE.,
(PID.TID 0000.0001) > KPPwriteState   = .TRUE.,
(PID.TID 0000.0001) > kpp_dumpFreq    = 21600.0
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001)  KPP_INIT: finished reading data.kpp
(PID.TID 0000.0001)  OBCS_READPARMS: opening data.obcs
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.obcs
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.obcs"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># Open-boundaries
(PID.TID 0000.0001) > &OBCS_PARM01
(PID.TID 0000.0001) > OB_Ieast=81*-1,
(PID.TID 0000.0001) > OB_Iwest=81*1,
(PID.TID 0000.0001) ># OB_Jnorth=120*0,
(PID.TID 0000.0001) > useOBCSsponge=.TRUE.,
(PID.TID 0000.0001) > useOBCSprescribe=.TRUE.,
(PID.TID 0000.0001) > useOBCSbalance=.FALSE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) >#closing northern boundary after mike williams advice
(PID.TID 0000.0001) ># OB_Jnorth=120*1,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#useOrlanskiNorth=.TRUE.,
(PID.TID 0000.0001) >#useOrlanskiSouth=.TRUE.,
(PID.TID 0000.0001) >#useOrlanskiEast=.TRUE.,
(PID.TID 0000.0001) >#useOrlanskiWest=.FALSE.,
(PID.TID 0000.0001) ># useOBCSprescribe = .TRUE.,
(PID.TID 0000.0001) ># OBNuFile = 'OBmeridU.bin',
(PID.TID 0000.0001) ># OBSuFile = 'OBmeridU.bin',
(PID.TID 0000.0001) > OBWuFile='OBW_uvel_jan2000ecco_sn_b.data',
(PID.TID 0000.0001) > OBWvFile='OBW_vvel_jan2000ecco_sn_b.data',
(PID.TID 0000.0001) > OBWwFile='OBW_wvel_jan2000ecco_sn_b.data',
(PID.TID 0000.0001) > OBWtFile='obw_theta_sn_flip.bin',
(PID.TID 0000.0001) > OBWsFile='obw_salinity_sn_flip.bin',
(PID.TID 0000.0001) > OBEuFile='OBE_uvel_jan2000ecco_sn_b.data',
(PID.TID 0000.0001) > OBEvFile='OBE_vvel_jan2000ecco_sn_b.data',
(PID.TID 0000.0001) > OBEwFile='OBE_wvel_jan2000ecco_sn_b.data',
(PID.TID 0000.0001) > OBEtFile='obe_theta_sn_flip.bin',
(PID.TID 0000.0001) > OBEsFile='obe_salinity_sn_flip.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Orlanski parameters
(PID.TID 0000.0001) ># &OBCS_PARM02
(PID.TID 0000.0001) >#Cmax=0.45,
(PID.TID 0000.0001) >#cVelTimeScale=1000.,
(PID.TID 0000.0001) ># &
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Boundary parameters
(PID.TID 0000.0001) > &OBCS_PARM03
(PID.TID 0000.0001) > spongeThickness=5,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
 OB Jn =           0           0           0           0           0           0           0           0
     0           0           0           0           0           0           0           0           0
   0           0           0           0           0           0           0           0           0
 0           0           0           0           0           0           0           0           0           0
           0           0           0           0           0           0           0           0           0
         0           0           0           0           0           0           0           0           0
       0           0           0           0           0           0           0           0           0
     0           0           0           0           0           0           0           0           0
   0           0           0           0           0           0           0           0           0
 0           0           0           0           0           0           0           0           0           0
           0           0           0           0           0           0           0           0           0
         0           0           0           0           0           0           0           0           0
       0           0           0           0           0           0           0           0           0
     0           0
 OB Js =           0           0           0           0           0           0           0           0
     0           0           0           0           0           0           0           0           0
   0           0           0           0           0           0           0           0           0
 0           0           0           0           0           0           0           0           0           0
           0           0           0           0           0           0           0           0           0
         0           0           0           0           0           0           0           0           0
       0           0           0           0           0           0           0           0           0
     0           0           0           0           0           0           0           0           0
   0           0           0           0           0           0           0           0           0
 0           0           0           0           0           0           0           0           0           0
           0           0           0           0           0           0           0           0           0
         0           0           0           0           0           0           0           0           0
       0           0           0           0           0           0           0           0           0
     0           0
 OB Ie =         120         120         120         120         120         120         120         120
   120         120         120         120         120         120         120         120         120
 120         120         120         120         120         120         120         120         120         1
20         120         120         120         120         120         120         120         120         120
         120         120         120         120         120         120         120         120         120
       120         120         120         120         120         120         120         120         120
     120         120         120         120         120         120         120         120         120
   120         120         120         120         120         120         120         120         120
 120         120         120         120         120         120         120         120         120
 0
 OB Iw =           1           1           1           1           1           1           1           1
     1           1           1           1           1           1           1           1           1
   1           1           1           1           1           1           1           1           1
 1           1           1           1           1           1           1           1           1           1
           1           1           1           1           1           1           1           1           1
         1           1           1           1           1           1           1           1           1
       1           1           1           1           1           1           1           1           1
     1           1           1           1           1           1           1           1           1
   1           1           1           1           1           1           1           1           1
 0
(PID.TID 0000.0001)  OBCS_READPARMS: finished reading data.obcs
(PID.TID 0000.0001) CAL_READPARMS: opening data.cal
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.cal
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.cal"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *******************
(PID.TID 0000.0001) ># Calendar parameters
(PID.TID 0000.0001) ># *******************
(PID.TID 0000.0001) > &CAL_NML
(PID.TID 0000.0001) > TheCalendar='gregorian',
(PID.TID 0000.0001) > startdate_1=19580101
(PID.TID 0000.0001) > startdate_2=000000,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001) CAL_READPARMS: finished reading data.cal
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Calendar configuration >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) Calendar version: 0.2.0
(PID.TID 0000.0001)
(PID.TID 0000.0001) startTime =   /* Start time of the model integration [s] */
(PID.TID 0000.0001)                 0.000000000000000E+00
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) endTime =   /* End time of the model integration [s] */
(PID.TID 0000.0001)                 6.480000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) deltatclock =   /* Time interval for a model forward step [s] */
(PID.TID 0000.0001)                 2.160000000000000E+04
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingGregorianCalendar =  /* Calendar Type: Gregorian Calendar */
(PID.TID 0000.0001)                   T
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingJulianCalendar =  /* Calendar Type: Julian Calendar */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingModelCalendar =  /* Calendar Type: Model Calendar */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) usingNoCalendar =  /* Calendar Type: No Calendar */
(PID.TID 0000.0001)                   F
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) modelstartdate (YYYYMMDD) =   /* Model start date YYYY-MM-DD  */
(PID.TID 0000.0001)                19580101
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) modelstartdate (HHMMSS) =   /* Model start date HH-MM-SS  */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) modelenddate (YYYYMMDD) =   /* Model end date YYYY-MM-DD  */
(PID.TID 0000.0001)                19580101
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) modelenddate (HHMMSS) =   /* Model end date HH-MM-SS  */
(PID.TID 0000.0001)                  180000
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) intyears =   /* Number of calendar years affected by the integration  */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) intmonths =   /* Number of calendar months affected by the integration  */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) intdays =   /* Number of calendar days affected by the integration  */
(PID.TID 0000.0001)                       1
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nIter0 =   /* Base timestep number  */
(PID.TID 0000.0001)                       0
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nEndIter =   /* Final timestep number  */
(PID.TID 0000.0001)                       3
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001) nTimeSteps =   /* Number of model timesteps  */
(PID.TID 0000.0001)                       3
(PID.TID 0000.0001)     ;
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Calendar configuration  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) EXF_READPARMS: opening data.exf
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.exf
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.exf"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># $Header: /ross_sea/exp1_OnlyBath/res1/input/data.exf begun 15.9.2006 jns
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># External Forcing Data
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &EXF_NML
(PID.TID 0000.0001) > useExfYearlyFields = .TRUE.,
(PID.TID 0000.0001) > useExfCheckRange = .FALSE.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > atempstartdate1 = 19580101,
(PID.TID 0000.0001) > atempstartdate2 = 00000,
(PID.TID 0000.0001) > atempperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > aqhstartdate1 = 19580101,
(PID.TID 0000.0001) > aqhstartdate2 = 00000,
(PID.TID 0000.0001) > atempperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > precipstartdate1 = 19580101,
(PID.TID 0000.0001) > precipstartdate2 = 00000,
(PID.TID 0000.0001) > precipperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > uwindstartdate1 = 19580101,
(PID.TID 0000.0001) > uwindstartdate2 = 00000,
(PID.TID 0000.0001) > uwindperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > vwindstartdate1 = 19580101,
(PID.TID 0000.0001) > vwindstartdate2 = 00000,
(PID.TID 0000.0001) > vwindperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > swfluxstartdate1 = 19580101,
(PID.TID 0000.0001) > swfluxstartdate2 = 00000,
(PID.TID 0000.0001) > swfluxperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > lwfluxstartdate1 = 19580101,
(PID.TID 0000.0001) > lwfluxstartdate2 = 00000,
(PID.TID 0000.0001) > lwfluxperiod = 21600,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > obcsEstartdate1 = 19580101,
(PID.TID 0000.0001) > obcsEstartdate2 = 00000,
(PID.TID 0000.0001) > obcsEPeriod      = 2635200.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > obcsWstartdate1 = 19580101,
(PID.TID 0000.0001) > obcsWstartdate2 = 00000,
(PID.TID 0000.0001) > obcsWPeriod      = 2635200.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > exf_iprec        = 32,
(PID.TID 0000.0001) ># exf_yftype        = 'RL'
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > atempfile='dummy_2mTemp_1yr_b.bin',
(PID.TID 0000.0001) > aqhfile='dummy_spechum_1yr_b.bin',
(PID.TID 0000.0001) > uwindfile='dummy_uwind_1yr_b.bin',
(PID.TID 0000.0001) > vwindfile='dummy_vwind_1yr_b.bin',
(PID.TID 0000.0001) > evapfile='dummy_evap_1yr_b.bin',
(PID.TID 0000.0001) > precipfile='dummy_precip_1yr_b.bin',
(PID.TID 0000.0001) > lwdownfile='dummy_stdr_1yr_b.bin',
(PID.TID 0000.0001) > swdownfile='dummy_ssdr_1yr_b.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001) EXF_READPARMS: finished reading data.exf
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // External forcing configuration >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) External forcing version: 0.2.2
(PID.TID 0000.0001) Uses Calendar version:    0.2.0
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ALLOW_ATM_TEMP:                 NOT defined
(PID.TID 0000.0001) // ALLOW_ATM_WIND:                 NOT defined
(PID.TID 0000.0001) // ALLOW_DOWNWARD_RADIATION:       NOT defined
(PID.TID 0000.0001) // ALLOW_BULKFORMULAE:             NOT defined
(PID.TID 0000.0001)
(PID.TID 0000.0001)    Zonal wind stress forcing starts at                   0.
(PID.TID 0000.0001)    Zonal wind stress forcing period is                   0.
(PID.TID 0000.0001)    Zonal wind stress forcing is read from file:
(PID.TID 0000.0001)    >>    <<
(PID.TID 0000.0001)
(PID.TID 0000.0001)    Meridional wind stress forcing starts at              0.
(PID.TID 0000.0001)    Meridional wind stress forcing period is              0.
(PID.TID 0000.0001)    Meridional wind stress forcing is read from file:
(PID.TID 0000.0001)    >>    <<
(PID.TID 0000.0001)
(PID.TID 0000.0001)    Heat flux forcing starts at                          0.
(PID.TID 0000.0001)    Heat flux forcing period is                           0.
(PID.TID 0000.0001)    Heat flux forcing is read from file:
(PID.TID 0000.0001)    >>    <<
(PID.TID 0000.0001)
(PID.TID 0000.0001)    Salt flux forcing starts at                           0.
(PID.TID 0000.0001)    Salt flux forcing period is                           0.
(PID.TID 0000.0001)    Salt flux forcing is read from file:
(PID.TID 0000.0001)    >>    <<
(PID.TID 0000.0001)
(PID.TID 0000.0001)    Net shortwave flux forcing starts at                0.
(PID.TID 0000.0001)    Net shortwave flux forcing period is            21600.
(PID.TID 0000.0001)    Net shortwave flux forcing is read from file:
(PID.TID 0000.0001)    >>    <<
(PID.TID 0000.0001) // EXF_READ_EVAP:                  NOT defined
(PID.TID 0000.0001) // ALLOW_RUNOFF:                   NOT defined
(PID.TID 0000.0001)
(PID.TID 0000.0001)    Atmospheric pressure forcing starts at                0.
(PID.TID 0000.0001)    Atmospheric pressure forcing period is                0.
(PID.TID 0000.0001)    Atmospheric pressureforcing is read from file:
(PID.TID 0000.0001)    >>    <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // External forcing configuration  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) EXF_CLIM_READPARMS: opening data.exf_clim
(PID.TID 0000.0001)  OPEN_COPY_DATA_FILE: opening file data.exf_clim
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.exf_clim"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># Added 15.9.2006 jns - appears to be required for exf to run.
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *********************************
(PID.TID 0000.0001) ># External Forcing climatology data
(PID.TID 0000.0001) ># *********************************
(PID.TID 0000.0001) > &EXF_CLIM_NML
(PID.TID 0000.0001) > climtempstartdate1 = 19580101,
(PID.TID 0000.0001) > climtempstartdate2 = 00000,
(PID.TID 0000.0001) > climtempperiod     = 2635200.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > climsaltstartdate1 = 19580101,
(PID.TID 0000.0001) > climsaltstartdate2 = 00000,
(PID.TID 0000.0001) > climsaltperiod     = 2635200.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > climtempfile       = 'dummy_theta',
(PID.TID 0000.0001) > climsaltfile       = 'dummy_salinity',
(PID.TID 0000.0001) > climsstfile        = ' ',
(PID.TID 0000.0001) > climsssfile        = ' ',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > exf_clim_iprec     = 32,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001) EXF_CLIM_READPARMS: finished reading data.exf_clim
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // External forcing climatology configuration >>> START <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) External forcing version: 0.2.2
(PID.TID 0000.0001) Uses Calendar version:    0.2.0
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ALLOW_CLIMTEMP_RELAXATION:      NOT defined
(PID.TID 0000.0001) // ALLOW_CLIMSALT_RELAXATION:      NOT defined
(PID.TID 0000.0001) // ALLOW_CLIMSST_RELAXATION:       NOT defined
(PID.TID 0000.0001) // ALLOW_CLIMSSS_RELAXATION:       NOT defined
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // External forcing climatology configuration  >>> END <<<
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) S/R LOAD_GRID_SPACING: delX loaded from file: DELX3.ross_1deg_sn_ciso
(PID.TID 0000.0001) S/R LOAD_GRID_SPACING: delY loaded from file: DELY3.ross_1deg_sn_ciso
(PID.TID 0000.0001)  MDSREADFIELD: opening global file: bathy_etgebda_sn_iso1deg_b.bin
(PID.TID 0000.0001) // ================



More information about the MITgcm-support mailing list