forrtl: severe (174): SIGSEGV, segmentation fault occurred

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
manhnt181
Posts: 33
Joined: Wed Feb 03, 2021 5:07 pm
Location: Institute of Geophysics

forrtl: severe (174): SIGSEGV, segmentation fault occurred

#1 Unread post by manhnt181 »

Hi everyone,
I got the error when I ran the model with WRF and ROMS option which is described in the below:

Code: Select all

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
coawstM            00000000038D5BFA  Unknown               Unknown  Unknown
libpthread-2.17.s  00002B4AF22F95D0  Unknown               Unknown  Unknown
coawstM            00000000039C10E3  Unknown               Unknown  Unknown
coawstM            000000000393937C  Unknown               Unknown  Unknown
coawstM            00000000004219BE  cdecode_line_             363  read_coawst_par.f90
coawstM            0000000000422897  read_model_inputs          59  read_model_inputs.f90
coawstM            000000000041A37D  MAIN__                     68  master.f90
coawstM            0000000000417AE2  Unknown               Unknown  Unknown
libc-2.17.so       00002B4AF2A2E495  __libc_start_main     Unknown  Unknown
The 363 line of the read_coawst_par.f90 program is:

Code: Select all

Cval(Nval)=TRIM(ADJUSTL(string))
The 59 line of the read_model_inputs.f90 program is:

Code: Select all

status=cdecode_line(line, KeyWord, Nval, Cval, Rval)
My couple.in file is:

Code: Select all

! Number of parallel nodes assigned to each model in the coupled system.
! Their sum must be equal to the total number of processors.

   NnodesATM =  70                 ! atmospheric model
   NnodesWAV =  0                 ! wave model
   NnodesOCN =  30                 ! ocean model
   NnodesHYD =  0                  ! hydrology model

! Time interval (seconds) between coupling of models.

  TI_ATM2WAV =   0.0d0              ! atmosphere to wave coupling interval
  TI_ATM2OCN =   1800.0d0              ! atmosphere to ocean coupling interval
  TI_WAV2ATM =   0.0d0              ! wave to atmosphere coupling interval
  TI_WAV2OCN =   0.0d0              ! wave to ocean coupling interval
  TI_OCN2WAV =   0.0d0              ! ocean to wave coupling interval
  TI_OCN2ATM =   1800.0d0              ! ocean to atmosphere coupling interval
  TI_OCN2HYD =   0.0d0              ! ocean to hydro coupling interval
  TI_HYD2OCN =   0.0d0              ! hydro to ocean coupling interval

! Enter names of Atm, Wav, and Ocn input files.
! The Wav program needs multiple input files, one for each grid.

   ATM_name = Projects/Doksuri/namelist.input                            ! atmospheric model
!  WAV_name = Projects/Doksuri/swan.in
!              Projects/Doksuri/swan_Doksuri_ref3.in         ! wave model
!  WAV_name = ww3_grid.inp
   OCN_name = Projects/Doksuri/ocean.in             ! ocean model
   HYD_name = hydro.namelist                            ! hydro model
This is my ROMS's namelist:

Code: Select all

       TITLE = Hurricane Doksuri
! C-preprocessing Flag.

    MyAppCPP = Doksuri
! Input variable information file name.  This file needs to be processed
! first so all information arrays can be initialized properly.

     VARNAME = ROMS/External/varinfo.dat

! Number of nested grids.

      Ngrids =  1

! Number of grid nesting layers.  This parameter is used to allow refinement
! and composite grid combinations.

  NestLayers =  1

! Number of grids in each nesting layer [1:NestLayers].

GridsInLayer =  1 1

! Grid dimension parameters. See notes below in the Glossary for how to set
! these parameters correctly.

          Lm == 475           ! Number of I-direction INTERIOR RHO-points -2
          Mm == 301           ! Number of J-direction INTERIOR RHO-points -2
           N == 16            ! Number of vertical levels
          ND ==  0            ! Number of wave directional bins

        Nbed =  0             ! Number of sediment bed layers
       NBAND =  30            ! Number of spectral irradiance bands 

         NAT =  2             ! Number of active tracers (usually, 2)
         NPT =  0             ! Number of inactive passive tracers
         NCS =  0             ! Number of cohesive (mud) sediment tracers
         NNS =  0             ! Number of non-cohesive (sand) sediment tracers

! Domain decomposition parameters for serial, distributed-memory or
! shared-memory configurations used to determine tile horizontal range
! indices (Istr,Iend) and (Jstr,Jend), [1:Ngrids].

      NtileI == 6                                ! I-direction partition
      NtileJ == 5                                ! J-direction partition

! Set horizontal and vertical advection schemes for active and inert
! tracers. A different advection scheme is allowed for each tracer.
! For example, a positive-definite (monotonic) algorithm can be activated
! for salinity and inert tracers, while a different one is set for
! temperature. [1:NAT+NPT,Ngrids] values are expected.
!
!   Keyword    Advection Algorithm
!
!   A4         4th-order Akima (horizontal/vertical)
!   C2         2nd-order centered differences (horizontal/vertical)
!   C4         4th-order centered differences (horizontal/vertical)
!   HSIMT      3th-order HSIMT-TVD (horizontal/vertical)
!   MPDATA     recursive flux corrected MPDATA (horizontal/vertical)
!   SPLINES    parabolic splines (only vertical)
!   SU3        split third-order upstream (horizontal/vertical)
!   U3         3rd-order upstream-biased (only horizontal)
!
! The user has the option of specifying the full Keyword or the first
! two letters, regardless if using uppercase or lowercase. If nested
! grids, specify values for each grid (see glossary below).

   Hadvection == U3   \  
                 U3   \                        ! temperature
                 U3   \ 
                 U3                            ! salinity

   Vadvection == C4   \
                 C4   \                        ! temperature
                 C4   \
                 C4                            ! salinity

! Adjoint-based algorithms can have different horizontal and schemes
! for active and inert tracers.

ad_Hadvection == U3       \                     ! temperature
                 U3                             ! salinity

ad_Vadvection == C4       \                     ! temperature
                 C4                             ! salinity

! Set lateral boundary conditions keyword. Notice that a value is expected
! for each boundary segment per nested grid for each state variable.
!
! Each tracer variable requires [1:4,1:NAT+NPT,Ngrids] values. Otherwise,
! [1:4,1:Ngrids] values are expected for other variables. The boundary
! order is: 1=west, 2=south, 3=east, and 4=north. That is, anticlockwise
! starting at the western boundary.
!
! The keyword is case insensitive and usually has three characters. However,
! it is possible to have compound keywords, if applicable. For example, the
! keyword "RadNud" implies radiation boundary condition with nudging. This
! combination is usually used in active/passive radiation conditions.
!
!   Keyword    Lateral Boundary Condition Type
!
!   Cha        Chapman_implicit (free-surface)
!   Che        Chapman_explicit (free-surface)
!   Cla        Clamped
!   Clo        Closed
!   Fla        Flather (2D momentum)                  _____N_____     j=Mm
!   Gra        Gradient                              |     4     |
!   Nes        Nested (refinement)                   |           |
!   Nud        Nudging                             1 W           E 3
!   Per        Periodic                              |           |
!   Rad        Radiation                             |_____S_____|
!   Red        Reduced Physics (2D momentum)               2          j=1
!   Shc        Shchepetkin (2D momentum)            i=1         i=Lm
!
!                   W       S       E       N
!                   e       o       a       o
!                   s       u       s       r
!                   t       t       t       t
!                           h               h
!
!                   1       2       3       4

   LBC(isFsur) ==   Clo    Cha    Cha    Cha        ! free-surface
   LBC(isUbar) ==   Clo    Fla    Fla    Fla       ! 2D U-momentum
   LBC(isVbar) ==   Clo    Fla    Fla    Fla        ! 2D V-momentum
   LBC(isUvel) ==   Clo RadNud RadNud RadNud        ! 3D U-momentum
   LBC(isVvel) ==   Clo RadNud RadNud RadNud        ! 3D V-momentum
   LBC(isMtke) ==   Clo     Gra    Gra    Gra        ! mixing TKE

   LBC(isTvar) ==   Clo RadNud RadNud RadNud \       ! temperature
                    Clo RadNud RadNud RadNud         ! salinity
! Wec boundary conditions

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

! InWave boundary conditions

!   LBC(isAC3d) ==   Gra     Clo     Gra     Cla         ! 3D wave action density
!   LBC(isCT3d) ==   Gra     Clo     Gra     Gra         ! 3D wave theta celerity
!   LBC(isCX3d) ==   Gra     Clo     Gra     Gra         ! 3D wave x-dir celerity
!   LBC(isCY3d) ==   Gra     Clo     Gra     Gra         ! 3D wave y-dir celerity

! Ice boundary conditions

   LBC(isAice) ==   Clo     Clo     Clo     Clo         ! ice concentration
   LBC(isHice) ==   Clo     Clo     Clo     Clo         ! ice thickness
   LBC(isHsno) ==   Clo     Clo     Clo     Clo         ! snow thickness
   LBC(isTice) ==   Clo     Clo     Clo     Clo         ! ice temperature
!   LBC(isApond)==   Clo     Clo     Clo     Clo         ! surface water
   LBC(isSfwat)==   Clo     Clo     Clo     Clo         ! surface water
   LBC(isSig11)==   Clo     Clo     Clo     Clo         ! sigma-11
   LBC(isSig12)==   Clo     Clo     Clo     Clo         ! sigma-12
   LBC(isSig22)==   Clo     Clo     Clo     Clo         ! sigma-22
   LBC(isUice) ==   Clo     Clo     Clo     Clo         ! ice U-momentum
   LBC(isVice) ==   Clo     Clo     Clo     Clo         ! ice V-momentum

! Adjoint-based algorithms can have different lateral boundary
! conditions keywords.

ad_LBC(isFsur) ==   Cha     Per     Clo     Per         ! free-surface
ad_LBC(isUbar) ==   Fla     Per     Clo     Per         ! 2D U-momentum
ad_LBC(isVbar) ==   Fla     Per     Clo     Per         ! 2D U-momentum
ad_LBC(isUvel) ==   Gra     Per     Clo     Per         ! 3D U-momentum
ad_LBC(isVvel) ==   Gra     Per     Clo     Per         ! 3D V-momentum
ad_LBC(isMtke) ==   Gra     Per     Clo     Per         ! mixing TKE

ad_LBC(isTvar) ==   Gra     Per     Clo     Per \       ! temperature
                    Gra     Per     Clo     Per         ! salinity

! Set lateral open boundary edge volume conservation switch for
! nonlinear model and adjoint-based algorithms. Usually activated
! with radiation boundary conditions to enforce global mass
! conservation, except if tidal forcing is enabled. [1:Ngrids].

   VolCons(west)  ==  F                            ! western  boundary
   VolCons(east)  ==  F                            ! eastern  boundary
   VolCons(south) ==  F                            ! southern boundary
   VolCons(north) ==  F                            ! northern boundary

ad_VolCons(west)  ==  F                            ! western  boundary
ad_VolCons(east)  ==  F                            ! eastern  boundary
ad_VolCons(south) ==  F                            ! southern boundary
ad_VolCons(north) ==  F                            ! northern boundary

! Time-Stepping parameters.

      NTIMES == 86400    1
          DT == 4.0d0  3.0d0
     NDTFAST == 28      28
! Number of timesteps for computing observation impacts during the
! analysis-forecast cycle.

!  NTIMES_ANA == 1800                               ! analysis interval
!  NTIMES_FCT == 1800                               ! forecast interval

! Model iteration loops parameters.

       ERstr =  1
       ERend =  1
      Nouter =  1
      Ninner =  1
  Nintervals =  1

! Number of eigenvalues (NEV) and eigenvectors (NCV) to compute for the
! Lanczos/Arnoldi problem in the Generalized Stability Theory (GST)
! analysis. NCV must be greater than NEV (see documentation below).

         NEV =  2                               ! Number of eigenvalues
         NCV =  10                              ! Number of eigenvectors

! Input/Output parameters.

       NRREC == 0    0
   LcycleRST == F    T
        NRST == 5400   120
        NSTA == 1    1
        NFLT == 1    1
       NINFO == 1    1

! Output history, average, diagnostic files parameters.

     LDEFOUT == T    T
        NHIS == 5400 
     NDEFHIS == 0 
        NQCK == 0
     NDEFQCK == 0
      NTSAVG == 1    1
        NAVG == 5400 
     NDEFAVG == 0 
      NTSDIA == 1    
        NDIA == 5400
     NDEFDIA == 0 

! Output tangent linear and adjoint models parameters.

   LcycleTLM == F
        NTLM == 720
     NDEFTLM == 0
   LcycleADJ == F
        NADJ == 720
     NDEFADJ == 0
        NSFF == 720
        NOBC == 720

! GST output and check pointing restart parameters.

   LmultiGST =  F                               ! one eigenvector per file
     LrstGST =  F                               ! GST restart switch
  MaxIterGST =  500                             ! maximum number of iterations
        NGST =  10                              ! check pointing interval

! Relative accuracy of the Ritz values computed in the GST analysis.

    Ritz_tol =  1.0d-15

! Harmonic/biharmonic horizontal diffusion of tracer for nonlinear model
! and adjoint-based algorithms: [1:NAT+NPT,Ngrids].

        TNU2 == 0.2d0  0.2d0 0.2d0  0.2d0       ! m2/s
        TNU4 == 0.0d0  0.0d0                    ! m4/s

     ad_TNU2 == 0.0d0  0.0d0                    ! m2/s
     ad_TNU4 == 0.0d0  0.0d0                    ! m4/s

! Harmonic/biharmonic, horizontal viscosity coefficient for nonlinear model
! and adjoint-based algorithms: [Ngrids].

       VISC2 == 0.10d0 0.10d0                   ! m2/s
       VISC4 == 0.0d0                           ! m4/s

    ad_VISC2 == 0.0d0                           ! m2/s
    ad_VISC4 == 0.0d0                           ! m4/s

! Logical switches (TRUE/FALSE) to increase/decrease horizontal viscosity
! and/or diffusivity in specific areas of the application domain (like
! sponge areas) for the desired application grid.

    LuvSponge == F                              ! horizontal momentum
LtracerSponge == F F                            ! temperature, salinity, inert

! Vertical mixing coefficients for tracers in nonlinear model and
! basic state scale factor in adjoint-based algorithms: [1:NAT+NPT,Ngrids]

     AKT_BAK == 1.0d-6  1.0d-6 1.0d-6  1.0d-6   ! m2/s

  ad_AKT_fac == 1.0d0  1.0d0                    ! nondimensional

! Vertical mixing coefficient for momentum for nonlinear model and
! basic state scale factor in adjoint-based algorithms: [Ngrids].

     AKV_BAK == 1.0d-5                    ! m2/s

  ad_AKV_fac == 1.0d0                           ! nondimensional

! Upper threshold values to limit vertical mixing coefficients computed
! from vertical mixing parameterizations. Although this is an engineering
! fix, the vertical mixing values inferred from ocean observations are
! rarely higher than this upper limit value.

!   AKT_LIMIT == 1.0d-3 1.0d-3                   ! m2/s

!   AKV_LIMIT == 1.0d-3 1.0d-3                   ! m2/s

! Turbulent closure parameters.

     AKK_BAK == 5.0d-6                          ! m2/s
     AKP_BAK == 5.0d-6                          ! m2/s
      TKENU2 == 0.0d0                           ! m2/s
      TKENU4 == 0.0d0                           ! m4/s

! Generic length-scale turbulence closure parameters.

       GLS_P == 3.0d0                           ! K-epsilon
       GLS_M == 1.5d0
       GLS_N == -1.0d0
    GLS_Kmin == 7.6d-6
    GLS_Pmin == 1.0d-12

    GLS_CMU0 == 0.5477d0
      GLS_C1 == 1.44d0
      GLS_C2 == 1.92d0
     GLS_C3M == -0.4d0
     GLS_C3P == 1.0d0
    GLS_SIGK == 1.0d0
    GLS_SIGP == 1.30d0

! Constants used in surface turbulent kinetic energy flux computation.

  CHARNOK_ALPHA == 1400.0d0                     ! Charnok surface roughness
 ZOS_HSIG_ALPHA == 0.5d0                        ! roughness from wave amplitude
       SZ_ALPHA == 0.25d0                       ! roughness from wave dissipation
      CRGBAN_CW == 100.0d0                      ! Craig and Banner wave breaking
      WEC_ALPHA == 0.0d0                        ! 0: all wave dissip goes to break and none to roller.
                                                ! 1: all wave dissip goes to roller and none to breaking.

! Constants used in momentum stress computation.

        RDRG == 3.0d-04                    ! m/s
       RDRG2 == 0.025d0                    ! nondimensional
         Zob == 0.02d0                     ! m
         Zos == 0.02d0                     ! m

! Height (m) of atmospheric measurements for Bulk fluxes parameterization.

      BLK_ZQ ==  2.0d0                     ! air humidity
      BLK_ZT ==  2.0d0                     ! air temperature
      BLK_ZW == 10.0d0                     ! winds

! Minimum depth for wetting and drying.

       DCRIT == 0.10d0                     ! m

! Various parameters.

       WTYPE == 1
     LEVSFRC == 15
     LEVBFRC == 1

! Set vertical, terrain-following coordinates transformation equation and
! stretching function (see below for details), [1:Ngrids].

  Vtransform == 2                          ! transformation equation
 Vstretching == 4                          ! stretching function

! Vertical S-coordinates parameters (see below for details), [1:Ngrids].

     THETA_S == 10.0d0                      ! surface stretching parameter
     THETA_B == 0.4d0                      ! bottom  stretching parameter
      TCLINE == 50.0d0                     ! critical depth (m)

! Mean Density and Brunt-Vaisala frequency.

        RHO0 =  1025.0d0                   ! kg/m3
     BVF_BAK =  1.0d-5                     ! 1/s2

! Time-stamp assigned for model initialization, reference time
! origin for tidal forcing, and model reference time for output
! NetCDF units attribute.

      DSTART =  58009.0d0                      ! days
  TIDE_START =  58009.0d0                      ! days
    TIME_REF =  18581117.0d0               ! yyyymmdd.dd

! Nudging/relaxation time scales, inverse scales will be computed
! internally, [1:Ngrids].

       TNUDG == 1.0d0  1.0d0               ! days
       ZNUDG == 0.0d0                      ! days
      M2NUDG == 0.0d0                      ! days
      M3NUDG == 1.0d0                      ! days

! Nudging/relaxation time scale for surface salinity nudging, inverse
! scales will be computed internally, [1:Ngrids].

       TNUDG_SSS == 90.0d0                 ! days

! Threshold to trigger SSS correction toward climatolgy,
! needs SCORRECTION and SSSC_THRESHOLD defined

       SSS_MISMATCH_THRESHOLD = 0.2d0

! Factor between passive (outflow) and active (inflow) open boundary
! conditions, [1:Ngrids]. If OBCFAC > 1, nudging on inflow is stronger
! than on outflow (recommended).

      OBCFAC == 0.0d0                      ! nondimensional

! Linear equation of State parameters:

          R0 == 1027.0d0                   ! kg/m3
          T0 == 10.0d0                     ! Celsius
          S0 == 30.0d0                     ! nondimensional
       TCOEF == 1.7d-4                     ! 1/Celsius
       SCOEF == 7.6d-4                     ! nondimensional

! Slipperiness parameter: 1.0 (free slip) or -1.0 (no slip)

      GAMMA2 == 1.0d0

! Logical switches (TRUE/FALSE) to activate horizontal momentum transport
! point Sources/Sinks (like river runoff transport) and mass point
! Sources/Sinks (like volume vertical influx), [1:Ngrids].

      LuvSrc == F                          ! horizontal momentum transport
       LwSrc == F                          ! volume vertical influx

! Logical switches (TRUE/FALSE) to activate tracers point Sources/Sinks
! (like river runoff) and to specify which tracer variables to consider:
! [1:NAT+NPT,Ngrids].  See glossary below for details.

  LtracerSrc == F F                        ! temperature, salinity, inert

! Logical switches (TRUE/FALSE) to read and process climatology fields.
! See glossary below for details.

     LsshCLM == F                          ! sea-surface height
      Lm2CLM == F                          ! 2D momentum
      Lm3CLM == F                          ! 3D momentum

  LtracerCLM == F F                  ! temperature, salinity, inert

! Logical switches (TRUE/FALSE) to nudge the desired climatology field(s).
! If not analytical climatology fields, users need to turn ON the logical
! switches above to process the fields from the climatology NetCDF file
! that are needed for nudging. See glossary below for details.
/Sar
 LnudgeM2CLM == F                          ! 2D momentum
 LnudgeM3CLM == F                          ! 3D momentum

  LnudgeTCLM == F F                 ! temperature, salinity, inert
! Logical switches for ice climatology and nudging

      LmiCLM == F
      LaiCLM == F
      LsiCLM == F
 LnudgeMICLM == F
 LnudgeAICLM == F
 LnudgeSICLM == F

! Starting (DstrS) and ending (DendS) day for adjoint sensitivity forcing.
! DstrS must be less or equal to DendS. If both values are zero, their
! values are reset internally to the full range of the adjoint integration.

       DstrS == 0.0d0                      ! starting day
       DendS == 0.0d0                      ! ending day

! Starting and ending vertical levels of the 3D adjoint state variables
! whose sensitivity is required.

       KstrS == 1                          ! starting level
       KendS == 1                          ! ending level

! Logical switches (TRUE/FALSE) to specify the adjoint state variables
! whose sensitivity is required.

Lstate(isFsur) == F                        ! free-surface
Lstate(isUbar) == F                        ! 2D U-momentum
Lstate(isVbar) == F                        ! 2D V-momentum
Lstate(isUvel) == F                        ! 3D U-momentum
Lstate(isVvel) == F                        ! 3D V-momentum
Lstate(isWvel) == F

Lstate(isTvar) == F F                      ! NT tracers

! Logical switches (TRUE/FALSE) to specify the state variables for
! which Forcing Singular Vectors or Stochastic Optimals is required.

Fstate(isFsur) == F                        ! free-surface
Fstate(isUbar) == F                        ! 2D U-momentum
Fstate(isVbar) == F                        ! 2D V-momentum
Fstate(isUvel) == F                        ! 3D U-momentum
Fstate(isVvel) == F                        ! 3D V-momentum
Fstate(isTvar) == F F                      ! NT tracers

Fstate(isUstr) == T                        ! surface U-stress
Fstate(isVstr) == T                        ! surface V-stress
Fstate(isTsur) == F F                      ! NT surface tracers flux

! Stochastic Optimals time decorrelation scale (days) assumed for
! red noise processes.

      SO_decay == 2.0d0                    ! days

! Stochastic Optimals surface forcing standard deviation for
! dimensionalization.

SO_sdev(isFsur) == 1.0d0                   ! free-surface
SO_sdev(isUbar) == 1.0d0                   ! 2D U-momentum
SO_sdev(isVbar) == 1.0d0                   ! 2D V-momentum
SO_sdev(isUvel) == 1.0d0                   ! 3D U-momentum
SO_sdev(isVvel) == 1.0d0                   ! 3D V-momentum
SO_sdev(isTvar) == 1.0d0 1.0d0             ! NT tracers

SO_sdev(isUstr) == 1.0d0                   ! surface U-stress
SO_sdev(isVstr) == 1.0d0                   ! surface V-stress
SO_sdev(isTsur) == 1.0d0 1.0d0             ! NT surface tracers flux

! Logical switches (TRUE/FALSE) to activate writing of fields into
! HISTORY output file.

Hout(idUvel) == T      ! u                  3D U-velocity
Hout(idVvel) == T      ! v                  3D V-velocity
Hout(idu3dE) == F      ! u_eastward         3D U-eastward  at RHO-points
Hout(idv3dN) == F      ! v_northward        3D V-northward at RHO-points
Hout(idWvel) == T      ! w                  3D W-velocity
Hout(idOvel) == T      ! omega              omega vertical velocity
Hout(idUbar) == T      ! ubar               2D U-velocity
Hout(idVbar) == T      ! vbar               2D V-velocity
Hout(idu2dE) == F      ! ubar_eastward      2D U-eastward  at RHO-points
Hout(idv2dN) == F      ! vbar_northward     2D V-northward at RHO-points
Hout(idFsur) == T      ! zeta               free-surface
Hout(idBath) == T      ! bath               time-dependent bathymetry

Hout(idTvar) == T T  ! temp, salt         temperature and salinity
Hout(idUair) == F     ! Uwind              surface U-wind
Hout(idVair) == F      ! Vwind              surface V-wind
Hout(idUairE) == F     ! Uwind_eastward     surface U-wind
Hout(idVairN) == F     ! Vwind_northward    surface V-wind
Hout(idUsms) == T      ! sustr              surface U-stress
Hout(idVsms) == T      ! svstr              surface V-stress
Hout(idUbms) == T      ! bustr              bottom U-stress
Hout(idVbms) == T      ! bvstr              bottom V-stress
Hout(idUbrs) == T      ! bustrc             bottom U-current stress
Hout(idVbrs) == T      ! bvstrc             bottom V-current stress
Hout(idUbws) == T      ! bustrw             bottom U-wave stress
Hout(idVbws) == T      ! bvstrw             bottom V-wave stress
Hout(idUbcs) == T      ! bustrcwmax         bottom max wave-current U-stress
Hout(idVbcs) == T      ! bvstrcwmax         bottom max wave-current V-stress
Hout(idUVwc) == T      ! bstrcwmax          bottom max wave-current stress magnitude
Hout(idUbot) == T      ! Ubot               bed wave orbital U-velocity
Hout(idVbot) == T      ! Vbot               bed wave orbital V-velocity													 
Hout(idUbur) == T      ! Ur                 bottom U-velocity above bed																			 
Hout(idVbvr) == T      ! Vr                 bottom V-velocity above bed
Hout(idW2xx) == F      ! Sxx_bar            WEC_Mellor 2D Sxx radiation stress
Hout(idW2xy) == F      ! Sxy_bar            WEC_Mellor 2D Sxy radiation stress
Hout(idW2yy) == F      ! Syy_bar            WEC_Mellor 2D Syy radiation stress
Hout(idW3xx) == F      ! Sxx                WEC_Mellor 3D Sxx radiation stress
Hout(idW3xy) == F      ! Sxy                WEC_Mellor 3D Sxy radiation stress
Hout(idW3yy) == F      ! Syy                WEC_Mellor 3D Syy radiation stress
Hout(idW3zx) == F      ! Szx                WEC_Mellor 3D Szx radiation stress
Hout(idW3zy) == F      ! Szy                WEC_Mellor 3D Szy radiation stress
Hout(idWztw) == T      ! zetaw              WEC_VF quasi-static sea level adjustment
Hout(idWqsp) == T      ! qsp                WEC_VF quasi-static pressure
Hout(idWbeh) == T      ! bh                 WEC_VF Bernoulli head
Hout(idU2rs) == F      ! ubar_Wecstress     WEC 2D U-stress
Hout(idV2rs) == F      ! vbar_Wecstress     WEC 2D V-stress
Hout(idU3rs) == F      ! u_Rstress          WEC 3D U-stress
Hout(idV3rs) == F      ! v_Rstress          WEC 3D V-stress
Hout(idU2Sd) == T      ! ubar_stokes        2D U-Stokes velocity
Hout(idV2Sd) == T      ! vbar_stokes        2D V-Stokes velocity
Hout(idU3Sd) == T      ! u_stokes           3D U-Stokes velocity
Hout(idV3Sd) == T      ! v_stokes           3D V-Stokes velocity
Hout(idW3Sd) == T      ! omega_stokes       3D Omega-Stokes velocity
Hout(idW3St) == T      ! w_stokes           3D W-Stokes velocity
Hout(idWamp) == T      ! Hwave              wave height
Hout(idWlen) == T      ! Lwave              wave length-mean
Hout(idWlep) == T      ! Lwavep             wave length-peak
Hout(idWdir) == T      ! Dwave              wave direction
Hout(idWptp) == T      ! Pwave_top          wave surface period
Hout(idWpbt) == T      ! Pwave_bot          wave bottom period
Hout(idWorb) == T      ! Uwave_rms          wave bottom orbital velocity
Hout(idWbrk) == F      ! Wave_break         wave breaking (percent)
Hout(idUwav) == T      ! uWave              wave-depth avgeraged U-velocity
Hout(idVwav) == T      ! vWave              wave-depth avgeraged V-velocity
Hout(idWdif) == F      ! Dissip_fric        wave dissipation due to bottom friction
Hout(idWdib) == F      ! Dissip_break       wave dissipation due to breaking
Hout(idWdiw) == F      ! Dissip_wcap        wave dissipation due to white capping
Hout(idWdis) == F      ! Dissip_roller      wave roller dissipation
Hout(idWrol) == F      ! rollA              wave roller action density
Hout(idRunoff) == F    ! Runoff             surface runoff from land
Hout(idPair) == F      ! Pair               surface air pressure
Hout(idUair) == F      ! Uair               surface U-wind component
Hout(idVair) == F      ! Vair               surface V-wind component
Hout(idTsur) == F F ! shflux, ssflux     surface net heat and salt flux
Hout(idLhea) == F      ! latent             latent heat flux
Hout(idShea) == F      ! sensible           sensible heat flux
Hout(idLrad) == F      ! lwrad              longwave radiation flux
Hout(idSrad) == F      ! swrad              shortwave radiation flux
Hout(idEmPf) == F      ! EminusP            E-P flux
Hout(idevap) == F      ! evaporation        evaporation rate
Hout(idrain) == F      ! rain               precipitation rate
Hout(idDano) == F      ! rho                density anomaly
Hout(idVvis) == T      ! AKv                vertical viscosity
Hout(idTdif) == F      ! AKt                vertical T-diffusion
Hout(idSdif) == F      ! AKs                vertical Salinity diffusion
Hout(idHsbl) == F      ! Hsbl               depth of surface boundary layer
Hout(idHbbl) == F      ! Hbbl               depth of bottom boundary layer
Hout(idMtke) == T      ! tke                turbulent kinetic energy
Hout(idMtls) == T      ! gls                turbulent length scale

! Logical switches (TRUE/FALSE) to activate writing of extra inert passive
! tracers other than biological and sediment tracers. An inert passive tracer
																		  
! is one that it is only advected and diffused. Other processes are ignored.
! These tracers include, for example, dyes, pollutants, oil spills, etc.
! NPT values are expected. However, these switches can be activated using
! compact parameter specification.

 Hout(inert) == T       ! dye_01, ...        inert passive tracers

! Logical switches (TRUE/FALSE) to activate writing of exposed sediment

Qout(idUvel) == F       ! u                  3D U-velocity
Qout(idVvel) == F       ! v                  3D V-velocity
Qout(idu3dE) == F       ! u_eastward         3D U-eastward  at RHO-points
Qout(idv3dN) == F       ! v_northward        3D V-northward at RHO-points
Qout(idWvel) == F       ! w                  3D W-velocity
Qout(idOvel) == F       ! omega              omega vertical velocity
Qout(idUbar) == T       ! ubar               2D U-velocity
Qout(idVbar) == T       ! vbar               2D V-velocity
Qout(idu2dE) == T       ! ubar_eastward      2D U-eastward  at RHO-points
Qout(idv2dN) == T       ! vbar_northward     2D V-northward at RHO-points
Qout(idFsur) == T       ! zeta               free-surface
Qout(idBath) == T       ! bath               time-dependent bathymetry

Qout(idTvar) == F F     ! temp, salt         temperature and salinity

Qout(idUsur) == T       ! u_sur              surface U-velocity
Qout(idVsur) == T       ! v_sur              surface V-velocity
Qout(idUsuE) == T       ! u_sur_eastward     surface U-eastward  velocity
Qout(idVsuN) == T       ! v_sur_northward    surface V-northward velocity

Qout(idsurT) == T T     ! temp_sur, salt_sur surface temperature and salinity

Qout(idpthR) == F       ! z_rho              time-varying depths of RHO-points
Qout(idpthU) == F       ! z_u                time-varying depths of U-points
Qout(idpthV) == F       ! z_v                time-varying depths of V-points
Qout(idpthW) == F       ! z_w                time-varying depths of W-points

Qout(idUsms) == F       ! sustr              surface U-stress
Qout(idVsms) == F       ! svstr              surface V-stress
Qout(idUbms) == F       ! bustr              bottom U-stress
Qout(idVbms) == F       ! bvstr              bottom V-stress

Qout(idUbrs) == F       ! bustrc             bottom U-current stress
Qout(idVbrs) == F       ! bvstrc             bottom V-current stress
Qout(idUbws) == F       ! bustrw             bottom U-wave stress
Qout(idVbws) == F       ! bvstrw             bottom V-wave stress
Qout(idUbcs) == F       ! bustrcwmax         bottom max wave-current U-stress
Qout(idVbcs) == F       ! bvstrcwmax         bottom max wave-current V-stress

Qout(idUbot) == F       ! Ubot               bed wave orbital U-velocity
Qout(idVbot) == F       ! Vbot               bed wave orbital V-velocity
Qout(idUbur) == F       ! Ur                 bottom U-velocity above bed
Qout(idVbvr) == F       ! Vr                 bottom V-velocity above bed

Qout(idW2xx) == F       ! Sxx_bar            2D radiation stress, Sxx component
Qout(idW2xy) == F       ! Sxy_bar            2D radiation stress, Sxy component
Qout(idW2yy) == F       ! Syy_bar            2D radiation stress, Syy component
Qout(idU2rs) == F       ! Ubar_Rstress       2D radiation U-stress
Qout(idV2rs) == F       ! Vbar_Rstress       2D radiation V-stress
Qout(idU2Sd) == F       ! ubar_stokes        2D U-Stokes velocity
Qout(idV2Sd) == F       ! vbar_stokes        2D V-Stokes velocity

Qout(idW3xx) == F       ! Sxx                3D radiation stress, Sxx component
Qout(idW3xy) == F       ! Sxy                3D radiation stress, Sxy component
Qout(idW3yy) == F       ! Syy                3D radiation stress, Syy component
Qout(idW3zx) == F       ! Szx                3D radiation stress, Szx component
Qout(idW3zy) == F       ! Szy                3D radiation stress, Szy component
Qout(idU3rs) == F       ! u_Rstress          3D U-radiation stress
Qout(idV3rs) == F       ! v_Rstress          3D V-radiation stress
Qout(idU3Sd) == F       ! u_stokes           3D U-Stokes velocity
Qout(idV3Sd) == F       ! v_stokes           3D V-Stokes velocity

Qout(idWamp) == F       ! Hwave              wave height
Qout(idWlen) == F       ! Lwave              wave length
Qout(idWdir) == F       ! Dwave              wave direction-mean
Qout(idWdip) == F       ! Dwavep             wave direction-peak
Qout(idWptp) == F       ! Pwave_top          wave surface period
Qout(idWpbt) == F       ! Pwave_bot          wave bottom period
Qout(idWorb) == F       ! Ub_swan            wave bottom orbital velocity
Qout(idWdis) == F       ! Wave_dissip        wave dissipation

Qout(idPair) == F       ! Pair               surface air pressure
Qout(idTair) == F       ! Tair               surface air temperature
Qout(idUair) == F       ! Uair               surface U-wind component
Qout(idVair) == F       ! Vair               surface V-wind component

Qout(idTsur) == F F     ! shflux, ssflux     surface net heat and salt flux
Qout(idLhea) == F       ! latent             latent heat flux
Qout(idShea) == F       ! sensible           sensible heat flux
Qout(idLrad) == F       ! lwrad              longwave radiation flux
Qout(idSrad) == F       ! swrad              shortwave radiation flux
Qout(idEmPf) == F       ! EminusP            E-P flux
Qout(idevap) == F       ! evaporation        evaporation rate
Qout(idrain) == F       ! rain               precipitation rate

Qout(idDano) == F       ! rho                density anomaly
Qout(idVvis) == F       ! AKv                vertical viscosity
Qout(idTdif) == F       ! AKt                vertical T-diffusion
Qout(idSdif) == F       ! AKs                vertical Salinity diffusion
Qout(idHsbl) == F       ! Hsbl               depth of surface boundary layer
Qout(idHbbl) == F       ! Hbbl               depth of bottom boundary layer
Qout(idMtke) == F       ! tke                turbulent kinetic energy
Qout(idMtls) == F       ! gls                turbulent length scale

! Logical switches (TRUE/FALSE) to activate writing of extra inert passive
! tracers other than biological and sediment tracers into the QUICKSAVE
 output file. An inert passive tracer is one that it is only advected and
! diffused. Other processes are ignored. These tracers include, for example,
! dyes, pollutants, oil spills, etc. NPT values are expected. However, these
! switches can be activated using compact parameter specification.

 Qout(inert) == F       ! dye_01, ...        inert passive tracers
 Qout(Snert) == F       ! dye_01_sur, ..     surface inert passive tracers
! layer properties into HISTORY output file.  Currently, MBOTP properties
! are expected for the bottom boundary layer and/or sediment models:
!
! idBott( 1=isd50)   grain_diameter          mean grain diameter
! idBott( 2=idens)   grain_density           mean grain density
! idBott( 3=iwsed)   settling_vel            mean settling velocity
! idBott( 4=itauc)   erosion_stress          critical erosion stress
! idBott( 5=irlen)   ripple_length           ripple length
! idBott( 6=irhgt)   ripple_height           ripple height
! idBott( 7=ibwav)   bed_wave_amp            wave excursion amplitude
! idBott( 8=izdef)   Zo_def                  default bottom roughness
! idBott( 9=izapp)   Zo_app                  apparent bottom roughness
! idBott(10=izNik)   Zo_Nik                  Nikuradse bottom roughness
! idBott(11=izbio)   Zo_bio                  biological bottom roughness
! idBott(12=izbfm)   Zo_bedform              bed form bottom roughness
! idBott(13=izbld)   Zo_bedload              bed load bottom roughness
! idBott(14=izwbl)   Zo_wbl                  wave bottom roughness
! idBott(15=iactv)   active_layer_thickness  active layer thickness
! idBott(16=ishgt)   saltation               saltation height
! idBott(17=imaxD)   dep_net                 maximum inundation depth
! idBott(18=idnet)   net erosion + dep       Erosion or deposition
! idBott(19=idoff)   tau critical offset     dmix erodibility profile offset
! idBott(20=idslp)   tau critical slope      dmix or erodibility slope
! idBott(21=idtim)   erodibility time scale  erodibility profile restore time
! idBott(22=idbmx)   diffusivity db_max      Bed biodifusivity maximum
! idBott(23=idbmm)   diffusivity db_m        Bed biodifusivity minimum
! idBott(24=idbzs)   diffusivity db_zs       Bed biodifusivity zs
! idBott(25=idbzm)   diffusivity db_zm       Bed biodifusivity zm
! idBott(26=idbzp)   diffusivity db_zphi     Bed biodifusivity phi
! idBott(27=idprp)   cohesive behavior       cohesive behavior
!
!                                 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2
!               1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7
!
Hout(idBott) == T T T T T T T T T F F F F F T F T T F F F F F F F F F
! Logical switches (TRUE/FALSE) to activate writing of time-averaged
! fields into AVERAGE output file.

Aout(idUvel) == F       ! u                  3D U-velocity
Aout(idVvel) == F       ! v                  3D V-velocity
Aout(idu3dE) == F       ! u_eastward         3D U-eastward  at RHO-points
Aout(idv3dN) == F       ! v_northward        3D V-northward at RHO-points
Aout(idWvel) == F       ! w                  3D W-velocity
Aout(idOvel) == F       ! omega              omega vertical velocity
Aout(idUbar) == F       ! ubar               2D U-velocity
Aout(idVbar) == F       ! vbar               2D V-velocity
Aout(idu2dE) == F       ! ubar_eastward      2D U-eastward  at RHO-points
Aout(idv2dN) == F       ! vbar_northward     2D V-northward at RHO-points
Aout(idFsur) == F       ! zeta               free-surface
Aout(idBath) == F       ! bath               time-dependent bathymetry
Aout(idTvar) == F F     ! temp, salt         temperature and salinity
Aout(idUair) == F       ! Uwind              surface U-wind
Aout(idVair) == F       ! Vwind              surface V-wind
Aout(idUairE) == F      ! Uwind_eastward     surface U-wind
Aout(idVairN) == F      ! Vwind_northward    surface V-wind
Aout(idUsms) == F       ! sustr              surface U-stress
Aout(idVsms) == F       ! svstr              surface V-stress
Aout(idUbms) == F       ! bustr              bottom U-stress
Aout(idVbms) == F       ! bvstr              bottom V-stress
Aout(idUbrs) == F       ! bustrc             bottom U-current stress
Aout(idVbrs) == F       ! bvstrc             bottom V-current stress
Aout(idUbws) == F       ! bustrw             bottom U-wave stress
Aout(idVbws) == F       ! bvstrw             bottom V-wave stress
Aout(idUbcs) == F       ! bustrcwmax         bottom max wave-current U-stress
Aout(idVbcs) == F       ! bvstrcwmax         bottom max wave-current V-stress
Aout(idUVwc) == F       ! bstrcwmax          bottom max wave-current stress magnitude
Aout(idUbot) == F       ! Ubot               bed wave orbital U-velocity
Aout(idVbot) == F       ! Vbot               bed wave orbital V-velocity
Aout(idUbur) == F       ! Ur                 bottom U-velocity above bed
Aout(idVbvr) == F       ! Vr                 bottom V-velocity above bed
Aout(idW2xx) == F       ! Sxx_bar            WEC_Mellor 2D Sxx radiation stress
Aout(idW2xy) == F       ! Sxy_bar            WEC_Mellor 2D Sxy radiation stress
Aout(idW2yy) == F       ! Syy_bar            WEC_Mellor 2D Syy radiation stress
Aout(idW3xx) == F       ! Sxx                WEC_Mellor 3D Sxx radiation stress
Aout(idW3xy) == F       ! Sxy                WEC_Mellor 3D Sxy radiation stress
Aout(idW3yy) == F       ! Syy                WEC_Mellor 3D Syy radiation stress
Aout(idW3zx) == F       ! Szx                WEC_Mellor 3D Szx radiation stress
Aout(idW3zy) == F       ! Szy                WEC_Mellor 3D Szy radiation stress
Aout(idWztw) == F       ! zetaw              WEC_VF quasi-static sea level adjustment
Aout(idWqsp) == F       ! qsp                WEC_VF quasi-static pressure
Aout(idWbeh) == F       ! bh                 WEC_VF Bernoulli head
Aout(idU2rs) == F       ! ubar_Wecstress     WEC 2D U-stress
Aout(idV2rs) == F       ! vbar_Wecstress     WEC 2D V-stress
Aout(idU3rs) == F       ! u_Rstress          WEC 3D U-stress
Aout(idV3rs) == F       ! v_Rstress          WEC 3D V-stress
Aout(idU2Sd) == F       ! ubar_stokes        2D U-Stokes velocity
Aout(idV2Sd) == F       ! vbar_stokes        2D V-Stokes velocity
Aout(idU3Sd) == F       ! u_stokes           3D U-Stokes velocity
Aout(idV3Sd) == F       ! v_stokes           3D V-Stokes velocity
Aout(idW3Sd) == F       ! omega_stokes       3D Omega-Stokes velocity
Aout(idW3St) == F       ! w_stokes           3D W-Stokes velocity
Aout(idWamp) == F       ! Hwave              wave height
Aout(idWlen) == F       ! Lwave              wave length-mean
Aout(idWlep) == F       ! Lwavep             wave length-peak
Aout(idWdir) == F       ! Dwave              wave direction
Aout(idWptp) == F       ! Pwave_top          wave surface period
Aout(idWpbt) == F       ! Pwave_bot          wave bottom period
Aout(idWorb) == F       ! Uwave_rms          wave bottom orbital velocity
Aout(idWbrk) == F       ! Wave_break         wave breaking (percent)
Aout(idUwav) == F       ! uWave              wave-depth avgeraged U-velocity
Aout(idVwav) == F       ! vWave              wave-depth avgeraged V-velocity
Aout(idWdif) == F       ! Dissip_fric        wave dissipation due to bottom friction
Aout(idWdib) == F       ! Dissip_break       wave dissipation due to breaking
Aout(idWdiw) == F       ! Dissip_wcap        wave dissipation due to white capping
Aout(idWdis) == F       ! Dissip_roller      wave roller dissipation
Aout(idWrol) == F       ! rollA              wave roller action density
Aout(idRunoff) == F     ! Runoff             surface runoff from land
Aout(idPair) == F       ! Pair               surface air pressure
Aout(idUair) == F       ! Uair               surface U-wind component
Aout(idVair) == F       ! Vair               surface V-wind component												 
Aout(idTsur) == F F     ! shflux, ssflux     surface net heat and salt flux
Aout(idLhea) == F       ! latent             latent heat flux
Aout(idShea) == F       ! sensible           sensible heat flux
Aout(idLrad) == F       ! lwrad              longwave radiation flux
Aout(idSrad) == F       ! swrad              shortwave radiation flux
Aout(idevap) == F       ! evaporation        evaporation rate
Aout(idrain) == F       ! rain               precipitation rate
Aout(idDano) == F       ! rho                density anomaly
Aout(idVvis) == F       ! AKv                vertical viscosity
Aout(idTdif) == F       ! AKt                vertical T-diffusion
Aout(idSdif) == T       ! AKs                vertical Salinity diffusion
Aout(idHsbl) == F       ! Hsbl               depth of surface boundary layer
Aout(idHbbl) == F       ! Hbbl               depth of bottom boundary layer
Aout(id2dRV) == F       ! pvorticity_bar     2D relative vorticity
Aout(id3dRV) == F       ! pvorticity         3D relative vorticity
Aout(id2dPV) == F       ! rvorticity_bar     2D potential vorticity
Aout(id3dPV) == F       ! rvorticity         3D potential vorticity
Aout(idu3dD) == F       ! u_detided          detided 3D U-velocity
Aout(idv3dD) == F       ! v_detided          detided 3D V-velocity
Aout(idu2dD) == F       ! ubar_detided       detided 2D U-velocity
Aout(idv2dD) == F       ! vbar_detided       detided 2D V-velocity
Aout(idFsuD) == F       ! zeta_detided       detided free-surface
Aout(idTrcD) == F F     ! temp_detided, ...  detided temperature and salinity
Aout(idHUav) == F       ! Huon               u-volume flux, Huon
Aout(idHVav) == F       ! Hvom               v-volume flux, Hvom
Aout(idUUav) == F       ! uu                 quadratic <u*u> term
Aout(idUVav) == F       ! uv                 quadratic <u*v> term
Aout(idVVav) == F       ! vv                 quadratic <v*v> term
Aout(idU2av) == F       ! ubar2              quadratic <ubar*ubar> term
Aout(idV2av) == F       ! vbar2              quadratic <vbar*vbar> term
Aout(idZZav) == F       ! zeta2              quadratic <zeta*zeta> term
Aout(idTTav) == F F     ! temp_2, ...        quadratic <t*t> tracer terms
Aout(idUTav) == F F     ! u_temp, ...        quadratic <u*t> tracer terms
Aout(idVTav) == F F     ! v_temp, ...        quadratic <v*t> tracer terms
Aout(iHUTav) == F F     ! Huon_temp, ...     tracer volume flux, <Huon*t>
Aout(iHVTav) == F F     ! Hvom_temp, ...     tracer volume flux, <Hvom*t>

! Logical switches (TRUE/FALSE) to activate writing of extra inert passive
! tracers other than biological and sediment tracers into the AVERAGE file.

 Aout(inert) == F       ! dye_01, ...        inert passive tracers

! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
! 2D momentum (ubar,vbar) diagnostic terms into DIAGNOSTIC output file.

Dout(M2rate) == F       ! ubar_accel, ...    acceleration
Dout(M2pgrd) == F       ! ubar_prsgrd, ...   pressure gradient
Dout(M2fcor) == F       ! ubar_cor, ...      Coriolis force
Dout(M2hadv) == F       ! ubar_hadv, ...     horizontal total advection
Dout(M2xadv) == F       ! ubar_xadv, ...     horizontal XI-advection
Dout(M2yadv) == F       ! ubar_yadv, ...     horizontal ETA-advection
Dout(M2hrad) == F       ! ubar_hrad, ...     horizontal total wec_mellor radiation stress
Dout(M2hvis) == F       ! ubar_hvisc, ...    horizontal total viscosity
Dout(M2xvis) == F       ! ubar_xvisc, ...    horizontal XI-viscosity
Dout(M2yvis) == F       ! ubar_yvisc, ...    horizontal ETA-viscosity
Dout(M2sstr) == F       ! ubar_sstr, ...     surface stress
Dout(M2bstr) == F       ! ubar_bstr, ...     bottom stress														  
Dout(M2fveg) == T 							
Dout(M2hjvf) == F       ! 2D wec_vf horizontal J vortex force
Dout(M2kvrf) == F       ! 2D wec_vf K vortex force
Dout(M2fsco) == F       ! 2D wec_vf coriolis-stokes
Dout(M2bstm) == F       ! 2D wec_vf bottom streaming
Dout(M2sstm) == F       ! 2D wec_vf surface streaming
Dout(M2wrol) == F       ! 2D wec_vf wave roller accel
Dout(M2wbrk) == F       ! 2D wec_vf wave breaking
Dout(M2zeta) == F       ! 2D wec_vf Eulerian sea level adjustment
Dout(M2zetw) == F       ! 2D wec_vf quasi-static sea level adjustment
Dout(M2zqsp) == F       ! 2D wec_vf quasi-static pressure
Dout(M2zbeh) == F       ! 2D wec_vf Bernoulli head

! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
! 3D momentum (u,v) diagnostic terms into DIAGNOSTIC output file.

Dout(M3rate) == F       ! u_accel, ...       acceleration
Dout(M3pgrd) == F       ! u_prsgrd, ...      pressure gradient
Dout(M3fcor) == F       ! u_cor, ...         Coriolis force
Dout(M3hadv) == F       ! u_hadv, ...        horizontal total advection
Dout(M3xadv) == F       ! u_xadv, ...        horizontal XI-advection
Dout(M3yadv) == F       ! u_yadv, ...        horizontal ETA-advection
Dout(M3vadv) == F       ! u_vadv, ...        vertical advection
Dout(M3hrad) == F       ! u_hrad, ...        horizontal total wec_mellor radiation stress
Dout(M3vrad) == F       ! v_hrad, ...        vertical total wec_mellor radiation stress
Dout(M3hvis) == F       ! u_hvisc, ...       horizontal total viscosity
Dout(M3xvis) == F       ! u_xvisc, ...       horizontal XI-viscosity
Dout(M3yvis) == F       ! u_yvisc, ...       horizontal ETA-viscosity
Dout(M3vvis) == F       ! u_vvisc, ...       vertical viscosity																  
Dout(M3fveg) == T 		! u_fveg, ...        vegetation drag force																		   
Dout(M3vjvf) == F       ! 3D wec_vf vertical J vortex force
Dout(M3hjvf) == F       ! 3D wec_vf horizontal J vortex force
Dout(M3kvrf) == F       ! 3D wec_vf K vortex force
Dout(M3fsco) == F       ! 3D wec_vf coriolis-stokes
Dout(M3bstm) == F       ! 3D wec_vf bottom streaming
Dout(M3sstm) == F       ! 3D wec_vf surface streaming
Dout(M3wrol) == F       ! 3D wec_vf wave roller accel
Dout(M3wbrk) == F       ! 3D wec_vf wave breaking

! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
! active (temperature and salinity) and passive (inert) tracer diagnostic
! terms into DIAGNOSTIC output file: [1:NAT+NPT,Ngrids].

Dout(iTrate) == F F     ! temp_rate, ...     time rate of change
Dout(iThadv) == F F     ! temp_hadv, ...     horizontal total advection
Dout(iTxadv) == F F     ! temp_xadv, ...     horizontal XI-advection
Dout(iTyadv) == F F     ! temp_yadv, ...     horizontal ETA-advection
Dout(iTvadv) == F F     ! temp_vadv, ...     vertical advection
Dout(iThdif) == F F     ! temp_hdiff, ...    horizontal total diffusion
Dout(iTxdif) == F F     ! temp_xdiff, ...    horizontal XI-diffusion
Dout(iTydif) == F F     ! temp_ydiff, ...    horizontal ETA-diffusion
Dout(iTsdif) == F F     ! temp_sdiff, ...    horizontal S-diffusion
Dout(iTvdif) == F F     ! temp_vdiff, ...    vertical diffusion

! Generic User parameters, [1:NUSER].

       NUSER =  0
        USER =  0.d0
! NetCDF-4/HDF5 compression parameters for output files.

  NC_SHUFFLE =  1                 ! if non-zero, turn on shuffle filter
  NC_DEFLATE =  1                 ! if non-zero, turn on deflate filter
   NC_DLEVEL =  1                 ! deflate level [0-9]

! Input NetCDF file names, [1:Ngrids].

     GRDNAME == Projects/Doksuri/ROMS_COAWST_grd1.nc

     ININAME == Projects/Doksuri/Doksuri_ini.nc 
     ITLNAME == ocean_itl.nc
     IRPNAME == ocean_irp.nc
     IADNAME == ocean_iad.nc
     FWDNAME == ocean_fwd.nc
     ADSNAME == ocean_ads.nc
     IWININAME == inwave_ini.nc
     IWSWNNAME == point1.spc2d
! Input adjoint forcing NetCDF filenames for computing observations
! impacts during the analysis-forecast cycle. If the forecast error
! metric is defined in state-space, then FOInameA and FOInameB should
! be regular adjoint forcing files just like ADSname. If the forecast
! error metric is defined in observation space (OBS_SPACE is activated)
! then the forecast is initialized OIFnameA and OIFnameB (specified in
! s4dvar.in input script) will have the structure of a 4D-Var observation
! file.

    FOInameA == dogbone_foi_a.nc
    FOInameB == dogbone_foi_b.nc

! Input NetCDF filenames for the forecasts initialized from the analysis
! of the current 4D-Var cycle (FCTnameA) and initialized from the analysis
! of the previous 4D-Var cycle (FCTnameB).

    FCTnameA == dogbone_fct_a.nc
    FCTnameB == dogbone_fct_b.nc

! Nesting grids connectivity data: contact points information. This
! NetCDF file is special and complex. It is currently generated using
! the script "matlab/grid/contact.m" from the Matlab repository.

     NGCNAME =  Projects/Doksuri/Doksuri_roms_contact.nc

! Input lateral boundary conditions and climatology file names.  The
! USER has the option to split input data time records into several
! NetCDF files (see prologue instructions above). If so, use a single
! line per entry with a vertical bar (|) symbol after each entry,
! except the last one.

     NCLMFILES == 1                        ! number of climate files
     CLMNAME == Projects/Doksuri/Doksuri_clm.nc
!	 CLMNAME == Projects/Sarika/coawst_clm_20161016.nc \ 
!                Projects/Sarika/coawst_clm_20161017.nc 
                

     NBCFILES == 1                         ! number of boundary files
     BRYNAME == Projects/Doksuri/Doksuri_bdy.nc
!     BRYNAME == Projects/Sarika4/coawst_bdy_20161016.nc \
!                Projects/Sarika4/coawst_bdy_20161017.nc 

! Input climatology nudging coefficients file name.

     NUDNAME == ocean_nud.nc

! Input Sources/Sinks forcing (like river runoff) file name.

     SSFNAME == ocean_rivers.nc

! Input tidal forcing file name.

    TIDENAME == Projects/Doksuri/tide_forc_Doksuri.nc
   
! Input forcing NetCDF file name(s).  The USER has the option to enter
! several file names for each nested grid.  For example, the USER may
! have different files for wind products, heat fluxes, tides, etc.
! The model will scan the file list and will read the needed data from
! the first file in the list containing the forcing field. Therefore,
! the order of the file names is very important. If using multiple forcing
! files per grid, first enter all the file names for grid 1, then grid 2,
! and so on.  It is also possible to split input data time records into
! several NetCDF files (see prologue instructions above). Use a single line
! per entry with a continuation (\) or vertical bar (|) symbol after each
! entry, except the last one.
! Input forcing NetCDF file name(s).
!
! The USER has the option to enter several sets of file names for each
! nested grid. For example, the USER may have different data for the
! wind products, heat fluxes, etc. Alternatively, if the all the forcing
! files are the same for nesting and the data is in its native resolution,
! we could enter only one set of files names and ROMS will replicate those
! files internally to the remaining grids using the plural KEYWORD protocol.
!
! The model will scan the files and will read the needed data from the first
! file in the list containing the forcing field. Therefore, the order of the
! filenames is critical. If using multiple forcing files per grid, first
! enter all the file names for grid one followed by two, and so on.  It is
! also possible to split input data time records into several NetCDF files
! (see Prolog instructions above). Use a single line per entry with a
! continuation (\) or a vertical bar (|) symbol after each entry, except
! the last one.

     NFFILES == 1
     FRCNAME == Projects/Doksuri/tide_forc_Doksuri.nc
                

! Output NetCDF file names, [1:Ngrids].

     GSTNAME == ocean_gst.nc
     RSTNAME == ocean_rst.nc 							
     HISNAME == ocean_hist.nc 		
     TLMNAME == ocean_tlm.nc
     TLFNAME == ocean_tlf.nc
     ADJNAME == ocean_adj.nc
     AVGNAME == ocean_avg.nc 
     DIANAME == ocean_dia.nc
     STANAME == ocean_sta.nc
     FLTNAME == ocean_flt.nc

! Input ASCII parameter filenames.

     APARNAM =  ROMS/External/s4dvar.in
     SPOSNAM =  ROMS/External/stations.in
     FPOSNAM =  ROMS/External/floats.in
     BPARNAM =  ROMS/External/bioFasham.in
     SPARNAM =  ROMS/External/sediment.in 
     USRNAME =  ROMS/External/MyFile.dat
Please help me to resolve this problem!
Any suggestions will be appreciated.

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

Re: forrtl: severe (174): SIGSEGV, segmentation fault occurred

#2 Unread post by jcwarner »

please post this as an issue here:
https://github.com/jcwarner-usgs/COAWST/issues
thanks,
j

Post Reply