(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: checkpoint68q (PID.TID 0000.0001) // Build user: jm_c (PID.TID 0000.0001) // Build host: node374 (PID.TID 0000.0001) // Build date: Tue Jul 25 22:52:12 EDT 2023 (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 namelist terminator (as shown here). (PID.TID 0000.0001) ># Other systems use a / character. (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 = 2 ; /* No. tiles in X per process */ (PID.TID 0000.0001) nSy = 2 ; /* No. tiles in Y per process */ (PID.TID 0000.0001) sNx = 45 ; /* Tile size in X */ (PID.TID 0000.0001) sNy = 20 ; /* 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 = 15 ; /* No. levels in the vertical */ (PID.TID 0000.0001) Nx = 90 ; /* Total domain size in X ( = nPx*nSx*sNx ) */ (PID.TID 0000.0001) Ny = 40 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */ (PID.TID 0000.0001) nTiles = 4 ; /* 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) useNest2W_parent = F ;/* Control 2-W Nesting comm */ (PID.TID 0000.0001) useNest2W_child = F ;/* Control 2-W Nesting comm */ (PID.TID 0000.0001) debugMode = F ; /* print debug msg. (sequence of S/R calls) */ (PID.TID 0000.0001) printMapIncludesZeros= F ; /* print zeros in Std.Output maps */ (PID.TID 0000.0001) maxLengthPrt1D= 65 /* maxLength of 1D array printed to StdOut */ (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: 2, 1: 2) (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. = 000000) (PID.TID 0000.0001) // WEST: Tile = 000002, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000001 (PID.TID 0000.0001) // EAST: Tile = 000002, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000001 (PID.TID 0000.0001) // SOUTH: Tile = 000003, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000002 (PID.TID 0000.0001) // NORTH: Tile = 000003, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000002 (PID.TID 0000.0001) // Tile number: 000002 (process no. = 000000) (PID.TID 0000.0001) // WEST: Tile = 000001, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000001 (PID.TID 0000.0001) // EAST: Tile = 000001, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000001 (PID.TID 0000.0001) // SOUTH: Tile = 000004, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000002 (PID.TID 0000.0001) // NORTH: Tile = 000004, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000002 (PID.TID 0000.0001) // Tile number: 000003 (process no. = 000000) (PID.TID 0000.0001) // WEST: Tile = 000004, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000002 (PID.TID 0000.0001) // EAST: Tile = 000004, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000002 (PID.TID 0000.0001) // SOUTH: Tile = 000001, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000001 (PID.TID 0000.0001) // NORTH: Tile = 000001, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000001 (PID.TID 0000.0001) // Tile number: 000004 (process no. = 000000) (PID.TID 0000.0001) // WEST: Tile = 000003, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000002 (PID.TID 0000.0001) // EAST: Tile = 000003, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000001, bj = 000002 (PID.TID 0000.0001) // SOUTH: Tile = 000002, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000001 (PID.TID 0000.0001) // NORTH: Tile = 000002, Process = 000000, Comm = put (PID.TID 0000.0001) // bi = 000002, bj = 000001 (PID.TID 0000.0001) (PID.TID 0000.0001) INI_PARMS: opening model parameter file "data" (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // 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) > &PARM01 (PID.TID 0000.0001) > tRef = 15*20., (PID.TID 0000.0001) > sRef = 15*35., (PID.TID 0000.0001) > viscAr=1.E-3, (PID.TID 0000.0001) > viscAh=5.E5, (PID.TID 0000.0001) > diffKhT=0.0, (PID.TID 0000.0001) > diffKrT=3.E-5, (PID.TID 0000.0001) > diffKhS=0.0, (PID.TID 0000.0001) > diffKrS=3.E-5, (PID.TID 0000.0001) > rhoConst=1035., (PID.TID 0000.0001) > rotationPeriod=86400., (PID.TID 0000.0001) > gravity=9.81, (PID.TID 0000.0001) > eosType = 'JMD95Z', (PID.TID 0000.0001) > ivdc_kappa=100., (PID.TID 0000.0001) > implicitDiffusion=.TRUE., (PID.TID 0000.0001) > allowFreezing=.TRUE., (PID.TID 0000.0001) > useRealFreshWaterFlux=.TRUE., (PID.TID 0000.0001) > useCDscheme=.TRUE., (PID.TID 0000.0001) ># turn on looped cells (PID.TID 0000.0001) > hFacMin=.05, (PID.TID 0000.0001) > hFacMindr=50., (PID.TID 0000.0001) ># set precision of data files (PID.TID 0000.0001) > readBinaryPrec=32, (PID.TID 0000.0001) > / (PID.TID 0000.0001) > (PID.TID 0000.0001) ># Elliptic solver parameters (PID.TID 0000.0001) > &PARM02 (PID.TID 0000.0001) > cg2dMaxIters=1000, (PID.TID 0000.0001) > cg2dTargetResidual=1.E-19, (PID.TID 0000.0001) > / (PID.TID 0000.0001) > (PID.TID 0000.0001) ># Time stepping parameters (PID.TID 0000.0001) > &PARM03 (PID.TID 0000.0001) > nIter0 = 0, (PID.TID 0000.0001) > nTimeSteps = 4, (PID.TID 0000.0001) ># 100 years of integration will yield a reasonable flow field (PID.TID 0000.0001) ># startTime = 0., (PID.TID 0000.0001) ># endTime = 3110400000., (PID.TID 0000.0001) > deltaTMom = 1200.0, (PID.TID 0000.0001) > tauCD = 321428., (PID.TID 0000.0001) > deltaTtracer= 43200.0, (PID.TID 0000.0001) > deltaTClock = 43200.0, (PID.TID 0000.0001) ># if you are using a version later than checkpoint45d on the main branch (PID.TID 0000.0001) ># you can uncomment the following line and increase the time step (PID.TID 0000.0001) ># deltaTtracer and deltaTClock to 172800.0 as well to speed up the (PID.TID 0000.0001) ># asynchronous time stepping (PID.TID 0000.0001) ># deltaTfreesurf = 172800.0, (PID.TID 0000.0001) > abEps = 0.1, (PID.TID 0000.0001) > pChkptFreq= 311040000., (PID.TID 0000.0001) > dumpFreq = 2592000., (PID.TID 0000.0001) > adjDumpFreq = 2592000., (PID.TID 0000.0001) > monitorFreq = 0., (PID.TID 0000.0001) > adjMonitorFreq = 0., (PID.TID 0000.0001) ># 2 months restoring timescale for temperature (PID.TID 0000.0001) > tauThetaClimRelax = 5184000.0, (PID.TID 0000.0001) ># 6 months restoring timescale for salinity (PID.TID 0000.0001) > tauSaltClimRelax = 15552000.0, (PID.TID 0000.0001) > periodicExternalForcing=.TRUE., (PID.TID 0000.0001) > externForcingPeriod=2592000., (PID.TID 0000.0001) > externForcingCycle=31104000., (PID.TID 0000.0001) > / (PID.TID 0000.0001) > (PID.TID 0000.0001) ># Gridding parameters (PID.TID 0000.0001) > &PARM04 (PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE., (PID.TID 0000.0001) > delR= 50., 70., 100., 140., 190., (PID.TID 0000.0001) > 240., 290., 340., 390., 440., (PID.TID 0000.0001) > 490., 540., 590., 640., 690., (PID.TID 0000.0001) > ygOrigin=-80., (PID.TID 0000.0001) > dySpacing=4., (PID.TID 0000.0001) > dxSpacing=4., (PID.TID 0000.0001) > / (PID.TID 0000.0001) > (PID.TID 0000.0001) ># Input datasets (PID.TID 0000.0001) > &PARM05 (PID.TID 0000.0001) > bathyFile= 'bathymetry.bin', (PID.TID 0000.0001) > hydrogThetaFile='lev_t.bin', (PID.TID 0000.0001) > hydrogSaltFile= 'lev_s.bin', (PID.TID 0000.0001) > zonalWindFile= 'trenberth_taux.bin', (PID.TID 0000.0001) > meridWindFile= 'trenberth_tauy.bin', (PID.TID 0000.0001) > thetaClimFile= 'lev_sst.bin', (PID.TID 0000.0001) > saltClimFile= 'lev_sss.bin', (PID.TID 0000.0001) > surfQnetFile= 'ncep_qnet.bin', (PID.TID 0000.0001) ># fresh water flux is turned off, uncomment next line to turn on (PID.TID 0000.0001) ># (not recommened together with surface salinity restoring) (PID.TID 0000.0001) ># EmPmRFile= 'ncep_emp.bin', (PID.TID 0000.0001) > / (PID.TID 0000.0001) (PID.TID 0000.0001) INI_PARMS ; starts to read PARM01 (PID.TID 0000.0001) INI_PARMS ; read PARM01 : OK (PID.TID 0000.0001) S/R INI_PARMS: No request for barotropic solver (PID.TID 0000.0001) S/R INI_PARMS: => Use implicitFreeSurface as default (PID.TID 0000.0001) INI_PARMS ; starts to read PARM02 (PID.TID 0000.0001) INI_PARMS ; read PARM02 : OK (PID.TID 0000.0001) INI_PARMS ; starts to read PARM03 (PID.TID 0000.0001) INI_PARMS ; read PARM03 : OK (PID.TID 0000.0001) INI_PARMS ; starts to read PARM04 (PID.TID 0000.0001) INI_PARMS ; read PARM04 : OK (PID.TID 0000.0001) INI_PARMS ; starts to read PARM05 (PID.TID 0000.0001) INI_PARMS ; read PARM05 : OK (PID.TID 0000.0001) INI_PARMS: finished reading file "data" (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) > useGGL90 = .FALSE., (PID.TID 0000.0001) > useKPP = .FALSE., (PID.TID 0000.0001) > useGrdchk = .TRUE., (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) PACKAGES_BOOT: On/Off package Summary -------- pkgs with a standard "usePKG" On/Off switch in "data.pkg": -------- pkg/ggl90 compiled but not used ( useGGL90 = F ) pkg/kpp compiled but not used ( useKPP = F ) pkg/gmredi compiled and used ( useGMRedi = T ) pkg/autodiff compiled and used ( useAUTODIFF = T ) pkg/grdchk compiled and used ( useGrdchk = T ) pkg/ctrl compiled and used ( useCTRL = T ) -------- pkgs without standard "usePKG" On/Off switch in "data.pkg": -------- pkg/generic_advdiff compiled and used ( useGAD = T ) pkg/mom_common compiled and used ( momStepping = T ) pkg/mom_vecinv compiled but not used ( +vectorInvariantMomentum = F ) pkg/mom_fluxform compiled and used ( & not vectorInvariantMom = T ) pkg/cd_code compiled and used ( useCDscheme = T ) pkg/debug compiled but not used ( debugMode = F ) pkg/rw compiled and used pkg/mdsio compiled and used pkg/autodiff compiled and used pkg/openad compiled and used pkg/cost compiled and used (PID.TID 0000.0001) PACKAGES_BOOT: End of package Summary (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-12, (PID.TID 0000.0001) > GM_slopeSqCutoff = 1.D+08, (PID.TID 0000.0001) > GM_AdvForm = .FALSE., (PID.TID 0000.0001) > GM_isopycK = 1.0D+3, (PID.TID 0000.0001) > GM_background_K = 1.0D+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.D0, (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) > (PID.TID 0000.0001) (PID.TID 0000.0001) GM_READPARMS: finished reading data.gmredi (PID.TID 0000.0001) AUTODIFF_READPARMS: opening data.autodiff (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.autodiff (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Parameter file "data.autodiff" (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) ># ========================= (PID.TID 0000.0001) ># pkg AUTODIFF parameters : (PID.TID 0000.0001) ># ========================= (PID.TID 0000.0001) ># inAdExact :: get an exact adjoint (no approximation) (def=.True.) (PID.TID 0000.0001) ># (PID.TID 0000.0001) > &AUTODIFF_PARM01 (PID.TID 0000.0001) ># inAdExact = .FALSE., (PID.TID 0000.0001) > / (PID.TID 0000.0001) (PID.TID 0000.0001) AUTODIFF_READPARMS: finished reading data.autodiff (PID.TID 0000.0001) // =================================== (PID.TID 0000.0001) // AUTODIFF parameters : (PID.TID 0000.0001) // =================================== (PID.TID 0000.0001) inAdExact = /* get an exact adjoint (no approximation) */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) useApproxAdvectionInAdMode = /* approximate AD-advection */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) cg2dFullAdjoint = /* use full hand written cg2d adjoint (no approximation) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useKPPinAdMode = /* use KPP in adjoint mode */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useGMRediInAdMode = /* use GMRedi in adjoint mode */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) useSEAICEinAdMode = /* use SEAICE in adjoint mode */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useGGL90inAdMode = /* use GGL90 in adjoint mode */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useSALT_PLUMEinAdMode = /* use SALT_PLUME in adjoint mode */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) dumpAdVarExch = /* control adexch before dumpinp */ (PID.TID 0000.0001) 2 (PID.TID 0000.0001) ; (PID.TID 0000.0001) mon_AdVarExch = /* control adexch before monitor */ (PID.TID 0000.0001) 2 (PID.TID 0000.0001) ; (PID.TID 0000.0001) viscFacInFw = /* viscosity factor for forward model */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) viscFacInAd = /* viscosity factor for adjoint */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) SIregFacInAd = /* sea ice factor for adjoint model */ (PID.TID 0000.0001) 1.234567000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) SIregFacInFw = /* sea ice factor for forward model */ (PID.TID 0000.0001) 1.234567000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) (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) 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) ># ECCO controlvariables (PID.TID 0000.0001) ># ********************* (PID.TID 0000.0001) > &CTRL_NML (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) > / (PID.TID 0000.0001) ># (PID.TID 0000.0001) ># ********************* (PID.TID 0000.0001) ># names for CTRL_GENARR, CTRL_GENTIM (PID.TID 0000.0001) ># ********************* (PID.TID 0000.0001) > &CTRL_NML_GENARR (PID.TID 0000.0001) ># the bounds are commented out to reproduce old results with (PID.TID 0000.0001) ># the CPP-flag DISABLE_CTRL_THETA_LIMIT defined (PID.TID 0000.0001) > xx_genarr3d_file(1) = 'xx_theta', (PID.TID 0000.0001) > xx_genarr3d_weight(1) = 'ones_64b.bin', (PID.TID 0000.0001) >#xx_genarr3d_bounds(1:5,1) = -2.,-1.9,39.,40.,5., (PID.TID 0000.0001) ># (PID.TID 0000.0001) > xx_genarr3d_file(2) = 'xx_salt', (PID.TID 0000.0001) > xx_genarr3d_weight(2) = 'ones_64b.bin', (PID.TID 0000.0001) >#xx_genarr3d_bounds(1:5,2) = 29.,29.5,40.5,41.,5., (PID.TID 0000.0001) ># (PID.TID 0000.0001) > xx_genarr2d_file(1) = 'xx_qnet', (PID.TID 0000.0001) > xx_genarr2d_weight(1) = 'ones_64b.bin', (PID.TID 0000.0001) ># (PID.TID 0000.0001) > xx_genarr2d_file(2) = 'xx_empmr', (PID.TID 0000.0001) > xx_genarr2d_weight(2) = 'ones_64b.bin', (PID.TID 0000.0001) ># (PID.TID 0000.0001) > xx_genarr2d_file(3) = 'xx_fu', (PID.TID 0000.0001) > xx_genarr2d_weight(3) = 'ones_64b.bin', (PID.TID 0000.0001) ># (PID.TID 0000.0001) > xx_genarr2d_file(4) = 'xx_fv', (PID.TID 0000.0001) > xx_genarr2d_weight(4) = 'ones_64b.bin', (PID.TID 0000.0001) ># (PID.TID 0000.0001) > / (PID.TID 0000.0001) (PID.TID 0000.0001) CTRL_READPARMS: finished reading data.ctrl (PID.TID 0000.0001) read-write ctrl files from current run directory (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) ># ECCO cost function (PID.TID 0000.0001) ># ****************** (PID.TID 0000.0001) > &COST_NML (PID.TID 0000.0001) ># (PID.TID 0000.0001) > mult_tracer = 1., (PID.TID 0000.0001) > mult_test = 1., (PID.TID 0000.0001) > mult_atl = 1., (PID.TID 0000.0001) > / (PID.TID 0000.0001) (PID.TID 0000.0001) COST_READPARMS: finished reading data.cost (PID.TID 0000.0001) GRDCHK_READPARMS: opening data.grdchk (PID.TID 0000.0001) OPEN_COPY_DATA_FILE: opening file data.grdchk (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Parameter file "data.grdchk" (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) > (PID.TID 0000.0001) ># ******************* (PID.TID 0000.0001) ># ECCO gradient check (PID.TID 0000.0001) ># ******************* (PID.TID 0000.0001) > &GRDCHK_NML (PID.TID 0000.0001) > grdchk_eps = 1.d-2, (PID.TID 0000.0001) > iGloTile = 2, (PID.TID 0000.0001) > jGloTile = 2, (PID.TID 0000.0001) > iGloPos = 26, (PID.TID 0000.0001) >#jGloPos = 19, (PID.TID 0000.0001) ># This original setting would start the gradient check for temperatures (PID.TID 0000.0001) ># that are below the freezing point in the inital conditions. Because of (PID.TID 0000.0001) ># allowFreezing=True in data&PARM01 and the unphysical nature of this (PID.TID 0000.0001) ># freezing parameterisation, the FD gradients depend very much on the (PID.TID 0000.0001) ># whether they are evaluated above and below freezing or not leading (PID.TID 0000.0001) ># to inconsistent results. Moving the starting point to jGloPos=17 (PID.TID 0000.0001) ># avoids this. This is similar to input_ad but here we use jGloTile = 2. (PID.TID 0000.0001) > jGloPos = 17, (PID.TID 0000.0001) > kGloPos = 1, (PID.TID 0000.0001) >### nbeg = 1, (PID.TID 0000.0001) > nstep = 1, (PID.TID 0000.0001) > nend = 7, (PID.TID 0000.0001) > grdchkvarindex = 201, (PID.TID 0000.0001) > / (PID.TID 0000.0001) (PID.TID 0000.0001) GRDCHK_READPARMS: finished reading data.grdchk (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Gradient check configuration >>> START <<< (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) grdchkvarindex : 201 (PID.TID 0000.0001) eps: 0.100E-01 (PID.TID 0000.0001) First location: 0 (PID.TID 0000.0001) Last location: 7 (PID.TID 0000.0001) Increment: 1 (PID.TID 0000.0001) grdchkWhichProc: 0 (PID.TID 0000.0001) iLocTile = 2 , jLocTile = 2 (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Gradient check configuration >>> END <<< (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) SET_PARMS: done (PID.TID 0000.0001) Enter INI_VERTICAL_GRID: setInterFDr= T ; setCenterDr= F (PID.TID 0000.0001) GAD_INIT_FIXED: GAD_OlMinSize= 1 0 1 (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) 2 (PID.TID 0000.0001) ; (PID.TID 0000.0001) tempVertAdvScheme = /* Temp. Vert. Advection scheme selector */ (PID.TID 0000.0001) 2 (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) T (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) 2 (PID.TID 0000.0001) ; (PID.TID 0000.0001) saltVertAdvScheme = /* Salt. Vert. Advection scheme selector */ (PID.TID 0000.0001) 2 (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) T (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) ctrl-wet 1: nvarlength = 67878 (PID.TID 0000.0001) ctrl-wet 2: surface wet C = 669 (PID.TID 0000.0001) ctrl-wet 3: surface wet W = 640 (PID.TID 0000.0001) ctrl-wet 4: surface wet S = 610 (PID.TID 0000.0001) ctrl-wet 4a:surface wet V = 0 (PID.TID 0000.0001) ctrl-wet 5: 3D wet points = 8651 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 1 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 2 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 3 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 4 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 5 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 6 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 7 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 8 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 9 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 10 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 11 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 12 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 13 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 14 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 15 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 16 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 17 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 18 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 19 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 20 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 21 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 22 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 23 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 24 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 25 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 26 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 27 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 28 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 29 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 30 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 31 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 32 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 33 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 34 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 35 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 36 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 37 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 38 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 39 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 40 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 41 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 42 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 43 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 44 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 45 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 46 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 47 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 48 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 49 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 50 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 51 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 52 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 53 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 54 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 55 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 56 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 57 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 58 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 59 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 60 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 61 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 62 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 63 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 64 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 65 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 66 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 67 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 68 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 69 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 70 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 71 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 72 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 73 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 74 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 75 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 76 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 77 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 78 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 79 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 80 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 81 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 82 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 83 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 84 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 85 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 86 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 87 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 88 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 89 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 90 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 91 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 92 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 93 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 94 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 95 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 96 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 97 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 98 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 99 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 100 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 101 1 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 102 1 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 103 1 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 104 1 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 105 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 106 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 107 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 108 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 109 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 110 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 111 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 112 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 113 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 114 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 115 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 116 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 117 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 118 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 119 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 120 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 121 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 122 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 123 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 124 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 125 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 126 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 127 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 128 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 129 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 130 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 131 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 132 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 133 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 134 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 135 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 136 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 137 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 138 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 139 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 140 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 141 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 142 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 143 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 144 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 145 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 146 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 147 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 148 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 149 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 150 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 151 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 152 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 153 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 154 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 155 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 156 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 157 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 158 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 159 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 160 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 161 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 162 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 163 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 164 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 165 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 166 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 167 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 168 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 169 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 170 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 171 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 172 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 173 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 174 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 175 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 176 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 177 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 178 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 179 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 180 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 181 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 182 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 183 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 184 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 185 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 186 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 187 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 188 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 189 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 190 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 191 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 192 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 193 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 194 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 195 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 196 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 197 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 198 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 199 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 200 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 201 1 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 202 1 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 203 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 204 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 205 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 206 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 207 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 208 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 209 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 210 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 211 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 212 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 213 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 214 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 215 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 216 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 217 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 218 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 219 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 220 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 221 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 222 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 223 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 224 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 225 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 226 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 227 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 228 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 229 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 230 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 231 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 232 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 233 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 234 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 235 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 236 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 237 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 238 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 239 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 240 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 241 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 242 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 243 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 244 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 245 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 246 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 247 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 248 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 249 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 250 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 251 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 252 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 253 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 254 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 255 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 256 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 257 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 258 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 259 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 260 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 261 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 262 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 263 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 264 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 265 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 266 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 267 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 268 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 269 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 270 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 271 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 272 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 273 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 274 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 275 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 276 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 277 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 278 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 279 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 280 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 281 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 282 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 283 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 284 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 285 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 286 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 287 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 288 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 289 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 290 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 291 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 292 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 293 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 294 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 295 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 296 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 297 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 298 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 299 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 300 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 301 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 302 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 303 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 304 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 305 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 306 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 307 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 308 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 309 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 310 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 311 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 312 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 313 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 314 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 315 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 316 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 317 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 318 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 319 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 320 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 321 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 322 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 323 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 324 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 325 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 326 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 327 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 328 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 329 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 330 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 331 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 332 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 333 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 334 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 335 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 336 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 337 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 338 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 339 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 340 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 341 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 342 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 343 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 344 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 345 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 346 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 347 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 348 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 349 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 350 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 351 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 352 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 353 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 354 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 355 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 356 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 357 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 358 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 359 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 360 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 361 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 362 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 363 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 364 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 365 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 366 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 367 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 368 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 369 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 370 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 371 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 372 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 373 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 374 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 375 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 376 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 377 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 378 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 379 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 380 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 381 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 382 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 383 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 384 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 385 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 386 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 387 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 388 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 389 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 390 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 391 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 392 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 393 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 394 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 395 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 396 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 397 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 398 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 399 0 (PID.TID 0000.0001) ctrl-wet 6: no recs for i = 400 0 (PID.TID 0000.0001) ctrl-wet 7: flux 17302 (PID.TID 0000.0001) ctrl-wet 8: atmos 17302 (PID.TID 0000.0001) ctrl-wet ------------------------------------------------- (PID.TID 0000.0001) ctrl-wet 13: global nvarlength for Nr = 15 67878 (PID.TID 0000.0001) ctrl-wet ------------------------------------------------- (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 1 2315 2149 2206 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 2 2315 2149 2206 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 3 2254 2102 2146 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 4 2215 2058 2104 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 5 2178 2027 2070 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 6 2142 1987 2029 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 7 2114 1959 2004 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 8 2076 1918 1959 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 9 2048 1887 1925 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 10 1999 1831 1869 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 11 1948 1771 1808 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 12 1850 1653 1705 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 13 1655 1404 1458 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 14 1372 1118 1164 0 (PID.TID 0000.0001) ctrl-wet 14: global nWet C/S/W/V k= 15 828 623 671 0 (PID.TID 0000.0001) ctrl-wet ------------------------------------------------- (PID.TID 0000.0001) ctrl-wet ------------------------------------------------- (PID.TID 0000.0001) ctrl-wet ------------------------------------------------- (PID.TID 0000.0001) ctrl_init_wet: no. of control variables: 6 (PID.TID 0000.0001) ctrl_init_wet: control vector length: 67878 (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // control vector configuration >>> START <<< (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) Total number of ocean points per tile: (PID.TID 0000.0001) -------------------------------------- (PID.TID 0000.0001) snx*sny*nr = 13500 (PID.TID 0000.0001) (PID.TID 0000.0001) Number of ocean points per tile: (PID.TID 0000.0001) -------------------------------- (PID.TID 0000.0001) bi,bj,#(c/s/w): 0001 0001 008651 007584 008016 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0002 0001 010278 009360 009834 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0001 0002 003504 003232 003049 (PID.TID 0000.0001) bi,bj,#(c/s/w): 0002 0002 006876 006460 006425 (PID.TID 0000.0001) (PID.TID 0000.0001) Settings of generic controls: (PID.TID 0000.0001) ----------------------------- (PID.TID 0000.0001) (PID.TID 0000.0001) -> 2D control, genarr2d no. 1 is in use (PID.TID 0000.0001) file = xx_qnet (PID.TID 0000.0001) weight = ones_64b.bin (PID.TID 0000.0001) index = 0101 (PID.TID 0000.0001) ncvarindex = 0201 (PID.TID 0000.0001) -> 2D control, genarr2d no. 2 is in use (PID.TID 0000.0001) file = xx_empmr (PID.TID 0000.0001) weight = ones_64b.bin (PID.TID 0000.0001) index = 0102 (PID.TID 0000.0001) ncvarindex = 0202 (PID.TID 0000.0001) -> 2D control, genarr2d no. 3 is in use (PID.TID 0000.0001) file = xx_fu (PID.TID 0000.0001) weight = ones_64b.bin (PID.TID 0000.0001) index = 0103 (PID.TID 0000.0001) ncvarindex = 0203 (PID.TID 0000.0001) -> 2D control, genarr2d no. 4 is in use (PID.TID 0000.0001) file = xx_fv (PID.TID 0000.0001) weight = ones_64b.bin (PID.TID 0000.0001) index = 0104 (PID.TID 0000.0001) ncvarindex = 0204 (PID.TID 0000.0001) -> 3d control, genarr3d no. 1 is in use (PID.TID 0000.0001) file = xx_theta (PID.TID 0000.0001) weight = ones_64b.bin (PID.TID 0000.0001) index = 0201 (PID.TID 0000.0001) ncvarindex = 0301 (PID.TID 0000.0001) -> 3d control, genarr3d no. 2 is in use (PID.TID 0000.0001) file = xx_salt (PID.TID 0000.0001) weight = ones_64b.bin (PID.TID 0000.0001) index = 0202 (PID.TID 0000.0001) ncvarindex = 0302 (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // control vector configuration >>> END <<< (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) INI_CG2D: CG2D normalisation factor = 6.5682677425711703E-05 (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model configuration (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // (PID.TID 0000.0001) // "Physical" paramters ( PARM01 in namelist ) (PID.TID 0000.0001) // (PID.TID 0000.0001) buoyancyRelation = /* Type of relation to get Buoyancy */ (PID.TID 0000.0001) 'OCEANIC' (PID.TID 0000.0001) ; (PID.TID 0000.0001) fluidIsAir = /* fluid major constituent is Air */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) fluidIsWater = /* fluid major constituent is Water */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) usingPCoords = /* use p (or p*) vertical coordinate */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) usingZCoords = /* use z (or z*) vertical coordinate */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) tRef = /* Reference temperature profile ( oC or K ) */ (PID.TID 0000.0001) 15 @ 2.000000000000000E+01 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) sRef = /* Reference salinity profile ( g/kg ) */ (PID.TID 0000.0001) 15 @ 3.500000000000000E+01 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rhoRef = /* Density vertical profile from (Ref,sRef)( kg/m^3 ) */ (PID.TID 0000.0001) 1.024872626184147E+03, /* K = 1 */ (PID.TID 0000.0001) 1.025135462285008E+03, /* K = 2 */ (PID.TID 0000.0001) 1.025507198938228E+03, /* K = 3 */ (PID.TID 0000.0001) 1.026030780760464E+03, /* K = 4 */ (PID.TID 0000.0001) 1.026748377776259E+03, /* K = 5 */ (PID.TID 0000.0001) 1.027679406285166E+03, /* K = 6 */ (PID.TID 0000.0001) 1.028820735595355E+03, /* K = 7 */ (PID.TID 0000.0001) 1.030168558073105E+03, /* K = 8 */ (PID.TID 0000.0001) 1.031718419899614E+03, /* K = 9 */ (PID.TID 0000.0001) 1.033465256541184E+03, /* K = 10 */ (PID.TID 0000.0001) 1.035403432414885E+03, /* K = 11 */ (PID.TID 0000.0001) 1.037526784183520E+03, /* K = 12 */ (PID.TID 0000.0001) 1.039828667078104E+03, /* K = 13 */ (PID.TID 0000.0001) 1.042302003623418E+03, /* K = 14 */ (PID.TID 0000.0001) 1.044939334132512E+03 /* K = 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dBdrRef = /* Vertical grad. of reference buoyancy [(m/s/r)^2] */ (PID.TID 0000.0001) 15 @ 0.000000000000000E+00 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) useStrainTensionVisc= /* Use StrainTension Form of Viscous Operator */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useVariableVisc = /* Use variable horizontal viscosity */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useHarmonicVisc = /* Use harmonic horizontal viscosity */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) useBiharmonicVisc= /* Use biharmonic horiz. viscosity */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useSmag3D = /* Use isotropic 3-D Smagorinsky viscosity */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) viscAh = /* Lateral harmonic viscosity ( m^2/s ) */ (PID.TID 0000.0001) 5.000000000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) viscA4 = /* Lateral biharmonic viscosity ( m^4/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) no_slip_sides = /* Viscous BCs: No-slip sides */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) sideDragFactor = /* side-drag scaling factor (non-dim) */ (PID.TID 0000.0001) 2.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) viscArNr = /* vertical profile of vertical viscosity ( m^2/s )*/ (PID.TID 0000.0001) 15 @ 1.000000000000000E-03 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) no_slip_bottom = /* Viscous BCs: No-slip bottom */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) bottomVisc_pCell = /* Partial-cell in bottom Visc. BC */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) bottomDragLinear = /* linear bottom-drag coefficient ( m/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) bottomDragQuadratic = /* quadratic bottom-drag coefficient (-) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectBotDragQuadr = /* select quadratic bottom drag options */ (PID.TID 0000.0001) -1 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKhT = /* Laplacian diffusion of heat laterally ( m^2/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffK4T = /* Biharmonic diffusion of heat laterally ( m^4/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKhS = /* Laplacian diffusion of salt laterally ( m^2/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffK4S = /* Biharmonic diffusion of salt laterally ( m^4/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKrNrT = /* vertical profile of vertical diffusion of Temp ( m^2/s )*/ (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKrNrS = /* vertical profile of vertical diffusion of Salt ( m^2/s )*/ (PID.TID 0000.0001) 15 @ 3.000000000000000E-05 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKrBL79surf = /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKrBL79deep = /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKrBL79scl = /* Depth scale for Bryan and Lewis 1979 ( m ) */ (PID.TID 0000.0001) 2.000000000000000E+02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) diffKrBL79Ho = /* Turning depth for Bryan and Lewis 1979 ( m ) */ (PID.TID 0000.0001) -2.000000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) ivdc_kappa = /* Implicit Vertical Diffusivity for Convection ( m^2/s) */ (PID.TID 0000.0001) 1.000000000000000E+02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) hMixCriteria= /* Criteria for mixed-layer diagnostic */ (PID.TID 0000.0001) -8.000000000000000E-01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) dRhoSmall = /* Parameter for mixed-layer diagnostic */ (PID.TID 0000.0001) 1.000000000000000E-06 (PID.TID 0000.0001) ; (PID.TID 0000.0001) hMixSmooth= /* Smoothing parameter for mixed-layer diagnostic */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) eosType = /* Type of Equation of State */ (PID.TID 0000.0001) 'JMD95Z' (PID.TID 0000.0001) ; (PID.TID 0000.0001) eosRefP0 = /* Reference atmospheric pressure for EOS ( Pa ) */ (PID.TID 0000.0001) 1.013250000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectP_inEOS_Zc = /* select pressure to use in EOS (0,1,2,3) */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) 0= -g*rhoConst*z ; 1= pRef (from tRef,sRef); 2= Hyd P ; 3= Hyd+NH P (PID.TID 0000.0001) ; (PID.TID 0000.0001) surf_pRef = /* Surface reference pressure ( Pa ) */ (PID.TID 0000.0001) 1.013250000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) HeatCapacity_Cp = /* Specific heat capacity ( J/kg/K ) */ (PID.TID 0000.0001) 3.994000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) celsius2K = /* 0 degree Celsius converted to Kelvin ( K ) */ (PID.TID 0000.0001) 2.731500000000000E+02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rhoConst = /* Reference density (Boussinesq) ( kg/m^3 ) */ (PID.TID 0000.0001) 1.035000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rhoFacC = /* normalized Reference density @ cell-Center (-) */ (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rhoFacF = /* normalized Reference density @ W-Interface (-) */ (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rhoConstFresh = /* Fresh-water reference density ( kg/m^3 ) */ (PID.TID 0000.0001) 1.035000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) gravity = /* Gravitational acceleration ( m/s^2 ) */ (PID.TID 0000.0001) 9.810000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) gBaro = /* Barotropic gravity ( m/s^2 ) */ (PID.TID 0000.0001) 9.810000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) gravFacC = /* gravity factor (vs surf.) @ cell-Center (-) */ (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) gravFacF = /* gravity factor (vs surf.) @ W-Interface (-) */ (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rotationPeriod = /* Rotation Period ( s ) */ (PID.TID 0000.0001) 8.640000000000000E+04 (PID.TID 0000.0001) ; (PID.TID 0000.0001) omega = /* Angular velocity ( rad/s ) */ (PID.TID 0000.0001) 7.272205216643040E-05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) f0 = /* Reference coriolis parameter ( 1/s ) */ (PID.TID 0000.0001) 1.000000000000000E-04 (PID.TID 0000.0001) ; (PID.TID 0000.0001) beta = /* Beta ( 1/(m.s) ) */ (PID.TID 0000.0001) 9.999999999999999E-12 (PID.TID 0000.0001) ; (PID.TID 0000.0001) fPrime = /* Second coriolis parameter ( 1/s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rigidLid = /* Rigid lid on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicitFreeSurface = /* Implicit free surface on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) freeSurfFac = /* Implicit free surface factor */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicSurfPress = /* Surface Pressure implicit factor (0-1) */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicDiv2DFlow = /* Barot. Flow Div. implicit factor (0-1) */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) uniformLin_PhiSurf = /* use uniform Bo_surf on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) uniformFreeSurfLev = /* free-surface level-index is uniform */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) sIceLoadFac = /* scale factor for sIceLoad (0-1) */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) hFacMin = /* minimum partial cell factor (hFac) */ (PID.TID 0000.0001) 5.000000000000000E-02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) hFacMinDr = /* minimum partial cell thickness ( m) */ (PID.TID 0000.0001) 5.000000000000000E+01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) exactConserv = /* Exact Volume Conservation on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) linFSConserveTr = /* Tracer correction for Lin Free Surface on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) nonlinFreeSurf = /* Non-linear Free Surf. options (-1,0,1,2,3)*/ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) -1,0= Off ; 1,2,3= On, 2=+rescale gU,gV, 3=+update cg2d solv. (PID.TID 0000.0001) ; (PID.TID 0000.0001) hFacInf = /* lower threshold for hFac (nonlinFreeSurf only)*/ (PID.TID 0000.0001) 2.000000000000000E-01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) hFacSup = /* upper threshold for hFac (nonlinFreeSurf only)*/ (PID.TID 0000.0001) 2.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) select_rStar = /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) useRealFreshWaterFlux = /* Real Fresh Water Flux on/off flag*/ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) temp_EvPrRn = /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/ (PID.TID 0000.0001) 1.234567000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) salt_EvPrRn = /* Salin. of Evap/Prec/R (UNSET=use local S)(g/kg)*/ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectAddFluid = /* option for mass source/sink of fluid (=0: off) */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) temp_addMass = /* Temp. of addMass array (UNSET=use local T)(oC)*/ (PID.TID 0000.0001) 1.234567000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) salt_addMass = /* Salin. of addMass array (UNSET=use local S)(g/kg)*/ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) convertFW2Salt = /* convert F.W. Flux to Salt Flux (-1=use local S)(g/kg)*/ (PID.TID 0000.0001) -1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) use3Dsolver = /* use 3-D pressure solver on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) nonHydrostatic = /* Non-Hydrostatic on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) nh_Am2 = /* Non-Hydrostatic terms scaling factor */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicitNHPress = /* Non-Hyd Pressure implicit factor (0-1)*/ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectNHfreeSurf = /* Non-Hyd (free-)Surface option */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) quasiHydrostatic = /* Quasi-Hydrostatic on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) calc_wVelocity = /* vertical velocity calculation on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) momStepping = /* Momentum equation on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) vectorInvariantMomentum= /* Vector-Invariant Momentum on/off */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) momAdvection = /* Momentum advection on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) momViscosity = /* Momentum viscosity on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) momImplVertAdv= /* Momentum implicit vert. advection on/off*/ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicitViscosity = /* Implicit viscosity on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectImplicitDrag= /* Implicit bot Drag options (0,1,2)*/ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) 0= Expl. ; 1= Impl. on provis. Vel ; 2= Fully Impl (with surf.P) (PID.TID 0000.0001) ; (PID.TID 0000.0001) metricTerms = /* metric-Terms on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) useNHMTerms = /* Non-Hydrostatic Metric-Terms on/off */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectCoriMap = /* Coriolis Map options (0,1,2,3)*/ (PID.TID 0000.0001) 2 (PID.TID 0000.0001) 0= f-Plane ; 1= Beta-Plane ; 2= Spherical ; 3= read from file (PID.TID 0000.0001) ; (PID.TID 0000.0001) use3dCoriolis = /* 3-D Coriolis on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useCoriolis = /* Coriolis on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) useCDscheme = /* CD scheme on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectCoriScheme= /* Scheme selector for Coriolis-Term */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) = 0 : original discretization (simple averaging, no hFac) (PID.TID 0000.0001) = 1 : Wet-point averaging (Jamar & Ozer 1986) (PID.TID 0000.0001) = 2 : energy conserving scheme (no hFac weight) (PID.TID 0000.0001) = 3 : energy conserving scheme using Wet-point averaging (PID.TID 0000.0001) ; (PID.TID 0000.0001) momForcing = /* Momentum forcing on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) momTidalForcing = /* Momentum Tidal forcing on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) momPressureForcing = /* Momentum pressure term on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicitIntGravWave= /* Implicit Internal Gravity Wave flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) staggerTimeStep = /* Stagger time stepping on/off flag */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) doResetHFactors = /* reset thickness factors @ each time-step */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) multiDimAdvection = /* enable/disable Multi-Dim Advection */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) useMultiDimAdvec = /* Multi-Dim Advection is/is-not used */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) implicitDiffusion = /* Implicit Diffusion on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) tempStepping = /* Temperature equation on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) tempAdvection = /* Temperature advection on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) tempImplVertAdv = /* Temp. implicit vert. advection on/off */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) tempForcing = /* Temperature forcing on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) doThetaClimRelax = /* apply SST relaxation on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) tempIsActiveTr = /* Temp. is a dynamically Active Tracer */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) saltStepping = /* Salinity equation on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) saltAdvection = /* Salinity advection on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) saltImplVertAdv = /* Sali. implicit vert. advection on/off */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) saltForcing = /* Salinity forcing on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) doSaltClimRelax = /* apply SSS relaxation on/off flag */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) saltIsActiveTr = /* Salt is a dynamically Active Tracer */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) readBinaryPrec = /* Precision used for reading binary files */ (PID.TID 0000.0001) 32 (PID.TID 0000.0001) ; (PID.TID 0000.0001) writeBinaryPrec = /* Precision used for writing binary files */ (PID.TID 0000.0001) 32 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rwSuffixType = /* select format of mds file suffix */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) = 0 : myIter (I10.10) ; = 1 : 100*myTime (100th sec) ; (PID.TID 0000.0001) = 2 : myTime (seconds); = 3 : myTime/360 (10th of hr); (PID.TID 0000.0001) = 4 : myTime/3600 (hours) (PID.TID 0000.0001) ; (PID.TID 0000.0001) globalFiles = /* write "global" (=not per tile) files */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useSingleCpuIO = /* only master MPI process does I/O */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useSingleCpuInput = /* only master process reads input */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) /* debLev[*] : level of debug & auxiliary message printing */ (PID.TID 0000.0001) debLevZero = 0 ; /* level of disabled aux. msg printing */ (PID.TID 0000.0001) debLevA = 1 ; /* level of minimum aux. msg printing */ (PID.TID 0000.0001) debLevB = 2 ; /* level of low aux. print (report read-file opening)*/ (PID.TID 0000.0001) debLevC = 3 ; /* level of moderate debug prt (most pkgs debug msg) */ (PID.TID 0000.0001) debLevD = 4 ; /* level of enhanced debug prt (add DEBUG_STATS prt) */ (PID.TID 0000.0001) debLevE = 5 ; /* level of extensive debug printing */ (PID.TID 0000.0001) debugLevel = /* select debug printing level */ (PID.TID 0000.0001) 1 (PID.TID 0000.0001) ; (PID.TID 0000.0001) plotLevel = /* select PLOT_FIELD printing level */ (PID.TID 0000.0001) 1 (PID.TID 0000.0001) ; (PID.TID 0000.0001) // (PID.TID 0000.0001) // Elliptic solver(s) paramters ( PARM02 in namelist ) (PID.TID 0000.0001) // (PID.TID 0000.0001) cg2dMaxIters = /* Upper limit on 2d con. grad iterations */ (PID.TID 0000.0001) 1000 (PID.TID 0000.0001) ; (PID.TID 0000.0001) cg2dMinItersNSA = /* Minimum number of iterations of 2d con. grad solver */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) cg2dUseMinResSol= /* use cg2d last-iter(=0) / min-resid.(=1) solution */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) cg2dTargetResidual = /* 2d con. grad target residual */ (PID.TID 0000.0001) 1.000000000000000E-19 (PID.TID 0000.0001) ; (PID.TID 0000.0001) cg2dTargetResWunit = /* CG2d target residual [W units] */ (PID.TID 0000.0001) -1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) cg2dPreCondFreq = /* Freq. for updating cg2d preconditioner */ (PID.TID 0000.0001) 1 (PID.TID 0000.0001) ; (PID.TID 0000.0001) useSRCGSolver = /* use single reduction CG solver(s) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useNSACGSolver = /* use not-self-adjoint CG solver */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) printResidualFreq = /* Freq. for printing CG residual */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) // (PID.TID 0000.0001) // Time stepping paramters ( PARM03 in namelist ) (PID.TID 0000.0001) // (PID.TID 0000.0001) deltaTMom = /* Momentum equation timestep ( s ) */ (PID.TID 0000.0001) 1.200000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) deltaTFreeSurf = /* FreeSurface equation timestep ( s ) */ (PID.TID 0000.0001) 1.200000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) dTtracerLev = /* Tracer equation timestep ( s ) */ (PID.TID 0000.0001) 15 @ 4.320000000000000E+04 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) deltaTClock = /* Model clock timestep ( s ) */ (PID.TID 0000.0001) 4.320000000000000E+04 (PID.TID 0000.0001) ; (PID.TID 0000.0001) cAdjFreq = /* Convective adjustment interval ( s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) momForcingOutAB = /* =1: take Momentum Forcing out of Adams-Bash. stepping */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) tracForcingOutAB = /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) momDissip_In_AB = /* put Dissipation Tendency in Adams-Bash. stepping */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) doAB_onGtGs = /* apply AB on Tendencies (rather than on T,S)*/ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) abEps = /* Adams-Bashforth-2 stabilizing weight */ (PID.TID 0000.0001) 1.000000000000000E-01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) applyExchUV_early = /* Apply EXCH to U,V earlier in time-step */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) tauCD = /* CD coupling time-scale ( s ) */ (PID.TID 0000.0001) 3.214280000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rCD = /* Normalised CD coupling parameter */ (PID.TID 0000.0001) 9.962666600296178E-01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) epsAB_CD = /* AB-2 stabilizing weight for CD-scheme*/ (PID.TID 0000.0001) 1.000000000000000E-01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) pickupStrictlyMatch= /* stop if pickup do not strictly match */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) nIter0 = /* Run starting timestep number */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) nTimeSteps = /* Number of timesteps */ (PID.TID 0000.0001) 4 (PID.TID 0000.0001) ; (PID.TID 0000.0001) nEndIter = /* Run ending timestep number */ (PID.TID 0000.0001) 4 (PID.TID 0000.0001) ; (PID.TID 0000.0001) baseTime = /* Model base time ( s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) startTime = /* Run start time ( s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) endTime = /* Integration ending time ( s ) */ (PID.TID 0000.0001) 1.728000000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) pChkPtFreq = /* Permanent restart/pickup file interval ( s ) */ (PID.TID 0000.0001) 3.110400000000000E+08 (PID.TID 0000.0001) ; (PID.TID 0000.0001) chkPtFreq = /* Rolling restart/pickup file interval ( s ) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) pickup_write_mdsio = /* Model IO flag. */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) pickup_read_mdsio = /* Model IO flag. */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) writePickupAtEnd = /* Model IO flag. */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) dumpFreq = /* Model state write out interval ( s ). */ (PID.TID 0000.0001) 2.592000000000000E+06 (PID.TID 0000.0001) ; (PID.TID 0000.0001) dumpInitAndLast= /* write out Initial & Last iter. model state */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) snapshot_mdsio = /* Model IO flag. */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) monitorFreq = /* Monitor output interval ( s ). */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) monitorSelect = /* select group of variables to monitor */ (PID.TID 0000.0001) 3 (PID.TID 0000.0001) ; (PID.TID 0000.0001) monitor_stdio = /* Model IO flag. */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) externForcingPeriod = /* forcing period (s) */ (PID.TID 0000.0001) 2.592000000000000E+06 (PID.TID 0000.0001) ; (PID.TID 0000.0001) externForcingCycle = /* period of the cyle (s). */ (PID.TID 0000.0001) 3.110400000000000E+07 (PID.TID 0000.0001) ; (PID.TID 0000.0001) tauThetaClimRelax = /* relaxation time scale (s) */ (PID.TID 0000.0001) 5.184000000000000E+06 (PID.TID 0000.0001) ; (PID.TID 0000.0001) tauSaltClimRelax = /* relaxation time scale (s) */ (PID.TID 0000.0001) 1.555200000000000E+07 (PID.TID 0000.0001) ; (PID.TID 0000.0001) latBandClimRelax = /* max. Lat. where relaxation */ (PID.TID 0000.0001) 1.800000000000000E+02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) // (PID.TID 0000.0001) // Gridding paramters ( PARM04 in namelist ) (PID.TID 0000.0001) // (PID.TID 0000.0001) usingCartesianGrid = /* Cartesian coordinates flag ( True/False ) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) usingCylindricalGrid = /* Cylindrical coordinates flag ( True/False ) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) usingSphericalPolarGrid = /* Spherical coordinates flag ( True/False ) */ (PID.TID 0000.0001) T (PID.TID 0000.0001) ; (PID.TID 0000.0001) usingCurvilinearGrid = /* Curvilinear coordinates flag ( True/False ) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) useMin4hFacEdges = /* set hFacW,S as minimum of adjacent hFacC factor */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) interViscAr_pCell = /* account for partial-cell in interior vert. viscosity */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) interDiffKr_pCell = /* account for partial-cell in interior vert. diffusion */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) pCellMix_select = /* option to enhance mixing near surface & bottom */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) selectSigmaCoord = /* Hybrid-Sigma Vert. Coordinate option */ (PID.TID 0000.0001) 0 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rSigmaBnd = /* r/sigma transition ( units of r == m ) */ (PID.TID 0000.0001) 1.234567000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rkSign = /* index orientation relative to vertical coordinate */ (PID.TID 0000.0001) -1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) gravitySign = /* gravity orientation relative to vertical coordinate */ (PID.TID 0000.0001) -1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) seaLev_Z = /* reference height of sea-level [m] */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) top_Pres = /* reference pressure at the top [Pa] */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) mass2rUnit = /* convert mass per unit area [kg/m2] to r-units [m] */ (PID.TID 0000.0001) 9.661835748792270E-04 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rUnit2mass = /* convert r-units [m] to mass per unit area [kg/m2] */ (PID.TID 0000.0001) 1.035000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) drC = /* C spacing ( units of r ) */ (PID.TID 0000.0001) 2.500000000000000E+01, /* K = 1 */ (PID.TID 0000.0001) 6.000000000000000E+01, /* K = 2 */ (PID.TID 0000.0001) 8.500000000000000E+01, /* K = 3 */ (PID.TID 0000.0001) 1.200000000000000E+02, /* K = 4 */ (PID.TID 0000.0001) 1.650000000000000E+02, /* K = 5 */ (PID.TID 0000.0001) 2.150000000000000E+02, /* K = 6 */ (PID.TID 0000.0001) 2.650000000000000E+02, /* K = 7 */ (PID.TID 0000.0001) 3.150000000000000E+02, /* K = 8 */ (PID.TID 0000.0001) 3.650000000000000E+02, /* K = 9 */ (PID.TID 0000.0001) 4.150000000000000E+02, /* K = 10 */ (PID.TID 0000.0001) 4.650000000000000E+02, /* K = 11 */ (PID.TID 0000.0001) 5.150000000000000E+02, /* K = 12 */ (PID.TID 0000.0001) 5.650000000000000E+02, /* K = 13 */ (PID.TID 0000.0001) 6.150000000000000E+02, /* K = 14 */ (PID.TID 0000.0001) 6.650000000000000E+02, /* K = 15 */ (PID.TID 0000.0001) 3.450000000000000E+02 /* K = 16 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) drF = /* W spacing ( units of r ) */ (PID.TID 0000.0001) 5.000000000000000E+01, /* K = 1 */ (PID.TID 0000.0001) 7.000000000000000E+01, /* K = 2 */ (PID.TID 0000.0001) 1.000000000000000E+02, /* K = 3 */ (PID.TID 0000.0001) 1.400000000000000E+02, /* K = 4 */ (PID.TID 0000.0001) 1.900000000000000E+02, /* K = 5 */ (PID.TID 0000.0001) 2.400000000000000E+02, /* K = 6 */ (PID.TID 0000.0001) 2.900000000000000E+02, /* K = 7 */ (PID.TID 0000.0001) 3.400000000000000E+02, /* K = 8 */ (PID.TID 0000.0001) 3.900000000000000E+02, /* K = 9 */ (PID.TID 0000.0001) 4.400000000000000E+02, /* K = 10 */ (PID.TID 0000.0001) 4.900000000000000E+02, /* K = 11 */ (PID.TID 0000.0001) 5.400000000000000E+02, /* K = 12 */ (PID.TID 0000.0001) 5.900000000000000E+02, /* K = 13 */ (PID.TID 0000.0001) 6.400000000000000E+02, /* K = 14 */ (PID.TID 0000.0001) 6.900000000000000E+02 /* K = 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) delX = /* U spacing ( m - cartesian, degrees - spherical ) */ (PID.TID 0000.0001) 90 @ 4.000000000000000E+00 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) delY = /* V spacing ( m - cartesian, degrees - spherical ) */ (PID.TID 0000.0001) 40 @ 4.000000000000000E+00 /* J = 1: 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) xgOrigin = /* X-axis origin of West edge (cartesian: m, lat-lon: deg) */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) ygOrigin = /* Y-axis origin of South edge (cartesian: m, lat-lon: deg) */ (PID.TID 0000.0001) -8.000000000000000E+01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) rSphere = /* Radius ( ignored - cartesian, m - spherical ) */ (PID.TID 0000.0001) 6.370000000000000E+06 (PID.TID 0000.0001) ; (PID.TID 0000.0001) deepAtmosphere = /* Deep/Shallow Atmosphere flag (True/False) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) xC = /* xC(:,1,:,1) : P-point X coord ( deg. or m if cartesian) */ (PID.TID 0000.0001) 2.000000000000000E+00, /* I = 1 */ (PID.TID 0000.0001) 6.000000000000000E+00, /* I = 2 */ (PID.TID 0000.0001) 1.000000000000000E+01, /* I = 3 */ (PID.TID 0000.0001) . . . (PID.TID 0000.0001) 8.200000000000000E+01, /* I = 21 */ (PID.TID 0000.0001) 8.600000000000000E+01, /* I = 22 */ (PID.TID 0000.0001) 9.000000000000000E+01, /* I = 23 */ (PID.TID 0000.0001) 9.400000000000000E+01, /* I = 24 */ (PID.TID 0000.0001) 9.800000000000000E+01, /* I = 25 */ (PID.TID 0000.0001) 1.020000000000000E+02, /* I = 26 */ (PID.TID 0000.0001) . . . (PID.TID 0000.0001) 1.700000000000000E+02, /* I = 43 */ (PID.TID 0000.0001) 1.740000000000000E+02, /* I = 44 */ (PID.TID 0000.0001) 1.780000000000000E+02, /* I = 45 */ (PID.TID 0000.0001) 1.820000000000000E+02, /* I = 46 */ (PID.TID 0000.0001) 1.860000000000000E+02, /* I = 47 */ (PID.TID 0000.0001) 1.900000000000000E+02, /* I = 48 */ (PID.TID 0000.0001) . . . (PID.TID 0000.0001) 2.580000000000000E+02, /* I = 65 */ (PID.TID 0000.0001) 2.620000000000000E+02, /* I = 66 */ (PID.TID 0000.0001) 2.660000000000000E+02, /* I = 67 */ (PID.TID 0000.0001) 2.700000000000000E+02, /* I = 68 */ (PID.TID 0000.0001) 2.740000000000000E+02, /* I = 69 */ (PID.TID 0000.0001) 2.780000000000000E+02, /* I = 70 */ (PID.TID 0000.0001) . . . (PID.TID 0000.0001) 3.500000000000000E+02, /* I = 88 */ (PID.TID 0000.0001) 3.540000000000000E+02, /* I = 89 */ (PID.TID 0000.0001) 3.580000000000000E+02 /* I = 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) yC = /* yC(1,:,1,:) : P-point Y coord ( deg. or m if cartesian) */ (PID.TID 0000.0001) -7.800000000000000E+01, /* J = 1 */ (PID.TID 0000.0001) -7.400000000000000E+01, /* J = 2 */ (PID.TID 0000.0001) -7.000000000000000E+01, /* J = 3 */ (PID.TID 0000.0001) -6.600000000000000E+01, /* J = 4 */ (PID.TID 0000.0001) -6.200000000000000E+01, /* J = 5 */ (PID.TID 0000.0001) -5.800000000000000E+01, /* J = 6 */ (PID.TID 0000.0001) -5.400000000000000E+01, /* J = 7 */ (PID.TID 0000.0001) -5.000000000000000E+01, /* J = 8 */ (PID.TID 0000.0001) -4.600000000000000E+01, /* J = 9 */ (PID.TID 0000.0001) -4.200000000000000E+01, /* J = 10 */ (PID.TID 0000.0001) -3.800000000000000E+01, /* J = 11 */ (PID.TID 0000.0001) -3.400000000000000E+01, /* J = 12 */ (PID.TID 0000.0001) -3.000000000000000E+01, /* J = 13 */ (PID.TID 0000.0001) -2.600000000000000E+01, /* J = 14 */ (PID.TID 0000.0001) -2.200000000000000E+01, /* J = 15 */ (PID.TID 0000.0001) -1.800000000000000E+01, /* J = 16 */ (PID.TID 0000.0001) -1.400000000000000E+01, /* J = 17 */ (PID.TID 0000.0001) -1.000000000000000E+01, /* J = 18 */ (PID.TID 0000.0001) -6.000000000000000E+00, /* J = 19 */ (PID.TID 0000.0001) -2.000000000000000E+00, /* J = 20 */ (PID.TID 0000.0001) 2.000000000000000E+00, /* J = 21 */ (PID.TID 0000.0001) 6.000000000000000E+00, /* J = 22 */ (PID.TID 0000.0001) 1.000000000000000E+01, /* J = 23 */ (PID.TID 0000.0001) 1.400000000000000E+01, /* J = 24 */ (PID.TID 0000.0001) 1.800000000000000E+01, /* J = 25 */ (PID.TID 0000.0001) 2.200000000000000E+01, /* J = 26 */ (PID.TID 0000.0001) 2.600000000000000E+01, /* J = 27 */ (PID.TID 0000.0001) 3.000000000000000E+01, /* J = 28 */ (PID.TID 0000.0001) 3.400000000000000E+01, /* J = 29 */ (PID.TID 0000.0001) 3.800000000000000E+01, /* J = 30 */ (PID.TID 0000.0001) 4.200000000000000E+01, /* J = 31 */ (PID.TID 0000.0001) 4.600000000000000E+01, /* J = 32 */ (PID.TID 0000.0001) 5.000000000000000E+01, /* J = 33 */ (PID.TID 0000.0001) 5.400000000000000E+01, /* J = 34 */ (PID.TID 0000.0001) 5.800000000000000E+01, /* J = 35 */ (PID.TID 0000.0001) 6.200000000000000E+01, /* J = 36 */ (PID.TID 0000.0001) 6.600000000000000E+01, /* J = 37 */ (PID.TID 0000.0001) 7.000000000000000E+01, /* J = 38 */ (PID.TID 0000.0001) 7.400000000000000E+01, /* J = 39 */ (PID.TID 0000.0001) 7.800000000000000E+01 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rcoord = /* P-point R coordinate ( units of r ) */ (PID.TID 0000.0001) -2.500000000000000E+01, /* K = 1 */ (PID.TID 0000.0001) -8.500000000000000E+01, /* K = 2 */ (PID.TID 0000.0001) -1.700000000000000E+02, /* K = 3 */ (PID.TID 0000.0001) -2.900000000000000E+02, /* K = 4 */ (PID.TID 0000.0001) -4.550000000000000E+02, /* K = 5 */ (PID.TID 0000.0001) -6.700000000000000E+02, /* K = 6 */ (PID.TID 0000.0001) -9.350000000000000E+02, /* K = 7 */ (PID.TID 0000.0001) -1.250000000000000E+03, /* K = 8 */ (PID.TID 0000.0001) -1.615000000000000E+03, /* K = 9 */ (PID.TID 0000.0001) -2.030000000000000E+03, /* K = 10 */ (PID.TID 0000.0001) -2.495000000000000E+03, /* K = 11 */ (PID.TID 0000.0001) -3.010000000000000E+03, /* K = 12 */ (PID.TID 0000.0001) -3.575000000000000E+03, /* K = 13 */ (PID.TID 0000.0001) -4.190000000000000E+03, /* K = 14 */ (PID.TID 0000.0001) -4.855000000000000E+03 /* K = 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rF = /* W-Interf. R coordinate ( units of r ) */ (PID.TID 0000.0001) 0.000000000000000E+00, /* K = 1 */ (PID.TID 0000.0001) -5.000000000000000E+01, /* K = 2 */ (PID.TID 0000.0001) -1.200000000000000E+02, /* K = 3 */ (PID.TID 0000.0001) -2.200000000000000E+02, /* K = 4 */ (PID.TID 0000.0001) -3.600000000000000E+02, /* K = 5 */ (PID.TID 0000.0001) -5.500000000000000E+02, /* K = 6 */ (PID.TID 0000.0001) -7.900000000000000E+02, /* K = 7 */ (PID.TID 0000.0001) -1.080000000000000E+03, /* K = 8 */ (PID.TID 0000.0001) -1.420000000000000E+03, /* K = 9 */ (PID.TID 0000.0001) -1.810000000000000E+03, /* K = 10 */ (PID.TID 0000.0001) -2.250000000000000E+03, /* K = 11 */ (PID.TID 0000.0001) -2.740000000000000E+03, /* K = 12 */ (PID.TID 0000.0001) -3.280000000000000E+03, /* K = 13 */ (PID.TID 0000.0001) -3.870000000000000E+03, /* K = 14 */ (PID.TID 0000.0001) -4.510000000000000E+03, /* K = 15 */ (PID.TID 0000.0001) -5.200000000000000E+03 /* K = 16 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) deepFacC = /* deep-model grid factor @ cell-Center (-) */ (PID.TID 0000.0001) 15 @ 1.000000000000000E+00 /* K = 1: 15 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) deepFacF = /* deep-model grid factor @ W-Interface (-) */ (PID.TID 0000.0001) 16 @ 1.000000000000000E+00 /* K = 1: 16 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rotateGrid = /* use rotated grid ( True/False ) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) phiEuler = /* Euler angle, rotation about original z-coordinate [rad] */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) thetaEuler = /* Euler angle, rotation about new x-coordinate [rad] */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) psiEuler = /* Euler angle, rotation about new z-coordinate [rad] */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxF = /* dxF(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxF = /* dxF(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */ (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */ (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */ (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */ (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */ (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */ (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */ (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */ (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */ (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */ (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */ (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */ (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */ (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */ (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */ (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */ (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */ (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */ (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */ (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */ (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */ (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */ (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */ (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */ (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */ (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */ (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */ (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */ (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */ (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */ (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */ (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */ (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */ (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */ (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */ (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */ (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */ (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */ (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyF = /* dyF(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyF = /* dyF(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxG = /* dxG(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxG = /* dxG(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */ (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */ (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */ (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */ (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */ (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */ (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */ (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */ (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */ (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */ (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */ (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */ (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */ (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */ (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */ (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */ (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */ (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */ (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */ (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */ (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */ (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */ (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */ (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */ (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */ (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */ (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */ (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */ (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */ (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */ (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */ (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */ (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */ (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */ (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */ (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */ (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */ (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */ (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */ (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyG = /* dyG(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyG = /* dyG(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxC = /* dxC(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 9.246038586187513E+04 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxC = /* dxC(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 9.246038586187513E+04, /* J = 1 */ (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 2 */ (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 3 */ (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 4 */ (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 5 */ (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 6 */ (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 7 */ (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 8 */ (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 9 */ (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 10 */ (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 11 */ (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 12 */ (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 13 */ (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 14 */ (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 15 */ (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 16 */ (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 17 */ (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 18 */ (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 19 */ (PID.TID 0000.0001) 2 @ 4.444389881567502E+05, /* J = 20: 21 */ (PID.TID 0000.0001) 4.422737260813012E+05, /* J = 22 */ (PID.TID 0000.0001) 4.379537508695838E+05, /* J = 23 */ (PID.TID 0000.0001) 4.315001090065870E+05, /* J = 24 */ (PID.TID 0000.0001) 4.229442419867491E+05, /* J = 25 */ (PID.TID 0000.0001) 4.123278331341585E+05, /* J = 26 */ (PID.TID 0000.0001) 3.997026045255871E+05, /* J = 27 */ (PID.TID 0000.0001) 3.851300650057323E+05, /* J = 28 */ (PID.TID 0000.0001) 3.686812105223082E+05, /* J = 29 */ (PID.TID 0000.0001) 3.504361782409254E+05, /* J = 30 */ (PID.TID 0000.0001) 3.304838561248741E+05, /* J = 31 */ (PID.TID 0000.0001) 3.089214498819034E+05, /* J = 32 */ (PID.TID 0000.0001) 2.858540093877838E+05, /* J = 33 */ (PID.TID 0000.0001) 2.613939168938713E+05, /* J = 34 */ (PID.TID 0000.0001) 2.356603395120765E+05, /* J = 35 */ (PID.TID 0000.0001) 2.087786486446736E+05, /* J = 36 */ (PID.TID 0000.0001) 1.808798091874300E+05, /* J = 37 */ (PID.TID 0000.0001) 1.520997414818001E+05, /* J = 38 */ (PID.TID 0000.0001) 1.225786591246834E+05, /* J = 39 */ (PID.TID 0000.0001) 9.246038586187513E+04 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyC = /* dyC(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyC = /* dyC(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxV = /* dxV(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 7.722306258078101E+04 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dxV = /* dxV(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 7.722306258078101E+04, /* J = 1 */ (PID.TID 0000.0001) 1.075850604052271E+05, /* J = 2 */ (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 3 */ (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 4 */ (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 5 */ (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 6 */ (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 7 */ (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 8 */ (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 9 */ (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 10 */ (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 11 */ (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 12 */ (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 13 */ (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 14 */ (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 15 */ (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 16 */ (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 17 */ (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 18 */ (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 19 */ (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 20 */ (PID.TID 0000.0001) 4.447098934081552E+05, /* J = 21 */ (PID.TID 0000.0001) 4.436266024588540E+05, /* J = 22 */ (PID.TID 0000.0001) 4.403820072955634E+05, /* J = 23 */ (PID.TID 0000.0001) 4.349919152597734E+05, /* J = 24 */ (PID.TID 0000.0001) 4.274825863380723E+05, /* J = 25 */ (PID.TID 0000.0001) 4.178906052261313E+05, /* J = 26 */ (PID.TID 0000.0001) 4.062627030917454E+05, /* J = 27 */ (PID.TID 0000.0001) 3.926555299052806E+05, /* J = 28 */ (PID.TID 0000.0001) 3.771353784467131E+05, /* J = 29 */ (PID.TID 0000.0001) 3.597778613338690E+05, /* J = 30 */ (PID.TID 0000.0001) 3.406675426453503E+05, /* J = 31 */ (PID.TID 0000.0001) 3.198975259328452E+05, /* J = 32 */ (PID.TID 0000.0001) 2.975690006299821E+05, /* J = 33 */ (PID.TID 0000.0001) 2.737907490675810E+05, /* J = 34 */ (PID.TID 0000.0001) 2.486786164970726E+05, /* J = 35 */ (PID.TID 0000.0001) 2.223549467040777E+05, /* J = 36 */ (PID.TID 0000.0001) 1.949479859617814E+05, /* J = 37 */ (PID.TID 0000.0001) 1.665912582279823E+05, /* J = 38 */ (PID.TID 0000.0001) 1.374229146297914E+05, /* J = 39 */ (PID.TID 0000.0001) 1.075850604052271E+05 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyU = /* dyU(:,1,:,1) ( units: m ) */ (PID.TID 0000.0001) 90 @ 4.447098934081552E+05 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) dyU = /* dyU(1,:,1,:) ( units: m ) */ (PID.TID 0000.0001) 40 @ 4.447098934081552E+05 /* J = 1: 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rA = /* rA (:,1,:,1) ( units: m^2 ) */ (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rA = /* rA (1,:,1,:) ( units: m^2 ) */ (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */ (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 2 */ (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */ (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 4 */ (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 5 */ (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */ (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */ (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */ (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 9 */ (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */ (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */ (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */ (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 13 */ (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */ (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */ (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */ (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */ (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 18 */ (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */ (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */ (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */ (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 23 */ (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */ (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */ (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */ (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */ (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 28 */ (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */ (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */ (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */ (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 32 */ (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */ (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */ (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */ (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 36 */ (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 37 */ (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */ (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 39 */ (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rAw = /* rAw(:,1,:,1) ( units: m^2 ) */ (PID.TID 0000.0001) 90 @ 4.110969866729047E+10 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rAw = /* rAw(1,:,1,:) ( units: m^2 ) */ (PID.TID 0000.0001) 4.110969866729047E+10, /* J = 1 */ (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 2 */ (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 3 */ (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 4 */ (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 5 */ (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 6 */ (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 7 */ (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 8 */ (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 9 */ (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 10 */ (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 11 */ (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 12 */ (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 13 */ (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 14 */ (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 15 */ (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 16 */ (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 17 */ (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 18 */ (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 19 */ (PID.TID 0000.0001) 2 @ 1.976062798008856E+11, /* J = 20: 21 */ (PID.TID 0000.0001) 1.966435618690094E+11, /* J = 22 */ (PID.TID 0000.0001) 1.947228162702493E+11, /* J = 23 */ (PID.TID 0000.0001) 1.918534006840893E+11, /* J = 24 */ (PID.TID 0000.0001) 1.880492946148330E+11, /* J = 25 */ (PID.TID 0000.0001) 1.833290312848623E+11, /* J = 26 */ (PID.TID 0000.0001) 1.777156073426415E+11, /* J = 27 */ (PID.TID 0000.0001) 1.712363708253574E+11, /* J = 28 */ (PID.TID 0000.0001) 1.639228879220326E+11, /* J = 29 */ (PID.TID 0000.0001) 1.558107891862220E+11, /* J = 30 */ (PID.TID 0000.0001) 1.469395959475426E+11, /* J = 31 */ (PID.TID 0000.0001) 1.373525277677230E+11, /* J = 32 */ (PID.TID 0000.0001) 1.270962918792468E+11, /* J = 33 */ (PID.TID 0000.0001) 1.162208556324093E+11, /* J = 34 */ (PID.TID 0000.0001) 1.047792030594103E+11, /* J = 35 */ (PID.TID 0000.0001) 9.282707674147525E+10, /* J = 36 */ (PID.TID 0000.0001) 8.042270623659996E+10, /* J = 37 */ (PID.TID 0000.0001) 6.762652439100235E+10, /* J = 38 */ (PID.TID 0000.0001) 5.450087291636666E+10, /* J = 39 */ (PID.TID 0000.0001) 4.110969866729047E+10 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rAs = /* rAs(:,1,:,1) ( units: m^2 ) */ (PID.TID 0000.0001) 90 @ 3.433488626798251E+10 /* I = 1: 90 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) rAs = /* rAs(1,:,1,:) ( units: m^2 ) */ (PID.TID 0000.0001) 3.433488626798251E+10, /* J = 1 */ (PID.TID 0000.0001) 4.783442523123625E+10, /* J = 2 */ (PID.TID 0000.0001) 6.110091968306415E+10, /* J = 3 */ (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 4 */ (PID.TID 0000.0001) 8.667769319778081E+10, /* J = 5 */ (PID.TID 0000.0001) 9.886336479107465E+10, /* J = 6 */ (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 7 */ (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 8 */ (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 9 */ (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 10 */ (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 11 */ (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 12 */ (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 13 */ (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 14 */ (PID.TID 0000.0001) 1.806323556642997E+11, /* J = 15 */ (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 16 */ (PID.TID 0000.0001) 1.900671314104744E+11, /* J = 17 */ (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 18 */ (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 19 */ (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 20 */ (PID.TID 0000.0001) 1.977267295821496E+11, /* J = 21 */ (PID.TID 0000.0001) 1.972450772065981E+11, /* J = 22 */ (PID.TID 0000.0001) 1.958024666419019E+11, /* J = 23 */ (PID.TID 0000.0001) 1.934059261417216E+11, /* J = 24 */ (PID.TID 0000.0001) 1.900671314104744E+11, /* J = 25 */ (PID.TID 0000.0001) 1.858023487204767E+11, /* J = 26 */ (PID.TID 0000.0001) 1.806323556642997E+11, /* J = 27 */ (PID.TID 0000.0001) 1.745823399284268E+11, /* J = 28 */ (PID.TID 0000.0001) 1.676817765813788E+11, /* J = 29 */ (PID.TID 0000.0001) 1.599642844741385E+11, /* J = 30 */ (PID.TID 0000.0001) 1.514674624524945E+11, /* J = 31 */ (PID.TID 0000.0001) 1.422327061792377E+11, /* J = 32 */ (PID.TID 0000.0001) 1.323050064586578E+11, /* J = 33 */ (PID.TID 0000.0001) 1.217327300458638E+11, /* J = 34 */ (PID.TID 0000.0001) 1.105673840088173E+11, /* J = 35 */ (PID.TID 0000.0001) 9.886336479107465E+10, /* J = 36 */ (PID.TID 0000.0001) 8.667769319778081E+10, /* J = 37 */ (PID.TID 0000.0001) 7.406973659603818E+10, /* J = 38 */ (PID.TID 0000.0001) 6.110091968306415E+10, /* J = 39 */ (PID.TID 0000.0001) 4.783442523123625E+10 /* J = 40 */ (PID.TID 0000.0001) ; (PID.TID 0000.0001) globalArea = /* Integrated horizontal Area (m^2) */ (PID.TID 0000.0001) 3.450614146649749E+14 (PID.TID 0000.0001) ; (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // End of Model config. summary (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) == Packages configuration : Check & print summary == (PID.TID 0000.0001) (PID.TID 0000.0001) GMREDI_CHECK: #define GMREDI (PID.TID 0000.0001) GM_AdvForm = /* if FALSE => use SkewFlux Form */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_InMomAsStress = /* if TRUE => apply as Eddy Stress */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_AdvSeparate = /* Calc Bolus & Euler Adv. separately */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_ExtraDiag = /* Tensor Extra Diag (line 1&2) non 0 */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_isopycK = /* Background Isopyc. Diffusivity [m^2/s] */ (PID.TID 0000.0001) 1.000000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_skewflx*K = /* Background GM_SkewFlx Diffusivity [m^2/s] */ (PID.TID 0000.0001) 1.000000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_advec*K = /* Backg. GM-Advec(=Bolus) Diffusivity [m^2/s]*/ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_Kmin_horiz = /* Minimum Horizontal Diffusivity [m^2/s] */ (PID.TID 0000.0001) 5.000000000000000E+01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_Visbeck_alpha = /* Visbeck alpha coeff. [-] */ (PID.TID 0000.0001) 0.000000000000000E+00 (PID.TID 0000.0001) ; OAD: TIMING: stamp 0: 1690340614.362698 (PID.TID 0000.0001) GM_Small_Number = /* epsilon used in slope calc */ (PID.TID 0000.0001) 1.000000000000000E-12 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_slopeSqCutoff = /* Slope^2 cut-off value */ (PID.TID 0000.0001) 1.000000000000000E+08 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_taper_scheme = /* Type of Tapering/Clipping scheme */ (PID.TID 0000.0001) 'dm95 ' (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_maxSlope = /* Maximum Slope (Tapering/Clipping) */ (PID.TID 0000.0001) 1.000000000000000E-02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_facTrL2dz = /* Minimum Trans.Layer Thick. (factor of dz) */ (PID.TID 0000.0001) 1.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_facTrL2ML = /* Max.Trans.Layer Thick. (factor of MxL Depth)*/ (PID.TID 0000.0001) 5.000000000000000E+00 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_maxTransLay = /* Maximum Transition Layer Thickness [m] */ (PID.TID 0000.0001) 5.000000000000000E+02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_UseBVP = /* if TRUE => use bvp a la Ferrari et al. (2010) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_BVP_ModeNumber = /* Vertical mode number for BVP wave speed */ (PID.TID 0000.0001) 1 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_BVP_cMin = /* Minimum wave speed for BVP [m/s] */ (PID.TID 0000.0001) 1.000000000000000E-01 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_useSubMeso = /* if TRUE => use Sub-Meso param. (B.Fox-Kemper) */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) subMeso_Ceff = /* efficiency coeff. of Mixed-Layer Eddies [-] */ (PID.TID 0000.0001) 7.000000000000001E-02 (PID.TID 0000.0001) ; (PID.TID 0000.0001) subMeso_invTau = /* inverse of Sub-Meso mixing time-scale [/s] */ (PID.TID 0000.0001) 2.000000000000000E-06 (PID.TID 0000.0001) ; (PID.TID 0000.0001) subMeso_LfMin = /* minimum length-scale "Lf" [m] */ (PID.TID 0000.0001) 1.000000000000000E+03 (PID.TID 0000.0001) ; (PID.TID 0000.0001) subMeso_Lmax = /* maximum grid-scale length [m] */ (PID.TID 0000.0001) 1.100000000000000E+05 (PID.TID 0000.0001) ; (PID.TID 0000.0001) GM_useLeithQG = /* if TRUE => add QG Leith viscosity to GMRedi tensor */ (PID.TID 0000.0001) F (PID.TID 0000.0001) ; (PID.TID 0000.0001) CTRL_CHECK: --> Starts to check CTRL set-up (PID.TID 0000.0001) CTRL_CHECK: <-- Ends Normally (PID.TID 0000.0001) (PID.TID 0000.0001) COST_CHECK: #define ALLOW_COST (PID.TID 0000.0001) GRDCHK_CHECK: grdchk package (PID.TID 0000.0001) GAD_CHECK: #define ALLOW_GENERIC_ADVDIFF (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Check Model config. (CONFIG_CHECK): (PID.TID 0000.0001) // CONFIG_CHECK : Normal End (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) OAD: IT+ 1048576 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -5.21804821573824E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -7.43849426498855E-15 6.29008403763954E+00 OAD: DT+ 1048576 OAD: IT+ 17825792 OAD: IT+ 34603008 OAD: DT+ 17825792 OAD: IT+ 51380224 OAD: IT+ 68157440 cg2d: Sum(rhs),rhsMax = -3.55271367880050E-15 6.70003306073569E+00 OAD: TIMING: stamp 1: 1690340618.338944 OAD: TIMING: delta stamps 1-0: 3.976246 OAD: TIMING: stamp 2: 1690340624.350184 OAD: TIMING: delta stamps 2-1: 6.011240 OAD: TIMING: delta stamps 2-0: 9.987486 OAD: TIMING: ratio stamps (2-1)/(1-0): 6.011240e+06/3.976246e+06=1.511788e+00 (PID.TID 0000.0001) %CHECKPOINT 4 ckptA (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003975202165D+06 (PID.TID 0000.0001) local fc = 0.804687744759400D+06 (PID.TID 0000.0001) global fc = 0.804687744759400D+06 cg2d: Sum(rhs),rhsMax = 0.00000000000000E+00 0.00000000000000E+00 cg2d: Sum(rhs),rhsMax = 7.88258347483861E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -3.10862446895044E-15 7.86983807726857E-04 cg2d: Sum(rhs),rhsMax = 1.05471187339390E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -5.32907051820075E-15 1.33029393421496E-03 EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = 1.16573417585641E-14 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 5.32907051820075E-15 1.40191055720658E-03 ph-pack: packing ecco_cost ph-pack: packing ecco_ctrl (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Gradient-check starts (grdchk_main) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) grdchk reference fc: fcref = 8.04687744759400E+05 grad-res ------------------------------- grad-res proc # i j k bi bj iobc fc ref fc + eps fc - eps grad-res proc # i j k bi bj iobc adj grad fd grad 1 - fd/adj grad-res exact position met: grad-res 0 526 26 17 1 2 2 (PID.TID 0000.0001) ====== Starts gradient-check number 1 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 526 29309 1 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 0 1 ph-grd -->hit<-- 26 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 26 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -9.10382880192628E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -2.72004641033163E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -2.55351295663786E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = 3.55271367880050E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003997193242D+06 (PID.TID 0000.0001) local fc = 0.804687766750477D+06 (PID.TID 0000.0001) global fc = 0.804687766750477D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687766750477E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -9.10382880192628E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -6.77236045021345E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 9.10382880192628E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -1.99840144432528E-14 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003953394813D+06 (PID.TID 0000.0001) local fc = 0.804687722952049D+06 (PID.TID 0000.0001) global fc = 0.804687722952049D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687722952049E+05 grad-res ------------------------------- grad-res 0 1 26 17 1 2 2 1 8.04687744759E+05 8.04687766750E+05 8.04687722952E+05 grad-res 0 1 1 526 0 2 2 1 2.18992150034E+00 2.18992143637E+00 2.92152290138E-08 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 2.18992150034475E+00 (PID.TID 0000.0001) ADM finite-diff_grad = 2.18992143636569E+00 (PID.TID 0000.0001) ====== End of gradient-check number 1 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 2 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 527 29309 2 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 526 2 ph-grd -->hit<-- 30 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 30 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -3.49720252756924E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 4.32986979603811E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = 2.66453525910038E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003984692549D+06 (PID.TID 0000.0001) local fc = 0.804687754249784D+06 (PID.TID 0000.0001) global fc = 0.804687754249784D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687754249784E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -9.10382880192628E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 7.10542735760100E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -7.66053886991358E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -5.77315972805081E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003965894237D+06 (PID.TID 0000.0001) local fc = 0.804687735451472D+06 (PID.TID 0000.0001) global fc = 0.804687735451472D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687735451472E+05 grad-res ------------------------------- grad-res 0 2 30 17 1 2 2 1 8.04687744759E+05 8.04687754250E+05 8.04687735451E+05 grad-res 0 2 2 527 0 2 2 1 9.39915722870E-01 9.39915596973E-01 1.33944963232E-07 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 9.39915722870128E-01 (PID.TID 0000.0001) ADM finite-diff_grad = 9.39915596973151E-01 (PID.TID 0000.0001) ====== End of gradient-check number 2 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 3 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 528 29309 3 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 527 3 ph-grd -->hit<-- 31 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 31 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -8.32667268468867E-16 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 6.43929354282591E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = 6.66133814775094E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003980663090D+06 (PID.TID 0000.0001) local fc = 0.804687750220325D+06 (PID.TID 0000.0001) global fc = 0.804687750220325D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687750220325E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.65973959207622E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 1.22124532708767E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 4.66293670342566E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -2.08721928629529E-14 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003969921301D+06 (PID.TID 0000.0001) local fc = 0.804687739478537D+06 (PID.TID 0000.0001) global fc = 0.804687739478537D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687739478537E+05 grad-res ------------------------------- grad-res 0 3 31 17 1 2 2 1 8.04687744759E+05 8.04687750220E+05 8.04687739479E+05 grad-res 0 3 3 528 0 2 2 1 5.37089397045E-01 5.37089433055E-01 -6.70464879171E-08 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 5.37089397045413E-01 (PID.TID 0000.0001) ADM finite-diff_grad = 5.37089433055371E-01 (PID.TID 0000.0001) ====== End of gradient-check number 3 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 4 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 529 29309 4 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 528 4 ph-grd -->hit<-- 32 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 32 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.65973959207622E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -6.49480469405717E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -2.99760216648792E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -2.66453525910038E-14 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003997321381D+06 (PID.TID 0000.0001) local fc = 0.804687766878617D+06 (PID.TID 0000.0001) global fc = 0.804687766878617D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687766878617E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 3.38618022510673E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 5.77315972805081E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -1.33226762955019E-14 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003953270627D+06 (PID.TID 0000.0001) local fc = 0.804687722827862D+06 (PID.TID 0000.0001) global fc = 0.804687722827862D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687722827862E+05 grad-res ------------------------------- grad-res 0 4 32 17 1 2 2 1 8.04687744759E+05 8.04687766879E+05 8.04687722828E+05 grad-res 0 4 4 529 0 2 2 1 2.20253900008E+00 2.20253771986E+00 5.81250153009E-07 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 2.20253900008494E+00 (PID.TID 0000.0001) ADM finite-diff_grad = 2.20253771985881E+00 (PID.TID 0000.0001) ====== End of gradient-check number 4 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 5 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 530 29309 5 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 529 5 ph-grd -->hit<-- 36 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 36 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.65973959207622E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 8.88178419700125E-16 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 5.10702591327572E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -6.66133814775094E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205004049861352D+06 (PID.TID 0000.0001) local fc = 0.804687819418588D+06 (PID.TID 0000.0001) global fc = 0.804687819418588D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687819418588E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 3.71924713249427E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -3.33066907387547E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = 4.44089209850063E-16 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089993D+06 --> objf_test(bi,bj) = 0.205003900726674D+06 (PID.TID 0000.0001) local fc = 0.804687670283909D+06 (PID.TID 0000.0001) global fc = 0.804687670283909D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687670283909E+05 grad-res ------------------------------- grad-res 0 5 36 17 1 2 2 1 8.04687744759E+05 8.04687819419E+05 8.04687670284E+05 grad-res 0 5 5 530 0 2 2 1 7.45673438915E+00 7.45673393249E+00 6.12419525092E-08 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 7.45673438915078E+00 (PID.TID 0000.0001) ADM finite-diff_grad = 7.45673393248580E+00 (PID.TID 0000.0001) ====== End of gradient-check number 5 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 6 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 531 29309 6 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 530 6 ph-grd -->hit<-- 37 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 37 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -1.22124532708767E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -2.10942374678780E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -1.68753899743024E-14 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089995D+06 --> objf_test(bi,bj) = 0.205004038769940D+06 (PID.TID 0000.0001) local fc = 0.804687808327177D+06 (PID.TID 0000.0001) global fc = 0.804687808327177D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687808327177E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.65973959207622E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -1.88737914186277E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -3.33066907387547E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = 8.88178419700125E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089991D+06 --> objf_test(bi,bj) = 0.205003911816880D+06 (PID.TID 0000.0001) local fc = 0.804687681374114D+06 (PID.TID 0000.0001) global fc = 0.804687681374114D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687681374114E+05 grad-res ------------------------------- grad-res 0 6 37 17 1 2 2 1 8.04687744759E+05 8.04687808327E+05 8.04687681374E+05 grad-res 0 6 6 531 0 2 2 1 6.34765354129E+00 6.34765317664E+00 5.74465573999E-08 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 6.34765354128635E+00 (PID.TID 0000.0001) ADM finite-diff_grad = 6.34765317663550E+00 (PID.TID 0000.0001) ====== End of gradient-check number 6 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 7 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 532 29309 7 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 531 7 ph-grd -->hit<-- 38 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 38 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -3.60822483003176E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = 1.33226762955019E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = 2.66453525910038E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089996D+06 --> objf_test(bi,bj) = 0.205004039670558D+06 (PID.TID 0000.0001) local fc = 0.804687809227796D+06 (PID.TID 0000.0001) global fc = 0.804687809227796D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687809227796E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = 6.66133814775094E-16 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -3.10862446895044E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -1.06581410364015E-14 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089990D+06 --> objf_test(bi,bj) = 0.205003910916544D+06 (PID.TID 0000.0001) local fc = 0.804687680473777D+06 (PID.TID 0000.0001) global fc = 0.804687680473777D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687680473777E+05 grad-res ------------------------------- grad-res 0 7 38 17 1 2 2 1 8.04687744759E+05 8.04687809228E+05 8.04687680474E+05 grad-res 0 7 7 532 0 2 2 1 6.43770134172E+00 6.43770095776E+00 5.96433897826E-08 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 6.43770134172468E+00 (PID.TID 0000.0001) ADM finite-diff_grad = 6.43770095775835E+00 (PID.TID 0000.0001) ====== End of gradient-check number 7 (ierr= 0) ======= (PID.TID 0000.0001) ====== Starts gradient-check number 8 (=ichknum) ======= ph-test icomp, ncvarcomp, ichknum 533 29309 8 ph-grd _loc: bi, bj, icomptest, ichknum 2 2 532 8 ph-grd -->hit<-- 39 17 1 1 (PID.TID 0000.0001) grdchk pos: i,j,k= 39 17 1 ; bi,bj= 2 2 ; iobc= 1 ; rec= 1 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -2.55351295663786E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -6.99440505513849E-15 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -5.77315972805081E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804090017D+06 --> objf_test(bi,bj) = 0.205004042215454D+06 (PID.TID 0000.0001) local fc = 0.804687811772713D+06 (PID.TID 0000.0001) global fc = 0.804687811772713D+06 (PID.TID 0000.0001) grdchk perturb(+)fc: fcpertplus = 8.04687811772713E+05 (PID.TID 0000.0001) Start initial hydrostatic pressure computation (PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Model current state (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) EXTERNAL_FIELDS_LOAD, it= 0 : Reading new data, i0,i1= 12 1 (prev= 12 0 ) (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F cg2d: Sum(rhs),rhsMax = -8.88178419700125E-15 2.53674886388737E+00 cg2d: Sum(rhs),rhsMax = -2.27595720048157E-15 5.04325900480939E+00 cg2d: Sum(rhs),rhsMax = -1.01030295240889E-14 6.29008403763954E+00 cg2d: Sum(rhs),rhsMax = -4.88498130835069E-15 6.70003306073569E+00 (PID.TID 0000.0001) Did not write pickup because writePickupAtEnd = FALSE (PID.TID 0000.0001) early fc = 0.000000000000000D+00 --> objf_test(bi,bj) = 0.221431048007422D+06 --> objf_test(bi,bj) = 0.241096917459820D+06 --> objf_test(bi,bj) = 0.137155804089968D+06 --> objf_test(bi,bj) = 0.205003908381863D+06 (PID.TID 0000.0001) local fc = 0.804687677939074D+06 (PID.TID 0000.0001) global fc = 0.804687677939074D+06 (PID.TID 0000.0001) grdchk perturb(-)fc: fcpertminus = 8.04687677939074E+05 grad-res ------------------------------- grad-res 0 8 39 17 1 2 2 1 8.04687744759E+05 8.04687811773E+05 8.04687677939E+05 grad-res 0 8 8 533 0 2 2 1 6.69168724574E+00 6.69168197201E+00 7.88101599980E-07 (PID.TID 0000.0001) ADM ref_cost_function = 8.04687744759400E+05 (PID.TID 0000.0001) ADM adjoint_gradient = 6.69168724573576E+00 (PID.TID 0000.0001) ADM finite-diff_grad = 6.69168197200634E+00 (PID.TID 0000.0001) ====== End of gradient-check number 8 (ierr= 0) ======= (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Gradient check results >>> START <<< (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) EPS = 1.000000E-02 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output h.p: Id Itile Jtile LAYER bi bj X(Id) X(Id)+/-EPS (PID.TID 0000.0001) grdchk output h.c: Id FC FC1 FC2 (PID.TID 0000.0001) grdchk output h.g: Id FC1-FC2/(2*EPS) ADJ GRAD(FC) 1-FDGRD/ADGRD (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 1 26 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 1 8.0468774475940E+05 8.0468776675048E+05 8.0468772295205E+05 (PID.TID 0000.0001) grdchk output (g): 1 2.1899214363657E+00 2.1899215003448E+00 2.9215229013779E-08 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 2 30 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 2 8.0468774475940E+05 8.0468775424978E+05 8.0468773545147E+05 (PID.TID 0000.0001) grdchk output (g): 2 9.3991559697315E-01 9.3991572287013E-01 1.3394496323205E-07 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 3 31 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 3 8.0468774475940E+05 8.0468775022033E+05 8.0468773947854E+05 (PID.TID 0000.0001) grdchk output (g): 3 5.3708943305537E-01 5.3708939704541E-01 -6.7046487917111E-08 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 4 32 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 4 8.0468774475940E+05 8.0468776687862E+05 8.0468772282786E+05 (PID.TID 0000.0001) grdchk output (g): 4 2.2025377198588E+00 2.2025390000849E+00 5.8125015300892E-07 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 5 36 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 5 8.0468774475940E+05 8.0468781941859E+05 8.0468767028391E+05 (PID.TID 0000.0001) grdchk output (g): 5 7.4567339324858E+00 7.4567343891508E+00 6.1241952509228E-08 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 6 37 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 6 8.0468774475940E+05 8.0468780832718E+05 8.0468768137411E+05 (PID.TID 0000.0001) grdchk output (g): 6 6.3476531766355E+00 6.3476535412863E+00 5.7446557399921E-08 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 7 38 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 7 8.0468774475940E+05 8.0468780922780E+05 8.0468768047378E+05 (PID.TID 0000.0001) grdchk output (g): 7 6.4377009577584E+00 6.4377013417247E+00 5.9643389782593E-08 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk output (p): 8 39 17 1 2 2 0.000000000E+00 -1.000000000E-02 (PID.TID 0000.0001) grdchk output (c): 8 8.0468774475940E+05 8.0468781177271E+05 8.0468767793907E+05 (PID.TID 0000.0001) grdchk output (g): 8 6.6916819720063E+00 6.6916872457358E+00 7.8810159998000E-07 (PID.TID 0000.0001) (PID.TID 0000.0001) grdchk summary : RMS of 8 ratios = 3.5228786244996E-07 (PID.TID 0000.0001) (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) // Gradient check results >>> END <<< (PID.TID 0000.0001) // ======================================================= (PID.TID 0000.0001) (PID.TID 0000.0001) Seconds in section "ALL [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 40.898254823172465 (PID.TID 0000.0001) System time: 0.72334497747942805 (PID.TID 0000.0001) Wall clock time: 42.798789978027344 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) Seconds in section "INITIALISE_FIXED [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 0.13449899875558913 (PID.TID 0000.0001) System time: 3.5695001017302275E-002 (PID.TID 0000.0001) Wall clock time: 0.60699796676635742 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) Seconds in section "THE_MAIN_LOOP (F) [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 2.9874020069837570 (PID.TID 0000.0001) System time: 0.58907900378108025 (PID.TID 0000.0001) Wall clock time: 3.9762380123138428 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) Seconds in section "INITIALISE_VARIA [THE_MAIN_LOOP]": (PID.TID 0000.0001) User time: 0.93090231716632843 (PID.TID 0000.0001) System time: 5.4338030517101288E-002 (PID.TID 0000.0001) Wall clock time: 1.1042797565460205 (PID.TID 0000.0001) No. starts: 17 (PID.TID 0000.0001) No. stops: 17 (PID.TID 0000.0001) Seconds in section "MAIN LOOP [THE_MAIN_LOOP]": (PID.TID 0000.0001) User time: 33.635853126645088 (PID.TID 0000.0001) System time: 0.56941493600606918 (PID.TID 0000.0001) Wall clock time: 34.546782255172729 (PID.TID 0000.0001) No. starts: 17 (PID.TID 0000.0001) No. stops: 17 (PID.TID 0000.0001) Seconds in section "MAIN_DO_LOOP [THE_MAIN_LOOP]": (PID.TID 0000.0001) User time: 33.633706837892532 (PID.TID 0000.0001) System time: 0.56936888396739960 (PID.TID 0000.0001) Wall clock time: 34.544605016708374 (PID.TID 0000.0001) No. starts: 17 (PID.TID 0000.0001) No. stops: 17 (PID.TID 0000.0001) Seconds in section "LOAD_FIELDS_DRIVER [FORWARD_STEP]": (PID.TID 0000.0001) User time: 6.9913178682327271E-002 (PID.TID 0000.0001) System time: 5.1540508866310120E-003 (PID.TID 0000.0001) Wall clock time: 8.8630676269531250E-002 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "EXTERNAL_FLDS_LOAD [LOAD_FLDS_DRIVER]": (PID.TID 0000.0001) User time: 6.8762034177780151E-002 (PID.TID 0000.0001) System time: 5.1210001111030579E-003 (PID.TID 0000.0001) Wall clock time: 8.7509632110595703E-002 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "CTRL_MAP_FORCING [FORWARD_STEP]": (PID.TID 0000.0001) User time: 5.1298737525939941E-004 (PID.TID 0000.0001) System time: 1.3023614883422852E-005 (PID.TID 0000.0001) Wall clock time: 5.2905082702636719E-004 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "DO_ATMOSPHERIC_PHYS [FORWARD_STEP]": (PID.TID 0000.0001) User time: 4.1427105665206909E-002 (PID.TID 0000.0001) System time: 1.3829916715621948E-003 (PID.TID 0000.0001) Wall clock time: 4.2854547500610352E-002 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "DO_OCEANIC_PHYS [FORWARD_STEP]": (PID.TID 0000.0001) User time: 6.6596977710723877 (PID.TID 0000.0001) System time: 6.6282019019126892E-002 (PID.TID 0000.0001) Wall clock time: 6.7264490127563477 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "THERMODYNAMICS [FORWARD_STEP]": (PID.TID 0000.0001) User time: 10.776220470666885 (PID.TID 0000.0001) System time: 0.14712599664926529 (PID.TID 0000.0001) Wall clock time: 10.925431013107300 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "DYNAMICS [FORWARD_STEP]": (PID.TID 0000.0001) User time: 13.725170254707336 (PID.TID 0000.0001) System time: 0.24813997000455856 (PID.TID 0000.0001) Wall clock time: 13.975260734558105 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "SOLVE_FOR_PRESSURE [FORWARD_STEP]": (PID.TID 0000.0001) User time: 2.4812054038047791 (PID.TID 0000.0001) System time: 9.9399685859680176E-004 (PID.TID 0000.0001) Wall clock time: 2.4823095798492432 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "MOM_CORRECTION_STEP [FORWARD_STEP]": (PID.TID 0000.0001) User time: 0.39686328172683716 (PID.TID 0000.0001) System time: 1.0880455374717712E-003 (PID.TID 0000.0001) Wall clock time: 0.39801645278930664 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "INTEGR_CONTINUITY [FORWARD_STEP]": (PID.TID 0000.0001) User time: 0.42757356166839600 (PID.TID 0000.0001) System time: 1.9810199737548828E-003 (PID.TID 0000.0001) Wall clock time: 0.42963600158691406 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "TRC_CORRECTION_STEP [FORWARD_STEP]": (PID.TID 0000.0001) User time: 5.7560205459594727E-004 (PID.TID 0000.0001) System time: 9.0003013610839844E-006 (PID.TID 0000.0001) Wall clock time: 5.8245658874511719E-004 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "BLOCKING_EXCHANGES [FORWARD_STEP]": (PID.TID 0000.0001) User time: 0.48233079910278320 (PID.TID 0000.0001) System time: 9.6899271011352539E-004 (PID.TID 0000.0001) Wall clock time: 0.48342084884643555 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "COST_TILE [FORWARD_STEP]": (PID.TID 0000.0001) User time: 0.77092266082763672 (PID.TID 0000.0001) System time: 0.0000000000000000 (PID.TID 0000.0001) Wall clock time: 0.77105093002319336 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "DO_THE_MODEL_IO [FORWARD_STEP]": (PID.TID 0000.0001) User time: 3.2958924770355225E-002 (PID.TID 0000.0001) System time: 5.9550404548645020E-003 (PID.TID 0000.0001) Wall clock time: 0.17225909233093262 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "DO_WRITE_PICKUP [FORWARD_STEP]": (PID.TID 0000.0001) User time: 4.3878376483917236E-002 (PID.TID 0000.0001) System time: 1.2197017669677734E-002 (PID.TID 0000.0001) Wall clock time: 0.10256338119506836 (PID.TID 0000.0001) No. starts: 71 (PID.TID 0000.0001) No. stops: 71 (PID.TID 0000.0001) Seconds in section "COST_FINAL [ADJOINT SPIN-DOWN]": (PID.TID 0000.0001) User time: 1.7390251159667969E-003 (PID.TID 0000.0001) System time: 3.6001205444335938E-005 (PID.TID 0000.0001) Wall clock time: 1.7826557159423828E-003 (PID.TID 0000.0001) No. starts: 17 (PID.TID 0000.0001) No. stops: 17 (PID.TID 0000.0001) Seconds in section "THE_MAIN_LOOP (A) [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 5.9712274074554443 (PID.TID 0000.0001) System time: 2.3562014102935791E-002 (PID.TID 0000.0001) Wall clock time: 6.0112929344177246 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 2.5646209716796875E-002 (PID.TID 0000.0001) System time: 4.7619938850402832E-003 (PID.TID 0000.0001) Wall clock time: 3.6576032638549805E-002 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) Seconds in section "CTRL_PACK [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 2.7779579162597656E-002 (PID.TID 0000.0001) System time: 2.1290183067321777E-003 (PID.TID 0000.0001) Wall clock time: 3.7117958068847656E-002 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) Seconds in section "GRDCHK_MAIN [THE_MODEL_MAIN]": (PID.TID 0000.0001) User time: 31.751627922058105 (PID.TID 0000.0001) System time: 6.8093955516815186E-002 (PID.TID 0000.0001) Wall clock time: 32.130488872528076 (PID.TID 0000.0001) No. starts: 1 (PID.TID 0000.0001) No. stops: 1 (PID.TID 0000.0001) // ====================================================== (PID.TID 0000.0001) // Tile <-> Tile communication statistics (PID.TID 0000.0001) // ====================================================== (PID.TID 0000.0001) // o Tile number: 000001 (PID.TID 0000.0001) // No. X exchanges = 0 (PID.TID 0000.0001) // Max. X spins = 0 (PID.TID 0000.0001) // Min. X spins = 1000000000 (PID.TID 0000.0001) // Total. X spins = 0 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00 (PID.TID 0000.0001) // No. Y exchanges = 0 (PID.TID 0000.0001) // Max. Y spins = 0 (PID.TID 0000.0001) // Min. Y spins = 1000000000 (PID.TID 0000.0001) // Total. Y spins = 0 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00 (PID.TID 0000.0001) // o Tile number: 000002 (PID.TID 0000.0001) // No. X exchanges = 0 (PID.TID 0000.0001) // Max. X spins = 0 (PID.TID 0000.0001) // Min. X spins = 1000000000 (PID.TID 0000.0001) // Total. X spins = 0 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00 (PID.TID 0000.0001) // No. Y exchanges = 0 (PID.TID 0000.0001) // Max. Y spins = 0 (PID.TID 0000.0001) // Min. Y spins = 1000000000 (PID.TID 0000.0001) // Total. Y spins = 0 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00 (PID.TID 0000.0001) // o Tile number: 000003 (PID.TID 0000.0001) // No. X exchanges = 0 (PID.TID 0000.0001) // Max. X spins = 0 (PID.TID 0000.0001) // Min. X spins = 1000000000 (PID.TID 0000.0001) // Total. X spins = 0 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00 (PID.TID 0000.0001) // No. Y exchanges = 0 (PID.TID 0000.0001) // Max. Y spins = 0 (PID.TID 0000.0001) // Min. Y spins = 1000000000 (PID.TID 0000.0001) // Total. Y spins = 0 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00 (PID.TID 0000.0001) // o Tile number: 000004 (PID.TID 0000.0001) // No. X exchanges = 0 (PID.TID 0000.0001) // Max. X spins = 0 (PID.TID 0000.0001) // Min. X spins = 1000000000 (PID.TID 0000.0001) // Total. X spins = 0 (PID.TID 0000.0001) // Avg. X spins = 0.00E+00 (PID.TID 0000.0001) // No. Y exchanges = 0 (PID.TID 0000.0001) // Max. Y spins = 0 (PID.TID 0000.0001) // Min. Y spins = 1000000000 (PID.TID 0000.0001) // Total. Y spins = 0 (PID.TID 0000.0001) // Avg. Y spins = 0.00E+00 (PID.TID 0000.0001) // o Thread number: 000001 (PID.TID 0000.0001) // No. barriers = 17686 (PID.TID 0000.0001) // Max. barrier spins = 1 (PID.TID 0000.0001) // Min. barrier spins = 1 (PID.TID 0000.0001) // Total barrier spins = 17686 (PID.TID 0000.0001) // Avg. barrier spins = 1.00E+00 PROGRAM MAIN: Execution ended Normally