CBLAST with SWAN_COUPLING

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

Moderators: arango, robertson

Post Reply
Message
Author
nganju
Posts: 82
Joined: Mon Aug 16, 2004 8:47 pm
Location: U.S. Geological Survey, Woods Hole
Contact:

CBLAST with SWAN_COUPLING

#1 Post by nganju » Tue Apr 29, 2008 6:53 pm

I successfully compiled and ran TEST_HEAD and INLET_TEST with waves-ocean coupling. I also successfully compiled and ran the CBLAST application, without coupling. I also ran CBLAST with SWAN alone, no problems.

I compiled CBLAST successfully with SWAN_COUPLING activated, but when I try to run the application, ROMS initializes and SWAN begins computations, but the wave fields never get sent (e.g. TEST_HEAD output like this isn't there...

Code: Select all

== SWAN sent wave fields and Myerror=            0
 Wave-Ocean models coupling, receive                           t=       0.0135
 Wave-Ocean models coupling, send,                             t=       0.0135
 == SWAN recvd ocean fields and Myerror=            0
Then it looks like ROMS is sitting there, patiently waiting for data, before it begins. See output below (I killed the process after a few minutes)...

Is it possibly a time-stamp problem? In ocean_cblast.in I changed DSTART to 12426 (timeref=-2) which should correspond to 5/31/2002 (I also changed timestamps in all forcing files). I have SWAN doing computations starting on 5/31/2002. I have tried other timespans for SWAN, with the same result. Any ideas are appreciated...

Code: Select all

 Coupled Input File name = Projects/cblast/external/coupling_cblast.in

 Waves-Ocean Models Coupling: 

       Ocean Model MPI nodes: 000 - 000

       Waves Model MPI nodes: 001 - 001
 Process Information:

 Node #  0 (pid=       0) is active.

 Model Input Parameters:  ROMS/TOMS version 3.0  
                          Tuesday - April 29, 2008 - 11:27:53 AM
 -----------------------------------------------------------------------------

 SWAN is preparing computation


 ROMS3.0 CBLAST2002 MY25 Winds Clim Tides

 Operating system : CYGWIN
 CPU/hardware     : i686
 Compiler system  : ifort
 Compiler command : ifort
 Compiler flags   : -O2 /align /G7 /MD /Ox -Ic:\\ROMS\\MPICH2\\include -Ic:\\ROMS\\MCT\\mct -Ic:\\ROMS\\MCT\\mpeu   -free

 Input Script  : Projects/cblast/external/ocean_cblast.in

 SVN Root URL  : https://www.myroms.org/svn/cstm/branches/nganju/cstm_trunk
 SVN Revision  : 1346M

 Local Root    : /cygdrive/c/ROMS/nganju_branch
 Header Dir    : ./ROMS/Include
 Header file   : cblast.h
 Analytical Dir: /cygdrive/c/ROMS/nganju_branch/ROMS/Functionals

 Resolution, Grid 01: 0160x0128x020,  Parallel Nodes:   1,  Tiling: 001x001


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

        100  ntimes          Number of timesteps for 3-D equations.
    120.000  dt              Timestep size (s) for 3-D equations.
         20  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
    120.000  TI_WAV_OCN      Time interval (s) between coupling WAV-OCN models.
          1  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.
        720  nRST            Number of timesteps between the writing of data
                               into restart fields.
          1  ninfo           Number of timesteps between print of information
                               to standard output.
          1  nSTA            Number of timesteps between the writing of data
                               the stations file.
          T  ldefout         Switch to create a new output NetCDF file(s).
         90  nHIS            Number of timesteps between the writing fields
                               into history file.
        720  ndefHIS         Number of timesteps between creation of new
                               history files.
          1  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
         90  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
        720  ndefAVG         Number of timesteps between creation of new
                               time-averaged file.
          1  ntsDIA          Starting timestep for the accumulation of output
                               time-averaged diagnostics data.
         90  nDIA            Number of timesteps between the writing of
                               time-averaged data into diagnostics file.
        720  ndefDIA         Number of timesteps between creation of new
                               diagnostic file.
 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.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2           Quadratic bottom drag coefficient.
 2.0000E-02  Zob             Bottom roughness (m).
 2.0000E+00  blk_ZQ          Height (m) of surface air humidity measurement.
 2.0000E+00  blk_ZT          Height (m) of surface air temperature measurement.
 1.0000E+01  blk_ZW          Height (m) of surface winds measurement.
          1  lmd_Jwt         Jerlov water type.
 3.0000E+00  theta_s         S-coordinate surface control parameter.
 4.0000E-01  theta_b         S-coordinate bottom  control parameter.
     50.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
  12426.000  dstart          Time-stamp assigned to model initialization (days).
  12426.000  tide_start      Reference time origin for tidal forcing (days).
      -2.00  time_ref        Reference time for units attribute (yyyymmdd.dd)
 1.0000E+01  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 1.0000E+01  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.
 1.0000E+01  M3nudg          Nudging/relaxation time scale (days)
                               for 3D momentum.
 2.0000E+01  obcfac          Factor between passive and active
                               open boundary conditions.
     14.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(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(idUbms)    Write out bottom U-momentum stress.
          T  Hout(idVbms)    Write out bottom V-momentum stress.
          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  Hout(idLrad)    Write out longwave radiation flux.
          T  Hout(idLhea)    Write out latent heat flux.
          T  Hout(idShea)    Write out sensible heat flux.
          T  Hout(idVvis)    Write out vertical viscosity coefficient.
          T  Hout(idTdif)    Write out vertical T-diffusion coefficient.

 Output/Input Files:

             Output Restart File:  ./rst_cblast2002_ex.nc
      Prefix for History Files:  ./his_cblast2002_ex
       Prefix for Averages Files:  ./avg_cblast2002_ex
    Prefix for Diagnostics Files:  ./dia_cblast2002_ex
            Output Stations File:  ./sta_cblast2002_ex.nc
        Physical parameters File:  Projects/cblast/external/ocean_cblast.in
                 Input Grid File:  ./Projects/cblast/grids/cblast_grid_rps.nc
    Input Nonlinear Initial File:  ./Projects/cblast/forcing/cblast_ini_clm612002.nc
           Input Forcing File 01:  ./Projects/cblast/forcing/cblast_frc_mooringA_2002_swrad.nc
           Input Forcing File 02:  ./Projects/cblast/forcing/cblast_frc_mooringA_2002_lwrad_down.nc
           Input Forcing File 03:  ./Projects/cblast/forcing/cblast_frc_tide_v9.nc
           Input Forcing File 04:  ./Projects/cblast/forcing/cblast_frc_summer02_coamps_wind.nc
           Input Forcing File 05:  ./Projects/cblast/forcing/cblast_frc_summer02_coamps_air.nc
             Input Boundary File:  ./Projects/cblast/forcing/cblast_bry_clm02.nc
          Station positions File:  ./Projects/cblast/external/cblast_stations.in

 Tile partition information for Grid 01:  0160x0128x0020  tiling: 001x001

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1      160        1      128   409600

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     338
               HaloSizeJ(1) =     274
                TileSide(1) =     163
                TileSize(1) =   21353



 Stations Parameters, Grid: 01
 =============================

         90  Nstation        Number of stations to write out into stations file.
          T  Sout(idFsur)    Write out free-surface.
          T  Sout(idUbar)    Write out 2D U-momentum component.
          T  Sout(idVbar)    Write out 2D V-momentum component.
          T  Sout(idUvel)    Write out 3D U-momentum component.
          T  Sout(idVvel)    Write out 3D V-momentum component.
          T  Sout(idWvel)    Write out W-momentum component.
          T  Hout(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt
          T  Sout(idUsms)    Write out surface U-momentum stress.
          T  Sout(idVsms)    Write out surface V-momentum stress.
          T  Sout(idUbms)    Write out bottom U-momentum stress.
          T  Sout(idVbms)    Write out bottom V-momentum stress.
          T  Sout(idTsur)    Write out surface net heat flux.
          T  Sout(idTsur)    Write out surface net salt flux.
          T  Sout(idSrad)    Write out shortwave radiation flux.
          T  Sout(idLrad)    Write out longwave radiation flux.
          T  Sout(idLhea)    Write out latent heat flux.
          T  Sout(idShea)    Write out sensible heat flux.
          T  Sout(idVvis)    Write out vertical viscosity coefficient.
          T  Sout(idTdif)    Write out vertical T-diffusion coefficient.
          T  Sout(idMtke)    Write out turbulent kinetic energy.
          T  Sout(idMtls)    Write out turbulent generic length-scale.
 
             Flag and positions for station 0001:  0    52.0000   79.0000
             Flag and positions for station 0002:  0    53.0000   79.0000
             Flag and positions for station 0003:  0    54.0000   79.0000
             Flag and positions for station 0004:  0    55.0000   79.0000
             Flag and positions for station 0005:  0    56.0000   79.0000
             Flag and positions for station 0006:  0    52.0000   80.0000
             Flag and positions for station 0007:  0    53.0000   80.0000
             Flag and positions for station 0008:  0    54.0000   80.0000
             Flag and positions for station 0009:  0    55.0000   80.0000
             Flag and positions for station 0010:  0    56.0000   80.0000
             Flag and positions for station 0011:  0    52.0000   81.0000
             Flag and positions for station 0012:  0    53.0000   81.0000
             Flag and positions for station 0013:  0    54.0000   81.0000
             Flag and positions for station 0014:  0    55.0000   81.0000
             Flag and positions for station 0015:  0    56.0000   81.0000
             Flag and positions for station 0016:  0    52.0000   82.0000
             Flag and positions for station 0017:  0    53.0000   82.0000
             Flag and positions for station 0018:  0    54.0000   82.0000
             Flag and positions for station 0019:  0    55.0000   82.0000
             Flag and positions for station 0020:  0    56.0000   82.0000
             Flag and positions for station 0021:  0    52.0000   83.0000
             Flag and positions for station 0022:  0    53.0000   83.0000
             Flag and positions for station 0023:  0    54.0000   83.0000
             Flag and positions for station 0024:  0    55.0000   83.0000
             Flag and positions for station 0025:  0    56.0000   83.0000
             Flag and positions for station 0026:  0    50.0000   48.0000
             Flag and positions for station 0027:  0    51.0000   48.0000
             Flag and positions for station 0028:  0    52.0000   48.0000
             Flag and positions for station 0029:  0    53.0000   48.0000
             Flag and positions for station 0030:  0    54.0000   48.0000
             Flag and positions for station 0031:  0    50.0000   49.0000
             Flag and positions for station 0032:  0    51.0000   49.0000
             Flag and positions for station 0033:  0    52.0000   49.0000
             Flag and positions for station 0034:  0    53.0000   49.0000
             Flag and positions for station 0035:  0    54.0000   49.0000
             Flag and positions for station 0036:  0    50.0000   50.0000
             Flag and positions for station 0037:  0    51.0000   50.0000
             Flag and positions for station 0038:  0    52.0000   50.0000
             Flag and positions for station 0039:  0    53.0000   50.0000
             Flag and positions for station 0040:  0    54.0000   50.0000
             Flag and positions for station 0041:  0    50.0000   51.0000
             Flag and positions for station 0042:  0    51.0000   51.0000
             Flag and positions for station 0043:  0    52.0000   51.0000
             Flag and positions for station 0044:  0    53.0000   51.0000
             Flag and positions for station 0045:  0    54.0000   51.0000
             Flag and positions for station 0046:  0    50.0000   52.0000
             Flag and positions for station 0047:  0    51.0000   52.0000
             Flag and positions for station 0048:  0    52.0000   52.0000
             Flag and positions for station 0049:  0    53.0000   52.0000
             Flag and positions for station 0050:  0    54.0000   52.0000
             Flag and positions for station 0051:  0    52.0000   30.0000
             Flag and positions for station 0052:  0    52.0000   35.0000
             Flag and positions for station 0053:  0    52.0000   40.0000
             Flag and positions for station 0054:  0    52.0000   45.0000
             Flag and positions for station 0055:  0    52.0000   50.0000
             Flag and positions for station 0056:  0    52.0000   55.0000
             Flag and positions for station 0057:  0    52.0000   60.0000
             Flag and positions for station 0058:  0    52.0000   65.0000
             Flag and positions for station 0059:  0    52.0000   70.0000
             Flag and positions for station 0060:  0    20.0000   40.0000
             Flag and positions for station 0061:  0    25.0000   40.0000
             Flag and positions for station 0062:  0    30.0000   40.0000
             Flag and positions for station 0063:  0    38.0000   38.0000
             Flag and positions for station 0064:  0    45.0000   35.0000
             Flag and positions for station 0065:  0    50.0000   30.0000
             Flag and positions for station 0066:  0    55.0000   27.0000
             Flag and positions for station 0067:  0    60.0000   23.0000
             Flag and positions for station 0068:  0    65.0000   20.0000
             Flag and positions for station 0069:  0    70.0000   15.0000
             Flag and positions for station 0070:  0    75.0000   15.0000
             Flag and positions for station 0071:  0    80.0000   10.0000
             Flag and positions for station 0072:  0    52.0000   50.0000
             Flag and positions for station 0073:  0    54.0000   81.0000
             Flag and positions for station 0074:  0    57.0000   85.0000
             Flag and positions for station 0075:  0    53.0000   76.0000
             Flag and positions for station 0076:  0    65.0000   83.0000
             Flag and positions for station 0077:  0    43.0000   82.0000
             Flag and positions for station 0078:  0    54.0000   82.0000
             Flag and positions for station 0079:  0    74.0000   76.0000
             Flag and positions for station 0080:  0    77.0000   34.0000
             Flag and positions for station 0081:  0     3.0000   60.0000
             Flag and positions for station 0082:  0    80.0000    3.0000
             Flag and positions for station 0083:  0   159.0000   70.0000
             Flag and positions for station 0084:  0   135.0000  127.0000
             Flag and positions for station 0085:  0    53.0000   50.0000
             Flag and positions for station 0086:  0    38.0000   55.0000
             Flag and positions for station 0087:  0    66.0000   45.0000
             Flag and positions for station 0088:  0    89.0000   20.0000
             Flag and positions for station 0089:  0    52.0000   66.0000
             Flag and positions for station 0090:  0    53.0000   76.0000

 Activated C-preprocessing Options:

  CBLAST             ROMS3.0 CBLAST2002 MY25 Winds Clim Tides
  ADD_FSOBC          Adding tidal elevation to proccesed OBC data.
  ADD_M2OBC          Adding tidal currents to proccesed OBC data.
  ANA_BSFLUX         Analytical kinematic bottom salinity flux.
  ANA_BTFLUX         Analytical kinematic bottom temperature flux.
  ANA_RAIN           Analytical rain fall rate.
  ANA_SSFLUX         Analytical kinematic surface salinity flux.
  ASSUMED_SHAPE      Using assumed-shape arrays.
  AVERAGES           Writing out time-averaged fields.
  AVERAGES_FLUXES    Writing out time-averaged surface fluxes.
  AVERAGES_QUADRATIC Writing out time-averaged quadratic terms.
  BULK_FLUXES        Surface bulk fluxes parametererization.
  CURVGRID           Orthogonal curvilinear grid.
  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.
  EAST_FSCHAPMAN     Eastern edge, free-surface, Chapman condition.
  EAST_M2FLATHER     Eastern edge, 2D momentum, Flather condition.
  EAST_M3NUDGING     Eastern edge, 3D momentum, passive/active outflow/inflow.
  EAST_M3RADIATION   Eastern edge, 3D momentum, radiation condition.
  EAST_TNUDGING      Eastern edge, tracers, passive/active outflow/inflow.
  EAST_TRADIATION    Eastern edge, tracers, radiation condition.
  KANTHA_CLAYSON     Kantha and Clayson stability function formulation.
  LONGWAVE_OUT       Compute outgoing longwave radiation internally.
  MASKING            Land/Sea masking.
  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_M3NUDGING    Northern edge, 3D momentum, passive/active outflow/inflow.
  NORTH_M3RADIATION  Northern edge, 3D momentum, radiation condition.
  NORTH_TNUDGING     Northern edge, tracers, passive/active outflow/inflow.
  NORTH_TRADIATION   Northern edge, tracers, radiation condition.
  N2S2_HORAVG        Horizontal smoothing of buoyancy and shear.
  POWER_LAW          Power-law shape time-averaging barotropic filter.
  PROFILE            Time profiling activated .
  K_GSCHEME          Third-order upstream bias advection of TKE fields.
  RADIATION_2D       Use tangential phase speed in radiation conditions.
  RAMP_TIDES         Ramping tidal forcing for one day.
  !RST_SINGLE        Double precision fields in restart NetCDF file.
  SALINITY           Using salinity.
  SOLAR_SOURCE       Solar Radiation Source Term.
  SOLVE3D            Solving 3D Primitive Equations.
  SOUTH_FSCHAPMAN    Southern edge, free-surface, Chapman condition.
  SOUTH_M2FLATHER    Southern edge, 2D momentum, Flather condition.
  SOUTH_M3NUDGING    Southern edge, 3D momentum, passive/active outflow/inflow.
  SOUTH_M3RADIATION  Southern edge, 3D momentum, radiation condition.
  SOUTH_TNUDGING     Southern edge, tracers, passive/active outflow/inflow.
  SOUTH_TRADIATION   Southern edge, tracers, radiation condition.
  SPLINES            Conservative parabolic spline reconstruction.
  SSH_TIDES          Add tidal elevation to SSH climatology.
  STATIONS           Writing out station data.
  SWAN_COUPLING      Two-way SWAN/ROMS coupling.
  TS_U3HADVECTION    Third-order upstream bias horizontal advection of tracers.
  TS_SVADVECTION     Parabolic splines vertical advection of tracers.
  UV_ADV             Advection of momentum.
  UV_COR             Coriolis term.
  UV_U3HADVECTION    Third-order upstream bias advection of momentum.
  UV_SADVECTION      Parabolic splines vertical advection of momentum.
  UV_QDRAG           Quadratic bottom stress.
  UV_TIDES           Add tidal currents to 2D momentum climatologies.
  VAR_RHO_2D         Variable density barotropic mode.
  WAVES_OCEAN        Two-way wave-ocean models coupling.
  WEST_FSCHAPMAN     Western edge, free-surface, Chapman condition.
  WEST_M2FLATHER     Western edge, 2D momentum, Flather condition.
  WEST_M3NUDGING     Western edge, 3D momentum, passive/active outflow/inflow.
  WEST_M3RADIATION   Western edge, 3D momentum, radiation condition.
  WEST_TNUDGING      Western edge, tracers, passive/active outflow/inflow.
  WEST_TRADIATION    Western edge, tracers, radiation condition.

 INITIAL: Configurating and initializing forward nonlinear model ...


 Vertical S-coordinate System: 

 level   S-coord     Cs-curve          at_hmin  over_slope     at_hmax

    20   0.0000000   0.0000000           0.000       0.000       0.000
    19  -0.0500000  -0.0158884          -0.100      -1.653      -3.206
    18  -0.1000000  -0.0340357          -0.200      -3.526      -6.853
    17  -0.1500000  -0.0551241          -0.300      -5.688     -11.075
    16  -0.2000000  -0.0798591          -0.400      -8.205     -16.010
    15  -0.2500000  -0.1089095          -0.500     -11.144     -21.789
    14  -0.3000000  -0.1428156          -0.600     -14.558     -28.516
    13  -0.3500000  -0.1818751          -0.700     -18.476     -36.251
    12  -0.4000000  -0.2260381          -0.800     -22.892     -44.984
    11  -0.4500000  -0.2748556          -0.900     -27.763     -54.626
    10  -0.5000000  -0.3275288          -1.000     -33.011     -65.022
     9  -0.5500000  -0.3830768          -1.100     -38.540     -75.980
     8  -0.6000000  -0.4405834          -1.200     -44.261     -87.321
     7  -0.6500000  -0.4994459          -1.300     -50.114     -98.927
     6  -0.7000000  -0.5595464          -1.400     -56.088    -110.775
     5  -0.7500000  -0.6213093          -1.500     -62.224    -122.948
     4  -0.8000000  -0.6856605          -1.600     -68.613    -135.627
     3  -0.8500000  -0.7539359          -1.700     -75.386    -149.073
     2  -0.9000000  -0.8277858          -1.800     -82.704    -163.608
     1  -0.9500000  -0.9091051          -1.900     -90.752    -179.604
     0  -1.0000000  -1.0000000          -2.000     -99.735    -197.471

 Time Splitting Weights: ndtfast =  20    nfast =  29

    Primary            Secondary            Accumulated to Current Step

  1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
  2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
  3-0.0011514592651645 0.0501156998674241-0.0034654566136460 0.1501639558342179
  4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649002
  5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792858
  6 0.0029141799764308 0.0501428056118168 0.0000580677400948 0.3006719858911026
  7 0.0054132615310267 0.0499970966129952 0.0054713292711215 0.3506690825040978
  8 0.0084687837865132 0.0497264335364439 0.0139401130576347 0.4003955160405417
  9 0.0120633394191050 0.0493029943471183 0.0260034524767397 0.4496985103876600
 10 0.0161716623600090 0.0486998273761630 0.0421751148367486 0.4983983377638230
 11 0.0207585511322367 0.0478912442581626 0.0629336659689853 0.5462895820219856
 12 0.0257765478740990 0.0468533167015507 0.0887102138430843 0.5931428987235363
 13 0.0311633730493853 0.0455644893078458 0.1198735868924696 0.6387073880313821
 14 0.0368391158442262 0.0440063206553765 0.1567127027366958 0.6827137086867585
 15 0.0427031802506397 0.0421643648631652 0.1994158829873354 0.7248780735499237
 16 0.0486309868367616 0.0400292058506332 0.2480468698240970 0.7649072794005569
 17 0.0544704302037591 0.0375976565087951 0.3025173000278562 0.8025049359093520
 18 0.0600380921294285 0.0348741349986072 0.3625553921572847 0.8373790709079592
 19 0.0651152103984763 0.0318722303921358 0.4276706025557610 0.8692513013000949
 20 0.0694434033194839 0.0286164698722119 0.4971140058752449 0.8978677711723068
 21 0.0727201499285569 0.0251442997062377 0.5698341558038018 0.9230120708785445
 22 0.0745940258796570 0.0215082922098099 0.6444281816834588 0.9445203630883544
 23 0.0746596950216180 0.0177785909158270 0.7190878767050768 0.9622989540041814
 24 0.0724526566618460 0.0140456061647461 0.7915405333669228 0.9763445601689276
 25 0.0674437485167025 0.0104229733316538 0.8589842818836253 0.9867675335005814
 26 0.0590334053485720 0.0070507859058187 0.9180176872321973 0.9938183194064002
 27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447904
 28 0.0292219798521905 0.0017718319739095 0.9937853403740005 0.9996892670186999
 29 0.0062146596259994 0.0003107329813000 0.9999999999999999 0.9999999999999999

 ndtfast, nfast =   20  29   nfast/ndtfast = 1.45000

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

    1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.14200

 Minimum X-grid spacing, DXmin =  1.14712071E+00 km
 Maximum X-grid spacing, DXmax =  1.17049278E+00 km
 Minimum Y-grid spacing, DYmin =  1.12894768E+00 km
 Maximum Y-grid spacing, DYmax =  1.15194948E+00 km
 Minimum Z-grid spacing, DZmin =  1.00000000E-01 m
 Maximum Z-grid spacing, DZmax =  1.78672997E+01 m

 Minimum barotropic Courant Number =  3.27753577E-02
 Maximum barotropic Courant Number =  3.27888340E-01
 Maximum Coriolis   Courant Number =  1.16604965E-02


 NLM: GET_STATE - Read state initial conditions,               t =   12426.0000
                   (File: cblast_ini_clm612002.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -1.65490489E-02 Max =  8.02945867E-02)
                - vertically integrated u-momentum component
                   (Min = -1.47108823E-01 Max =  1.41498417E-01)
                - vertically integrated v-momentum component
                   (Min = -1.78467423E-01 Max =  1.30967095E-01)
+time 20020601.000000   , step    1; iteration    1; sweep 1
                - u-momentum component
                   (Min = -1.89053997E-01 Max =  1.60559401E-01)
                - v-momentum component
                   (Min = -2.45919555E-01 Max =  1.71909615E-01)
                - potential temperature
                   (Min =  4.20529461E+00 Max =  1.38932610E+01)
                - salinity
                   (Min =  3.08298759E+01 Max =  3.35688896E+01)
    GET_NGFLD   - tidal period
                   (Min =  4.30920000E+04 Max =  9.67320000E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Min =  1.01909846E-04 Max =  8.79483093E-01)
    GET_2DFLD   - tidal elevation phase angle
                   (Min = -3.21134579E+00 Max =  6.28318399E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Min =  0.00000000E+00 Max =  3.14158876E+00)
    GET_2DFLD   - tidal current phase angle
                   (Min =  2.31585156E-05 Max =  6.28313351E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Min =  3.06969177E-05 Max =  1.51886454E+00)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Min = -6.22084939E-01 Max =  1.78097595E-01)
    GET_2DFLD   - surface u-wind component,                    t =   12426.0000
                   (File: cblast_frc_summer02_coamps_wind.nc, Rec=0001, Index=1)
                   (Tmin=      12426.0000 Tmax=      12524.2080)
                   (Min =  8.26364934E-01 Max =  4.98556662E+00)
    GET_2DFLD   - surface v-wind component,                    t =   12426.0000
                   (File: cblast_frc_summer02_coamps_wind.nc, Rec=0001, Index=1)
                   (Tmin=      12426.0000 Tmax=      12524.2080)
                   (Min =  4.33963299E+00 Max =  1.10936642E+01)
    GET_2DFLD   - solar shortwave radiation flux,              t =   12426.0000
                   (File: cblast_frc_mooringA_2002_swrad.nc, Rec=0144, Index=1)
                   (Tmin=      12424.0139 Tmax=      12521.0621)
                   (Min =  3.00000012E-01 Max =  3.00000012E-01)
    GET_2DFLD   - downwelling longwave radiation flux,         t =   12426.0000
                   (File: cblast_frc_mooringA_2002_lwrad_down.nc, Rec=0144, Index=1)
                   (Tmin=      12424.0139 Tmax=      12521.0621)
                   (Min =  3.49000000E+02 Max =  3.49000000E+02)
    GET_2DFLD   - surface air pressure,                        t =   12426.0000
                   (File: cblast_frc_summer02_coamps_air.nc, Rec=0001, Index=1)
                   (Tmin=      12426.0000 Tmax=      12524.2080)
                   (Min =  1.00156512E+03 Max =  1.00581506E+03)
    GET_2DFLD   - surface air temperature,                     t =   12426.0000
                   (File: cblast_frc_summer02_coamps_air.nc, Rec=0001, Index=1)
                   (Tmin=      12426.0000 Tmax=      12524.2080)
                   (Min =  1.47910137E+01 Max =  1.79712734E+01)
    GET_2DFLD   - surface air relative humidity,               t =   12426.0000
                   (File: cblast_frc_summer02_coamps_air.nc, Rec=0001, Index=1)
                   (Tmin=      12426.0000 Tmax=      12524.2080)
                   (Min =  9.26330490E-01 Max =  9.90294571E-01)
    GET_NGFLD   - free-surface western boundary condition,     t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  1.54311322E-02 Max =  7.32154027E-02)
    GET_NGFLD   - free-surface eastern boundary condition,     t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.35922851E-03 Max =  1.84328966E-02)
    GET_NGFLD   - free-surface southern boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.65490489E-02 Max =  1.73498932E-02)
    GET_NGFLD   - free-surface northern boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  3.09711578E-03 Max =  7.05213025E-02)
    GET_NGFLD   - 2D u-momentum western boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -6.13183714E-02 Max =  2.18863264E-02)
    GET_NGFLD   - 2D v-momentum western boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.19962662E-01 Max =  3.96904647E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -8.43802765E-02 Max =  2.51779184E-02)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -5.00205792E-02 Max = -6.57607848E-03)
    GET_NGFLD   - 2D u-momentum southern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.47108823E-01 Max =  1.81566961E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.19962662E-01 Max =  4.80175205E-02)
    GET_NGFLD   - 2D u-momentum northern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -6.46762550E-02 Max =  7.97862336E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -9.46933553E-02 Max =  1.30967095E-01)
    GET_NGFLD   - 3D u-momentum western boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -7.93634728E-02 Max =  4.00631279E-02)
    GET_NGFLD   - 3D v-momentum western boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.65154204E-01 Max =  4.99015003E-02)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.26414552E-01 Max =  8.33319873E-02)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,    t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.17083892E-01 Max =  2.91435998E-02)
    GET_NGFLD   - 3D u-momentum southern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.89053997E-01 Max =  3.97326760E-02)
    GET_NGFLD   - 3D v-momentum southern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.87171161E-01 Max =  6.43106401E-02)
    GET_NGFLD   - 3D u-momentum northern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -7.05740526E-02 Max =  1.08818837E-01)
    GET_NGFLD   - 3D v-momentum northern boundary condition,   t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min = -1.38612866E-01 Max =  1.71909615E-01)
    GET_NGFLD   - temperature western boundary condition,      t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  5.49383068E+00 Max =  1.31007767E+01)
    GET_NGFLD   - salinity western boundary condition,         t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  3.08595428E+01 Max =  3.27931786E+01)
    GET_NGFLD   - temperature eastern boundary condition,      t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  4.22192049E+00 Max =  1.09779930E+01)
    GET_NGFLD   - salinity eastern boundary condition,         t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  3.16230736E+01 Max =  3.34656715E+01)
    GET_NGFLD   - temperature southern boundary condition,     t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  5.55179501E+00 Max =  1.26785898E+01)
    GET_NGFLD   - salinity southern boundary condition,        t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  3.19673405E+01 Max =  3.31789131E+01)
    GET_NGFLD   - temperature northern boundary condition,     t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  4.25002718E+00 Max =  1.28432846E+01)
    GET_NGFLD   - salinity northern boundary condition,        t =   12425.0000
                   (File: cblast_bry_clm02.nc, Rec=0003, Index=2)
                   (Tmin=      12305.0000 Tmax=      12605.0000)
                   (Min =  3.09957561E+01 Max =  3.35249901E+01)

 Maximum grid stiffness ratios:  rx0 =   4.387035E-01 (Beckmann and Haidvogel)
                                 rx1 =   1.070178E+01 (Haney)


 Initial basin volumes: TotVolume =  1.14859438805428E+12 m3
                        MinVolume =  1.30119903092678E+05 m3
                        MaxVolume =  2.31828546578703E+07 m3
                          Max/Min =  1.78165323727288E+02

+time 20020601.000000   , step    1; iteration    1; sweep 2
+time 20020601.000000   , step    1; iteration    1; sweep 3
+time 20020601.000000   , step    1; iteration    1; sweep 4


User avatar
arango
Site Admin
Posts: 1084
Joined: Wed Feb 26, 2003 4:41 pm
Location: IMCS, Rutgers University
Contact:

#2 Post by arango » Tue Apr 29, 2008 7:32 pm

This is the old version of the coupling. Notice that ROMS has been updated to standardize the coupling with the MCT library, :arrow: ticket 144 and the SWAN model was updated to the latest version 40.51, :arrow: ticket 146.

I recommend you to update your code to see if you still have the same problem.

nganju
Posts: 82
Joined: Mon Aug 16, 2004 8:47 pm
Location: U.S. Geological Survey, Woods Hole
Contact:

#3 Post by nganju » Tue Apr 29, 2008 8:44 pm

I updated to release 177, and now compilation fails...

Code: Select all

cd Build; ifort -c -O2 /align /G7 /MD /Ox -I/cygdrive/c/ROMS/MPICH2/include init
ial.f90 /object:initial.o
Intel(R) Fortran Compiler for 32-bit applications, Version 9.0    Build 20050430
Z Package ID: W_FC_P_9.0.018
Copyright (C) 1985-2005 Intel Corporation.  All rights reserved.

initial.f90(31) : Error: Name in only-list does not exist.   [OCN2WAV_COUPLING]
      USE ocean_coupler_mod, ONLY : ocn2wav_coupling
------------------------------------^
compilation aborted for initial.f90 (code 1)
make: *** [Build/initial.o] Error 1
I have the latest MCT release, installed outside of the ROMS directory. Is it critical to use the MCT release distributed with ROMS?

User avatar
arango
Site Admin
Posts: 1084
Joined: Wed Feb 26, 2003 4:41 pm
Location: IMCS, Rutgers University
Contact:

#4 Post by arango » Tue Apr 29, 2008 9:50 pm

What application are you running?

This a good example of the compilation errors that users can get when using ROMS. Users need to be curious, read carefully the compiling errors, and track it:
initial.f90(31) : Error: Name in only-list does not exist. [OCN2WAV_COUPLING]
USE ocean_coupler_mod, ONLY : ocn2wav_coupling
------------------------------------^
compilation aborted for initial.f90 (code 1)
make: *** [Build/initial.o] Error 1
It tells you that the problem is in initial.F. The compiler cannot find ocn2wav_coupling in the module ocean_coupler_mod. This module is located in file Mater/ocean_coupler.F. If you don't know where that file is located, just do a grep in the Build directory:

Code: Select all

cd Build
grep ocean_coupler *.f90
ocean_coupler.f90:      MODULE ocean_coupler_mod
ocean_coupler.f90:      END MODULE ocean_coupler_mod
Then, go and find ocean_coupler.F in the source code directory tree. If you check this file, you see something like:

Code: Select all

# ifdef SWAN_COUPLING
#  include "mct_roms_swan.h"
# endif
# ifdef REFDIF_COUPLING
#  include "mct_roms_refdif.h"
# endif
# ifdef WRF_COUPLING
#  include "mct_roms_wrf.h"
# endif
You will see that the routine ocn2wav_coupling is included in file mct_roms_swan.h. If the compiler cannot find this routine is because it was not included during C-preprocessing. Probably, there is nothing at the bottom of ocean_coupler.f90. Did you activated SWAN_COUPLING :?:

Notice that in a previous :arrow: post. I mentioned that you need to have the following in your header file:

Code: Select all

#define SWAN_COUPLING
#ifdef SWAN_COUPLING
# define MCT_LIB
#endif
Mostly all the time, we see a print out of the compiler errors posted in this forum. The errors sometimes are very clear. Users need to be a detective and track them to see what they may be doing something wrong. We do our best to test all the new features that we add to ROMS.

nganju
Posts: 82
Joined: Mon Aug 16, 2004 8:47 pm
Location: U.S. Geological Survey, Woods Hole
Contact:

#5 Post by nganju » Tue Apr 29, 2008 10:20 pm

Hello Hernan-
Thanks for the information; I try my best to troubleshoot before posting to the board. I don't expect to have my problems solved by others, I only want to share the information in case others duplicate the problem.

In the original post, what is odd is that TEST_HEAD and INLET_TEST run fine, but the CBLAST application did not (which I ran successfully as separate ROMS and SWAN runs). This was CSTM revision 1498, with MCT v2.5.

With the revision 177, I get "47 unresolved externals" when compiling TEST_HEAD, INLET_TEST, and CBLAST (with #define MCT_LIB added). They are all MPI or MCT-related. I had a previous post about this same problem, which was resovled when I went to the version mentioned in the above paragraph.

I am now totally confused about how to proceed, because other posts suggest that r166 was stable with SWAN, but the current revision is not (e.g. viewtopic.php?p=2815)

User avatar
arango
Site Admin
Posts: 1084
Joined: Wed Feb 26, 2003 4:41 pm
Location: IMCS, Rutgers University
Contact:

#6 Post by arango » Wed Apr 30, 2008 12:01 am

I think that the problem is yours. If you are having problems with the MPI and MCT libraries, this is very trivial. The compiler is not finding those libraries :!: It is as simple as that. This doesn't have anything to do with ROMS source code but your compiling configuration file (Compilers/*.mk) and the full path for those libraries. Recall that this are user adjusted file. There is plenty of information about this in :arrow: wikiROMS.

The compiler always gives you the library paths during compilation. Check. You can check if that library path actually exist in your computer.

I ran all these coupling test cases with no problems. I did a correction but that was related to a compiler bug. I made the code more robust to about this error with IFORT.

nganju
Posts: 82
Joined: Mon Aug 16, 2004 8:47 pm
Location: U.S. Geological Survey, Woods Hole
Contact:

Re: CBLAST with SWAN_COUPLING

#7 Post by nganju » Mon Jul 14, 2008 8:53 pm

Just to update, my problem (for both versions of the code) was that I had commented out the matlab file output commands (BLOCK) in the SWAN input file. Apparently those output commands must be specified for SWAN to pass information to ROMS.

One additional note: SWAN's begin time is not considered when the models are coupled, they both begin running and coupling at the same time. All that is considered is the time step and coupling interval.

One more note: Regardless of your SWAN convention for directions (nautical vs. cartesian), wind vector components (as with water currents) are always cartesian: the direction the vector is pointing.

Post Reply