updated test problem

[pooran]

Base on your suggestion to use new version of COAWST, I have updated my coupled test and then each model (ROMS and SWAN) ran ,separately without any problem . But their coupling faced with blow up very very soon, in the first step and therefore all data in rst and his files has their initial values, yet. Roms, swan and coupled's time steps (90s,900s,900s,respectively) and almost all options in header file were same with my coupled test in COAWST old version (except lateral wall condition which defined in ocean.in as 'clo' for all LBC flags ). It makes me confused why same coupled test has ran in COAWST's old version correctly but in new version faced with blow up
#undef ROLLER_RENIERS and decreasing time step can not solve the problem.
Would you please help me . I really need you help to solve it
I have send my set ups. I ashamed to this request but I become so appreciate if you look at it.
Many thanks in advance for any reply
cheers
pooran

Code: Select all

pooran@fk1:~/program_install/COAWST/my_test_Project/Rip_current1$ mpirun -np 4 ./coawstG coupling_rip_current.in 
 Coupled Input File name = coupling_rip_current.in 

 Model Coupling: 


       Ocean Model MPI nodes: 000 - 001 

       Waves Model MPI nodes: 002 - 003 

 Model Input Parameters:  ROMS/TOMS version 3.7  
                          Tuesday - September 23, 2014 -  2:39:04 PM 
 ----------------------------------------------------------------------------- 

SWAN grid   1 is preparing computation 


 Rip Current 

 Operating system : Linux 
 CPU/hardware     : i686 
 Compiler system  : gfortran 
 Compiler command : /home/pooran/program_install/mpich-3.0.4/bin/mpif90 
 Compiler flags   : -frepack-arrays -g -fbounds-check -I/home/pooran/program_install/mct/include -I/usr/include -ffree-form -ffree-line-length-none 

 Input Script  : 

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

 Local Root    : /home/pooran/program_install/COAWST 
 Header Dir    : /home/pooran/program_install/COAWST/my_test_Project/Rip_current1 
 Header file   : rip_current.h 
 Analytical Dir: /home/pooran/program_install/COAWST/my_test_Project/Rip_current1 

 Resolution, Grid 01: 0194x0344x032,  Parallel Nodes:   2,  Tiling: 001x002 


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

      57600  ntimes            Number of timesteps for 3-D equations. 
     90.000  dt                Timestep size (s) for 3-D equations. 
         30  ndtfast           Number of timesteps for 2-D equations between 
                                 each 3D timestep. 
          1  ERstr             Starting ensemble/perturbation run number. 
          1  ERend             Ending ensemble/perturbation run number. 
          0  nrrec             Number of restart records to read from disk. 
          T  LcycleRST         Switch to recycle time-records in restart file. 
          1  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). 
          1  nHIS              Number of timesteps between the writing fields 
                                 into history file. 
          1  ntsAVG            Starting timestep for the accumulation of output 
                                 time-averaged data. 
       9600  nAVG              Number of timesteps between the writing of 
                                 time-averaged data into averages file. 
          1  ntsDIA            Starting timestep for the accumulation of output 
                                 time-averaged diagnostics data. 
      28800  nDIA              Number of timesteps between the writing of 
                                 time-averaged data into diagnostics file. 
 1.0000E+02  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient 
                                 (m2/s) for tracer 01: temp 
 1.0000E+02  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient 
                                 (m2/s) for tracer 02: salt 
 1.0000E+01  nl_visc2          NLM Horizontal, harmonic mixing coefficient 
                                 (m2/s) for momentum. 
          F  LuvSponge         Turning OFF sponge on horizontal momentum. 
          F  LtracerSponge(01) Turning OFF sponge on tracer 01: temp 
          F  LtracerSponge(02) Turning OFF sponge on tracer 02: salt 
 5.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s) 
                                 for tracer 01: temp 
 5.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.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. 
   1400.000  charnok_alpha     Charnok factor for Zos calculation. 
      0.500  zos_hsig_alpha    Factor for Zos calculation using Hsig(Awave). 
      0.020  sz_alpha          Factor for Wave dissipation surface tke flux . 
    100.000  crgban_cw         Factor for Craig/Banner surface tke flux. 
      0.500  wec_alpha         WEC factor for roller/breaking energy distribution. 
 2.0000E-04  rdrg              Linear bottom drag coefficient (m/s). 
 3.3000E-03  rdrg2             Quadratic bottom drag coefficient. 
 4.0000E-03  Zob               Bottom roughness (m). 
 4.0000E-03  Zos               Surface roughness (m). 
          1  lmd_Jwt           Jerlov water type. 
          2  Vtransform        S-coordinate transformation equation. 
          4  Vstretching       S-coordinate stretching function. 
 8.0000E+00  theta_s           S-coordinate surface control parameter. 
 4.0000E-01  theta_b           S-coordinate bottom  control parameter. 
      6.000  Tcline            S-coordinate surface/bottom layer width (m) used 
                                 in vertical coordinate stretching. 
   1025.000  rho0              Mean density (kg/m3) for Boussinesq approximation. 
    304.000  dstart            Time-stamp assigned to model initialization (days). 
20060101.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. 
          F  VolCons(1)        NLM western  edge boundary volume conservation. 
          F  VolCons(2)        NLM southern edge boundary volume conservation. 
          F  VolCons(3)        NLM eastern  edge boundary volume conservation. 
          F  VolCons(4)        NLM northern edge boundary volume conservation. 
     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  LuvSrc            Turning ON  momentum point Sources/Sinks. 
          F  LwSrc             Turning OFF volume influx point Sources/Sinks. 
          T  LtracerSrc(01)    Turning ON  point Sources/Sinks on tracer 01: temp 
          T  LtracerSrc(02)    Turning ON  point Sources/Sinks on tracer 02: salt 
          F  LsshCLM           Turning OFF processing of SSH climatology. 
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology. 
          F  Lm3CLM            Turning OFF processing of 3D momentum climatology. 
          F  LtracerCLM(01)    Turning OFF processing of climatology tracer 01: temp 
          F  LtracerCLM(02)    Turning OFF processing of climatology tracer 02: salt 
          F  LnudgeM2CLM       Turning OFF nudging of 2D momentum climatology. 
          F  LnudgeM3CLM       Turning OFF nudging of 3D momentum climatology. 
          F  LnudgeTCLM(01)    Turning OFF nudging of climatology tracer 01: temp 
          F  LnudgeTCLM(02)    Turning OFF nudging of climatology tracer 02: salt 
          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(idUsms)      Write out surface U-momentum stress. 
          T  Hout(idVsms)      Write out surface V-momentum stress. 
          T  Hout(idU2Sd)      Write out 2D u-momentum stokes velocity. 
          T  Hout(idV2Sd)      Write out 2D v-momentum stokes velocity. 
          T  Hout(idU3Sd)      Write out 3D u-momentum stokes velocity. 
          T  Hout(idV3Sd)      Write out 3D v-momentum stokes velocity. 
          T  Hout(idW3Sd)      Write out 3D omega-momentum stokes velocity. 
          T  Hout(idW3St)      Write out 3D w-momentum stokes velocity. 
          T  Hout(idWamp)      Write out wave height. 
          T  Hout(idWlen)      Write out wavelength. 
          T  Hout(idWdir)      Write out wave direction. 
          T  Hout(idWptp)      Write out wave surface period. 
          T  Hout(idWdif)      Write out wave dissipation due to bottom friction. 
          T  Hout(idUwav)      Wave-avg surface u-velocity. 
          T  Hout(idVwav)      Wave-avg surface v-velocity. 
          T  Hout(idTsur)      Write out surface net heat flux. 
          T  Hout(idTsur)      Write out surface net salt flux. 
          T  Hout(idSrad)      Write out shortwave radiation flux. 

          T  Aout(idUbar)      Write out averaged 2D U-momentum component. 
          T  Aout(idVbar)      Write out averaged 2D V-momentum component. 
          T  Aout(idUvel)      Write out averaged 3D U-momentum component. 
          T  Aout(idVvel)      Write out averaged 3D V-momentum component. 
          T  Aout(idWvel)      Write out averaged W-momentum component. 
          T  Aout(idTvar)      Write out averaged tracer 01: temp 
          T  Aout(idTvar)      Write out averaged tracer 02: salt 
          T  Aout(idU2Sd)      Write out averaged 2D u-momentum stokes velocity. 
          T  Aout(idV2Sd)      Write out averaged 2D v-momentum stokes velocity. 
          T  Aout(idTsur)      Write out averaged surface net heat flux. 
          T  Aout(idTsur)      Write out averaged surface net salt flux. 
          T  Aout(idSrad)      Write out averaged shortwave radiation flux. 
          T  Aout(idSdif)      Write out averaged vertical diffusion: AKt(isalt). 




 Output/Input Files: 
 
             Output Restart File:  ocean_rip_current_rst.nc 
             Output History File:  ocean_rip_current_his.nc 
            Output Averages File:  ocean_rip_current_avg.nc 
         Output Diagnostics File:  ocean_rip_current_dia.nc 
        Physical parameters File:  
                 Input Grid File:  /home/pooran/my_test_data/grid/test_grd_3.nc 
    Input Nonlinear Initial File:  /home/pooran/my_test_data/clim-init/CS_initial_Coh_01nov06.nc 
        Input Sources/Sinks File:  /home/pooran/my_test_data/river/rivers_force_06.nc 
           Input Forcing File 01:  /home/pooran/my_test_data/Forcing/CS_shflux2_era_06.nc 
           Input Forcing File 02:  /home/pooran/my_test_data/Forcing/wind2stress/CS_wind2stress_06.nc 
           Input Forcing File 03:  /home/pooran/my_test_data/Forcing/CS_swflux2_era_06.nc 
           Input Forcing File 04:  /home/pooran/my_test_data/Forcing/CS_swrad2_era_06.nc 

 Tile partition information for Grid 01:  0194x0344x0032  tiling: 001x002 

     tile     Istr     Iend     Jstr     Jend     Npts 

 Number of tracers:            2 
        0        1      194        1      172  1067776 
        1        1      194      173      344  1067776 

 Tile minimum and maximum fractional coordinates for Grid 01: 
   (interior points only) 

     tile     Xmin     Xmax     Ymin     Ymax     grid 

        0     0.50   194.50     0.50   172.50  RHO-points 
        1     0.50   194.50   172.50   344.50  RHO-points 

        0     1.00   194.00     0.50   172.50    U-points 
        1     1.00   194.00   172.50   344.50    U-points 

        0     0.50   194.50     1.00   172.50    V-points 
        1     0.50   194.50   172.50   344.00    V-points 

 Maximum halo size in XI and ETA directions: 

               HaloSizeI(1) =     410 
               HaloSizeJ(1) =     368 
                TileSide(1) =     199 
                TileSize(1) =   35422 


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

 Variable               Grid  West Edge    South Edge   East Edge    North Edge 
 ---------              ----  ----------   ----------   ----------   ---------- 
 
 zeta                     1   Closed       Closed       Closed       Closed 

 ubar                     1   Closed       Closed       Closed       Closed 

 vbar                     1   Closed       Closed       Closed       Closed 

 u                        1   Closed       Closed       Closed       Closed 

 v                        1   Closed       Closed       Closed       Closed 

 temp                     1   Closed       Closed       Closed       Closed 

 salt                     1   Closed       Closed       Closed       Closed 

 tke                      1   Closed       Closed       Closed       Closed 

 ubar_stokes              1   Closed       Closed       Closed       Closed 

 vbar_stokes              1   Closed       Closed       Closed       Closed 

 u_stokes                 1   Closed       Closed       Closed       Closed 

 v_stokes                 1   Closed       Closed       Closed       Closed 

 Activated C-preprocessing Options: 

 RIP_CURRENT         Rip Current from F.Komijani 
 ANA_BSFLUX          Analytical kinematic bottom salinity flux. 
 ANA_BTFLUX          Analytical kinematic bottom temperature flux. 
 ASSUMED_SHAPE       Using assumed-shape arrays. 
 AVERAGES            Writing out time-averaged nonlinear model fields. 
 DIAGNOSTICS_TS      Computing and writing tracer diagnostic terms. 
 DIAGNOSTICS_UV      Computing and writing momentum diagnostic terms. 
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002). 
 DOUBLE_PRECISION    Double precision arithmetic. 
 GLS_MIXING          Generic Length-Scale turbulence closure. 
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation. 
 MASKING             Land/Sea masking. 
 MCT_LIB             Using Model Coupling Toolkit library. 
 MIX_S_TS            Mixing of tracers along constant S-surfaces. 
 MIX_S_UV            Mixing of momentum along constant S-surfaces. 
 MPI                 MPI distributed-memory configuration. 
 WEC_VF              Vortex Force wave current interaction- UNDER DEVELOPMENT!!!!. 
 WDISS_WAVEMOD       Wave energy dissipation acquired from coupled wave model. 
 ROLLER_RENIERS      Wave energy roller based on Reniers 2004. 
 NONLINEAR           Nonlinear Model. 
 NONLIN_EOS          Nonlinear Equation of State for seawater. 
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear. 
 OUT_DOUBLE          Double precision output fields in NetCDF files. 
 PERFECT_RESTART     Processing perfect restart variables. 
 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. 
 SPLINES             Conservative parabolic spline reconstruction. 
 SWAN_COUPLING       SWAN model coupling. 
 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_KIRBY            Compute uwave and vwave Kirby avg velocities. 
 UV_QDRAG            Quadratic bottom stress. 
 UV_VIS2             Harmonic mixing of momentum. 
 VAR_RHO_2D          Variable density barotropic mode. 
 WAVES_OCEAN         Two-way wave-ocean models coupling. 

 Process Information: 

 Node #  0 (pid=    7379) is active. 
 Node #  1 (pid=    7380) is active. 

 INITIAL: Configuring and initializing forward nonlinear model ... 
 ******* 

 Vertical S-coordinate System, Grid 01: 

 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.0000256       -0.101      -0.094      -0.198      -0.211 
    30  -0.0625000  -0.0001040       -0.202      -0.188      -0.421      -0.474 
    29  -0.0937500  -0.0002400       -0.304      -0.282      -0.672      -0.792 
    28  -0.1250000  -0.0004423       -0.406      -0.376      -0.956      -1.175 
    27  -0.1562500  -0.0007236       -0.508      -0.471      -1.278      -1.635 
    26  -0.1875000  -0.0011014       -0.610      -0.566      -1.647      -2.188 
    25  -0.2187500  -0.0015996       -0.713      -0.661      -2.074      -2.858 
    24  -0.2500000  -0.0022492       -0.816      -0.757      -2.575      -3.675 
    23  -0.2812500  -0.0030910       -0.920      -0.853      -3.170      -4.679 
    22  -0.3125000  -0.0041778       -1.025      -0.950      -3.884      -5.921 
    21  -0.3437500  -0.0055774       -1.132      -1.048      -4.749      -7.467 
    20  -0.3750000  -0.0073773       -1.239      -1.147      -5.810      -9.402 
    19  -0.4062500  -0.0096898       -1.349      -1.248      -7.119     -11.834 
    18  -0.4375000  -0.0126588       -1.461      -1.350      -8.748     -14.904 
    17  -0.4687500  -0.0164685       -1.576      -1.456     -10.786     -18.791 
    16  -0.5000000  -0.0213546       -1.696      -1.564     -13.346     -23.723 
    15  -0.5312500  -0.0276179       -1.820      -1.677     -16.577     -29.993 
    14  -0.5625000  -0.0356419       -1.952      -1.794     -20.663     -37.973 
    13  -0.5937500  -0.0459141       -2.091      -1.919     -25.843     -48.136 
    12  -0.6250000  -0.0590529       -2.242      -2.052     -32.416     -61.084 
    11  -0.6562500  -0.0758390       -2.406      -2.196     -40.763     -77.575 
    10  -0.6875000  -0.0972542       -2.588      -2.354     -51.360     -98.561 
     9  -0.7187500  -0.1245249       -2.791      -2.530     -64.804    -125.235 
     8  -0.7500000  -0.1591703       -3.023      -2.728     -81.834    -159.072 
     7  -0.7812500  -0.2030516       -3.289      -2.953    -103.354    -201.879 
     6  -0.8125000  -0.2584153       -3.599      -3.213    -130.457    -255.839 
     5  -0.8437500  -0.3279186       -3.962      -3.515    -164.434    -323.533 
     4  -0.8750000  -0.4146160       -4.390      -3.869    -206.771    -407.927 
     3  -0.9062500  -0.5218760       -4.895      -4.284    -259.106    -512.293 
     2  -0.9375000  -0.6531834       -5.491      -4.772    -323.132    -640.016 
     1  -0.9687500  -0.8117684       -6.190      -5.342    -400.420    -794.233 
     0  -1.0000000  -1.0000000       -7.000      -6.000    -492.123    -977.246 

 Time Splitting Weights for Grid 01:    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.2343461315221332 
  8-0.0005032751608632 0.0336903762267453-0.0112145619632232 0.2680365077488785 
  9 0.0003887272597151 0.0337071520654408-0.0108258347035082 0.3017436598143193 
 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.4026337758325832 
 13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801077 
 14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905599 
 15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029189 
 16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429149 
 17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276979 
 18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780106 
 19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489551 
 20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976139 
 21 0.0256737146312911 0.0297153720153352 0.1342125541712343 0.6900486906129492 
 22 0.0286498597812017 0.0288595815276255 0.1628624139524359 0.7189082721405747 
 23 0.0316613792205220 0.0279045862015855 0.1945237931729579 0.7468128583421602 
 24 0.0346736416507075 0.0268492068942347 0.2291974348236654 0.7736620652363949 
 25 0.0376471948657328 0.0256934188392112 0.2668446296893983 0.7993554840756061 
 26 0.0405373376992233 0.0244385123436867 0.3073819673886216 0.8237939964192929 
 27 0.0432936737565711 0.0230872677537126 0.3506756411451927 0.8468812641730055 
 28 0.0458596469320356 0.0216441452951603 0.3965352880772283 0.8685254094681657 
 29 0.0481720587108285 0.0201154903974257 0.4447073467880568 0.8886408998655915 
 30 0.0501605672561820 0.0185097551070648 0.4948679140442388 0.9071506549726562 
 31 0.0517471682814031 0.0168377361985254 0.5466150823256418 0.9239883911711816 
 32 0.0528456577069106 0.0151128305891453 0.5994607400325525 0.9391012217603268 
 33 0.0533610761022577 0.0133513086655816 0.6528218161348102 0.9524525304259084 
 34 0.0531891349131380 0.0115726061288397 0.7060109510479481 0.9640251365547481 
 35 0.0522156244733761 0.0097996349650684 0.7582265755213241 0.9738247715198165 
 36 0.0503158038019030 0.0080591141492892 0.8085423793232271 0.9818838856691057 
 37 0.0473537721847153 0.0063819206892258 0.8558961515079425 0.9882658063583314 
 38 0.0431818225418188 0.0048034616164019 0.8990779740497613 0.9930692679747333 
 39 0.0376397765791563 0.0033640675316746 0.9367177506289176 0.9964333355064079 
 40 0.0305543017255205 0.0021094083123694 0.9672720523544381 0.9985427438187773 
 41 0.0217382098544503 0.0010909315881854 0.9890102622088884 0.9996336754069627 
 42 0.0109897377911117 0.0003663245930371 1.0000000000000000 0.9999999999999998 

 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 

 Metrics information for Grid 01: 
 =============================== 

 Minimum X-grid spacing, DXmin =  3.33514538E+00 km 
 Maximum X-grid spacing, DXmax =  3.98831267E+00 km 
 Minimum Y-grid spacing, DYmin =  3.33800882E+00 km 
 Maximum Y-grid spacing, DYmax =  3.98739982E+00 km 
 Minimum Z-grid spacing, DZmin =  1.01057987E-01 m 
 Maximum Z-grid spacing, DZmax =  1.83012378E+02 m 

 Minimum barotropic Courant Number =  8.90418085E-03 
 Maximum barotropic Courant Number =  1.08159505E-01 
 Maximum Coriolis   Courant Number =  9.58513554E-03 


 NLM: GET_STATE - Read state initial conditions,             t =   304 00:00:00 
                   (Grid 01, File: CS_initial_Coh_01nov06.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.24582748E-06 Max =  2.29034667E+01) 
                - salinity 
                   (Min =  0.00000000E+00 Max =  1.29654990E+01) 
    GET_NGFLD   - river runoff  XI-positions at RHO-points 
                   (Min =  1.80000000E+01 Max =  1.11000000E+02) 
    GET_NGFLD   - river runoff ETA-positions at RHO-points 
                   (Min =  8.50000000E+01 Max =  3.31000000E+02) 
    GET_NGFLD   - river runoff direction 
                   (Min =  0.00000000E+00 Max =  1.00000000E+00) 
    GET_NGFLD   - river runoff mass transport vertical profile 
                   (Min =  0.00000000E+00 Max =  7.63425333E-02) 
    GET_NGFLD   - river runoff mass transport,               t =   289 03:45:00 
                   (Rec=0000010, Index=2, File: rivers_force_06.nc) 
                   (Tmin=         15.2188 Tmax=        350.0313) 
                   (Min =  1.95000000E+01 Max =  1.34300000E+02) 
    GET_NGFLD   - river runoff potential temperature,        t =   289 03:45:00 
                   (Rec=0000010, Index=2, File: rivers_force_06.nc) 
                   (Tmin=         15.2188 Tmax=        350.0313) 
                   (Min =  2.00000000E+00 Max =  1.78000000E+01) 
    GET_NGFLD   - river runoff salinity,                     t =   289 03:45:00 
                   (Rec=0000010, Index=2, File: rivers_force_06.nc) 
                   (Tmin=         15.2188 Tmax=        350.0313) 
                   (Min =  0.00000000E+00 Max =  0.00000000E+00) 
    GET_2DFLD   - surface u-momentum stress,                 t =   304 00:00:00 
                   (Rec=0002432, Index=1, File: CS_wind2stress_06.nc) 
                   (Tmin=          0.1250 Tmax=        365.0000) 
                   (Min = -8.32389989E-06 Max =  7.00611541E-06) 
    GET_2DFLD   - surface v-momentum stress,                 t =   304 00:00:00 
                   (Rec=0002432, Index=1, File: CS_wind2stress_06.nc) 
                   (Tmin=          0.1250 Tmax=        365.0000) 
                   (Min = -2.51533043E-06 Max =  9.34370765E-06) 
    GET_2DFLD   - solar shortwave radiation flux,            t =   303 22:30:00 
                   (Rec=0002432, Index=1, File: CS_swrad2_era_06.nc) 
                   (Tmin=          0.0625 Tmax=        364.9375) 
                   (Min =  0.00000000E+00 Max =  0.00000000E+00) 
    GET_2DFLD   - surface net heat flux,                     t =   303 22:30:00 
                   (Rec=0002432, Index=1, File: CS_shflux2_era_06.nc) 
                   (Tmin=          0.0625 Tmax=        364.9375) 
                   (Min = -3.27425753E-05 Max = -1.92829103E-06) 
    GET_2DFLD   - surface net freswater flux, (E-P),         t =   303 22:30:00 
                   (Rec=0002432, Index=1, File: CS_swflux2_era_06.nc) 
                   (Tmin=          0.0625 Tmax=        364.9375) 
                   (Min = -1.48378471E-07 Max =  2.19167402E-08) 

 Basin information for Grid 01: 

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

 Initial basin volumes: TotVolume =  7.6973949537E+13 m3 
                        MinVolume =  1.1438083862E+06 m3 
                        MaxVolume =  2.6995814696E+09 m3 
                          Max/Min =  2.3601693274E+03 
 == SWAN grid  1 sent wave data to ROMS grid  1 
 ** ROMS grid  1 recv data from SWAN grid  1 


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



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

    GET_NGFLD   - river runoff mass transport,               t =   319 14:15:00 
                   (Rec=0000011, Index=1, File: rivers_force_06.nc) 
                   (Tmin=         15.2188 Tmax=        350.0313) 
                   (Min =  1.70000000E+01 Max =  1.34300000E+02) 
    GET_NGFLD   - river runoff potential temperature,        t =   319 14:15:00 
                   (Rec=0000011, Index=1, File: rivers_force_06.nc) 
                   (Tmin=         15.2188 Tmax=        350.0313) 
                   (Min =  0.00000000E+00 Max =  9.00000000E+00) 
    GET_NGFLD   - river runoff salinity,                     t =   319 14:15:00 
                   (Rec=0000011, Index=1, File: rivers_force_06.nc) 
                   (Tmin=         15.2188 Tmax=        350.0313) 
                   (Min =  0.00000000E+00 Max =  0.00000000E+00) 
    GET_2DFLD   - surface u-momentum stress,                 t =   304 03:00:00 
                   (Rec=0002433, Index=2, File: CS_wind2stress_06.nc) 
                   (Tmin=          0.1250 Tmax=        365.0000) 
                   (Min = -8.18458388E-06 Max =  6.08474362E-06) 
    GET_2DFLD   - surface v-momentum stress,                 t =   304 03:00:00 
                   (Rec=0002433, Index=2, File: CS_wind2stress_06.nc) 
                   (Tmin=          0.1250 Tmax=        365.0000) 
                   (Min = -4.47947393E-06 Max =  7.40512407E-06) 
    GET_2DFLD   - solar shortwave radiation flux,            t =   304 01:30:00 
                   (Rec=0002433, Index=2, File: CS_swrad2_era_06.nc) 
                   (Tmin=          0.0625 Tmax=        364.9375) 
                   (Min =  0.00000000E+00 Max =  1.63846996E-08) 
    GET_2DFLD   - surface net heat flux,                     t =   304 01:30:00 
                   (Rec=0002433, Index=2, File: CS_shflux2_era_06.nc) 
                   (Tmin=          0.0625 Tmax=        364.9375) 
                   (Min = -3.82565374E-05 Max = -2.69276225E-06) 
+time 20061101.001500   , step      1; iteration    1; sweep 1 grid   1 
    GET_2DFLD   - surface net freswater flux, (E-P),         t =   304 01:30:00 
                   (Rec=0002433, Index=2, File: CS_swflux2_era_06.nc) 
                   (Tmin=          0.0625 Tmax=        364.9375) 
                   (Min = -1.21615651E-07 Max =  3.29965440E-08) 

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

      0   304 00:00:00  0.000000E+00  2.675339E+03  2.675339E+03  8.067723E+13 
          (000,000,00)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00 
      DEF_HIS   - creating history file, Grid 01: ocean_rip_current_his.nc 
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001 
      DEF_AVG   - creating average file, Grid 01: ocean_rip_current_avg.nc 
      DEF_DIAGS - creating diagnostics file, Grid 01: ocean_rip_current_dia.nc 
+time 20061101.001500   , step      1; iteration    1; sweep 2 grid   1 
+time 20061101.001500   , step      1; iteration    1; sweep 3 grid   1 
+time 20061101.001500   , step      1; iteration    1; sweep 4 grid   1 
 1   304 00:01:30           NaN           NaN           NaN           NaN 
          (145,323,26)  2.828234E-06  2.260028E-06  3.989284E-03  4.544180E-03 

 Blowing-up: Saving latest model state into  RESTART file 

(145,323)is boundary cell.
mask_h(145,323)=1
mask_u(145,323)=1
mask_v(145,323)=0

[jcwarner]

We have updated the rip_current test case and i just ran it and it works fine. If you want to do an svn update you can get the applications with a working set up. We are about to announce that all the test cases have been updated to the new roms input files etc. sorry for the delay.

I looked through your input files and have a few suggestions. for the ocean.in
- this must be a different grids than we distribute
- 32 layers for a shallow water rip current is a lot, maybe 16 is a good number?
- we have eastern edge as gradient, not all 4 closed.
- 90 sec time step might be too high. check your courant number.
- visc2 and 4 is your input file are too high. see what we used as a suggestion.
- the thetas b etc params are stretched rather tight in your setup. might want to relax them a bit.

for swan input you have
CIRcle 24 0.045 .65
suggest you include the tail
CIRcle 24 0.045 1.0

[pooran]

thanks for rapid reply.
My domain in Caspian sea (closed basin)which i want to do nesting (4 times) in coupled model to investigate rip current near the coast.
1-

this must be a different grids than we distribute

yes, i used
Lm == 194 ! Number of I-direction INTERIOR RHO-points
Mm == 344 ! Number of J-direction INTERIOR RHO-points

2-we have eastern edge as gradient, not all 4 closed.

As i told My domain is close basin which all boundary are closed. I think you choosed gradient for east edge because you want to inter waves from this edge to your domain. But in my test i have not define that wave come from all edge, instead of it i define wind force in swan model for whole domain Therefore i guess that i should choose closed for Wec boundary conditions. Is it correct?
3- also by choosing

Code: Select all

LBC(isU2Sd) ==   Gra     Gra     Gra     Gra         ! 2D U-stokes
  LBC(isV2Sd) ==   Gra     Gra     Gra     Gra         ! 2D V-stokes
   LBC(isU3Sd) ==   Gra     Gra     Gra     Gra         ! 3D U-stokes
   LBC(isV3Sd) ==   Gra     Gra     Gra     Gra         ! 3D V-stokes

problem does not solved.
3-

32 layers for a shallow water rip current is a lot, maybe 16 is a good number

the thetas b etc params are stretched rather tight in your setup. might want to relax them a bit
but my test in COAWST old version works with this kind of set up?
4-

90 sec time step might be too high. check your courant number.

I have changes 90s to 30s but it cant help to solve blow up. Courant numbers for 90s are:

Code: Select all

Minimum barotropic Courant Number =  8.90418085E-03 
 Maximum barotropic Courant Number =  1.08159505E-01 
 Maximum Coriolis   Courant Number =  9.58513554E-03 

5-

visc2 and 4 is your input file are too high. see what we used as a suggestion.

I used what you suggest in rip_current test. But cant solve problem.
It is strange for me that by all mentioned changes, blow up happen in same step (step 1)
please help me with your experience

[jcwarner]

well you need to dig in a little here. there is no magic fix.
is it blowing up on the boundary?
is it blowing up at a wet/dry point?
is it blowing up in a region of strong bathy gradients? ....
you need to dig into the rst file, maybe turn off the roller and other options and just get it to work with a general set up first. then start to add other options.

[pooran]

Code: Select all

1   304 00:01:30           NaN           NaN           NaN           NaN 
          (145,323,26)  2.828234E-06  2.260028E-06  3.989284E-03  4.544180E-03 

 Blowing-up: Saving latest model state into  RESTART file 

as i told there are just u and v for initial time step (equal to zero) in rst file (because model has not been solved for next steps). if base on model result, blow up happened on(145,323,26), masking in mask_h(145,323)=1 and for points near it are:
mask_h(145,323)=1,
mask_h(144,323)=1,
mask_h(146,323)=1,
mask_h(147,323)=0;
mask_h(145,324)=0
mask_h(145,322)=1
i don't know this point is on the boundary or is wet/dry point .
this point's depth are 7m which is equal to other points near it.
as i told, turning off the roller cant solved problem

[kate]

Actually, you don't know that it blew up at that point. That point is returned from the MPI reduction to find the largest Courant number, which is probably ignoring the NaNs. You need to look at your restart file with ncview or some other viewer (panoply, perhaps, haven't tried it yet).

[pooran]

Dear Kate and Warner
coupled model has been ran for different scenarios:
a) decreasing dt of ROMS from 90s to 1s, b) adding WEC_VF, WDISS_WAVEMOD,ROLLER_RENIERS,UV_KIRBY to header file step by step,c) changing of Wec boundary conditions types (Clo to Gra), d) using ana_initai instead of created initial file. But model faced with blow up exactly in first step for all scenarios,.therefore there is just one time step in restart file which has zero values for u and v. In the restart file i saw some strange things:
1- rho ,rubar, rvbar, rzeta, u_stokes,v_stokes, ubar, vbar have NAN in most parts of domain (as you can sea in attached fig)
2- AKS,AKT,AKV and gls have fix values in land cells as well as sea cells, 5e-6,5e-6,1.2e-5 and 1e-12, respectively.
To be honest, i can not find what exactly cause blow up. i hope to receive your help to solve model.
I become so appreciate if you guide me to run model successfully.

[jcwarner]

how about trying to run it as just roms by itself, with no wave forcing.

[pooran]

roms itself and swan itselp have ran without problem

[jcwarner]

if you undefine all these:
#undef WEC_VF
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

then there should be no affect of waves to roms. does that work?

[pooran]

exactly
it works with out any swan output

[jcwarner]

ok. so if you couple to swan, and
#undef WEC_VF
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

what do the Hwave, Dwave, Lwave variables look like?

[pooran]

thay are in attached fig

[jcwarner]

what if you try
#define WEC_VF
#define WDISS_THORGUZA
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

[pooran]

I have tried
#define WEC_VF
#define WDISS_THORGUZA
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

and

#define WEC_VF
#define WDISS_CHURTHOR
#undef WDISS_THORGUZA
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY

but same blow up appeared.

[pooran]

Dear
I have tried another states which in all of them GRA chosen for all lateral WEC boundary condition (just defined options have been written and others are undef)
a) #define WEC_VF and #define ROLLER_RENIERS
b) #define WEC_VF and #define ROLLER_SVENDSEN
c)#define WEC_VF , #define WAVE_MIXING
d)#define WEC_VF , #define WAVE_MIXING and # define WDISS_WAVEMOD
f)#define WEC_VF , #define WAVE_MIXING , # define WDISS_WAVEMOD and #define ROLLER_RENIERS

for a and b states, blow up appears in the first step and in the same position as seventh post . But in remaining states ( c, d and f), model faces with error in the first step

Code: Select all

+time 20061101.001500   , step      1; iteration    1; sweep 1 grid   1
   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
          C => (i,j,k)       Cu            Cv            Cw         Max Speed
      0   304 00:00:00  0.000000E+00  2.675339E+03  2.675339E+03  8.067723E+13
          (000,000,00)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
      DEF_HIS   - creating history file, Grid 01: ocean_rip_current_his.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
+time 20061101.001500   , step      1; iteration    1; sweep 2 grid   1
+time 20061101.001500   , step      1; iteration    1; sweep 3 grid   1
+time 20061101.001500   , step      1; iteration    1; sweep 4 grid   1
At line 202 of file wec_wave_mix.f90
Fortran runtime error: Index '33' of dimension 3 of array 'z_r' above upper bound of 32
At line 202 of file wec_wave_mix.f90
Fortran runtime error: Index '33' of dimension 3 of array 'z_r' above upper bound of 32
=========================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   EXIT CODE: 2
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES

in wec_wave_mix.f90 we have:

Code: Select all

! 
! Compute contribution of wave breaking induced mixing to Akv
!
        fac1=0.5_r8*SQRT(2.0_r8)
line 202        fac2=(z_r(i,j,k)-z_w(i,j,0))*oDstp(i,j)
        DO j=Jstr,Jend
          DO i=Istr,Iend
            cff2=fac2*gamr(i,j)
            cff3=COSH(2.0_r8*pi*cff2)
            cff1=Dissip_wcap(i,j)+Dissip_break(i,j)
            cff4=((1.0_r8-wec_alpha(ng))*cff1)**(1.0_r8/3.0_r8)
            cff4=cff4+(Dissip_roller(i,j))**(1.0_r8/3.0_r8)
            Akv(i,j,k)=Akv(i,j,k)+                                      &
     &                 Cb*cff4*fac1*Hwave(i,j)*                         &
     &                 Dstp(i,j)*cff3*oroller(i,j)
          END DO
        END DO

what is the meaning of this error? i know that z_r is Actual depths of variables at ρ-points. i have 32 layers in grid file.
i become grateful to receive your help

[jcwarner]

so that was a bug in wec_wave_mix and i just corrected it. It had to deal with a k loop issue. thanks.
But in the bigger picutre here are some suggestions:
1) - this is a coawst issue, not roms. so lets move this to
https://coawstmodel-trac.sourcerepo.com ... el_COAWST/
and post a ticket there.
2) undef all the wec stuff and try
#define SSW_BBL
#define SSSW_CALC_ZNOT
and see if that works. If it does, then it is probably the wec stuff. if it does not, then it is the way the wave stuff is coming into roms.
-john

[pooran]

Dear warner
sorry for sending post in roms forum again. I prefer to continue this subject here but for another one I will use COAWST forum.
I used
#undef WEC_VF
#undef WDISS_WAVEMOD
#undef ROLLER_RENIERS
#undef UV_KIRBY
#define SSW_BBL
#define SSSW_CALC_ZNOT
#undef UV_QDRAG
model could work and past step 1.
But what is the meaning of it? Do those 2 suggested options works as wave effect on current options? What should I do now?

[jcwarner]

i just sent you an email with a new version of the wec_stokes.F. We made some changes to the way the stokes velocities were computed. can you try this new version
-john

[pooran]

WOW, model is running with new wec_stokes.F
Sorry but I am worried to continue using COAWST recent version for next approach. For example nesting. I am worried that their sources need update, too. And it gives long time
What is your recommendation?

[jcwarner]

that is good to know about the stokes.
I am not sure what you mean about the nesting "sources need upate"? coawst nesting is consistent with the rutgers version. I am doing some tests with realistic applications now. and if there are problems i will correct it.