[MITgcm-support] open boundary condition cost
Suneet Dwivedi
suneetdwivedi at gmail.com
Mon Mar 17 19:52:00 EDT 2008
Hi Patrick,
Thanks for your reply. As per your suggestion, I tried to find out the
reason for the error by looking into the data files but I didn't able
to figure it out. Please find attached herewith the output.txt file
for your perusal. Can you please suggest me what may be the error?
Without obcs defined in ECCO_CPPOPTIONS.h, the model is running fine.
The problem occurs as soon as I try to define either/both of
ALLOW_OBCSN_CONTROL and ALLOW_OBCSN_COST_CONTRIBUTION. Please help.
Regards,
Suneet
On Fri, Mar 14, 2008 at 4:15 AM, Patrick Heimbach <heimbach at mit.edu> wrote:
>
> Hi Suneet,
>
> what's happening is that either during the read of the
> OB control variable (xx_obc...) or the read of the
> OB cost function, you are trying to read a non-existing record
> (not clear from the output provided where the model stops).
> This could have three reasons,
> 1. not enough records in your cost file
> 2. not enough records in your control file (less likely)
> 3. mismatch in file precision for cost file (e.g.
> you provide 32bit file, but try to read 64bit file)
> Have you checked and properly set open boundary values in
> data.ctrl (for OB control) and data.ecco (for OB cost)?
> You could try to first only do
> #define ALLOW_OBCSN_CONTROL
> #undef ALLOW_OBCSN_COST_CONTRIBUTION
> i.e. don't use the cost function penalty term.
> I think tough the code wants both (at least it used to),
> and complains if you don't have the cost term, so you need to
> temporarily comment the "complaint".
> I (embarrassingly) forgot which routine,
> but should be easy to find once you get to that point.
>
> -Patrick
>
>
>
>
>
> On Mar 13, 2008, at 2:10 PM, Suneet Dwivedi wrote:
>
> > Hi Everyone,
> > When I am trying to calculate the model cost corresponding to open
> > boundary conditions by defining them in ECCO_CPPOPTIONS.h (define
> > ALLOW_OBCSN_COST_CONTRIBUTION; define ALLOW_OBCSN_CONTROL), I am
> > getting the following errors:
> > ----------------------------------------------------------------------
> > ------------------------------------------------
> > forrtl: severe (24): end-of-file during read, unit 11, file fortqtPIW1
> > Image PC Routine Line
> > Source
> > libintlc.so.5 0000002A95D7A05E Unknown Unknown
> > Unknown
> > forrtl: severe (24): end-of-file during read, unit 11, file fortC2Du11
> > Image PC Routine Line
> > Source
> > libintlc.so.5 0000002A95D7A05E Unknown Unknown
> > Unknown
> > ----------------------------------------------------------------------
> > ------------------------------------------------
> >
> > If I don't make these changes in ECCO_CPPOPTIONS.h, then the model
> > runs fine and open boundary conditions show an effect on cost of
> > others variables but as soon as I define them this problem occurs. Can
> > anyone of you please tell me what is the meaning of this error and how
> > to tackle this problem.
> > Regards,
> > Suneet
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> ---
> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
(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: checkpoint59g
(PID.TID 0000.0001) // Build user: mitgcm
(PID.TID 0000.0001) // Build host: dolphin.eps.jhu.edu
(PID.TID 0000.0001) // Build date: Mon Mar 17 19:22:59 EDT 2008
(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 = 20 ; /* Tile size in X */
(PID.TID 0000.0001) sNy = 10 ; /* Tile size in Y */
(PID.TID 0000.0001) OLx = 4 ; /* Tile overlap distance in X */
(PID.TID 0000.0001) OLy = 4 ; /* 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 = 16 ; /* No. levels in the vertical */
(PID.TID 0000.0001) nX = 20 ; /* Total domain size in X ( = nPx*nSx*sNx ) */
(PID.TID 0000.0001) nY = 10 ; /* 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) ># ******************************
(PID.TID 0000.0001) ># Model parameters
(PID.TID 0000.0001) ># Continuous equation parameters
(PID.TID 0000.0001) ># ******************************
(PID.TID 0000.0001) > &PARM01
(PID.TID 0000.0001) ># tRef= 24.0 , 23.0 , 22.0 , 21.0 , 20.0 ,
(PID.TID 0000.0001) ># 19.0 , 18.0 , 17.0 , 16.0 , 15.0 ,
(PID.TID 0000.0001) ># 14.0 , 13.0 , 12.0 , 11.0 , 10.0 ,
(PID.TID 0000.0001) ># 9.0 , 8.0 , 7.0 , 6.0, 5.0 ,
(PID.TID 0000.0001) ># sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
(PID.TID 0000.0001) ># 34.90, 34.86, 34.78, 34.69, 34.60,
(PID.TID 0000.0001) ># 34.58, 34.62, 34.68, 34.72, 34.73,
(PID.TID 0000.0001) ># 34.74, 34.73, 34.73, 34.72, 34.72,
(PID.TID 0000.0001) > viscAz =1.E-3,
(PID.TID 0000.0001) > viscA4 =3.125E9,
(PID.TID 0000.0001) > diffK4T=3.125E9,
(PID.TID 0000.0001) > diffKzT=1.E-5,
(PID.TID 0000.0001) > diffK4S=3.125E9,
(PID.TID 0000.0001) > diffKzS=1.E-5,
(PID.TID 0000.0001) > no_slip_sides=.FALSE.,
(PID.TID 0000.0001) > no_slip_bottom=.TRUE.,
(PID.TID 0000.0001) > bottomDragQuadratic = 0.001,
(PID.TID 0000.0001) > rigidLid=.FALSE.,
(PID.TID 0000.0001) > implicitFreeSurface=.TRUE.,
(PID.TID 0000.0001) > eosType='JMD95Z',
(PID.TID 0000.0001) > readBinaryPrec=64,
(PID.TID 0000.0001) > writeBinaryPrec=64,
(PID.TID 0000.0001) > useSingleCpuIo=.FALSE.,
(PID.TID 0000.0001) > hFacMin=0.1,
(PID.TID 0000.0001) > hFacMinDz=5.,
(PID.TID 0000.0001) > allowFreezing=.FALSE.,
(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) > implicitDiffusion=.TRUE.,
(PID.TID 0000.0001) > implicitViscosity=.TRUE.,
(PID.TID 0000.0001) > vectorInvariantMomentum=.TRUE.,
(PID.TID 0000.0001) > staggerTimeStep=.TRUE.,
(PID.TID 0000.0001) > tempAdvScheme=30,
(PID.TID 0000.0001) > saltAdvScheme=30,
(PID.TID 0000.0001) > MultiDimAdvection=.FALSE.,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># **************************
(PID.TID 0000.0001) ># Elliptic solver parameters
(PID.TID 0000.0001) ># **************************
(PID.TID 0000.0001) > &PARM02
(PID.TID 0000.0001) > cg2dMaxIters=1000,
(PID.TID 0000.0001) > cg2dTargetResidual=1.E-13,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># ************************
(PID.TID 0000.0001) ># Time stepping parameters
(PID.TID 0000.0001) ># ************************
(PID.TID 0000.0001) > &PARM03
(PID.TID 0000.0001) > nIter0=0,
(PID.TID 0000.0001) ># endTime=72000.,
(PID.TID 0000.0001) > endTime=14400.,
(PID.TID 0000.0001) ># endTime=28800.,
(PID.TID 0000.0001) > deltaT=600.,
(PID.TID 0000.0001) > tauCD=0.,
(PID.TID 0000.0001) > cAdjFreq=0.,
(PID.TID 0000.0001) > abEps=0.1,
(PID.TID 0000.0001) > forcing_In_AB=.FALSE.,
(PID.TID 0000.0001) > pChkptFreq=2592000.,
(PID.TID 0000.0001) > dumpFreq=43200.,
(PID.TID 0000.0001) > diagFreq=86400.,
(PID.TID 0000.0001) > monitorFreq=12000.,
(PID.TID 0000.0001) > tauThetaClimRelax=0.,
(PID.TID 0000.0001) > tauSaltClimRelax=0.,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># *******************
(PID.TID 0000.0001) ># Gridding parameters
(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) > delZ=50,50,50,50,50,50,50,50,50,50,50,50,100,100,100,100,
(PID.TID 0000.0001) > phiMin=-57,
(PID.TID 0000.0001) > thetaMin=-50,
(PID.TID 0000.0001) > delY=10*1.00e+00,
(PID.TID 0000.0001) > delX=20*1.00e+00,
(PID.TID 0000.0001) > rSphere = 6371.D3,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># **********
(PID.TID 0000.0001) ># Data Files
(PID.TID 0000.0001) ># **********
(PID.TID 0000.0001) > &PARM05
(PID.TID 0000.0001) > bathyFile= 'so_bathymetry.bin',
(PID.TID 0000.0001) > hydrogThetaFile= 'Clim_TEMP_lev.bin',
(PID.TID 0000.0001) > hydrogSaltFile = 'Clim_SALT_lev.bin',
(PID.TID 0000.0001) > uVelInitFile = 'ncep_6hrly_uwind_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > vVelInitFile = 'ncep_6hrly_vwind_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) ># pSurfInitFile,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># Open boundaries
(PID.TID 0000.0001) > &PARM06
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># KPP Mixing scheme
(PID.TID 0000.0001) > &PARM07
(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) >#
(PID.TID 0000.0001) ># ********
(PID.TID 0000.0001) ># Packages
(PID.TID 0000.0001) ># ********
(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) > useEXF = .TRUE.,
(PID.TID 0000.0001) > useECCO = .TRUE.,
(PID.TID 0000.0001) > usePTRACERS = .FALSE.,
(PID.TID 0000.0001) > useOBCS = .TRUE.,
(PID.TID 0000.0001) > useGrdchk = .FALSE.,
(PID.TID 0000.0001) > useDiagnostics = .FALSE.,
(PID.TID 0000.0001) > useMNC = .FALSE.,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
(PID.TID 0000.0001) PACKAGES_BOOT: finished reading data.pkg
(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) >#
(PID.TID 0000.0001) > &CAL_NML
(PID.TID 0000.0001) > TheCalendar = 'gregorian',
(PID.TID 0000.0001) > startDate_1 = 20030601,
(PID.TID 0000.0001) > startDate_2 = 00000,
(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) 1.440000000000000E+04
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) deltatclock = /* Time interval for a model forward step [s] */
(PID.TID 0000.0001) 6.000000000000000E+02
(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) 20030601
(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) 20030601
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) modelenddate (HHMMSS) = /* Model end date HH-MM-SS */
(PID.TID 0000.0001) 40000
(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) 24
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) nTimeSteps = /* Number of model timesteps */
(PID.TID 0000.0001) 24
(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) 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) ># GM+Redi package parameters:
(PID.TID 0000.0001) ># GM_Small_Number :: epsilon used in computing the slope
(PID.TID 0000.0001) ># GM_slopeSqCutoff :: slope^2 cut-off value
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#-from MOM :
(PID.TID 0000.0001) ># GM_background_K: G & Mc.W diffusion coefficient
(PID.TID 0000.0001) ># GM_maxSlope : max slope of isopycnals
(PID.TID 0000.0001) ># GM_Scrit : transition for scaling diffusion coefficient
(PID.TID 0000.0001) ># GM_Sd : half width scaling for diffusion coefficient
(PID.TID 0000.0001) ># GM_taper_scheme: slope clipping or one of the tapering schemes
(PID.TID 0000.0001) ># GM_Kmin_horiz : horizontal diffusion minimum value
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#-Option parameters (needs to "define" options in GMREDI_OPTIONS.h")
(PID.TID 0000.0001) ># GM_isopycK : isopycnal diffusion coefficient (default=GM_background_K)
(PID.TID 0000.0001) ># GM_AdvForm : turn on GM Advective form (default=Skew flux form)
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &GM_PARM01
(PID.TID 0000.0001) > GM_Small_Number = 1.D-20,
(PID.TID 0000.0001) > GM_slopeSqCutoff = 1.D+08,
(PID.TID 0000.0001) > GM_AdvForm = .FALSE.,
(PID.TID 0000.0001) ># GM_isopycK = 1.1D+3,
(PID.TID 0000.0001) ># GM_background_K = 0.9D+3,
(PID.TID 0000.0001) > GM_background_K = 1.D+3,
(PID.TID 0000.0001) > GM_taper_scheme = 'dm95',
(PID.TID 0000.0001) > GM_maxSlope = 1.D-2,
(PID.TID 0000.0001) > GM_Kmin_horiz = 50.,
(PID.TID 0000.0001) > GM_Scrit = 4.D-3,
(PID.TID 0000.0001) > GM_Sd = 1.D-3,
(PID.TID 0000.0001) ># GM_Visbeck_alpha = 1.5D-2,
(PID.TID 0000.0001) > GM_Visbeck_alpha = 0.,
(PID.TID 0000.0001) > GM_Visbeck_length = 2.D+5,
(PID.TID 0000.0001) > GM_Visbeck_depth = 1.D+3,
(PID.TID 0000.0001) > GM_Visbeck_maxval_K= 2.5D+3,
(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) ># ***************
(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 = 864000.,
(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) >#
(PID.TID 0000.0001) ># ***************
(PID.TID 0000.0001) ># Open boundaries
(PID.TID 0000.0001) ># ***************
(PID.TID 0000.0001) > &OBCS_PARM01
(PID.TID 0000.0001) > OB_Jnorth=20*10,
(PID.TID 0000.0001) > OB_Jsouth=20*10,
(PID.TID 0000.0001) > OB_Ieast=10*20,
(PID.TID 0000.0001) > OB_Iwest=10*20,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > useOrlanskiNorth=.FALSE.
(PID.TID 0000.0001) > useOrlanskiSouth=.FALSE.
(PID.TID 0000.0001) > useOrlanskiEast =.FALSE.
(PID.TID 0000.0001) > useOrlanskiWest =.FALSE.
(PID.TID 0000.0001) > useOBCSsponge =.FALSE.
(PID.TID 0000.0001) > useOBCSbalance =.FALSE.
(PID.TID 0000.0001) > useOBCSprescribe=.TRUE.
(PID.TID 0000.0001) > OBCSprintDiags =.FALSE.
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > OBNuFile='U_N.bin',
(PID.TID 0000.0001) > OBNvFile='U_N.bin',
(PID.TID 0000.0001) > OBNtFile='T_N.bin',
(PID.TID 0000.0001) > OBNsFile='S_N.bin',
(PID.TID 0000.0001) > OBSuFile='U_N.bin',
(PID.TID 0000.0001) > OBSvFile='U_N.bin',
(PID.TID 0000.0001) > OBStFile='T_N.bin',
(PID.TID 0000.0001) > OBSsFile='S_N.bin',
(PID.TID 0000.0001) > OBEuFile='U_E.bin',
(PID.TID 0000.0001) > OBEvFile='U_E.bin',
(PID.TID 0000.0001) > OBEtFile='T_E.bin',
(PID.TID 0000.0001) > OBEsFile='S_E.bin',
(PID.TID 0000.0001) > OBWuFile='U_E.bin',
(PID.TID 0000.0001) > OBWvFile='U_E.bin',
(PID.TID 0000.0001) > OBWtFile='T_E.bin',
(PID.TID 0000.0001) > OBWsFile='S_E.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *****************************************
(PID.TID 0000.0001) ># Orlanski Boundary Condition Parameters.
(PID.TID 0000.0001) ># *****************************************
(PID.TID 0000.0001) ># &OBCS_PARM02
(PID.TID 0000.0001) ># cvelTimeScale = 2000.,
(PID.TID 0000.0001) ># *****************************************
(PID.TID 0000.0001) ># Sponge Layer Parameters.
(PID.TID 0000.0001) ># *****************************************
(PID.TID 0000.0001) ># &OBCS_PARM03
(PID.TID 0000.0001) ># Urelaxobcsinner=7776000.E0,
(PID.TID 0000.0001) ># Urelaxobcsbound=2592000.E0,
(PID.TID 0000.0001) ># Vrelaxobcsinner=7776000.E0,
(PID.TID 0000.0001) ># Vrelaxobcsbound=2592000.E0,
(PID.TID 0000.0001) ># spongeThickness=10,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
OB Jn = 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
OB Js = 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
OB Ie = 20 20 20 20 20 20 20 20 20 20
OB Iw = 20 20 20 20 20 20 20 20 20 20
(PID.TID 0000.0001) OBCS_READPARMS: finished reading data.obcs
(PID.TID 0000.0001) OPTIM_READPARMS: opening data.optim
(PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.optim
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.optim"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ********************************
(PID.TID 0000.0001) ># Off-line optimization parameters
(PID.TID 0000.0001) ># ********************************
(PID.TID 0000.0001) > &OPTIM
(PID.TID 0000.0001) > optimcycle=0,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &ECCO_OPTIM
(PID.TID 0000.0001) > optimcycle=0,
(PID.TID 0000.0001) > numiter=1,
(PID.TID 0000.0001) > nfunc=5,
(PID.TID 0000.0001) > fmin=18300.,
(PID.TID 0000.0001) > iprint=10,
(PID.TID 0000.0001) > nupdate=8,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
(PID.TID 0000.0001) OPTIM_READPARMS: finished reading data.optim
(PID.TID 0000.0001) CTRL_READPARMS: opening data.ctrl
(PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.ctrl
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.ctrl"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># ECCO controlvariables
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &CTRL_NML
(PID.TID 0000.0001) > xx_theta_file ='xx_theta',
(PID.TID 0000.0001) > xx_salt_file ='xx_salt',
(PID.TID 0000.0001) > xx_uvel_file ='xx_uvel',
(PID.TID 0000.0001) > xx_vvel_file ='xx_vvel',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># xx_hfluxstartdate1 = 20030601,
(PID.TID 0000.0001) ># xx_hfluxstartdate2 = 220000,
(PID.TID 0000.0001) ># xx_hfluxperiod = 864000.0,
(PID.TID 0000.0001) ># xx_hflux_file = 'xx_hfl',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># xx_sfluxstartdate1 = 20030601,
(PID.TID 0000.0001) ># xx_sfluxstartdate2 = 220000,
(PID.TID 0000.0001) ># xx_sfluxperiod = 864000.0,
(PID.TID 0000.0001) ># xx_sflux_file = 'xx_sfl',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># xx_tauustartdate1 = 20030601,
(PID.TID 0000.0001) ># xx_tauustartdate2 = 220000,
(PID.TID 0000.0001) ># xx_tauuperiod = 864000.0,
(PID.TID 0000.0001) ># xx_tauu_file = 'xx_tauu',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># xx_tauvstartdate1 = 20030601,
(PID.TID 0000.0001) ># xx_tauvstartdate2 = 220000,
(PID.TID 0000.0001) ># xx_tauvperiod = 864000.0,
(PID.TID 0000.0001) ># xx_tauv_file = 'xx_tauv',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_atempstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_atempstartdate2 = 220000,
(PID.TID 0000.0001) > xx_atempperiod = 864000.0,
(PID.TID 0000.0001) > xx_atemp_file = 'xx_atemp',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_aqhstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_aqhstartdate2 = 220000,
(PID.TID 0000.0001) > xx_aqhperiod = 864000.0,
(PID.TID 0000.0001) > xx_aqh_file = 'xx_aqh',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_precipstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_precipstartdate2 = 220000,
(PID.TID 0000.0001) > xx_precipperiod = 864000.0,
(PID.TID 0000.0001) > xx_precip_file = 'xx_precip',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_swfluxstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_swfluxstartdate2 = 220000,
(PID.TID 0000.0001) > xx_swfluxperiod = 864000.0,
(PID.TID 0000.0001) > xx_swflux_file = 'xx_swflux',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># xx_swdownstartdate1 = 20030601,
(PID.TID 0000.0001) ># xx_swdownstartdate2 = 220000,
(PID.TID 0000.0001) ># xx_swdownperiod = 864000.0,
(PID.TID 0000.0001) ># xx_swdown_file = 'xx_swdown',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_uwindstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_uwindstartdate2 = 220000,
(PID.TID 0000.0001) > xx_uwindperiod = 864000.0,
(PID.TID 0000.0001) > xx_uwind_file = 'xx_uwind',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_vwindstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_vwindstartdate2 = 220000,
(PID.TID 0000.0001) > xx_vwindperiod = 864000.0,
(PID.TID 0000.0001) > xx_vwind_file = 'xx_vwind',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_obcsnstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_obcsnstartdate2 = 220000,
(PID.TID 0000.0001) > xx_obcsnperiod = 864000.0,
(PID.TID 0000.0001) > xx_obcsn_file = 'xx_obcsn'
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_obcssstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_obcssstartdate2 = 220000,
(PID.TID 0000.0001) > xx_obcssperiod = 864000.0,
(PID.TID 0000.0001) > xx_obcss_file = 'xx_obcss'
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_obcswstartdate1 = 20030601,
(PID.TID 0000.0001) > xx_obcswstartdate2 = 220000,
(PID.TID 0000.0001) > xx_obcswperiod = 864000.0,
(PID.TID 0000.0001) > xx_obcsw_file = 'xx_obcsw'
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > xx_obcsestartdate1 = 20030601,
(PID.TID 0000.0001) > xx_obcsestartdate2 = 220000,
(PID.TID 0000.0001) > xx_obcseperiod = 864000.0,
(PID.TID 0000.0001) > xx_obcse_file = 'xx_obcse'
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) ># names for ctrl_pack/unpack
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &CTRL_PACKNAMES
(PID.TID 0000.0001) > yadmark = 'ad',
(PID.TID 0000.0001) > ctrlname = 'JHU_ctrl',
(PID.TID 0000.0001) > costname = 'JHU_cost',
(PID.TID 0000.0001) > scalname = 'JHU_scal',
(PID.TID 0000.0001) > maskname = 'JHU_mask',
(PID.TID 0000.0001) > metaname = 'JHU_meta',
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
(PID.TID 0000.0001) CTRL_READPARMS: finished reading data.ctrl
(PID.TID 0000.0001) COST_READPARMS: opening data.cost
(PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.cost
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.cost"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ******************
(PID.TID 0000.0001) ># cost function
(PID.TID 0000.0001) ># ******************
(PID.TID 0000.0001) > &COST_NML
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
(PID.TID 0000.0001) COST_READPARMS: finished reading data.cost
(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: /u/gcmpack/MITgcm/verification/lab_sea/input/data.exf,v 1.11 2007/04/18 19:57:21 heimbach Exp $
(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_01
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > useExfCheckRange = .FALSE.,
(PID.TID 0000.0001) > repeatPeriod = 31622400.0,
(PID.TID 0000.0001) > exf_iprec = 64,
(PID.TID 0000.0001) > exf_yftype = 'RL',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) >&EXF_NML_02
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > hfluxstartdate1 = 20030601,
(PID.TID 0000.0001) > hfluxstartdate2 = 000000,
(PID.TID 0000.0001) > hfluxperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > atempstartdate1 = 20030601,
(PID.TID 0000.0001) > atempstartdate2 = 000000,
(PID.TID 0000.0001) > atempperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > aqhstartdate1 = 20030601,
(PID.TID 0000.0001) > aqhstartdate2 = 000000,
(PID.TID 0000.0001) > aqhperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > evapstartdate1 = 20030601,
(PID.TID 0000.0001) > evapstartdate2 = 000000,
(PID.TID 0000.0001) > evapperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > precipstartdate1 = 20030601,
(PID.TID 0000.0001) > precipstartdate2 = 000000,
(PID.TID 0000.0001) > precipperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > uwindstartdate1 = 20030601,
(PID.TID 0000.0001) > uwindstartdate2 = 000000,
(PID.TID 0000.0001) > uwindperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > vwindstartdate1 = 20030601,
(PID.TID 0000.0001) > vwindstartdate2 = 000000,
(PID.TID 0000.0001) > vwindperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ustressstartdate1 = 20030601,
(PID.TID 0000.0001) ># ustressstartdate2 = 000000,
(PID.TID 0000.0001) ># ustressperiod = 21600.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># vstressstartdate1 = 20030601,
(PID.TID 0000.0001) ># vstressstartdate2 = 000000,
(PID.TID 0000.0001) ># vstressperiod = 21600.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > swdownstartdate1 = 20030601,
(PID.TID 0000.0001) > swdownstartdate2 = 000000,
(PID.TID 0000.0001) > swdownperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > lwdownstartdate1 = 20030601,
(PID.TID 0000.0001) > lwdownstartdate2 = 000000,
(PID.TID 0000.0001) > lwdownperiod = 21600.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># hfluxfile = ' ',
(PID.TID 0000.0001) > atempfile = 'ncep_6hrly_atemp_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > aqhfile = 'ncep_6hrly_rhum_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > uwindfile = 'SeaWinds_6hrly_uwind_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > vwindfile = 'SeaWinds_6hrly_vwind_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) ># ustressfile = 'ncep_6hrly_tauu_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) ># vstressfile = 'ncep_6hrly_tauv_1Jun03_to_Aug03.stretched.bin ',
(PID.TID 0000.0001) > evapfile = 'ncep_6hrly_evap_latenthf_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > precipfile = 'ncep_6hrly_precip_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > lwdownfile = 'ncep_6hrly_dw_lw2_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) > swdownfile = 'ncep_6hrly_dw_sw2_1Jun03_to_Aug03.stretched.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &EXF_NML_03
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &EXF_NML_04
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) ># *********************
(PID.TID 0000.0001) > &EXF_NML_OBCS
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > obcsNstartdate1 = 20030601,
(PID.TID 0000.0001) > obcsNstartdate2 = 00000,
(PID.TID 0000.0001) > obcsNperiod = 21600.0,
(PID.TID 0000.0001) ># obcsNperiod = 5637600.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > obcsSstartdate1 = 20030601,
(PID.TID 0000.0001) > obcsSstartdate2 = 00000,
(PID.TID 0000.0001) > obcsSperiod = 21600.0,
(PID.TID 0000.0001) ># obcsSperiod = 5637600.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > obcsEstartdate1 = 20030601,
(PID.TID 0000.0001) > obcsEstartdate2 = 00000,
(PID.TID 0000.0001) > obcsEperiod = 21600.0,
(PID.TID 0000.0001) ># obcsEperiod = 5637600.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > obcsWstartdate1 = 20030601,
(PID.TID 0000.0001) > obcsWstartdate2 = 00000,
(PID.TID 0000.0001) > obcsWperiod = 21600.0,
(PID.TID 0000.0001) ># obcsWperiod = 5637600.0,
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
(PID.TID 0000.0001) EXF_READPARMS: reading EXF_NML_01
(PID.TID 0000.0001) EXF_READPARMS: reading EXF_NML_02
(PID.TID 0000.0001) EXF_READPARMS: reading EXF_NML_03
(PID.TID 0000.0001) EXF_READPARMS: reading EXF_NML_04
(PID.TID 0000.0001) EXF_READPARMS: reading EXF_NML_OBCS
(PID.TID 0000.0001) EXF_READPARMS: finished reading data.exf
(PID.TID 0000.0001) ECCO_READPARMS: opening data.ecco
(PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.ecco
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.ecco"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ******************
(PID.TID 0000.0001) ># ECCO cost function
(PID.TID 0000.0001) ># ******************
(PID.TID 0000.0001) > &ECCO_COST_NML
(PID.TID 0000.0001) > data_errfile = 'data.err',
(PID.TID 0000.0001) > tbarfile = 'tbar',
(PID.TID 0000.0001) > sbarfile = 'sbar',
(PID.TID 0000.0001) > psbarfile = 'psbar',
(PID.TID 0000.0001) > ubarfile = 'ubar',
(PID.TID 0000.0001) > vbarfile = 'vbar',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > using_topex = .TRUE.,
(PID.TID 0000.0001) > topexstartdate1 = 20030601,
(PID.TID 0000.0001) > topexstartdate2 = 00000,
(PID.TID 0000.0001) > topexperiod = 86400.0,
(PID.TID 0000.0001) ># topexfile = 'SSH_TPJ.bin',
(PID.TID 0000.0001) > topexfile = 'SSH_ERS.bin',
(PID.TID 0000.0001) ># topexmeanfile = 'mit9yrmn_TP_mean_s',
(PID.TID 0000.0001) ># geoid_errfile ???? Don't know how this file is generated !!!!
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > using_ers = .TRUE.,
(PID.TID 0000.0001) > ersstartdate1 = 20030601,
(PID.TID 0000.0001) > ersstartdate2 = 00000,
(PID.TID 0000.0001) > ersperiod = 86400.0,
(PID.TID 0000.0001) > ersfile = 'SSH_ERS.bin',
(PID.TID 0000.0001) > ssh_errfile = 'ssh_error.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > sststartdate1 = 20030601,
(PID.TID 0000.0001) > sststartdate2 = 00000,
(PID.TID 0000.0001) > sstdatfile = 'SST.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># tmistartdate1 = 20030601,
(PID.TID 0000.0001) ># tmistartdate2 = 00000,
(PID.TID 0000.0001) ># tmidatfile = 'SST_TMI.bin',
(PID.TID 0000.0001) ># tmidatfile = ' ',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > tdatfile = 'PTMP_lev.bin',
(PID.TID 0000.0001) > sdatfile = 'SALT_lev.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># ctdtfile = 'PTMP_CTD.bin',
(PID.TID 0000.0001) ># ctdsfile = 'SALT_CTD.bin',
(PID.TID 0000.0001) ># xbtfile = 'PTMP_XBT.bin',
(PID.TID 0000.0001) ># argotfile = 'PTMP_ARGO.bin',
(PID.TID 0000.0001) ># argosfile = 'SALT_ARGO.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > temp0errfile = 'theta_error.bin',
(PID.TID 0000.0001) > salt0errfile = 'salt_error.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># temperrfile = ' ',
(PID.TID 0000.0001) ># salterrfile = ' ',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) >## hflux_errfile = 'NCEP_R6_turbflux_sd_1992-2003',
(PID.TID 0000.0001) >## sflux_errfile = 'NCEP_R6_emp_sd_1992-2003',
(PID.TID 0000.0001) >## tauu_errfile = 'NCEP_R6_ustr_anom_sd_1992-2003',
(PID.TID 0000.0001) >## tauv_errfile = 'NCEP_R6_vstr_anom_sd_1992-2003',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># geoid_errfile = ' ',
(PID.TID 0000.0001) ># geoid_covariancefile = ' ',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># atemp_errfile = ' ',
(PID.TID 0000.0001) ># aqh_errfile = ' ',
(PID.TID 0000.0001) ># precip_errfile = ' ',
(PID.TID 0000.0001) ># swflux_errfile = ' ',
(PID.TID 0000.0001) > uwind_errfile = 'uv_error.bin',
(PID.TID 0000.0001) > vwind_errfile = 'uv_error.bin',
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > whflux0 = 20.,
(PID.TID 0000.0001) ># whflux0 = 1.E-8,
(PID.TID 0000.0001) > wsflux0 = 3.0E-8,
(PID.TID 0000.0001) > wtau0 = 2.0E-2,
(PID.TID 0000.0001) > watemp0 = 0.5,
(PID.TID 0000.0001) > waqh0 = 5.E-4,
(PID.TID 0000.0001) > wprecip0 = 1.E-8,
(PID.TID 0000.0001) > wswflux0 = 20.,
(PID.TID 0000.0001) > wwind0 = 1.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > wmean_hflux = 30.,
(PID.TID 0000.0001) ># wmean_hflux = 1.E-8,
(PID.TID 0000.0001) > wmean_sflux = 1.6E-8,
(PID.TID 0000.0001) > wmean_tau = 0.1,
(PID.TID 0000.0001) > wmean_atemp = 1.,
(PID.TID 0000.0001) > wmean_aqh = 1.E-3,
(PID.TID 0000.0001) > wmean_precip = 1.5E-8,
(PID.TID 0000.0001) > wmean_swflux = 20.,
(PID.TID 0000.0001) > wmean_wind = 2.0,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># mult_hmean = 1.,
(PID.TID 0000.0001) > mult_h = 1.,
(PID.TID 0000.0001) > mult_tp = 1.,
(PID.TID 0000.0001) > mult_ers = 1.,
(PID.TID 0000.0001) ># mult_gfo = 1.,
(PID.TID 0000.0001) > mult_temp = 1.,
(PID.TID 0000.0001) > mult_salt = 1.,
(PID.TID 0000.0001) > mult_temp0 = 1.,
(PID.TID 0000.0001) > mult_salt0 = 1.,
(PID.TID 0000.0001) > mult_u0 = 1.,
(PID.TID 0000.0001) > mult_v0 = 1.,
(PID.TID 0000.0001) > mult_sst = 1.,
(PID.TID 0000.0001) > mult_hflux = 1.,
(PID.TID 0000.0001) > mult_sflux = 1.,
(PID.TID 0000.0001) > mult_tauu = 1.,
(PID.TID 0000.0001) > mult_tauv = 1.,
(PID.TID 0000.0001) > mult_atemp = 1.,
(PID.TID 0000.0001) > mult_aqh = 1.,
(PID.TID 0000.0001) > mult_precip= 1.,
(PID.TID 0000.0001) > mult_swflux= 1.,
(PID.TID 0000.0001) > mult_uwind = 1.,
(PID.TID 0000.0001) > mult_vwind = 1.,
(PID.TID 0000.0001) > mult_obcsn = 1.,
(PID.TID 0000.0001) > mult_obcss = 1.,
(PID.TID 0000.0001) > mult_obcse = 1.,
(PID.TID 0000.0001) > mult_obcsw = 1.,
(PID.TID 0000.0001) > mult_obcsvol=1.,
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > cost_iprec = 64,
(PID.TID 0000.0001) > cost_yftype = 'RL',
(PID.TID 0000.0001) > &
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) > &ECCO_PARMS
(PID.TID 0000.0001) > expId = 'MIT_CE_000',
(PID.TID 0000.0001) > &
(PID.TID 0000.0001)
(PID.TID 0000.0001) ECCO_READPARMS: finished reading data.ecco
(PID.TID 0000.0001) ECCO_READPARMS: start assigning cost dates
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) ECCO_READPARMS: end assigning cost dates
(PID.TID 0000.0001) SET_PARMS: done
(PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr= T ; setCenterDr= F
(PID.TID 0000.0001) %MON XC_max = -3.0500000000000E+01
(PID.TID 0000.0001) %MON XC_min = -4.9500000000000E+01
(PID.TID 0000.0001) %MON XC_mean = -4.0000000000000E+01
(PID.TID 0000.0001) %MON XC_sd = 5.7662812973354E+00
(PID.TID 0000.0001) %MON XG_max = -3.1000000000000E+01
(PID.TID 0000.0001) %MON XG_min = -5.0000000000000E+01
(PID.TID 0000.0001) %MON XG_mean = -4.0500000000000E+01
(PID.TID 0000.0001) %MON XG_sd = 5.7662812973354E+00
(PID.TID 0000.0001) %MON DXC_max = 7.5122203577606E+04
(PID.TID 0000.0001) %MON DXC_min = 6.1372592594193E+04
(PID.TID 0000.0001) %MON DXC_mean = 6.8372443343038E+04
(PID.TID 0000.0001) %MON DXC_sd = 4.3899851261152E+03
(PID.TID 0000.0001) %MON DXF_max = 7.5122203577606E+04
(PID.TID 0000.0001) %MON DXF_min = 6.1372592594193E+04
(PID.TID 0000.0001) %MON DXF_mean = 6.8372443343038E+04
(PID.TID 0000.0001) %MON DXF_sd = 4.3899851261152E+03
(PID.TID 0000.0001) %MON DXG_max = 7.4403928689702E+04
(PID.TID 0000.0001) %MON DXG_min = 6.0561097546259E+04
(PID.TID 0000.0001) %MON DXG_mean = 6.7606156926271E+04
(PID.TID 0000.0001) %MON DXG_sd = 4.4197252361186E+03
(PID.TID 0000.0001) %MON DXV_max = 7.4403928689702E+04
(PID.TID 0000.0001) %MON DXV_min = 6.0561097546259E+04
(PID.TID 0000.0001) %MON DXV_mean = 6.7606156926271E+04
(PID.TID 0000.0001) %MON DXV_sd = 4.4197252361186E+03
(PID.TID 0000.0001) %MON YC_max = -4.7500000000000E+01
(PID.TID 0000.0001) %MON YC_min = -5.6500000000000E+01
(PID.TID 0000.0001) %MON YC_mean = -5.2000000000000E+01
(PID.TID 0000.0001) %MON YC_sd = 2.8722813232690E+00
(PID.TID 0000.0001) %MON YG_max = -4.8000000000000E+01
(PID.TID 0000.0001) %MON YG_min = -5.7000000000000E+01
(PID.TID 0000.0001) %MON YG_mean = -5.2500000000000E+01
(PID.TID 0000.0001) %MON YG_sd = 2.8722813232690E+00
(PID.TID 0000.0001) %MON DYC_max = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYC_min = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYC_mean = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYC_sd = 8.7311491370201E-11
(PID.TID 0000.0001) %MON DYF_max = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYF_min = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYF_mean = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYF_sd = 8.7311491370201E-11
(PID.TID 0000.0001) %MON DYG_max = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYG_min = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYG_mean = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYG_sd = 8.7311491370201E-11
(PID.TID 0000.0001) %MON DYU_max = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYU_min = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYU_mean = 1.1119492664456E+05
(PID.TID 0000.0001) %MON DYU_sd = 8.7311491370201E-11
(PID.TID 0000.0001) %MON RA_max = 8.3531018943998E+09
(PID.TID 0000.0001) %MON RA_min = 6.8242343148678E+09
(PID.TID 0000.0001) %MON RA_mean = 7.6025723263491E+09
(PID.TID 0000.0001) %MON RA_sd = 4.8813787837650E+08
(PID.TID 0000.0001) %MON RAW_max = 8.3531018943998E+09
(PID.TID 0000.0001) %MON RAW_min = 6.8242343148678E+09
(PID.TID 0000.0001) %MON RAW_mean = 7.6025723263491E+09
(PID.TID 0000.0001) %MON RAW_sd = 4.8813787837650E+08
(PID.TID 0000.0001) %MON RAS_max = 8.2732343846474E+09
(PID.TID 0000.0001) %MON RAS_min = 6.7340013278230E+09
(PID.TID 0000.0001) %MON RAS_mean = 7.5173662459266E+09
(PID.TID 0000.0001) %MON RAS_sd = 4.9144478575377E+08
(PID.TID 0000.0001) %MON RAZ_max = 8.2732343846474E+09
(PID.TID 0000.0001) %MON RAZ_min = 6.7340013278230E+09
(PID.TID 0000.0001) %MON RAZ_mean = 7.5173662459266E+09
(PID.TID 0000.0001) %MON RAZ_sd = 4.9144478575377E+08
(PID.TID 0000.0001) %MON AngleCS_max = 1.0000000000000E+00
(PID.TID 0000.0001) %MON AngleCS_min = 1.0000000000000E+00
(PID.TID 0000.0001) %MON AngleCS_mean = 1.0000000000000E+00
(PID.TID 0000.0001) %MON AngleCS_sd = 0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_max = 0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_min = 0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_mean = 0.0000000000000E+00
(PID.TID 0000.0001) %MON AngleSN_sd = 0.0000000000000E+00
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Field Model R_low (ini_masks_etc) at iteration 1
(PID.TID 0000.0001) // CMIN = -1.000000000000000E+03
(PID.TID 0000.0001) // CMAX = -1.364635703720620E+02
(PID.TID 0000.0001) // CINT = 3.198283072696067E+01
(PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
(PID.TID 0000.0001) // 0.0: .
(PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):( -3: 24: 1)
(PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):( 14: -3: -1)
(PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):( 1: 1: 1)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) K = 1
(PID.TID 0000.0001) // I=6 I=16
(PID.TID 0000.0001) // |--J--|321012345|789012345|78901234
(PID.TID 0000.0001) // 14 -----------sp----c----------
(PID.TID 0000.0001) // 13 ---------------w+.w---------
(PID.TID 0000.0001) // 12 ------------------+---------
(PID.TID 0000.0001) // 11 ----------------------------
(PID.TID 0000.0001) // 10 ----------------------------
(PID.TID 0000.0001) // 9 ----------------------------
(PID.TID 0000.0001) // 8 ----------------------------
(PID.TID 0000.0001) // 7 ----------------------------
(PID.TID 0000.0001) // 6 ----------------------------
(PID.TID 0000.0001) // 5 ----------------------------
(PID.TID 0000.0001) // 4 -----------sp----c----------
(PID.TID 0000.0001) // 3 ---------------w+.w---------
(PID.TID 0000.0001) // 2 ------------------+---------
(PID.TID 0000.0001) // 1 ----------------------------
(PID.TID 0000.0001) // 0 ----------------------------
(PID.TID 0000.0001) // -1 ----------------------------
(PID.TID 0000.0001) // -2 ----------------------------
(PID.TID 0000.0001) // -3 ----------------------------
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // END OF FIELD =
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Field Model Ro_surf (ini_masks_etc) at iteration 1
(PID.TID 0000.0001) // CMIN = 1.000000000000000E+32
(PID.TID 0000.0001) // CMAX = -1.000000000000000E+32
(PID.TID 0000.0001) // CINT = 0.000000000000000E+00
(PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
(PID.TID 0000.0001) // 0.0: .
(PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):( -3: 24: 1)
(PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):( 14: -3: -1)
(PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):( 1: 1: 1)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // END OF FIELD =
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Field hFacC at iteration 1
(PID.TID 0000.0001) // CMIN = 1.000000000000000E+00
(PID.TID 0000.0001) // CMAX = 1.000000000000000E+00
(PID.TID 0000.0001) // CINT = 0.000000000000000E+00
(PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
(PID.TID 0000.0001) // 0.0: .
(PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):( -3: 24: 1)
(PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):( 14: -3: -1)
(PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):( 1: 1: 1)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // END OF FIELD =
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Field hFacW at iteration 1
(PID.TID 0000.0001) // CMIN = 1.000000000000000E+00
(PID.TID 0000.0001) // CMAX = 1.000000000000000E+00
(PID.TID 0000.0001) // CINT = 0.000000000000000E+00
(PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
(PID.TID 0000.0001) // 0.0: .
(PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):( -3: 24: 1)
(PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):( 14: -3: -1)
(PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):( 1: 1: 1)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // END OF FIELD =
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Field hFacS at iteration 1
(PID.TID 0000.0001) // CMIN = 1.000000000000000E+00
(PID.TID 0000.0001) // CMAX = 1.000000000000000E+00
(PID.TID 0000.0001) // CINT = 0.000000000000000E+00
(PID.TID 0000.0001) // SYMBOLS (CMIN->CMAX): -abcdefghijklmnopqrstuvwxyz+
(PID.TID 0000.0001) // 0.0: .
(PID.TID 0000.0001) // RANGE I (Lo:Hi:Step):( -3: 24: 1)
(PID.TID 0000.0001) // RANGE J (Lo:Hi:Step):( 14: -3: -1)
(PID.TID 0000.0001) // RANGE K (Lo:Hi:Step):( 1: 1: 1)
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // END OF FIELD =
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001)
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) // GAD parameters :
(PID.TID 0000.0001) // ===================================
(PID.TID 0000.0001) tempAdvScheme = /* Temp. Horiz.Advection scheme selector */
(PID.TID 0000.0001) 30
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) tempVertAdvScheme = /* Temp. Vert. Advection scheme selector */
(PID.TID 0000.0001) 30
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) tempMultiDimAdvec = /* use Muti-Dim Advec method for Temp */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) tempSOM_Advection = /* use 2nd Order Moment Advection for Temp */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) AdamsBashforthGt = /* apply Adams-Bashforth extrapolation on Gt */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) AdamsBashforth_T = /* apply Adams-Bashforth extrapolation on Temp */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) saltAdvScheme = /* Salt. Horiz.advection scheme selector */
(PID.TID 0000.0001) 30
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) saltVertAdvScheme = /* Salt. Vert. Advection scheme selector */
(PID.TID 0000.0001) 30
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) saltMultiDimAdvec = /* use Muti-Dim Advec method for Salt */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) saltSOM_Advection = /* use 2nd Order Moment Advection for Salt */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) AdamsBashforthGs = /* apply Adams-Bashforth extrapolation on Gs */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) AdamsBashforth_S = /* apply Adams-Bashforth extrapolation on Salt */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) // ===================================
(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) EXF general parameters:
(PID.TID 0000.0001)
(PID.TID 0000.0001) exf_yftype = /* ? */
(PID.TID 0000.0001) 'RL'
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_iprec = /* exf file precision */
(PID.TID 0000.0001) 64
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) useExfYearlyFields = /* add extension _YEAR to input file names */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) twoDigitYear = /* use 2-digit year extension */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) useExfCheckRange = /* check for fields range */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_monFreq = /* EXF monitor frequency [ s ] */
(PID.TID 0000.0001) 1.200000000000000E+04
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) repeatPeriod = /* period for cycling forcing dataset [ s ] */
(PID.TID 0000.0001) 3.162240000000000E+07
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) climTempFreeze= /* Minimum climatological temperature [deg.C] */
(PID.TID 0000.0001) -1.900000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) windStressMax = /* Maximum absolute windstress [ Pa ] */
(PID.TID 0000.0001) 2.000000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) stressIsOnCgrid = /* set u,v_stress on Arakawa C-grid */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cen2kel = /* conversion of deg. Centigrade to Kelvin [K] */
(PID.TID 0000.0001) 2.731500000000000E+02
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) gravity_mks= /* gravitational acceleration [m/s^2] */
(PID.TID 0000.0001) 9.810000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) atmrho = /* mean atmospheric density [kg/m^3] */
(PID.TID 0000.0001) 1.200000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) atmcp = /* mean atmospheric specific heat [J/kg/K] */
(PID.TID 0000.0001) 1.005000000000000E+03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) flamb = /* latent heat of evaporation [J/kg] */
(PID.TID 0000.0001) 2.500000000000000E+06
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) flami = /* latent heat of pure-ice melting [J/kg] */
(PID.TID 0000.0001) 3.340000000000000E+05
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cvapor_fac = /* const. for Saturation calculation [?] */
(PID.TID 0000.0001) 6.403800000000000E+05
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cvapor_exp = /* const. for Saturation calculation [?] */
(PID.TID 0000.0001) 5.107400000000000E+03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cvapor_fac_ice= /* const. for Saturation calculation [?] */
(PID.TID 0000.0001) 1.163780000000000E+07
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cvapor_exp_ice= /* const. for Saturation calculation [?] */
(PID.TID 0000.0001) 5.897800000000000E+03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) humid_fac = /* humidity coef. in virtual temp. [(kg/kg)^-1] */
(PID.TID 0000.0001) 6.060000000000000E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) gamma_blk = /* adiabatic lapse rate [?] */
(PID.TID 0000.0001) 1.000000000000000E-02
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) saltsat = /* reduction of Qsat over salty water [-] */
(PID.TID 0000.0001) 9.800000000000000E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cDrag_1 = /* coef used in drag calculation [?] */
(PID.TID 0000.0001) 2.700000000000000E-03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cDrag_2 = /* coef used in drag calculation [?] */
(PID.TID 0000.0001) 1.420000000000000E-04
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cDrag_3 = /* coef used in drag calculation [?] */
(PID.TID 0000.0001) 7.640000000000000E-05
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cStanton_1 = /* coef used in Stanton number calculation [?] */
(PID.TID 0000.0001) 3.270000000000000E-02
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cStanton_2 = /* coef used in Stanton number calculation [?] */
(PID.TID 0000.0001) 1.800000000000000E-02
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) cDalton = /* coef used in Dalton number calculation [?] */
(PID.TID 0000.0001) 3.460000000000000E-02
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_scal_BulkCdn= /* Drag coefficient scaling factor [-] */
(PID.TID 0000.0001) 1.000000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) zolmin = /* minimum stability parameter [?] */
(PID.TID 0000.0001) -1.000000000000000E+02
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) psim_fac = /* coef used in turbulent fluxes calculation [-] */
(PID.TID 0000.0001) 5.000000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) zref = /* reference height [ m ] */
(PID.TID 0000.0001) 1.000000000000000E+01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) hu = /* height of mean wind [ m ] */
(PID.TID 0000.0001) 1.000000000000000E+01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) ht = /* height of mean temperature [ m ] */
(PID.TID 0000.0001) 2.000000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) hq = /* height of mean spec.humidity [ m ] */
(PID.TID 0000.0001) 2.000000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) uMin = /* minimum wind speed [m/s] */
(PID.TID 0000.0001) 5.000000000000000E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) useStabilityFct_overIce= /* transfert Coeffs over sea-ice depend on stability */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_iceCd = /* drag coefficient over sea-ice (fixed) [-] */
(PID.TID 0000.0001) 1.630000000000000E-03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_iceCe = /* transfert coeff. over sea-ice, for Evap (fixed) [-] */
(PID.TID 0000.0001) 1.630000000000000E-03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_iceCh = /* transfert coeff. over sea-ice, Sens.Heat.(fixed)[-] */
(PID.TID 0000.0001) 1.630000000000000E-03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) exf_albedo = /* Sea-water albedo [-] */
(PID.TID 0000.0001) 1.000000000000000E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) ocean_emissivity = /* longwave ocean-surface emissivity [-] */
(PID.TID 0000.0001) 9.700176366843034E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) ice_emissivity = /* longwave seaice emissivity [-] */
(PID.TID 0000.0001) 9.500000000000000E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) snow_emissivity = /* longwave snow emissivity [-] */
(PID.TID 0000.0001) 9.500000000000000E-01
(PID.TID 0000.0001) ;
(PID.TID 0000.0001)
(PID.TID 0000.0001) EXF main CPP flags:
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ALLOW_ATM_TEMP: defined
(PID.TID 0000.0001) // ALLOW_ATM_WIND: defined
(PID.TID 0000.0001) // ALLOW_DOWNWARD_RADIATION: NOT defined
(PID.TID 0000.0001) // ALLOW_BULKFORMULAE: defined
(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 0.
(PID.TID 0000.0001) Net shortwave flux forcing is read from file:
(PID.TID 0000.0001) >> <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) Zonal wind forcing starts at 0.
(PID.TID 0000.0001) Zonal wind forcing period is 21600.
(PID.TID 0000.0001) Zonal wind forcing is read from file:
(PID.TID 0000.0001) >> SeaWinds_6hrly_uwind_1Jun03_to_Aug03.stretched.bin <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) Meridional wind forcing starts at 0.
(PID.TID 0000.0001) Meridional wind forcing period is 21600.
(PID.TID 0000.0001) Meridional wind forcing is read from file:
(PID.TID 0000.0001) >> SeaWinds_6hrly_vwind_1Jun03_to_Aug03.stretched.bin <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) Atmospheric temperature starts at 0.
(PID.TID 0000.0001) Atmospheric temperature period is 21600.
(PID.TID 0000.0001) Atmospheric temperature is read from file:
(PID.TID 0000.0001) >> ncep_6hrly_atemp_1Jun03_to_Aug03.stretched.bin <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) Atmospheric specific humidity starts at 0.
(PID.TID 0000.0001) Atmospheric specific humidity period is 21600.
(PID.TID 0000.0001) Atmospheric specific humidity is read from file:
(PID.TID 0000.0001) >> ncep_6hrly_rhum_1Jun03_to_Aug03.stretched.bin <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) Net longwave flux forcing starts at 0.
(PID.TID 0000.0001) Net longwave flux forcing period is 0.
(PID.TID 0000.0001) Net longwave flux forcing is read from file:
(PID.TID 0000.0001) >> <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) Precipitation data set starts at 0.
(PID.TID 0000.0001) Precipitation data period is 21600.
(PID.TID 0000.0001) Precipitation data is read from file:
(PID.TID 0000.0001) >> ncep_6hrly_precip_1Jun03_to_Aug03.stretched.bin <<
(PID.TID 0000.0001)
(PID.TID 0000.0001) // EXF_READ_EVAP: NOT defined
(PID.TID 0000.0001)
(PID.TID 0000.0001) // ALLOW_RUNOFF: NOT defined
(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)
(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) // 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) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) cal_CheckDate: Calendar date before predef. reference date
(PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.err
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Parameter file "data.err"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) > 0.25
(PID.TID 0000.0001) > 0.5201 0.2676 0.5 0.5
(PID.TID 0000.0001) > 0.5199 0.2224 0.5 0.5
(PID.TID 0000.0001) > 0.5201 0.1942 0.5 0.5
(PID.TID 0000.0001) > 0.5142 0.1751 0.5 0.5
(PID.TID 0000.0001) > 0.4917 0.1452 0.5 0.5
(PID.TID 0000.0001) > 0.4707 0.1223 0.5 0.5
(PID.TID 0000.0001) > 0.4324 0.1125 0.5 0.5
(PID.TID 0000.0001) > 0.3782 0.1078 0.5 0.5
(PID.TID 0000.0001) > 0.3103 0.0884 0.5 0.5
(PID.TID 0000.0001) > 0.2435 0.0785 0.5 0.5
(PID.TID 0000.0001) > 0.1994 0.0777 0.5 0.5
(PID.TID 0000.0001) > 0.1582 0.0702 0.5 0.5
(PID.TID 0000.0001) > 0.1144 0.0710 0.5 0.5
(PID.TID 0000.0001) > 0.0905 0.0599 0.5 0.5
(PID.TID 0000.0001) > 0.0659 0.0510 0.5 0.5
(PID.TID 0000.0001) > 0.0602 0.0408 0.5 0.5
(PID.TID 0000.0001)
More information about the MITgcm-support
mailing list