Need help regarding free surface in South China Sea

[often]

I am setting up a test case of South China Sea. The model currently run for 3 year forced by 2004 NCEP reanalysis 1 and SODA_2.2.4, with 10 sponge/nudging layer near boundary. Below are the std output and result in the 4th year (ssh and temperature at -230m). In January, free surface was too deep in the middle of South China Sea and too high to the west of Taiwan when comparing to SODA (see figure_WPAC_zeta.png). Temperature at -230m shows similar pattern (see figure_WPAC_temp(-230).png). I also check on this link http://www.myroms.org/Workshops/ROMS200 ... petkin.pdf. It seems that this problem also existed in that Pacific Ocean case. So my question is: Is there something I can try to further calibrate the model (free surface or thermoline), targeting South China Sea?

Code: Select all

 Model Input Parameters:  ROMS/TOMS version 3.6  
                          Saturday - February 16, 2013 -  1:54:19 AM
 -----------------------------------------------------------------------------

 South China Sea Test Case

 Operating system : Linux
 CPU/hardware     : i686
 Compiler system  : gfortran
 Compiler command : /usr/bin/mpif90
 Compiler flags   : -frepack-arrays -O3 -ffast-math -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none

 Input Script  : ./External/ocean_SCS.in

 SVN Root URL  : https://www.myroms.org/svn/src/trunk
 SVN Revision  : 631M

 Local Root    : /home/c606/Downloads/MyDir/trunk
 Header Dir    : /home/c606/Downloads/MyDir/MyProject/ljp_test01/Include
 Header file   : scs.h
 Analytical Dir: /home/c606/Downloads/MyDir/MyProject/ljp_test01/Functionals

 Resolution, Grid 01: 0166x0166x030,  Parallel Nodes:   4,  Tiling: 002x002


 Physical Parameters, Grid: 01
 =============================

     526320  ntimes          Number of timesteps for 3-D equations.
    600.000  dt              Timestep size (s) for 3-D equations.
         30  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
          1  ERstr           Starting ensemble/perturbation run number.
          1  ERend           Ending ensemble/perturbation run number.
          0  nrrec           Number of restart records to read from disk.
          T  LcycleRST       Switch to recycle time-records in restart file.
       4386  nRST            Number of timesteps between the writing of data
                               into restart fields.
       4386  ninfo           Number of timesteps between print of information
                               to standard output.
          T  ldefout         Switch to create a new output NetCDF file(s).
       4386  nHIS            Number of timesteps between the writing fields
                               into history file.
          0  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
       4386  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
 1.0000E+01  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 1.0000E+01  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 1.0000E+01  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 1.0000E-06  Akt_bak(01)     Background vertical mixing coefficient (m2/s)
                               for tracer 01: temp
 1.0000E-06  Akt_bak(02)     Background vertical mixing coefficient (m2/s)
                               for tracer 02: salt
 1.0000E-05  Akv_bak         Background vertical mixing coefficient (m2/s)
                               for momentum.
 3.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2           Quadratic bottom drag coefficient.
 2.0000E-02  Zob             Bottom roughness (m).
 4.2000E+01  blk_ZQ          Height (m) of surface air humidity measurement.
 2.0000E+00  blk_ZT          Height (m) of surface air temperature measurement.
 1.0000E+01  blk_ZW          Height (m) of surface winds measurement.
 1.0000E-01  Dcrit           Minimum depth for wetting and drying (m).
          1  lmd_Jwt         Jerlov water type.
          2  Vtransform      S-coordinate transformation equation.
          4  Vstretching     S-coordinate stretching function.
 5.0000E+00  theta_s         S-coordinate surface control parameter.
 4.0000E-01  theta_b         S-coordinate bottom  control parameter.
     10.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
      0.000  dstart          Time-stamp assigned to model initialization (days).
       0.00  time_ref        Reference time for units attribute (yyyymmdd.dd)
 3.6000E+02  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 3.6000E+02  Tnudg(02)       Nudging/relaxation time scale (days)
                               for tracer 02: salt
 3.6000E+02  Znudg           Nudging/relaxation time scale (days)
                               for free-surface.
 3.6000E+02  M2nudg          Nudging/relaxation time scale (days)
                               for 2D momentum.
 3.6000E+02  M3nudg          Nudging/relaxation time scale (days)
                               for 3D momentum.
 3.6000E+02  obcfac          Factor between passive and active
                               open boundary conditions.
          F  VolCons(1)      NLM western  edge boundary volume conservation.
          F  VolCons(2)      NLM southern edge boundary volume conservation.
          F  VolCons(3)      NLM eastern  edge boundary volume conservation.
          F  VolCons(4)      NLM northern edge boundary volume conservation.
      4.000  T0              Background potential temperature (C) constant.
     35.000  S0              Background salinity (PSU) constant.
      1.000  gamma2          Slipperiness variable: free-slip (1.0) or 
                                                    no-slip (-1.0).
          T  Hout(idFsur)    Write out free-surface.
          T  Hout(idu3dE)    Write out 3D U-wastward  component at RHO-points.
          T  Hout(idv3dN)    Write out 3D V-northward component at RHO-points.
          T  Hout(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt

          T  Aout(idFsur)    Write out averaged free-surface.
          T  Aout(idu3dE)    Write out averaged 3D U-eastward  at RHO-points.
          T  Aout(idv3dN)    Write out averaged 3D V-northward at RHO-points.
          T  Aout(idTvar)    Write out averaged tracer 01: temp
          T  Aout(idTvar)    Write out averaged tracer 02: salt

 Output/Input Files:

             Output Restart File:  ./Output/WPAC_rst.nc
             Output History File:  ./Output/WPAC_his.nc
            Output Averages File:  ./Output/WPAC_avg.nc
                 Input Grid File:  ./Data/WPAC_grd.nc
    Input Nonlinear Initial File:  ./Data/WPAC_ini.nc
           Input Forcing File 01:  ./Data/WPAC_wind_ncep_4xday.nc
           Input Forcing File 02:  ./Data/WPAC_lwrad_ncep_4xday.nc
           Input Forcing File 03:  ./Data/WPAC_swrad_ncep_1xday.nc
           Input Forcing File 04:  ./Data/WPAC_rain_ncep_4xday.nc
           Input Forcing File 05:  ./Data/WPAC_Tair_ncep_4xday.nc
           Input Forcing File 06:  ./Data/WPAC_Qair_ncep_4xday.nc
           Input Forcing File 07:  ./Data/WPAC_Pair_ncep_4xday.nc
          Input Climatology File:  ./Data/WPAC_clm.nc
             Input Boundary File:  ./Data/WPAC_obc.nc
          Input/Output USER File:  ROMS/External/MyFile.dat

 Generic User Parameters:

 1.0000E+01  user(01)        User parameter 01.

 Tile partition information for Grid 01:  0166x0166x0030  tiling: 002x002

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       83        1       83   206670
        1       84      166        1       83   206670
        2        1       83       84      166   206670
        3       84      166       84      166   206670

 Tile minimum and maximum fractional grid coordinates:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    83.50     0.50    83.50  RHO-points
        1    83.50   166.50     0.50    83.50  RHO-points
        2     0.50    83.50    83.50   166.50  RHO-points
        3    83.50   166.50    83.50   166.50  RHO-points

        0     1.00    83.50     0.50    83.50    U-points
        1    83.50   166.00     0.50    83.50    U-points
        2     1.00    83.50    83.50   166.50    U-points
        3    83.50   166.00    83.50   166.50    U-points

        0     0.50    83.50     1.00    83.50    V-points
        1    83.50   166.50     1.00    83.50    V-points
        2     0.50    83.50    83.50   166.00    V-points
        3    83.50   166.50    83.50   166.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     190
               HaloSizeJ(1) =     190
                TileSide(1) =      89
                TileSize(1) =    7921


 Lateral Boundary Conditions: NLM
 ============================

 Variable               Grid    West Edge   South Edge  East Edge   North Edge
 ---------              ----    ----------  ----------  ----------  ----------

 zeta                     1     Chapman     Chapman     Chapman     Chapman

 ubar                     1     Flather     Flather     Flather     Flather

 vbar                     1     Flather     Flather     Flather     Flather

 u                        1     Rad + Nud   Rad + Nud   Rad + Nud   Rad + Nud

 v                        1     Rad + Nud   Rad + Nud   Rad + Nud   Rad + Nud

 temp                     1     Rad + Nud   Rad + Nud   Rad + Nud   Rad + Nud

 salt                     1     Rad + Nud   Rad + Nud   Rad + Nud   Rad + Nud

 Activated C-preprocessing Options:

 SCS                 South China Sea Test Case
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 AVERAGES            Writing out time-averaged nonlinear model fields.
 BULK_FLUXES         Surface bulk fluxes parameterization.
 CURVGRID            Orthogonal curvilinear grid.
 DIURNAL_SRFLUX      Modulate shortwave radiation by the local diurnal cycle.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 EMINUSP             Compute Salt Flux using E-P.
 LMD_CONVEC          LMD convective mixing due to shear instability.
 LMD_MIXING          Large/McWilliams/Doney interior mixing.
 LMD_NONLOCAL        LMD convective nonlocal transport.
 LMD_RIMIX           LMD diffusivity due to shear instability.
 LMD_SKPP            KPP surface boundary layer mixing.
 LONGWAVE_OUT        Compute outgoing longwave radiation internally.
 MASKING             Land/Sea masking.
 MIX_GEO_TS          Mixing of tracers along geopotential surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 SPLINES             Conservative parabolic spline reconstruction.
 SPHERICAL           Spherical grid configuration.
 SPONGE              Enhanced horizontal mixing in the sponge areas.
 TCLIMATOLOGY        Processing tracer climatology data.
 TCLM_NUDGING        Nudging toward tracer climatology.
 TS_U3HADVECTION     Third-order upstream horizontal advection of tracers.
 TS_C4VADVECTION     Fourth-order centered vertical advection of tracers.
 TS_DIF2             Harmonic mixing of tracers.
 UV_ADV              Advection of momentum.
 UV_COR              Coriolis term.
 UV_U3HADVECTION     Third-order upstream horizontal advection of 3D momentum.
 UV_C4VADVECTION     Fourth-order centered vertical advection of momentum.
 UV_QDRAG            Quadratic bottom stress.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 WET_DRY             Wetting and drying activated.

 Process Information:

 Node #  0 (pid=   29719) is active.
 Node #  1 (pid=   29720) is active.
 Node #  2 (pid=   29721) is active.
 Node #  3 (pid=   29722) is active.

 INITIAL: Configuring and initializing forward nonlinear model ...


 Vertical S-coordinate System: 

 level   S-coord     Cs-curve   Z   at hmin       at hc    half way     at hmax

    30   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    29  -0.0333333  -0.0002307       -0.168      -0.168      -1.069      -1.806
    28  -0.0666667  -0.0009291       -0.338      -0.338      -3.631      -6.597
    27  -0.1000000  -0.0021144       -0.511      -0.511      -7.746     -14.497
    26  -0.1333333  -0.0038189       -0.686      -0.686     -13.520     -25.712
    25  -0.1666667  -0.0060892       -0.864      -0.864     -21.099     -40.539
    24  -0.2000000  -0.0089871       -1.045      -1.045     -30.682     -59.373
    23  -0.2333333  -0.0125913       -1.230      -1.230     -42.520     -82.716
    22  -0.2666667  -0.0169996       -1.418      -1.418     -56.924    -111.192
    21  -0.3000000  -0.0223307       -1.612      -1.612     -74.274    -145.561
    20  -0.3333333  -0.0287278       -1.810      -1.810     -95.027    -186.734
    19  -0.3666667  -0.0363614       -2.015      -2.015    -119.727    -235.802
    18  -0.4000000  -0.0454338       -2.227      -2.227    -149.020    -294.056
    17  -0.4333333  -0.0561830       -2.448      -2.448    -183.665    -363.014
    16  -0.4666667  -0.0688882       -2.678      -2.678    -224.554    -444.461
    15  -0.5000000  -0.0838752       -2.919      -2.919    -272.728    -540.475
    14  -0.5333333  -0.1015226       -3.174      -3.174    -329.393    -653.474
    13  -0.5666667  -0.1222681       -3.445      -3.445    -395.949    -786.252
    12  -0.6000000  -0.1466152       -3.733      -3.733    -474.001    -942.023
    11  -0.6333333  -0.1751391       -4.042      -4.042    -565.386   -1124.460
    10  -0.6666667  -0.2084930       -4.376      -4.376    -672.190   -1337.734
     9  -0.7000000  -0.2474113       -4.737      -4.737    -796.755   -1586.532
     8  -0.7333333  -0.2927113       -5.130      -5.130    -941.692   -1876.072
     7  -0.7666667  -0.3452901       -5.560      -5.560   -1109.865   -2212.083
     6  -0.8000000  -0.4061157       -6.031      -6.031   -1304.362   -2600.743
     5  -0.8333333  -0.4762095       -6.548      -6.548   -1528.445   -3048.574
     4  -0.8666667  -0.5566161       -7.116      -7.116   -1785.449   -3562.245
     3  -0.9000000  -0.6483588       -7.742      -7.742   -2078.639   -4148.289
     2  -0.9333333  -0.7523753       -8.429      -8.429   -2411.009   -4812.692
     1  -0.9666667  -0.8694298       -9.180      -9.180   -2784.998   -5560.333
     0  -1.0000000  -1.0000000      -10.000     -10.000   -3202.131   -6394.263

 Time Splitting Weights: ndtfast =  30    nfast =  42

    Primary            Secondary            Accumulated to Current Step

  1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
  2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
  3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
  4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
  5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
  6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
  7-0.0011849025289723 0.0336508794757796-0.0107112868023601 0.2343461315221331
  8-0.0005032751608632 0.0336903762267453-0.0112145619632232 0.2680365077488784
  9 0.0003887272597151 0.0337071520654408-0.0108258347035082 0.3017436598143192
 10 0.0014892209965583 0.0336941944901169-0.0093366137069499 0.3354378543044362
 11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
 12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325831
 13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801077
 14 0.0079087469427945 0.0332075372104522 0.0116826306292292 0.4692492050905598
 15 0.0099910761708919 0.0329439123123590 0.0216737068001212 0.5021931174029188
 16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
 17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
 18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
 19 0.0199444086685725 0.0311389577709445 0.0857756755402366 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399430 0.6603333185976138
 21 0.0256737146312910 0.0297153720153352 0.1342125541712340 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524357 0.7189082721405746
 23 0.0316613792205219 0.0279045862015855 0.1945237931729576 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236651 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893979 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886211 0.8237939964192927
 27 0.0432936737565710 0.0230872677537126 0.3506756411451921 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772277 0.8685254094681656
 29 0.0481720587108284 0.0201154903974257 0.4447073467880561 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442381 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256412 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325518 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348095 0.9524525304259084
 34 0.0531891349131379 0.0115726061288397 0.7060109510479474 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213234 0.9738247715198165
 36 0.0503158038019030 0.0080591141492892 0.8085423793232265 0.9818838856691057
 37 0.0473537721847154 0.0063819206892258 0.8558961515079418 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497606 0.9930692679747335
 39 0.0376397765791564 0.0033640675316746 0.9367177506289170 0.9964333355064081
 40 0.0305543017255207 0.0021094083123694 0.9672720523544377 0.9985427438187775
 41 0.0217382098544505 0.0010909315881854 0.9890102622088881 0.9996336754069629
 42 0.0109897377911118 0.0003663245930371 0.9999999999999999 1.0000000000000000

 ndtfast, nfast =   30  42   nfast/ndtfast =  1.40000

 Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

    1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

 Minimum X-grid spacing, DXmin =  2.71525249E+01 km
 Maximum X-grid spacing, DXmax =  3.30952409E+01 km
 Minimum Y-grid spacing, DYmin =  2.85531544E+01 km
 Maximum Y-grid spacing, DYmax =  3.48024618E+01 km
 Minimum Z-grid spacing, DZmin =  1.67820182E-01 m
 Maximum Z-grid spacing, DZmax =  8.33929317E+02 m

 Minimum barotropic Courant Number =  8.68661890E-03
 Maximum barotropic Courant Number =  2.52566146E-01
 Maximum Coriolis   Courant Number =  5.00304087E-02


 NLM: GET_STATE - Read state initial conditions,             t =     0 00:00:00
                   (File: WPAC_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  1.22439120E+00)
                - vertically integrated u-momentum component
                   (Min = -2.23804461E-01 Max =  2.50266579E-01)
                - vertically integrated v-momentum component
                   (Min = -3.23598644E-01 Max =  2.07418007E-01)
                - u-momentum component
                   (Min = -8.71930988E-01 Max =  1.17035636E+00)
                - v-momentum component
                   (Min = -1.07161324E+00 Max =  8.26374012E-01)
                - potential temperature
                   (Min =  0.00000000E+00 Max =  3.11600403E+01)
                - salinity
                   (Min =  0.00000000E+00 Max =  3.56537157E+01)
    GET_2DFLD   - surface u-wind component,                  t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_wind_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -8.17192311E+00 Max =  8.98236117E+00)
    GET_2DFLD   - surface v-wind component,                  t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_wind_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -1.21344838E+01 Max =  6.13346964E+00)
    GET_2DFLD   - surface air pressure,                      t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_Pair_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  1.00483394E+03 Max =  1.04817782E+03)
    GET_2DFLD   - solar shortwave radiation flux,            t =   362 24:00:00
                   (Rec=0362, Index=1, File: WPAC_swrad_ncep_1xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  1.34111520E-05 Max =  7.40225670E-05)
    GET_2DFLD   - downwelling longwave radiation flux,       t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_lwrad_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  1.91454537E-05 Max =  1.10117997E-04)
    GET_2DFLD   - surface air temperature,                   t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_Tair_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -2.20085964E+01 Max =  3.45231036E+01)
    GET_2DFLD   - surface air relative humidity,             t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_Qair_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  2.96213914E-01 Max =  1.00000000E+00)
    GET_2DFLD   - rain fall rate,                            t =   362 24:00:00
                   (Rec=1445, Index=1, File: WPAC_rain_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -2.32830644E-10 Max =  5.92064589E-04)
    GET_NGFLD   - free-surface western boundary condition,   t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  6.25749824E-01)
    GET_NGFLD   - free-surface eastern boundary condition,   t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  1.32274256E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  6.16081601E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  1.36994433E+00)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.37479047E-02 Max =  2.93913857E-02)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.30086418E-02 Max =  2.61413055E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.48327145E-01 Max =  1.11527912E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.54204850E-01 Max =  1.18102554E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -7.93827916E-03 Max =  2.25392292E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -7.28950773E-03 Max =  2.50784507E-02)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.34305693E-02 Max =  5.42625654E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.36480849E-02 Max =  5.52313674E-02)
    GET_NGFLD   - 3D u-momentum western boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -3.25224559E-01 Max =  2.58872651E-01)
    GET_NGFLD   - 3D v-momentum western boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -3.28389876E-01 Max =  2.59596502E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -5.51423810E-01 Max =  5.38648028E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -5.07877585E-01 Max =  5.17713060E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.08242782E-01 Max =  2.51820315E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.97628494E-01 Max =  2.58631065E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.62608113E-01 Max =  4.78339571E-01)
    GET_NGFLD   - 3D v-momentum northern boundary condition, t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.54034128E-01 Max =  4.32484092E-01)
    GET_NGFLD   - temperature western boundary condition,    t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.04730306E+01)
    GET_NGFLD   - salinity western boundary condition,       t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.54758940E+01)
    GET_NGFLD   - temperature eastern boundary condition,    t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.01480682E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.56205944E+01)
    GET_NGFLD   - temperature southern boundary condition,   t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.20708415E+01)
    GET_NGFLD   - salinity southern boundary condition,      t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.60100836E+01)
    GET_NGFLD   - temperature northern boundary condition,   t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  2.25741814E+01)
    GET_NGFLD   - salinity northern boundary condition,      t =   349 24:00:00
                   (Rec=0012, Index=2, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.49039353E+01)
    GET_3DFLD   - potential temperature climatology,         t =   349 24:00:00
                   (Rec=0012, Index=1, File: WPAC_clm.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.11600403E+01)
    GET_3DFLD   - salinity climatology,                      t =   349 24:00:00
                   (Rec=0012, Index=1, File: WPAC_clm.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.56537157E+01)

 Maximum grid stiffness ratios:  rx0 =   3.409140E-01 (Beckmann and Haidvogel)
                                 rx1 =   5.327743E+00 (Haney)


 Initial basin volumes: TotVolume =  5.5948638503E+16 m3
                        MinVolume =  1.3389253082E+08 m3
                        MaxVolume =  8.7151967084E+11 m3
                          Max/Min =  6.5090984948E+03


 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00526320)

    GET_2DFLD   - surface u-wind component,                  t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_wind_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -8.69136515E+00 Max =  1.22447927E+01)
    GET_2DFLD   - surface v-wind component,                  t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_wind_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -8.52942774E+00 Max =  8.06503930E+00)
    GET_2DFLD   - surface air pressure,                      t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_Pair_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  1.00428609E+03 Max =  1.04486133E+03)
    GET_2DFLD   - solar shortwave radiation flux,            t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_swrad_ncep_1xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  1.70221840E-05 Max =  7.10536291E-05)
    GET_2DFLD   - downwelling longwave radiation flux,       t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_lwrad_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  2.87331803E-05 Max =  1.09158491E-04)
    GET_2DFLD   - surface air temperature,                   t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_Tair_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min = -1.62258109E+01 Max =  3.20900977E+01)
    GET_2DFLD   - surface air relative humidity,             t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_Qair_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  2.55303849E-01 Max =  1.00000000E+00)
    GET_2DFLD   - rain fall rate,                            t =     1 24:00:00
                   (Rec=0001, Index=2, File: WPAC_rain_ncep_4xday.nc)
                   (Tmin=          2.0000 Tmax=        363.0000)
                   (Min =  9.98843461E-08 Max =  1.05174169E-03)
    GET_NGFLD   - free-surface western boundary condition,   t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  5.48380650E-01)
    GET_NGFLD   - free-surface eastern boundary condition,   t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  1.03682614E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  7.74583883E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  8.54006664E-01)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.30951235E-02 Max =  1.68044290E-02)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.32421957E-02 Max =  1.63933506E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.22875168E-01 Max =  1.20368084E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.20937758E-01 Max =  1.17347394E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.81692396E-02 Max =  1.36521768E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.22651651E-02 Max =  1.51419809E-02)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.85275264E-02 Max =  4.89064052E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -1.68151211E-02 Max =  4.35439429E-02)
    GET_NGFLD   - 3D u-momentum western boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.61556517E-01 Max =  2.23654487E-01)
    GET_NGFLD   - 3D v-momentum western boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.37013401E-01 Max =  2.18353432E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -7.59196824E-01 Max =  9.02618851E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -6.84316447E-01 Max =  8.76143222E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.90052285E-01 Max =  1.05781450E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -2.58984573E-01 Max =  1.12168103E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -5.39053258E-02 Max =  4.29254405E-01)
    GET_NGFLD   - 3D v-momentum northern boundary condition, t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min = -4.21430259E-02 Max =  4.38323941E-01)
    GET_NGFLD   - temperature western boundary condition,    t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  2.95494769E+01)
    GET_NGFLD   - salinity western boundary condition,       t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.54916284E+01)
    GET_NGFLD   - temperature eastern boundary condition,    t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.02461958E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.55536569E+01)
    GET_NGFLD   - temperature southern boundary condition,   t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.11600403E+01)
    GET_NGFLD   - salinity southern boundary condition,      t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.56537157E+01)
    GET_NGFLD   - temperature northern boundary condition,   t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  1.96588363E+01)
    GET_NGFLD   - salinity northern boundary condition,      t =    14 24:00:00
                   (Rec=0001, Index=1, File: WPAC_obc.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.47379970E+01)
    GET_3DFLD   - potential temperature climatology,         t =    14 24:00:00
                   (Rec=0001, Index=2, File: WPAC_clm.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.11600403E+01)
    GET_3DFLD   - salinity climatology,                      t =    14 24:00:00
                   (Rec=0001, Index=2, File: WPAC_clm.nc)
                   (Tmin=         15.0000 Tmax=        350.0000)
                   (Min =  0.00000000E+00 Max =  3.56537157E+01)

[mengqingjun]

I am wandering how you set up your grid file?
Using Seagrid,or other else?
Could you kindly tell me ?
Thank you!

[often]

it was seagrid, and i get those mex files compiled. But do not really know if the mex files really work.

[LiuZHQ]

Hi often,
How about you compare the solutions of zeta_east, zeta_south,zeta_north in your boundary file with SODA? it there big difference?
I think you included the tidal forcing in your case. Did the model you want to nest with included tide also?
Or, according to your description, you nested with SODA?

ZQ

[often]

Hi ZQ,
The solution did not forced with any tide. The zeta_xx in boundary file was generated by interpolating SODA monthly mean ssh into ROMS grid boundary, so zeta_xx should match with SODA ssh. But I have discard that solution for a while, so can not actully verify if zeta_xx==ssh.
I am not sure if the involvement of tide could make a different to my model's monthly mean zeta output. Thanks for the reminder to try that.
Beside that, what I can come up with is to increase the resolution of air-sea forcing(particularly wind), bathmetry and the grid, or to calibrate the wind measure height, if I can access to those higher resulution data and have some leisure time to play with ROMS(kind of busy.. right now).

Thanks for your comment!