Blowing up problem

Discussion on computers, ROMS installation and compiling

Moderators: arango, robertson

Post Reply
Message
Author
dhd_1979
Posts: 5
Joined: Sun Dec 26, 2010 9:38 pm
Location: Marine Blue Australia

Blowing up problem

#1 Unread post by dhd_1979 »

I have a succesful modelling of an area using a grid of 116x156x32, but when I changes the grid size of the same area to 0378x0358x032, It blows up just at the beginging, Could any one help me? the following is the output the model:

=====================================================================================================

Code: Select all

 Model Input Parameters:  ROMS/TOMS version 3.5  
                          Thursday - August 23, 2012 -  8:05:16 AM
 -----------------------------------------------------------------------------

 South_Sea

 Operating system : Darwin
 CPU/hardware     : i386
 Compiler system  : ifort
 Compiler command : /usr/local/mpich-intel/bin//mpif90
 Compiler flags   : -heap-arrays -fp-model precise -ip -O3 -axP -free

 Input Script  : para.ini

 SVN Root URL  : https:://myroms.org/svn/src
 SVN Revision  : 564M

 Local Root    : /Users/hd/roms/myroms34/trunk
 Header Dir    : /Users/hd/roms/myroms34/southsea
 Header file   : southsea.h
 Analytical Dir: /Users/hd/roms/myroms34/southsea

 Resolution, Grid 01: 0378x0358x032,  Parallel Nodes:   2,  Tiling: 001x002


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

      72000  ntimes          Number of timesteps for 3-D equations.
     90.000  dt              Timestep size (s) for 3-D equations.
          9  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.
          1  nrrec           Number of restart records to read from disk.
          T  LcycleRST       Switch to recycle time-records in restart file.
        240  nRST            Number of timesteps between the writing of data
                               into restart fields.
          1  ninfo           Number of timesteps between print of information
                               to standard output.
          T  ldefout         Switch to create a new output NetCDF file(s).
        600  nHIS            Number of timesteps between the writing fields
                               into history file.
 0.0000E+00  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 0.0000E+00  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 0.0000E+00  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 1.0000E-05  Akt_bak(01)     Background vertical mixing coefficient (m2/s)
                               for tracer 01: temp
 1.0000E-05  Akt_bak(02)     Background vertical mixing coefficient (m2/s)
                               for tracer 02: salt
 1.0000E-04  Akv_bak         Background vertical mixing coefficient (m2/s)
                               for momentum.
 5.0000E-06  Akk_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent energy.
 5.0000E-06  Akp_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent generic statistical field.
      3.000  gls_p           GLS stability exponent.
      1.500  gls_m           GLS turbulent kinetic energy exponent.
     -1.000  gls_n           GLS turbulent length scale exponent.
 7.6000E-06  gls_Kmin        GLS minimum value of turbulent kinetic energy.
 1.0000E-12  gls_Pmin        GLS minimum value of dissipation.
 5.4770E-01  gls_cmu0        GLS stability coefficient.
 1.4400E+00  gls_c1          GLS shear production coefficient.
 1.9200E+00  gls_c2          GLS dissipation coefficient.
-4.0000E-01  gls_c3m         GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p         GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk        GLS constant Schmidt number for TKE.
 1.3000E+00  gls_sigp        GLS constant Schmidt number for PSI.
      0.000  charnok_alpha   Charnok factor for Zos calculation.
      0.000  zos_hsig_alpha  Factor for Zos calculation using Hsig(Awave).
      0.000  sz_alpha        Factor for Wave dissipation surface tke flux .
      0.000  crgban_cw       Factor for Craig/Banner surface tke flux.
 3.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2           Quadratic bottom drag coefficient.
 0.0000E+00  Zob             Bottom roughness (m).
 0.0000E+00  Zos             Surface roughness (m).
 1.0000E+01  blk_ZQ          Height (m) of surface air humidity measurement.
 1.0000E+01  blk_ZT          Height (m) of surface air temperature measurement.
 1.0000E+01  blk_ZW          Height (m) of surface winds measurement.
          1  lmd_Jwt         Jerlov water type.
          2  Vtransform      S-coordinate transformation equation.
          1  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.
   1000.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)
 0.0000E+00  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 0.0000E+00  Tnudg(02)       Nudging/relaxation time scale (days)
                               for tracer 02: salt
 0.0000E+00  Znudg           Nudging/relaxation time scale (days)
                               for free-surface.
 0.0000E+00  M2nudg          Nudging/relaxation time scale (days)
                               for 2D momentum.
 0.0000E+00  M3nudg          Nudging/relaxation time scale (days)
                               for 3D momentum.
 0.0000E+00  obcfac          Factor between passive and active
                               open boundary conditions.
      0.000  T0              Background potential temperature (C) constant.
      0.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(idUbar)    Write out 2D U-momentum component.
          T  Hout(idVbar)    Write out 2D V-momentum component.
          T  Hout(idUvel)    Write out 3D U-momentum component.
          T  Hout(idVvel)    Write out 3D V-momentum component.
          T  Hout(idWvel)    Write out W-momentum component.
          T  Hout(idOvel)    Write out omega vertical velocity.
          T  Hout(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt
          T  Hout(idDano)    Write out density anomaly.
          T  Hout(idVvis)    Write out vertical viscosity: AKv.

 Output/Input Files:

             Output Restart File:  ocean_rst2.nc
             Output History File:  ocean_his2.nc
                 Input Grid File:  southsea2012_grd.nc
    Input Nonlinear Initial File:  southsea2012_ini.nc
           Input Forcing File 01:  southsea2012_frc.nc
             Input Boundary File:  rnt/southsea2012_bry.nc

 Tile partition information for Grid 01:  0378x0358x0032  tiling: 001x002

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1      378        1      179  2165184
        1        1      378      180      358  2165184

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50   378.50     0.50   179.50  RHO-points
        1     0.50   378.50   179.50   358.50  RHO-points

        0     1.00   378.00     0.50   179.50    U-points
        1     1.00   378.00   179.50   358.50    U-points

        0     0.50   378.50     1.00   179.50    V-points
        1     0.50   378.50   179.50   358.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     778
               HaloSizeJ(1) =     382
                TileSide(1) =     383
                TileSize(1) =   70855


 Activated C-preprocessing Options:

 southSEA            South_Sea
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 BULK_FLUXES         Surface bulk fluxes parametererization.
 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.
 EAST_FSCHAPMAN      Eastern edge, free-surface, Chapman condition.
 EAST_M2FLATHER      Eastern edge, 2D momentum, Flather condition.
 EAST_M3CLAMPED      Eastern edge, 3D momentum, Clamped condition.
 EAST_TCLAMPED       Eastern edge, tracers, Clamped condition.
 EMINUSP             Compute Salt Flux using E-P.
 GLS_MIXING          Generic Length-Scale turbulence closure.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 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.
 NORTH_FSCHAPMAN     Northern edge, free-surface, Chapman condition.
 NORTH_M2FLATHER     Northern edge, 2D momentum, Flather condition.
 NORTH_M3CLAMPED     Northern edge, 3D momentum, Clamped condition.
 NORTH_TCLAMPED      Northern edge, tracers, Clamped condition.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 SOUTH_FSCHAPMAN     Southern edge, free-surface, Chapman condition.
 SOUTH_M2FLATHER     Southern edge, 2D momentum, Flather condition.
 SOUTH_M3CLAMPED     Southern edge, 3D momentum, Clamped condition.
 SOUTH_TCLAMPED      Southern edge, tracers, Clamped condition.
 SPLINES             Conservative parabolic spline reconstruction.
 SPHERICAL           Spherical grid configuration.
 SPONGE              Enhanced horizontal mixing in the sponge areas.
 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.

 Process Information:

 Node #  0 (pid=   22488) is active.
 Node #  1 (pid=   19760) 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

    32   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    31  -0.0312500  -0.0022549       -0.168      -0.168      -6.350     -12.388
    30  -0.0625000  -0.0048948       -0.337      -0.337     -13.730     -26.839
    29  -0.0937500  -0.0080749       -0.509      -0.509     -22.557     -44.182
    28  -0.1250000  -0.0119914       -0.685      -0.685     -33.356     -65.469
    27  -0.1562500  -0.0168927       -0.866      -0.866     -46.793     -92.031
    26  -0.1875000  -0.0230908       -1.053      -1.053     -63.702    -125.538
    25  -0.2187500  -0.0309708       -1.249      -1.249     -85.114    -168.052
    24  -0.2500000  -0.0409944       -1.455      -1.455    -112.268    -222.047
    23  -0.2812500  -0.0536924       -1.675      -1.675    -146.582    -290.365
    22  -0.3125000  -0.0696342       -1.911      -1.911    -189.583    -376.055
    21  -0.3437500  -0.0893672       -2.166      -2.166    -242.736    -482.049
    20  -0.3750000  -0.1133185       -2.442      -2.442    -307.186    -610.636
    19  -0.4062500  -0.1416675       -2.740      -2.740    -383.411    -762.775
    18  -0.4375000  -0.1742175       -3.059      -3.059    -470.887    -937.413
    17  -0.4687500  -0.2103218       -3.395      -3.395    -567.880   -1131.086
    16  -0.5000000  -0.2489214       -3.745      -3.745    -671.555   -1338.123
    15  -0.5312500  -0.2887177       -4.100      -4.100    -778.435   -1551.570
    14  -0.5625000  -0.3284417       -4.455      -4.455    -885.121   -1764.630
    13  -0.5937500  -0.3671228       -4.804      -4.804    -989.014   -1972.103
    12  -0.6250000  -0.4042644       -5.146      -5.146   -1088.785   -2171.332
    11  -0.6562500  -0.4398848       -5.481      -5.481   -1184.482   -2362.414
    10  -0.6875000  -0.4744488       -5.810      -5.810   -1277.351   -2547.838
     9  -0.7187500  -0.5087464       -6.137      -6.137   -1369.506   -2731.836
     8  -0.7500000  -0.5437741       -6.469      -6.469   -1463.616   -2919.743
     7  -0.7812500  -0.5806477       -6.809      -6.809   -1562.669   -3117.537
     6  -0.8125000  -0.6205544       -7.165      -7.165   -1669.845   -3331.575
     5  -0.8437500  -0.6647399       -7.542      -7.542   -1788.478   -3568.528
     4  -0.8750000  -0.7145184       -7.948      -7.948   -1922.088   -3835.436
     3  -0.9062500  -0.7712996       -8.388      -8.388   -2074.451   -4139.847
     2  -0.9375000  -0.8366254       -8.871      -8.871   -2249.695   -4490.021
     1  -0.9687500  -0.9122123       -9.405      -9.405   -2452.416   -4895.149
     0  -1.0000000  -1.0000000      -10.000     -10.000   -2687.810   -5365.620

 Time Splitting Weights: ndtfast =   9    nfast =  13

    Primary            Secondary            Accumulated to Current Step

  1 0.0020364543700987 0.1111111111111111 0.0020364543700987 0.1111111111111111
  2 0.0081423516807323 0.1108848384033223 0.0101788060508311 0.2219959495144335
  3 0.0182864997349408 0.1099801326610188 0.0284653057857719 0.3319760821754522
  4 0.0323475955445464 0.1079482993571365 0.0608129013303183 0.4399243815325887
  5 0.0500102513402867 0.1043541220744091 0.1108231526706050 0.5442785036069977
  6 0.0706194309860016 0.0987974274810439 0.1814425836566066 0.6430759310880416
  7 0.0929932967968730 0.0909508240381548 0.2744358804534797 0.7340267551261964
  8 0.1151944667617182 0.0806182355051689 0.3896303472151979 0.8146449906313653
  9 0.1342596821693365 0.0678188503094224 0.5238900293845344 0.8824638409407878
 10 0.1458878856389088 0.0529011078461628 0.6697779150234432 0.9353649487869506
 11 0.1440867095544515 0.0366913427751730 0.8138646245778947 0.9720562915621236
 12 0.1207773749033226 0.0206817083802339 0.9346419994812172 0.9927379999423575
 13 0.0653580005187826 0.0072620000576425 0.9999999999999998 1.0000000000000000

 ndtfast, nfast =    9  13   nfast/ndtfast =  1.44444

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

    1.000000000000 1.083631825329 0.541815912665 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.00000

 Minimum X-grid spacing, DXmin =  2.38424540E+00 km
 Maximum X-grid spacing, DXmax =  8.53876089E+00 km
 Minimum Y-grid spacing, DYmin =  3.08795992E+00 km
 Maximum Y-grid spacing, DYmax =  9.94775766E+00 km
 Minimum Z-grid spacing, DZmin =  1.67524434E-01 m
 Maximum Z-grid spacing, DZmax =  4.70470903E+02 m

 Minimum barotropic Courant Number =  1.59192139E-02
 Maximum barotropic Courant Number =  4.27050946E-01
 Maximum Coriolis   Courant Number =  5.16361271E-03


 NLM: GET_STATE - Read state initial conditions,             t =     0 00:00:00
                   (File: southsea2012_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min =  1.00998491E+00 Max =  2.96956050E+01)
                - salinity
                   (Min =  2.94033954E+01 Max =  3.50299926E+01)
    GET_2DFLD   - surface u-wind component,                  t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.51582279E+00 Max =  3.39370625E+00)
    GET_2DFLD   - surface v-wind component,                  t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.12874444E+00 Max = -1.51403108E-01)
    GET_2DFLD   - surface air pressure,                      t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.00918056E+03 Max =  1.01960408E+03)
    GET_2DFLD   - solar shortwave radiation flux,            t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  2.64019335E-05 Max =  4.98503556E-05)
    GET_2DFLD   - downwelling longwave radiation flux,       t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  8.13027358E-05 Max =  1.04959257E-04)
    GET_2DFLD   - surface air temperature,                   t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.81269569E+01 Max =  2.82614824E+01)
    GET_2DFLD   - surface air relative humidity,             t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  7.37835999E-01 Max =  8.49440016E-01)
    GET_2DFLD   - rain fall rate,                            t =   345 00:00:00
                   (Rec=0012, Index=1, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.65831263E-06 Max =  1.19896346E-04)
    GET_NGFLD   - free-surface eastern boundary condition,   t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.04938970E-01 Max =  5.39443795E-02)
    GET_NGFLD   - free-surface southern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  0.00000000E+00 Max =  3.45140437E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -8.46820972E-02 Max =  2.37409635E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.05745179E-01 Max =  1.12963062E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.64431079E-01 Max =  1.32986371E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -4.25512775E-02 Max =  1.02525889E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -6.61152090E-02 Max =  5.38647064E-02)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  0.00000000E+00 Max =  1.20633476E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.54876184E-03 Max =  1.80039973E-02)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.36025823E-01 Max =  1.93699978E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.84188818E-01 Max =  4.79846932E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -7.41717477E-02 Max =  2.05983589E-02)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.40678430E-01 Max =  1.19201540E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.04374539E-03 Max =  3.01733123E-02)
    GET_NGFLD   - 3D v-momentum northern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -5.39660018E-03 Max =  4.06704033E-02)
    GET_NGFLD   - temperature eastern boundary condition,    t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  2.26234973E+00 Max =  2.92821597E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.31653714E+01 Max =  3.47828366E+01)
    GET_NGFLD   - temperature southern boundary condition,   t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.58342220E+00 Max =  2.93760120E+01)
    GET_NGFLD   - salinity southern boundary condition,      t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.03248198E+01 Max =  3.48744274E+01)
    GET_NGFLD   - temperature northern boundary condition,   t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  2.08878426E+01 Max =  2.42159190E+01)
    GET_NGFLD   - salinity northern boundary condition,      t =   345 00:00:00
                   (Rec=0012, Index=2, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  2.68756965E+01 Max =  3.47764573E+01)

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


 Initial basin volumes: TotVolume =  4.1073132888E+15 m3
                        MinVolume =  6.4173169552E+06 m3
                        MaxVolume =  2.5091605958E+10 m3
                          Max/Min =  3.9099838972E+03


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

    GET_2DFLD   - surface u-wind component,                  t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -8.47976738E+00 Max =  5.14942244E+00)
    GET_2DFLD   - surface v-wind component,                  t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -8.60785679E+00 Max =  1.70728176E-01)
    GET_2DFLD   - surface air pressure,                      t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.00913328E+03 Max =  1.01991374E+03)
    GET_2DFLD   - solar shortwave radiation flux,            t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  2.74973133E-05 Max =  5.20150240E-05)
    GET_2DFLD   - downwelling longwave radiation flux,       t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  8.29663275E-05 Max =  1.03522925E-04)
    GET_2DFLD   - surface air temperature,                   t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.60999184E+01 Max =  2.80809454E+01)
    GET_2DFLD   - surface air relative humidity,             t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  7.40901947E-01 Max =  8.58145370E-01)
    GET_2DFLD   - rain fall rate,                            t =    15 00:00:00
                   (Rec=0001, Index=2, File: southsea2012_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.55106197E-06 Max =  1.38071025E-04)
    GET_NGFLD   - free-surface eastern boundary condition,   t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.96120652E-01 Max =  5.08092590E-02)
    GET_NGFLD   - free-surface southern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  0.00000000E+00 Max =  4.25729693E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.79222626E-02 Max =  3.13438820E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.77855057E-01 Max =  2.26948792E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -4.29227736E-01 Max =  1.60314516E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.89490618E-02 Max =  1.88929858E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -8.35157011E-02 Max =  8.53220221E-02)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  0.00000000E+00 Max =  2.13904517E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.48704332E-03 Max =  2.96492477E-02)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.16814478E-01 Max =  4.92955030E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -7.32225337E-01 Max =  5.33770588E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -7.74111031E-02 Max =  2.98046285E-02)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.78695627E-01 Max =  1.84854131E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.21232575E-03 Max =  5.21252952E-02)
    GET_NGFLD   - 3D v-momentum northern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.54543610E-03 Max =  6.74086822E-02)
    GET_NGFLD   - temperature eastern boundary condition,    t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  2.26078187E+00 Max =  2.93209212E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.33929114E+01 Max =  3.48432600E+01)
    GET_NGFLD   - temperature southern boundary condition,   t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.59805859E+00 Max =  2.93209212E+01)
    GET_NGFLD   - salinity southern boundary condition,      t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.10507369E+01 Max =  3.47537850E+01)
    GET_NGFLD   - temperature northern boundary condition,   t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.78687857E+01 Max =  2.32269094E+01)
    GET_NGFLD   - salinity northern boundary condition,      t =    15 00:00:00
                   (Rec=0001, Index=1, File: southsea2012_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.19560246E+01 Max =  3.45651214E+01)
 volume =   4.924104978081719E+015 avgpe =   6.267835317610791E+019
 volume =   4.924104978081719E+015 avgpe =   6.267835317610791E+019

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
          C => (i,j,k)       Cu            Cv            Cw         Max Speed

      0     0 00:00:00  6.711673E-07  1.272888E+04  1.272888E+04  4.924105E+15
          (378,303,06)  2.006515E-03  0.000000E+00  0.000000E+00  2.028955E-01
      DEF_HIS   - creating history file: ocean_his2.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
 volume =   4.924105294627492E+015 avgpe =  NaN                    
 volume =   4.924105294627492E+015 avgpe =  NaN                    
      1     0 00:01:30  1.245361E-06           NaN           NaN  4.924105E+15
          (378,322,11)  9.248940E-04  2.067024E-04  1.258541E-02  2.703812E-01

 Blowing-up: Saving latest model state into  RESTART file

      WRT_RST   - wrote re-start fields (Index=2,2) into time record = 0000001

 Elapsed CPU time (seconds):

 Node   #  0 CPU:      20.938
 Node   #  1 CPU:      17.828
 Total:                38.766

 Nonlinear model elapsed time profile:

  Allocation and array initialization ..............         1.672  ( 4.3128 %)
  Ocean state initialization .......................         0.391  ( 1.0077 %)
  Reading of input data ............................         0.281  ( 0.7255 %)
  Processing of input data .........................         0.062  ( 0.1612 %)
  Computation of global information integrals ......         9.484  (24.4659 %)
  Writing of output data ...........................         3.266  ( 8.4240 %)
  Model 2D kernel ..................................         2.531  ( 6.5296 %)
  2D/3D coupling, vertical metrics .................         0.734  ( 1.8944 %)
  Omega vertical velocity ..........................         0.219  ( 0.5643 %)
  Equation of state for seawater ...................         1.547  ( 3.9903 %)
  Atmosphere-Ocean bulk flux parameterization ......         0.375  ( 0.9674 %)
  GLS vertical mixing parameterization .............         5.188  (13.3817 %)
  3D equations right-side terms ....................         1.250  ( 3.2245 %)
  3D equations predictor step ......................         1.375  ( 3.5470 %)
  Pressure gradient ................................         0.734  ( 1.8944 %)
  Harmonic mixing of tracers, geopotentials ........         1.641  ( 4.2322 %)
  Harmonic stress tensor, S-surfaces ...............         0.672  ( 1.7332 %)
  Corrector time-step for 3D momentum ..............         0.766  ( 1.9750 %)
  Corrector time-step for tracers ..................         1.031  ( 2.6602 %)
                                              Total:        33.219   85.6913

 Nonlinear model message Passage profile:

  Message Passage: 2D halo exchanges ...............         0.047  ( 0.1209 %)
  Message Passage: 3D halo exchanges ...............         0.047  ( 0.1209 %)
  Message Passage: 4D halo exchanges ...............         0.016  ( 0.0403 %)
  Message Passage: data broadcast ..................         0.188  ( 0.4837 %)
  Message Passage: data reduction ..................         0.016  ( 0.0403 %)
  Message Passage: data gathering ..................         2.750  ( 7.0939 %)
  Message Passage: data scattering..................         1.094  ( 2.8214 %)
                                              Total:         4.156   10.7215

 All percentages are with respect to total time =           38.766

 ROMS/TOMS - Output NetCDF summary for Grid 01:
             number of time records written in HISTORY file = 00000001
             number of time records written in RESTART file = 00000001

 Analytical header files used:

     ROMS/Functionals/ana_btflux.h
     ROMS/Functionals/ana_hmixcoef.h
     ROMS/Functionals/ana_srflux.h

 ROMS/TOMS: DONE... Thursday - August 23, 2012 -  8:06:29 AM
=====================================================================================================
Last edited by dhd_1979 on Fri Aug 24, 2012 3:21 pm, edited 1 time in total.

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blowing up problem

#2 Unread post by kate »

Did you change the timestep when you changed the resolution? How about the bathymetry steepness factors - did they change? You can also adjust the horizontal viscosity and diffusion coefficients.

dhd_1979
Posts: 5
Joined: Sun Dec 26, 2010 9:38 pm
Location: Marine Blue Australia

Re: Blowing up problem

#3 Unread post by dhd_1979 »

Thanks Kate!
I have change the time step from 720s to 90s, and the Maximum grid stiffness ratios: rx0 = 3.069334E-01 (Beckmann and Haidvogel) rx1 = 9.416266E+00 (Haney).

The SST and speed at the beginging is
屏幕快照 2012-08-24 下午10.58.48.png
屏幕快照 2012-08-24 下午10.58.48.png (26.52 KiB) Viewed 5032 times
at time record 2 (the total time = 2*90s)is:
屏幕快照 2012-08-24 下午10.55.17.png
屏幕快照 2012-08-24 下午10.55.17.png (29.58 KiB) Viewed 5032 times
at time record 3 (the total time = 4*90s)is:
屏幕快照 2012-08-24 下午10.59.02.png
屏幕快照 2012-08-24 下午10.59.02.png (11.2 KiB) Viewed 5032 times
.

It seems that the NaN propagates from the east to west, so I have test the east boundary for both close and open conditions, but the result not changed.

I also change the depth of the bathymetry to a constant 800m, the result still not changed much.

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blowing up problem

#4 Unread post by kate »

You also have trouble moving west to east. The trouble is moving as a barotropic wave, in the depth-integrated equations. At least it is happening early, which makes it easier to investigate. I don't know what it is - perhaps you can let us know when you fix it.

LiuZHQ
Posts: 46
Joined: Mon Apr 06, 2009 5:12 pm
Location: The Hong Kong University of Sci&Tech

Re: Blowing up problem

#5 Unread post by LiuZHQ »

Did you smooth your bathymetry?

dhd_1979
Posts: 5
Joined: Sun Dec 26, 2010 9:38 pm
Location: Marine Blue Australia

Re: Blowing up problem

#6 Unread post by dhd_1979 »

LiuZHQ wrote:Did you smooth your bathymetry?
Thanks LiuZHQ!

I smoothed the bathymetry with rfac=0.35(maximum slope of topography allowed in smoothing), and I have a test with the bathymetry to be a constant 800m, it is no useful.

I changed the initial state from levitas climate data to hycom exp08, but still no useful.

Post Reply