Model blowing up using elevation at the boundary

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
paim_oceano
Posts: 10
Joined: Mon Jun 16, 2014 4:33 pm
Location: University of Delaware

Model blowing up using elevation at the boundary

#1 Unread post by paim_oceano »

Dear all,

I'm trying to run COAWST (coupling ROMS and SWAN) using free-surface information from HYCOM at the boundary.
In my first test, I've created a file with information just in the East boundary and it worked.
For a second test, I've created a file with information also in the North and South boundary, but in this case the model blows immediately during the first time step. I've set the boundaries to be Chapman and Flather.

I've tried to decrease the time step, but even using 1 second it still blows in the first time step. I've looked at the results for this first time step, trying to figure what is causing the blow, but I can't seem to find anything out of normal.

Do you guys have any idea about what could be causing it or at least how can I try to figure it out?
Below is the log generated

Thanks a million!

Code: Select all

 =============================================================
 Coupled Input File name = coupling_inlet_test.in

 Model Coupling: 


       Ocean Model MPI nodes: 000 - 063

       Waves Model MPI nodes: 064 - 087
 Process Information:

 Node #  0 (pid=   18176) is active.
 Node #  4 (pid=   18180) is active.
 Node #  7 (pid=   18183) is active.
 Node # 10 (pid=   18186) is active.
 Node # 11 (pid=   18187) is active.
 Node # 13 (pid=   18189) is active.
 Node # 14 (pid=   18190) is active.
 Node # 15 (pid=   18191) is active.
 Node # 16 (pid=   18192) is active.
 Node # 17 (pid=   18193) is active.
 Node # 18 (pid=   18194) is active.
 Node # 19 (pid=   18195) is active.
 Node # 20 (pid=   18196) is active.
 Node # 21 (pid=   18197) is active.
 Node # 22 (pid=   18198) is active.
 Node # 23 (pid=   18199) is active.
 Node #  1 (pid=   18177) is active.
 Node #  2 (pid=   18178) is active.
 Node #  3 (pid=   18179) is active.
 Node #  5 (pid=   18181) is active.
 Node #  6 (pid=   18182) is active.
 Node #  8 (pid=   18184) is active.
 Node #  9 (pid=   18185) is active.
 Node # 12 (pid=   18188) is active.
 Node # 48 (pid=   10489) is active.
 Node # 49 (pid=   10490) is active.
 Node # 50 (pid=   10491) is active.
 Node # 51 (pid=   10492) is active.
 Node # 52 (pid=   10493) is active.
 Node # 53 (pid=   10494) is active.
 Node # 54 (pid=   10495) is active.
 Node # 55 (pid=   10496) is active.
 Node # 56 (pid=   10497) is active.
 Node # 57 (pid=   10498) is active.
 Node # 58 (pid=   10499) is active.
 Node # 59 (pid=   10500) is active.
 Node # 60 (pid=   10501) is active.
 Node # 61 (pid=   10502) is active.
 Node # 62 (pid=   10503) is active.
 Node # 63 (pid=   10504) is active.
 Node # 24 (pid=   22731) is active.
 Node # 25 (pid=   22732) is active.
 Node # 26 (pid=   22733) is active.
 Node # 27 (pid=   22734) is active.
 Node # 28 (pid=   22735) is active.
 Node # 29 (pid=   22736) is active.
 Node # 30 (pid=   22737) is active.
 Node # 31 (pid=   22738) is active.
 Node # 32 (pid=   22739) is active.
 Node # 33 (pid=   22740) is active.
 Node # 34 (pid=   22741) is active.
 Node # 35 (pid=   22742) is active.
 Node # 36 (pid=   22743) is active.
 Node # 37 (pid=   22744) is active.
 Node # 38 (pid=   22745) is active.
 Node # 39 (pid=   22746) is active.
 Node # 40 (pid=   22747) is active.
 Node # 41 (pid=   22748) is active.
 Node # 42 (pid=   22749) is active.
 Node # 43 (pid=   22750) is active.
 Node # 44 (pid=   22751) is active.
 Node # 45 (pid=   22752) is active.
 Node # 46 (pid=   22753) is active.
 Node # 47 (pid=   22754) is active.

 Model Input Parameters:  ROMS/TOMS version 3.4  
                          Thursday - June 18, 2015 -  6:25:44 PM
 -----------------------------------------------------------------------------

SWAN grid   1 is preparing computation


 delbayCC2012

 Operating system : Linux
 CPU/hardware     : x86_64
 Compiler system  : pgi
 Compiler command : /opt/shared/openmpi/1.4.4-pgi/bin/mpif90
 Compiler flags   :  -fastsse -Mipa=fast -tp k8-64 -Mfree

 Input Script  : ocean.in

 SVN Root URL  : https://www.myroms.org/svn/omlab/branches/jcwarner
 SVN Revision  : exported

 Local Root    : /lustre/work/kukulka_lab/androd/COAWST/Rob_Delbay/RobCOAWST
 Header Dir    : /lustre/work/kukulka_lab/androd/COAWST/Rob_Delbay/RobCOAWST/Projects/delbayCC201
 Header file   : delbaycc2012.h
 Analytical Dir: /lustre/work/kukulka_lab/androd/COAWST/Rob_Delbay/RobCOAWST/Projects/delbayCC201

 Resolution, Grid 01: 0148x0298x010,  Parallel Nodes:  64,  Tiling: 008x008


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

     192961  ntimes          Number of timesteps for 3-D equations.
    150.000  dt              Timestep size (s) for 3-D equations.
         30  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
   1800.000  TI_WAV_OCN      Time interval (s) between coupling WAV-OCN models.
         12  nOCN_WAV        Number of OCN timesteps between coupling to WAV.
          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.
        576  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).
         24  nHIS            Number of timesteps between the writing fields
                               into history file.
      17520  ndefHIS         Number of timesteps between creation of new
                               history files.
 5.0000E+00  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 5.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-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.
 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.5000E-04  rdrg            Linear bottom drag coefficient (m/s).
 5.0000E-03  rdrg2           Quadratic bottom drag coefficient.
 0.0000E+00  Zob             Bottom roughness (m).
 1.0000E+00  Dcrit           Minimum depth for wetting and drying (m).
          2  Vtransform      S-coordinate transformation equation.
          4  Vstretching     S-coordinate stretching function.
 3.0000E+00  theta_s         S-coordinate surface control parameter.
 4.0000E-01  theta_b         S-coordinate bottom  control parameter.
      0.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).
  -2191.042  tide_start      Reference time origin for tidal forcing (days).
20120101.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.
     10.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(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(idTvar)    Write out tracer 02: salt

 Output/Input Files:

             Output Restart File:  Output/ocean_rst.nc
        Prefix for History Files:  Output/ocean_his
        Physical parameters File:  ocean.in
                 Input Grid File:  delbay_grd.nc
    Input Nonlinear Initial File:  delbayini2012CCNAM.nc
             Input Boundary File:  delbay_bry_hycom2012_ENSfsm2.nc

 Tile partition information for Grid 01:  0148x0298x0010  tiling: 008x008

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       17        1       35     5950
        1       18       36        1       35     6650
        2       37       55        1       35     6650
        3       56       74        1       35     6650
        4       75       93        1       35     6650
        5       94      112        1       35     6650
        6      113      131        1       35     6650
        7      132      148        1       35     5950
        8        1       17       36       73     6460
        9       18       36       36       73     7220
       10       37       55       36       73     7220
       11       56       74       36       73     7220
       12       75       93       36       73     7220
       13       94      112       36       73     7220
       14      113      131       36       73     7220
       15      132      148       36       73     6460
       16        1       17       74      111     6460
       17       18       36       74      111     7220
       18       37       55       74      111     7220
       19       56       74       74      111     7220
       20       75       93       74      111     7220
       21       94      112       74      111     7220
       22      113      131       74      111     7220
       23      132      148       74      111     6460
       24        1       17      112      149     6460
       25       18       36      112      149     7220
       26       37       55      112      149     7220
       27       56       74      112      149     7220
       28       75       93      112      149     7220
       29       94      112      112      149     7220
       30      113      131      112      149     7220
       31      132      148      112      149     6460
       32        1       17      150      187     6460
       33       18       36      150      187     7220
       34       37       55      150      187     7220
       35       56       74      150      187     7220
       36       75       93      150      187     7220
       37       94      112      150      187     7220
       38      113      131      150      187     7220
       39      132      148      150      187     6460
       40        1       17      188      225     6460
       41       18       36      188      225     7220
       42       37       55      188      225     7220
       43       56       74      188      225     7220
       44       75       93      188      225     7220
       45       94      112      188      225     7220
       46      113      131      188      225     7220
       47      132      148      188      225     6460
       48        1       17      226      263     6460
       49       18       36      226      263     7220
       50       37       55      226      263     7220
       51       56       74      226      263     7220
       52       75       93      226      263     7220
       53       94      112      226      263     7220
       54      113      131      226      263     7220
       55      132      148      226      263     6460
       56        1       17      264      298     5950
       57       18       36      264      298     6650
       58       37       55      264      298     6650
       59       56       74      264      298     6650
       60       75       93      264      298     6650
       61       94      112      264      298     6650
       62      113      131      264      298     6650
       63      132      148      264      298     5950

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    17.50     0.50    35.50  RHO-points
        1    17.50    36.50     0.50    35.50  RHO-points
        2    36.50    55.50     0.50    35.50  RHO-points
        3    55.50    74.50     0.50    35.50  RHO-points
        4    74.50    93.50     0.50    35.50  RHO-points
        5    93.50   112.50     0.50    35.50  RHO-points
        6   112.50   131.50     0.50    35.50  RHO-points
        7   131.50   148.50     0.50    35.50  RHO-points
        8     0.50    17.50    35.50    73.50  RHO-points
        9    17.50    36.50    35.50    73.50  RHO-points
       10    36.50    55.50    35.50    73.50  RHO-points
       11    55.50    74.50    35.50    73.50  RHO-points
       12    74.50    93.50    35.50    73.50  RHO-points
       13    93.50   112.50    35.50    73.50  RHO-points
       14   112.50   131.50    35.50    73.50  RHO-points
       15   131.50   148.50    35.50    73.50  RHO-points
       16     0.50    17.50    73.50   111.50  RHO-points
       17    17.50    36.50    73.50   111.50  RHO-points
       18    36.50    55.50    73.50   111.50  RHO-points
       19    55.50    74.50    73.50   111.50  RHO-points
       20    74.50    93.50    73.50   111.50  RHO-points
       21    93.50   112.50    73.50   111.50  RHO-points
       22   112.50   131.50    73.50   111.50  RHO-points
       23   131.50   148.50    73.50   111.50  RHO-points
       24     0.50    17.50   111.50   149.50  RHO-points
       25    17.50    36.50   111.50   149.50  RHO-points
       26    36.50    55.50   111.50   149.50  RHO-points
       27    55.50    74.50   111.50   149.50  RHO-points
       28    74.50    93.50   111.50   149.50  RHO-points
       29    93.50   112.50   111.50   149.50  RHO-points
       30   112.50   131.50   111.50   149.50  RHO-points
       31   131.50   148.50   111.50   149.50  RHO-points
       32     0.50    17.50   149.50   187.50  RHO-points
       33    17.50    36.50   149.50   187.50  RHO-points
       34    36.50    55.50   149.50   187.50  RHO-points
       35    55.50    74.50   149.50   187.50  RHO-points
       36    74.50    93.50   149.50   187.50  RHO-points
       37    93.50   112.50   149.50   187.50  RHO-points
       38   112.50   131.50   149.50   187.50  RHO-points
       39   131.50   148.50   149.50   187.50  RHO-points
       40     0.50    17.50   187.50   225.50  RHO-points
       41    17.50    36.50   187.50   225.50  RHO-points
       42    36.50    55.50   187.50   225.50  RHO-points
       43    55.50    74.50   187.50   225.50  RHO-points
       44    74.50    93.50   187.50   225.50  RHO-points
       45    93.50   112.50   187.50   225.50  RHO-points
       46   112.50   131.50   187.50   225.50  RHO-points
       47   131.50   148.50   187.50   225.50  RHO-points
       48     0.50    17.50   225.50   263.50  RHO-points
       49    17.50    36.50   225.50   263.50  RHO-points
       50    36.50    55.50   225.50   263.50  RHO-points
       51    55.50    74.50   225.50   263.50  RHO-points
       52    74.50    93.50   225.50   263.50  RHO-points
       53    93.50   112.50   225.50   263.50  RHO-points
       54   112.50   131.50   225.50   263.50  RHO-points
       55   131.50   148.50   225.50   263.50  RHO-points
       56     0.50    17.50   263.50   298.50  RHO-points
       57    17.50    36.50   263.50   298.50  RHO-points
       58    36.50    55.50   263.50   298.50  RHO-points
       59    55.50    74.50   263.50   298.50  RHO-points
       60    74.50    93.50   263.50   298.50  RHO-points
       61    93.50   112.50   263.50   298.50  RHO-points
       62   112.50   131.50   263.50   298.50  RHO-points
       63   131.50   148.50   263.50   298.50  RHO-points

        0     1.00    17.50     0.50    35.50    U-points
        1    17.50    36.50     0.50    35.50    U-points
        2    36.50    55.50     0.50    35.50    U-points
        3    55.50    74.50     0.50    35.50    U-points
        4    74.50    93.50     0.50    35.50    U-points
        5    93.50   112.50     0.50    35.50    U-points
        6   112.50   131.50     0.50    35.50    U-points
        7   131.50   148.00     0.50    35.50    U-points
        8     1.00    17.50    35.50    73.50    U-points
        9    17.50    36.50    35.50    73.50    U-points
       10    36.50    55.50    35.50    73.50    U-points
       11    55.50    74.50    35.50    73.50    U-points
       12    74.50    93.50    35.50    73.50    U-points
       13    93.50   112.50    35.50    73.50    U-points
       14   112.50   131.50    35.50    73.50    U-points
       15   131.50   148.00    35.50    73.50    U-points
       16     1.00    17.50    73.50   111.50    U-points
       17    17.50    36.50    73.50   111.50    U-points
       18    36.50    55.50    73.50   111.50    U-points
       19    55.50    74.50    73.50   111.50    U-points
       20    74.50    93.50    73.50   111.50    U-points
       21    93.50   112.50    73.50   111.50    U-points
       22   112.50   131.50    73.50   111.50    U-points
       23   131.50   148.00    73.50   111.50    U-points
       24     1.00    17.50   111.50   149.50    U-points
       25    17.50    36.50   111.50   149.50    U-points
       26    36.50    55.50   111.50   149.50    U-points
       27    55.50    74.50   111.50   149.50    U-points
       28    74.50    93.50   111.50   149.50    U-points
       29    93.50   112.50   111.50   149.50    U-points
       30   112.50   131.50   111.50   149.50    U-points
       31   131.50   148.00   111.50   149.50    U-points
       32     1.00    17.50   149.50   187.50    U-points
       33    17.50    36.50   149.50   187.50    U-points
       34    36.50    55.50   149.50   187.50    U-points
       35    55.50    74.50   149.50   187.50    U-points
       36    74.50    93.50   149.50   187.50    U-points
       37    93.50   112.50   149.50   187.50    U-points
       38   112.50   131.50   149.50   187.50    U-points
       39   131.50   148.00   149.50   187.50    U-points
       40     1.00    17.50   187.50   225.50    U-points
       41    17.50    36.50   187.50   225.50    U-points
       42    36.50    55.50   187.50   225.50    U-points
       43    55.50    74.50   187.50   225.50    U-points
       44    74.50    93.50   187.50   225.50    U-points
       45    93.50   112.50   187.50   225.50    U-points
       46   112.50   131.50   187.50   225.50    U-points
       47   131.50   148.00   187.50   225.50    U-points
       48     1.00    17.50   225.50   263.50    U-points
       49    17.50    36.50   225.50   263.50    U-points
       50    36.50    55.50   225.50   263.50    U-points
       51    55.50    74.50   225.50   263.50    U-points
       52    74.50    93.50   225.50   263.50    U-points
       53    93.50   112.50   225.50   263.50    U-points
       54   112.50   131.50   225.50   263.50    U-points
       55   131.50   148.00   225.50   263.50    U-points
       56     1.00    17.50   263.50   298.50    U-points
       57    17.50    36.50   263.50   298.50    U-points
       58    36.50    55.50   263.50   298.50    U-points
       59    55.50    74.50   263.50   298.50    U-points
       60    74.50    93.50   263.50   298.50    U-points
       61    93.50   112.50   263.50   298.50    U-points
       62   112.50   131.50   263.50   298.50    U-points
       63   131.50   148.00   263.50   298.50    U-points

        0     0.50    17.50     1.00    35.50    V-points
        1    17.50    36.50     1.00    35.50    V-points
        2    36.50    55.50     1.00    35.50    V-points
        3    55.50    74.50     1.00    35.50    V-points
        4    74.50    93.50     1.00    35.50    V-points
        5    93.50   112.50     1.00    35.50    V-points
        6   112.50   131.50     1.00    35.50    V-points
        7   131.50   148.50     1.00    35.50    V-points
        8     0.50    17.50    35.50    73.50    V-points
        9    17.50    36.50    35.50    73.50    V-points
       10    36.50    55.50    35.50    73.50    V-points
       11    55.50    74.50    35.50    73.50    V-points
       12    74.50    93.50    35.50    73.50    V-points
       13    93.50   112.50    35.50    73.50    V-points
       14   112.50   131.50    35.50    73.50    V-points
       15   131.50   148.50    35.50    73.50    V-points
       16     0.50    17.50    73.50   111.50    V-points
       17    17.50    36.50    73.50   111.50    V-points
       18    36.50    55.50    73.50   111.50    V-points
       19    55.50    74.50    73.50   111.50    V-points
       20    74.50    93.50    73.50   111.50    V-points
       21    93.50   112.50    73.50   111.50    V-points
       22   112.50   131.50    73.50   111.50    V-points
       23   131.50   148.50    73.50   111.50    V-points
       24     0.50    17.50   111.50   149.50    V-points
       25    17.50    36.50   111.50   149.50    V-points
       26    36.50    55.50   111.50   149.50    V-points
       27    55.50    74.50   111.50   149.50    V-points
       28    74.50    93.50   111.50   149.50    V-points
       29    93.50   112.50   111.50   149.50    V-points
       30   112.50   131.50   111.50   149.50    V-points
       31   131.50   148.50   111.50   149.50    V-points
       32     0.50    17.50   149.50   187.50    V-points
       33    17.50    36.50   149.50   187.50    V-points
       34    36.50    55.50   149.50   187.50    V-points
       35    55.50    74.50   149.50   187.50    V-points
       36    74.50    93.50   149.50   187.50    V-points
       37    93.50   112.50   149.50   187.50    V-points
       38   112.50   131.50   149.50   187.50    V-points
       39   131.50   148.50   149.50   187.50    V-points
       40     0.50    17.50   187.50   225.50    V-points
       41    17.50    36.50   187.50   225.50    V-points
       42    36.50    55.50   187.50   225.50    V-points
       43    55.50    74.50   187.50   225.50    V-points
       44    74.50    93.50   187.50   225.50    V-points
       45    93.50   112.50   187.50   225.50    V-points
       46   112.50   131.50   187.50   225.50    V-points
       47   131.50   148.50   187.50   225.50    V-points
       48     0.50    17.50   225.50   263.50    V-points
       49    17.50    36.50   225.50   263.50    V-points
       50    36.50    55.50   225.50   263.50    V-points
       51    55.50    74.50   225.50   263.50    V-points
       52    74.50    93.50   225.50   263.50    V-points
       53    93.50   112.50   225.50   263.50    V-points
       54   112.50   131.50   225.50   263.50    V-points
       55   131.50   148.50   225.50   263.50    V-points
       56     0.50    17.50   263.50   298.00    V-points
       57    17.50    36.50   263.50   298.00    V-points
       58    36.50    55.50   263.50   298.00    V-points
       59    55.50    74.50   263.50   298.00    V-points
       60    74.50    93.50   263.50   298.00    V-points
       61    93.50   112.50   263.50   298.00    V-points
       62   112.50   131.50   263.50   298.00    V-points
       63   131.50   148.50   263.50   298.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =      62
               HaloSizeJ(1) =     100
                TileSide(1) =      44
                TileSize(1) =    1100


 Activated C-preprocessing Options:

 DELBAYCC2012        delbayCC2012
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_SMFLUX          Analytical kinematic surface momentum flux.
 ANA_SSFLUX          Analytical kinematic surface salinity flux.
 ANA_STFLUX          Analytical kinematic surface temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 CURVGRID            Orthogonal curvilinear grid.
 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_M3RADIATION    Eastern edge, 3D momentum, radiation condition.
 EAST_TRADIATION     Eastern edge, tracers, radiation condition.
 FSOBC_REDUCED       Using free-surface data in reduced physics condtions
 INLINE_2DIO         Processing 3D IO level by level to reduce memory needs.
 MASKING             Land/Sea masking.
 MCT_LIB             Using Model Coupling Toolkit library.
 MIX_GEO_TS          Mixing of tracers along geopotential surfaces.
 MIX_GEO_UV          Mixing of momentum along geopotential surfaces.
 MPI                 MPI distributed-memory configuration.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 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_M3RADIATION   Northern edge, 3D momentum, radiation condition.
 NORTH_TRADIATION    Northern edge, tracers, radiation condition.
 OUT_DOUBLE          Double precision output fields in NetCDF files.
 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.
 SOLVE3D             Solving 3D Primitive Equations.
 SOUTH_FSCHAPMAN     Southern edge, free-surface, Chapman condition.
 SOUTH_M2FLATHER     Southern edge, 2D momentum, Flather condition.
 SOUTH_M3RADIATION   Southern edge, 3D momentum, radiation condition.
 SOUTH_TRADIATION    Southern edge, tracers, radiation condition.
 SPLINES             Conservative parabolic spline reconstruction.
 SWAN_COUPLING       SWAN model coupling.
 TS_A4HADVECTION     Fourth-order Akima horizontal advection of tracers.
 TS_A4VADVECTION     Fourth-order Akima 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_TIDES            Add tidal currents to 2D momentum climatologies.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 WAVES_OCEAN         Two-way wave-ocean models coupling.
 WET_DRY             Wetting and drying activated.

 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

    10   0.0000000   0.0000000        0.000       0.000       0.000       0.000
     9  -0.1000000  -0.0060605       -0.018       0.000      -0.312      -0.606
     8  -0.2000000  -0.0247149       -0.074       0.000      -1.273      -2.471
     7  -0.3000000  -0.0573991       -0.172       0.000      -2.956      -5.740
     6  -0.4000000  -0.1065532       -0.320       0.000      -5.487     -10.655
     5  -0.5000000  -0.1756670       -0.527       0.000      -9.047     -17.567
     4  -0.6000000  -0.2692794       -0.808       0.000     -13.868     -26.928
     3  -0.7000000  -0.3928500       -1.179       0.000     -20.232     -39.285
     2  -0.8000000  -0.5523651       -1.657       0.000     -28.447     -55.237
     1  -0.9000000  -0.7534850       -2.260       0.000     -38.804     -75.349
     0  -1.0000000  -1.0000000       -3.000       0.000     -51.500    -100.000

 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.0093366137069498 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.0116826306292293 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.0857756755402367 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399431 0.6603333185976138
 21 0.0256737146312911 0.0297153720153352 0.1342125541712341 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524358 0.7189082721405746
 23 0.0316613792205220 0.0279045862015855 0.1945237931729577 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236652 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893980 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
 27 0.0432936737565710 0.0230872677537126 0.3506756411451923 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772279 0.8685254094681656
 29 0.0481720587108285 0.0201154903974257 0.4447073467880564 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256414 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325520 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348097 0.9524525304259084
 34 0.0531891349131380 0.0115726061288397 0.7060109510479476 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213237 0.9738247715198165
 36 0.0503158038019030 0.0080591141492892 0.8085423793232267 0.9818838856691057
 37 0.0473537721847153 0.0063819206892258 0.8558961515079421 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497609 0.9930692679747334
 39 0.0376397765791564 0.0033640675316746 0.9367177506289173 0.9964333355064080
 40 0.0305543017255206 0.0021094083123694 0.9672720523544379 0.9985427438187774
 41 0.0217382098544504 0.0010909315881854 0.9890102622088883 0.9996336754069628
 42 0.0109897377911118 0.0003663245930371 1.0000000000000002 0.9999999999999999

 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 =  6.83656115E-01 km
 Maximum X-grid spacing, DXmax =  7.98995075E+00 km
 Minimum Y-grid spacing, DYmin =  5.84204633E-01 km
 Maximum Y-grid spacing, DYmax =  4.61661121E+00 km
 Minimum Z-grid spacing, DZmin =  1.81813631E-02 m
 Maximum Z-grid spacing, DZmax =  2.46514996E+01 m

 Minimum barotropic Courant Number =  3.10445284E-02
 Maximum barotropic Courant Number =  2.70482781E-01
 Maximum Coriolis   Courant Number =  1.39805988E-02


 NLM: GET_STATE - Read state initial conditions,             t =     0 00:00:00
                   (File: delbayini2012CCNAM.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -2.00000000E+00 Max =  1.17945957E+00)
                - vertically integrated u-momentum component
                   (Min = -6.78342130E-01 Max =  6.58236907E-01)
                - vertically integrated v-momentum component
                   (Min = -5.32789570E-01 Max =  1.49376311E-01)
                - u-momentum component
                   (Min = -8.47913088E-01 Max =  1.19578593E+00)
                - v-momentum component
                   (Min = -9.65962735E-01 Max =  2.50923286E-01)
                - potential temperature
                   (Min = -8.73054692E-02 Max =  9.84094524E+00)
                - salinity
                   (Min = -3.05489670E-01 Max =  3.44343500E+01)
    GET_NGFLD   - tidal period
                   (Min =  1.49040000E+04 Max =  9.67320000E+04)
    GET_2DFLD   - tidal current inclination angle
                   (Min =  0.00000000E+00 Max =  6.28318531E+00)
    GET_2DFLD   - tidal current phase angle
                   (Min =  0.00000000E+00 Max =  6.28317864E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Min =  0.00000000E+00 Max =  1.91539911E+00)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Min = -3.28818765E-01 Max =  2.26390060E-01)
    GET_NGFLD   - free-surface eastern boundary condition,   t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -6.24522797E-01 Max = -5.56181716E-01)
    GET_NGFLD   - free-surface southern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -6.64299099E-01 Max = -6.13464345E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -5.67338138E-01 Max = -5.10653361E-01)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -1.78650608E-01 Max =  2.38413242E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -4.05335075E-01 Max =  2.60545953E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -1.56434706E-01 Max =  3.81790173E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -2.16383088E-01 Max =  1.65037623E-01)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -2.25385881E-02 Max =  1.06998015E-01)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -2.31242745E-01 Max = -2.67002905E-02)

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


 Initial basin volumes: TotVolume =  1.5638747369E+12 m3
                        MinVolume =  9.6644887677E+03 m3
                        MaxVolume =  8.6750340884E+08 m3
                          Max/Min =  8.9761955308E+04

 == SWAN grid  1 sent wave data to ROMS
 ** ROMS grid  1 recv data from SWAN
 ** ROMS grid  1 sent data to SWAN


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00192961)
    GET_NGFLD   - free-surface eastern boundary condition,   t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -6.27662354E-01 Max = -5.51600357E-01)
    GET_NGFLD   - free-surface southern boundary condition,  t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -6.69472724E-01 Max = -6.10971062E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -5.59730947E-01 Max = -5.06058428E-01)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -1.90737257E-01 Max =  1.25957661E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -3.49767328E-01 Max =  1.98288440E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -2.00448120E-01 Max = -1.47098131E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -1.93069435E-01 Max =  2.10398759E-01)
 == SWAN grid  1 recvd data from ROMS
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -5.45434644E-02 Max =  7.34075623E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =     0 03:00:00
                   (Rec=0002, Index=1, File: delbay_bry_hycom2012_ENSfsm2.nc)
                   (Tmin=          0.0000 Tmax=        365.8750)
                   (Min = -1.79090935E-01 Max = -1.54667225E-02)
+time 20120101.001500   , step      1; iteration    1; sweep 1 grid   1

  GRID  STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME

   1       0     0 00:00:00  4.965835E-03  2.557462E+02  2.557512E+02  1.645769E+12
      DEF_HIS   - creating history file: Output/ocean_his_0001.nc
+time 20120101.001500   , step      1; iteration    1; sweep 2 grid   1
+time 20120101.001500   , step      1; iteration    1; sweep 3 grid   1
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
+time 20120101.001500   , step      1; iteration    1; sweep 4 grid   1
   1       1     0 00:02:30           NaN           NaN           NaN           NaN

 ROMS/TOMS - Blows up ................ exit_flag:   1


 Blowing-up: Saving latest model state into  RESTART file

+time 20120101.003000   , step      2; iteration    1; sweep 1 grid   1
+time 20120101.003000   , step      2; iteration    1; sweep 2 grid   1
 Node   # 23 CPU:       5.568
 Node   #  1 CPU:       5.568
 Node   #  2 CPU:       5.568
 Node   # 12 CPU:       5.568
 Node   #  3 CPU:       5.568
 Node   #  5 CPU:       5.568
 Node   #  4 CPU:       5.568
 Node   #  6 CPU:       5.568
 Node   #  8 CPU:       5.568
 Node   #  7 CPU:       5.568
 Node   # 10 CPU:       5.568
 Node   #  9 CPU:       5.568
 Node   # 11 CPU:       5.568
 Node   # 13 CPU:       5.568
 Node   # 14 CPU:       5.568
 Node   # 15 CPU:       5.568
 Node   # 16 CPU:       5.569
 Node   # 17 CPU:       5.568
 Node   # 18 CPU:       5.568
 Node   # 19 CPU:       5.568
 Node   # 20 CPU:       5.568
 Node   # 21 CPU:       5.568
 Node   # 22 CPU:       5.568
      WRT_RST   - wrote re-start fields (Index=1,2) into time record = 0000001

 Elapsed CPU time (seconds):

 Node   #  0 CPU:       5.569
 Node   # 24 CPU:       5.568
 Node   # 25 CPU:       5.568
 Node   # 48 CPU:       5.569
 Node   # 26 CPU:       5.568
 Node   # 49 CPU:       5.569
 Node   # 27 CPU:       5.568
 Node   # 50 CPU:       5.569
 Node   # 28 CPU:       5.568
 Node   # 51 CPU:       5.569
 Node   # 29 CPU:       5.568
 Node   # 52 CPU:       5.569
 Node   # 30 CPU:       5.568
 Node   # 53 CPU:       5.569
 Node   # 31 CPU:       5.568
 Node   # 54 CPU:       5.569
 Node   # 32 CPU:       5.568
 Node   # 55 CPU:       5.569
 Node   # 33 CPU:       5.568
 Node   # 56 CPU:       5.569
 Node   # 34 CPU:       5.568
 Node   # 57 CPU:       5.569
 Node   # 35 CPU:       5.568
 Node   # 58 CPU:       5.569
 Node   # 36 CPU:       5.568
 Node   # 59 CPU:       5.569
 Node   # 37 CPU:       5.568
 Node   # 60 CPU:       5.569
 Node   # 38 CPU:       5.568
 Node   # 61 CPU:       5.569
 Node   # 39 CPU:       5.568
 Node   # 62 CPU:       5.569
 Node   # 40 CPU:       5.568
 Node   # 63 CPU:       5.569
 Node   # 41 CPU:       5.568
 Node   # 42 CPU:       5.569
 Node   # 43 CPU:       5.568
 Node   # 44 CPU:       5.568
 Node   # 45 CPU:       5.568
 Node   # 46 CPU:       5.568
 Node   # 47 CPU:       5.568
 Total:               356.370

 Nonlinear model elapsed time profile:

  Initialization ...................................       163.414  (45.8552 %)
  OI data assimilation .............................         2.222  ( 0.6236 %)
  Reading of input data ............................        95.292  (26.7396 %)
  Processing of input data .........................         0.983  ( 0.2759 %)
  Computation of vertical boundary conditions ......         0.215  ( 0.0604 %)
  Computation of global information integrals ......         0.760  ( 0.2132 %)
  Writing of output data ...........................        53.810  (15.0994 %)
  Model 2D kernel ..................................        13.850  ( 3.8863 %)
  Tidal forcing ....................................         0.277  ( 0.0776 %)
  2D/3D coupling, vertical metrics .................         0.129  ( 0.0363 %)
  Omega vertical velocity ..........................         0.338  ( 0.0949 %)
  Equation of state for seawater ...................         0.403  ( 0.1131 %)
  My2.5 vertical mixing parameterization ...........         1.133  ( 0.3180 %)
  3D equations right-side terms ....................         0.320  ( 0.0897 %)
  3D equations predictor step ......................         2.074  ( 0.5819 %)
  Pressure gradient ................................         1.360  ( 0.3818 %)
  Harmonic mixing of tracers, geopotentials ........         3.138  ( 0.8805 %)
  Harmonic stress tensor, geopotentials ............         3.393  ( 0.9520 %)
  Corrector time-step for 3D momentum ..............         0.781  ( 0.2192 %)
  Corrector time-step for tracers ..................         0.424  ( 0.1189 %)
                                              Total:       344.316   96.6174

 Nonlinear model message Passage profile:

  Message Passage: 2D halo exchanges ...............         8.907  ( 2.4995 %)
  Message Passage: 3D halo exchanges ...............         1.453  ( 0.4078 %)
  Message Passage: 4D halo exchanges ...............         0.513  ( 0.1439 %)
  Message Passage: data broadcast ..................       156.046  (43.7875 %)
  Message Passage: data reduction ..................         0.148  ( 0.0416 %)
  Message Passage: data gathering ..................        18.610  ( 5.2220 %)
  Message Passage: data scattering..................        18.562  ( 5.2085 %)
  Message Passage: boundary data gathering .........         0.239  ( 0.0671 %)
  Message Passage: multi-model coupling ............        22.428  ( 6.2935 %)
                                              Total:       226.906   63.6714

 All percentages are with respect to total time =          356.370

 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_smflux.h
     ROMS/Functionals/ana_stflux.h

 ROMS/TOMS - Blows up ................ exit_flag:   1


 MAIN: Abnormal termination: BLOWUP.
+time 20120101.003000   , step      2; iteration    1; sweep 3 grid   1
+time 20120101.003000   , step      2; iteration    1; sweep 4 grid   1
 == SWAN grid  1 sent wave data to ROMS

jcwarner
Posts: 1171
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Model blowing up using elevation at the boundary

#2 Unread post by jcwarner »

i am not sure that using free surface from HYCOM is the best choice. it really depends on your domain. i would think adcirc/osu tides would be a better place to start.

for the error, it is difficult to tell. You could try to run the modeling system with just roms, and that would limit where to look. If it runs for a few time steps, then output the his file every time step and see where it blows up. you can try to reduce the physics options, etc. you will need to dig around and try a few things.
-j

paim_oceano
Posts: 10
Joined: Mon Jun 16, 2014 4:33 pm
Location: University of Delaware

Re: Model blowing up using elevation at the boundary

#3 Unread post by paim_oceano »

Hey there John, thanks for your answer.

I thought about using Hycom because I am trying to simulate the surge caused by hurricane Sandy on Delaware Bay. When I applied just the winds and tides, the simulation fails to achieve that (in fact, instead of having an elevation in the sea surface, it generates a depression).

Do you think adcirc could work in this case? And if not, would you suggest something else for me to try?

About running without SWAN, I will try that.

Thanks again

User avatar
wilkin
Posts: 872
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: Model blowing up using elevation at the boundary

#4 Unread post by wilkin »

The real-time assimilation system we run for MARACOOS has a reasonable realization of the surge associated with Sandy along the NJ/Delaware coast. The Version 2 output at http://tds.marine.rutgers.edu:8080/thre ... talog.html would be the version you want.

In this version we never reconciled sea level datum with NAVD88 so you will find you need to add about 40 cm to make the values consistent with Mean Sea Level referenced tide gauge data, but you can be the judge of the best approach there for your purposes.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

jcwarner
Posts: 1171
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Model blowing up using elevation at the boundary

#5 Unread post by jcwarner »

you can also use our COAWST forecast to get open bc's for roms. We run romw+swan coupled daily, and have output available on a thredds server
model_url='http://geoport.whoi.edu/thredds/dodsC/c ... e_best.ncd'

One way to access this is to download nctoolbox from here:
https://github.com/nctoolbox/nctoolbox

then you can run this m file attached.

We had rather good tides and waves compared to Sandy observations.

-john

nacholibre
Posts: 81
Joined: Thu Dec 07, 2006 3:14 pm
Location: USGS
Contact:

Re: Model blowing up using elevation at the boundary

#6 Unread post by nacholibre »

For a second test, I've created a file with information also in the North and South boundary, but in this case the model blows immediately during the first time step.
I would check the masking and bathymetry near the coastline at the southern and northern boundaries. You might have bad values in your forcing files as a result of extrapolation at the mismatched coastline boundaries between the two models or you might have unmasked cells with negative bathymetry near the coastline.
Zafer

paim_oceano
Posts: 10
Joined: Mon Jun 16, 2014 4:33 pm
Location: University of Delaware

Re: Model blowing up using elevation at the boundary

#7 Unread post by paim_oceano »

Jcwarner and Wilkin - Thanks a lot for the suggestion guys. I will look into both dataset to see which one can be more suitable to my case (and thanks for the heads up on the mean sea level Wilkin) :)

Nacholibre - When you talk about bad values, do you mean NaN (not a number) values?
At some point I noticed that the southern boundary in my grid was set to water between the points 100 - 150 while my bryfile just had values for the points between 103-150. However, even after adding values to these three extra points I am running into the same problem.

And also your comment about negative bathymetry brought me to another question: is it possible to simulate flood using ROMS? The aim of my project is to simulate flood caused by storm surge. If I cannot use negative values in the bathymetry, is there another way through which I could do it?

Thanks a lot

jcwarner
Posts: 1171
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Model blowing up using elevation at the boundary

#8 Unread post by jcwarner »

Warner, J.C., Haas, K.H ., Defne, Z., and Arango, H.G.(2013). A wetting and drying scheme for ROMS, Computers and Geosciences, 58, 54-61.

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

Re: Model blowing up using elevation at the boundary

#9 Unread post by kate »

You can have negative values for bathymetry, you just can't have zero values. In my experience, if the depth goes negative, ROMS will keep running - unless you try to grow ice there, in which case the ice_frazil routine will complain about negative ice growth.

Or as John says, you can turn on the WET_DRY option.

paim_oceano
Posts: 10
Joined: Mon Jun 16, 2014 4:33 pm
Location: University of Delaware

Re: Model blowing up using elevation at the boundary

#10 Unread post by paim_oceano »

Kate and John: WET_DRY is definitely the path I am thinking to follow. I was just worried about it not being possible after Zafer comment.
Thanks once more.

poppai
Posts: 9
Joined: Wed Sep 30, 2009 6:05 pm
Location: IMHEN
Contact:

Re: Model blowing up using elevation at the boundary

#11 Unread post by poppai »

jcwarner wrote:you can also use our COAWST forecast to get open bc's for roms. We run romw+swan coupled daily, and have output available on a thredds server
model_url='http://geoport.whoi.edu/thredds/dodsC/c ... e_best.ncd'

One way to access this is to download nctoolbox from here:
https://github.com/nctoolbox/nctoolbox

then you can run this m file attached.

We had rather good tides and waves compared to Sandy observations.

-john
Dear John.
How to get BC condition from your COAWST forecast? Have you a m file for this work? I using roms_master_climatology_coawst_mw, but seem SST and zeta have problem (ex: SST obtain 35oC). I run for South china sea domain.

jcwarner
Posts: 1171
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Model blowing up using elevation at the boundary

#12 Unread post by jcwarner »

our coawst forecast is for the US east coast. so I suggest you use the HYCOMM data as a start, and set the correct hycomm url in roms_master_climatology_coawst_mw.
-j

hbzong
Posts: 36
Joined: Thu Oct 04, 2007 4:14 am
Location: Fathom Science
Contact:

Re: Model blowing up using elevation at the boundary

#13 Unread post by hbzong »

I think that this error was caused by NaN values in South/North boundary.

Post Reply