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

[lixinwen000]

I have errors when running the upwelling case, the version is arango
The model can run successfully with oceanG and oceanO.


The following error occurs with the executable oceanS:

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00000288)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed

0 0 00:00:00 0.000000E+00 6.585677E+02 6.585677E+02 3.884376E+11
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: ocean_avg.nc
DEF_DIAGS - creating diagnostics file, Grid 01: ocean_dia.nc
1 0 00:05:00 6.479416E-13 6.585677E+02 6.585677E+02 3.884376E+11
(01,01,01) 1.068554E-09 1.950202E-06 0.000000E+00 9.216771E-06
2 0 00:10:00 9.311375E-12 6.585677E+02 6.585677E+02 3.884376E+11
(01,02,08) 6.735275E-08 1.793582E-06 1.433645E-05 3.661905E-05
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanM 00000000004B5F15 Unknown Unknown Unknown
oceanM 00000000004AAB23 Unknown Unknown Unknown
oceanM 000000000048BD57 Unknown Unknown Unknown
oceanM 0000000000404BEC Unknown Unknown Unknown
oceanM 0000000000404A0B Unknown Unknown Unknown
oceanM 000000000040485C Unknown Unknown Unknown
libc.so.6 0000003C2FC1D8B4 Unknown Unknown Unknown
oceanM 0000000000404769 Unknown Unknown Unknown
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanM 00000000004B5F15 Unknown Unknown Unknown
oceanM 00000000004AAB23 Unknown Unknown Unknown
oceanM 000000000048BD57 Unknown Unknown Unknown
oceanM 0000000000404BEC Unknown Unknown Unknown
oceanM 0000000000404A0B Unknown Unknown Unknown
oceanM 000000000040485C Unknown Unknown Unknown
libc.so.6 0000003C2FC1D8B4 Unknown Unknown Unknown
oceanM 0000000000404769 Unknown Unknown Unknown
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanM 00000000004B5F15 Unknown Unknown Unknown
oceanM 00000000004AAB23 Unknown Unknown Unknown
oceanM 000000000048BD57 Unknown Unknown Unknown
oceanM 0000000000404BEC Unknown Unknown Unknown
oceanM 0000000000404A0B Unknown Unknown Unknown
oceanM 000000000040485C Unknown Unknown Unknown
libc.so.6 0000003C2FC1D8B4 Unknown Unknown Unknown
oceanM 0000000000404769 Unknown Unknown Unknown
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanM 00000000004B5F15 Unknown Unknown Unknown
oceanM 00000000004AAB23 Unknown Unknown Unknown
oceanM 000000000048BD57 Unknown Unknown Unknown
oceanM 0000000000404BEC Unknown Unknown Unknown
oceanM 0000000000404A0B Unknown Unknown Unknown
oceanM 000000000040485C Unknown Unknown Unknown
libc.so.6 0000003C2FC1D8B4 Unknown Unknown Unknown
oceanM 0000000000404769 Unknown Unknown Unknown
rank 2 in job 9 hksbs-s10.com_37297 caused collective abort of all ranks
exit status of rank 2: return code 174
rank 1 in job 9 hksbs-s10.com_37297 caused collective abort of all ranks
exit status of rank 1: return code 174
rank 0 in job 9 hksbs-s10.com_37297 caused collective abort of all ranks
exit status of rank 0: return code 174


The following error occurs with the executable oceanM:

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


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

0 0 00:00:00 0.000000E+00 6.585677E+02 6.585677E+02 3.884376E+11
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: ocean_avg.nc
DEF_DIAGS - creating diagnostics file, Grid 01: ocean_dia.nc
1 0 00:05:00 6.479416E-13 6.585677E+02 6.585677E+02 3.884376E+11
(01,01,01) 1.068554E-09 1.950202E-06 0.000000E+00 9.216771E-06
2 0 00:10:00 9.311373E-12 6.585677E+02 6.585677E+02 3.884376E+11
(01,02,08) 6.735275E-08 1.793580E-06 1.433645E-05 3.661904E-05
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
oceanS 000000000049E275 Unknown Unknown Unknown
oceanS 0000000000492D23 Unknown Unknown Unknown
oceanS 0000000000474667 Unknown Unknown Unknown
oceanS 0000000000403B3C Unknown Unknown Unknown
oceanS 000000000040395D Unknown Unknown Unknown
oceanS 000000000040386C Unknown Unknown Unknown
libc.so.6 0000003C2FC1D8B4 Unknown Unknown Unknown
oceanS 0000000000403779 Unknown Unknown Unknown

[jcwarner]

it looks like it took 2 time steps

1 0 00:05:00 6.479416E-13 6.585677E+02 6.585677E+02 3.884376E+11
(01,01,01) 1.068554E-09 1.950202E-06 0.000000E+00 9.216771E-06
2 0 00:10:00 9.311375E-12 6.585677E+02 6.585677E+02 3.884376E+11
(01,02,08) 6.735275E-08 1.793582E-06 1.433645E-05 3.661905E-05
forrtl: severe (174): SIGSEGV, segmentation fault occurred

so save the history every time step and look at the output his file. Where is it having trouble?
grid location 01,02 looks like a boundary problem.

-j

[kate]

I would also recompile with bounds checking turned on (often the default with USE_DEBUG). One way to get a seg fault is through an array bounds violation.

[matias.dinapoli]

Hi, I have the same issue, but I'm with oceanS and setting a grid netcdf file. When I run it:

log file:

Code: Select all

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


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

      0     0 00:00:00  0.000000E+00  0.000000E+00  0.000000E+00  4.584548E+13
              (00,000)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00

and in the console:

Code: Select all

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
oceanS             000000000044BE3D  Unknown               Unknown  Unknown
oceanS             00000000004423B2  Unknown               Unknown  Unknown
oceanS             000000000040375A  Unknown               Unknown  Unknown
oceanS             00000000004035A2  Unknown               Unknown  Unknown
oceanS             000000000040342C  Unknown               Unknown  Unknown
libc.so.6          00007FD8F0AA3D5D  Unknown               Unknown  Unknown
oceanS             0000000000403329  Unknown               Unknown  Unknown

I've used the option ulimit -c unlimited, but this go on. Is important the order of the parameters in the grid.nc? I made this grid file with easygrid.m (contains the bath and mask points)

Matias.

[kate]

This is what I know of seg faults. We really don't know what caused yours from what you posted. The model read the grid file, so it's probably not that (some grid file problems can cause another sort of blowing up).

[mathieu]

As Kate said, you have to use USE_DEBUG and maybe edit the corresponding file in Compiler directory for changing the debug compiler flags. The "Unknown" will then be replaced by a precise statement of the nature of the error.

In truth, it is more likely that there is an error in your configuration and/or blow-up than a true error in the ROMS model. Unfortunately, it is hard if not impossible to make detailed error trackback when a problem happens.

[matias.dinapoli]

I made the build.bahs with USE_DEBUG and I got this

Code: Select all

forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
oceanG             00000000004BDC68  set_vbc_mod_mp_se         180  set_vbc.f90
oceanG             00000000004B9C8C  set_vbc_mod_mp_se          63  set_vbc.f90
oceanG             0000000000497F48  main2d_                   122  main2d.f90
oceanG             0000000000403FA1  ocean_control_mod         151  ocean_control.f90
oceanG             0000000000403645  MAIN__                     86  master.f90
oceanG             000000000040344C  Unknown               Unknown  Unknown
libc.so.6          00007FEE7A6C6D5D  Unknown               Unknown  Unknown
oceanG             0000000000403349  Unknown               Unknown  Unknown

So I have somethings specific. And Kate your advice helped me to figere out some mistakes in my grid file.

[kate]

Code: Select all

Image              PC                Routine            Line        Source
oceanG             00000000004BDC68  set_vbc_mod_mp_se         180  set_vbc.f90

Now the next step is to look at line 180 of your $SCRATCH_DIR/set_vbc.f90 and see what's there.

[matias.dinapoli]

I solved my issue, it was a mix of grid setting and the forcing.

Thanks!!

[nmichelet]

Hi,
I have the same kind of problem, could you explain me how you solve it.
I am currently trying to make a two-way nest simulation, so I followed the example of the Lake Jersey AB example. The program don't even start and it gives me the log file I attached.
Do I miss something in the example?
Thanks for your time

[kate]

Did you read the posts higher in the thread? Try USE_DEBUG and see if it gives you useful line numbers at least. Then you can look at those files to see what is at the line number.

[nmichelet]

Deer Kate,
I read the entire post but it is not exactly the same problem.
I tried the USE_DEBUG and here the information it gives me

Code: Select all

 Model Input Parameters:  ROMS/TOMS version 3.5
                          Wednesday - December 2, 2015 -  3:32:04 PM
 -----------------------------------------------------------------------------
forrtl: severe (408): fort: (2): Subscript #1 of the array FILES has value 1 which is greater than the upper bound of 0

Image              PC                Routine            Line        Source
oceanG             0000000000F8DD4D  Unknown               Unknown  Unknown
oceanG             0000000000F8C855  Unknown               Unknown  Unknown
oceanG             0000000000F33A80  Unknown               Unknown  Unknown
oceanG             0000000000EF1FBF  Unknown               Unknown  Unknown
oceanG             0000000000EF23C2  Unknown               Unknown  Unknown
oceanG             0000000000B204F2  load_s1d_                2508  inp_par.f90
oceanG             0000000000AFC9A0  read_phypar_             1390  inp_par.f90
oceanG             0000000000AA521A  inp_par_                   87  inp_par.f90
oceanG             000000000040492C  ocean_control_mod          84  ocean_control.f90
oceanG             00000000004045E9  MAIN__                     95  master.f90
oceanG             00000000004044CC  Unknown               Unknown  Unknown
libc.so.6          00002B914712DBC6  Unknown               Unknown  Unknown
oceanG             00000000004043C9  Unknown               Unknown  Unknown

I checked what's in line 2508 of inp_par.f90

S(ng)%name=TRIM(S(ng)%files(1)) ! load first file

If you have some advice to solve it it would be great
Thanks

[kate]

How about what is here?

Code: Select all

oceanG             0000000000AFC9A0  read_phypar_             1390  inp_par.f90

Inp_par is reading the ocean.in file, presumably the names of the forcing files. What do you have in that section of the ocean.in file? Show us all of:

Code: Select all

     NFFILES == 1                          ! number of forcing files

     FRCNAME == ocean_frc.nc               ! forcing file 1, grid 1

[zhzhlin]

I have got similar errors a few weeks ago and it was due to stackspace exhaustion. I added -heap-arrays to the fortran compiler and it solved the issue. I'm not sure if this is your case, but if you don't see anything wrong with the code, it may be something worth trying. There is a post online explaining this error and I found it quite helpful: https://software.intel.com/en-us/articl ... us-errors/
Good luck! I know how much pain it is to see this error.

[nmichelet]

Hi,
I have checked and the -heap-array option is already activated
I have looked at the lines 1390 and 87 in the inp_par.f90, but I don't really understand what it means

line 1390 : Npts=load_s1d(Nval, Cval, line, label, igrid, Nfiles, ITL)

line 87 : CALL read_PhyPar (model, inp, out, Lwrite)

Finally in the ocean.in I have

Code: Select all

NCFILES == 2                           ! number of unique forcing files

FRCNAMES == ./frc_parent.nc \          ! forcing file 1, grid 1
            ./frc_child.nc             ! forcing file 1, grid 2

[kate]

Hmm, it's not the forcing file:

Code: Select all

            CASE ('ITLNAME')
              label='ITL - tangent linear model initial conditions'
              Npts=load_s1d(Nval, Cval, Cdim, line, label, igrid,       &
     &                      Nfiles, ITL)

Do you have:

Code: Select all

ITLNAME == ocean_itl.nc

in your ocean.in?

As for the forcing file, try removing the comments after the "\" and after the filenames. ROMS will ignore any line starting with "!" but an ending "!" with comment is more problematic.

[nmichelet]

I have this line in my ocean.in but no file named ocean_itl.nc, what kind of file is this?
I also have tried to remove the comment line but same thing, no change

[kate]

Maybe show us that whole chunk of your ocean.in file. Did you ask it to read an initial conditions file? Do you have a comment on the ININAME line? It's reading that just before the ITLNAME.

[lindandan]

Maybe show us that whole chunk of your ocean.in file. Did you ask it to read an initial conditions file? Do you have a comment on the ININAME line? It's reading that just before the ITLNAME.

Hello Kate
I may have a problem when the model read the focring/bulk file. But I can't find the point,could you help? Thanks very much.I will run the model one year form 2009.01.01.00:00:00 to 2009.12.31.00:00:00
the output used oceanG as follows:

Code: Select all

Model Input Parameters:  ROMS/TOMS version 3.7  
                          Wednesday - June 29, 2016 -  4:57:58 PM
 -----------------------------------------------------------------------------

  Internal tidal System, 1/12 degree resolution

 Operating system : Linux
 CPU/hardware     : x86_64
 Compiler system  : ifort
 Compiler command : /usr/local/mpi/bin/mpif90
 Compiler flags   : -heap-arrays -fp-model precise -I/usr/local/mpi/include -g -check bounds -traceback -check uninit -warn interfaces,nouncalled -gen-interfaces -free -free -free

 Input Script  : 

 SVN Root URL  : https://www.myroms.org/svn/src/trunk
 SVN Revision  : 783M

 Local Root    : /vol/home/student/lshw/roms/trunk
 Header Dir    : /vol/home/student/lshw/roms/MYPROJECT/rSSC12
 Header file   : ssc12.h
 Analytical Dir: /vol/home/student/lshw/roms/MYPROJECT/Functionals

 Resolution, Grid 01: 0275x0275x032,  Parallel Nodes:  16,  Tiling: 004x004


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

      87600  ntimes            Number of timesteps for 3-D equations.
    360.000  dt                Timestep size (s) for 3-D equations.
         72  ndtfast           Number of timesteps for 2-D equations between
                                 each 3D timestep.
          1  ERstr             Starting ensemble/perturbation run number.
          1  ERend             Ending ensemble/perturbation run number.
          0  nrrec             Number of restart records to read from disk.
          T  LcycleRST         Switch to recycle time-records in restart file.
         90  nRST              Number of timesteps between the writing of data
                                 into restart fields.
          1  ninfo             Number of timesteps between print of information
                                 to standard output.
          T  ldefout           Switch to create a new output NetCDF file(s).
       1200  nHIS              Number of timesteps between the writing fields
                                 into history file.
      14400  ndefHIS           Number of timesteps between creation of new
                                 history files.
     864000  ntsAVG            Starting timestep for the accumulation of output
                                 time-averaged data.
     864000  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
     864000  ndefAVG           Number of timesteps between creation of new
                                 time-averaged file.
     864000  ntsDIA            Starting timestep for the accumulation of output
                                 time-averaged diagnostics data.
     864000  nDIA              Number of timesteps between the writing of
                                 time-averaged data into diagnostics file.
     864000  ndefDIA           Number of timesteps between creation of new
                                 diagnostic file.
 2.5000E+01  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 2.5000E+01  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 2.5000E+01  nl_visc2          NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for momentum.
          T  LuvSponge         Turning ON  sponge on horizontal momentum.
          T  LtracerSponge(01) Turning ON  sponge on tracer 01: temp
          T  LtracerSponge(02) Turning ON  sponge on tracer 02: salt
 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).
 2.5000E-03  rdrg2             Quadratic bottom drag coefficient.
 2.0000E-02  Zob               Bottom roughness (m).
 1.0000E+01  blk_ZQ            Height (m) of surface air humidity measurement.
 1.0000E+01  blk_ZT            Height (m) of surface air temperature measurement.
 1.0000E+01  blk_ZW            Height (m) of surface winds measurement.
          1  lmd_Jwt           Jerlov water type.
          2  Vtransform        S-coordinate transformation equation.
          4  Vstretching       S-coordinate stretching function.
 6.0000E+00  theta_s           S-coordinate surface control parameter.
 5.0000E-01  theta_b           S-coordinate bottom  control parameter.
     15.000  Tcline            S-coordinate surface/bottom layer width (m) used
                                 in vertical coordinate stretching.
   1025.000  rho0              Mean density (kg/m3) for Boussinesq approximation.
      0.000  dstart            Time-stamp assigned to model initialization (days).
      0.000  tide_start        Reference time origin for tidal forcing (days).
20090101.00  time_ref          Reference time for units attribute (yyyymmdd.dd)
 1.2000E+02  Tnudg(01)         Nudging/relaxation time scale (days)
                                 for tracer 01: temp
 1.2000E+02  Tnudg(02)         Nudging/relaxation time scale (days)
                                 for tracer 02: salt
 1.2000E+02  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 1.2000E+02  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 1.2000E+02  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 1.2000E+02  obcfac            Factor between passive and active
                                 open boundary conditions.
          F  VolCons(1)        NLM western  edge boundary volume conservation.
          F  VolCons(2)        NLM southern edge boundary volume conservation.
          F  VolCons(3)        NLM eastern  edge boundary volume conservation.
          F  VolCons(4)        NLM northern edge boundary volume conservation.
      5.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).
          F  LuvSrc            Turning OFF momentum point Sources/Sinks.
          F  LwSrc             Turning OFF volume influx point Sources/Sinks.
          F  LtracerSrc(01)    Turning OFF point Sources/Sinks on tracer 01: temp
          F  LtracerSrc(02)    Turning OFF point Sources/Sinks on tracer 02: salt
          F  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          F  Lm3CLM            Turning OFF processing of 3D momentum climatology.
          F  LtracerCLM(01)    Turning OFF processing of climatology tracer 01: temp
          F  LtracerCLM(02)    Turning OFF processing of climatology tracer 02: salt
          F  LnudgeM2CLM       Turning OFF nudging of 2D momentum climatology.
          F  LnudgeM3CLM       Turning OFF nudging of 3D momentum climatology.
          F  LnudgeTCLM(01)    Turning OFF nudging of climatology tracer 01: temp
          F  LnudgeTCLM(02)    Turning OFF nudging of climatology tracer 02: salt
          T  Hout(idFsur)      Write out free-surface.
          T  Hout(idUbar)      Write out 2D U-momentum component.
          T  Hout(idVbar)      Write out 2D V-momentum component.
          T  Hout(idUvel)      Write out 3D U-momentum component.
          T  Hout(idVvel)      Write out 3D V-momentum component.
          T  Hout(idWvel)      Write out W-momentum component.
          T  Hout(idOvel)      Write out omega vertical velocity.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt
          T  Hout(idUsms)      Write out surface U-momentum stress.
          T  Hout(idVsms)      Write out surface V-momentum stress.
          T  Hout(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(idEmPf)      Write out E-P flux.
          T  Hout(idevap)      Write out evaporation rate.
          T  Hout(idrain)      Write out rain rate.
          T  Hout(idDano)      Write out density anomaly.
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idTdif)      Write out vertical diffusion: AKt(itemp).
          T  Hout(idSdif)      Write out vertical diffusion: AKt(isalt).
          T  Hout(idMtke)      Write out turbulent kinetic energy.
          T  Hout(idMtls)      Write out turbulent generic length-scale.

          T  Aout(idFsur)      Write out averaged free-surface.
          T  Aout(idUbar)      Write out averaged 2D U-momentum component.
          T  Aout(idVbar)      Write out averaged 2D V-momentum component.
          T  Aout(idUvel)      Write out averaged 3D U-momentum component.
          T  Aout(idVvel)      Write out averaged 3D V-momentum component.
          T  Aout(idWvel)      Write out averaged W-momentum component.
          T  Aout(idOvel)      Write out averaged omega vertical velocity.
          T  Aout(idTvar)      Write out averaged tracer 01: temp
          T  Aout(idTvar)      Write out averaged tracer 02: salt

          T  Dout(M2rate)      Write out 2D momentum acceleration.
          T  Dout(M2pgrd)      Write out 2D momentum pressure gradient.

          T  Dout(M3rate)      Write out 3D momentum acceleration.
          T  Dout(M3pgrd)      Write out 3D momentum pressure gradient.


 Output/Input Files:

             Output Restart File:  ssc12_rst.nc
        Prefix for History Files:  ssc12_his
       Prefix for Averages Files:  ssc12_avg
    Prefix for Diagnostics Files:  ssc12_dia
                 Input Grid File:  ./Data/ssc12_grd.nc
    Input Nonlinear Initial File:  ./Data/ssc12_ini.nc
           Input Forcing File 01:  ./Data/ssc12_blkm1.nc
                                   ./Data/ssc12_blkm2.nc
                                   ./Data/ssc12_blkm3.nc
                                   ./Data/ssc12_blkm4.nc
                                   ./Data/ssc12_blkm5.nc
                                   ./Data/ssc12_blkm6.nc
                                   ./Data/ssc12_blkm7.nc
                                   ./Data/ssc12_blkm8.nc
                                   ./Data/ssc12_blkm9.nc
                                   ./Data/ssc12_blkm10.nc
                                   ./Data/ssc12_blkm11.nc
                                   ./Data/ssc12_blkm12.nc
             Input Boundary File:  ./Data/ssc12_brym1.nc
                                   ./Data/ssc12_brym2.nc
                                   ./Data/ssc12_brym3.nc
                                   ./Data/ssc12_brym4.nc
                                   ./Data/ssc12_brym5.nc
                                   ./Data/ssc12_brym6.nc
                                   ./Data/ssc12_brym7.nc
                                   ./Data/ssc12_brym8.nc
                                   ./Data/ssc12_brym9.nc
                                   ./Data/ssc12_brym10.nc
                                   ./Data/ssc12_brym11.nc
                                   ./Data/ssc12_brym12.nc
          Input/Output USER File:  MyFile.dat

 Generic User Parameters:

 2.0000E+01  user(01)          User parameter 01.

 Tile partition information for Grid 01:  0275x0275x0032  tiling: 004x004

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       69        1       69   152352
        1       70      138        1       69   152352
        2      139      207        1       69   152352
        3      208      275        1       69   150144
        4        1       69       70      138   152352
        5       70      138       70      138   152352
        6      139      207       70      138   152352
        7      208      275       70      138   150144
        8        1       69      139      207   152352
        9       70      138      139      207   152352
       10      139      207      139      207   152352
       11      208      275      139      207   150144
       12        1       69      208      275   150144
       13       70      138      208      275   150144
       14      139      207      208      275   150144
       15      208      275      208      275   147968

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    69.50     0.50    69.50  RHO-points
        1    69.50   138.50     0.50    69.50  RHO-points
        2   138.50   207.50     0.50    69.50  RHO-points
        3   207.50   275.50     0.50    69.50  RHO-points
        4     0.50    69.50    69.50   138.50  RHO-points
        5    69.50   138.50    69.50   138.50  RHO-points
        6   138.50   207.50    69.50   138.50  RHO-points
        7   207.50   275.50    69.50   138.50  RHO-points
        8     0.50    69.50   138.50   207.50  RHO-points
        9    69.50   138.50   138.50   207.50  RHO-points
       10   138.50   207.50   138.50   207.50  RHO-points
       11   207.50   275.50   138.50   207.50  RHO-points
       12     0.50    69.50   207.50   275.50  RHO-points
       13    69.50   138.50   207.50   275.50  RHO-points
       14   138.50   207.50   207.50   275.50  RHO-points
       15   207.50   275.50   207.50   275.50  RHO-points

        0     1.00    69.50     0.50    69.50    U-points
        1    69.50   138.50     0.50    69.50    U-points
        2   138.50   207.50     0.50    69.50    U-points
        3   207.50   275.00     0.50    69.50    U-points
        4     1.00    69.50    69.50   138.50    U-points
        5    69.50   138.50    69.50   138.50    U-points
        6   138.50   207.50    69.50   138.50    U-points
        7   207.50   275.00    69.50   138.50    U-points
        8     1.00    69.50   138.50   207.50    U-points
        9    69.50   138.50   138.50   207.50    U-points
       10   138.50   207.50   138.50   207.50    U-points
       11   207.50   275.00   138.50   207.50    U-points
       12     1.00    69.50   207.50   275.50    U-points
       13    69.50   138.50   207.50   275.50    U-points
       14   138.50   207.50   207.50   275.50    U-points
       15   207.50   275.00   207.50   275.50    U-points

        0     0.50    69.50     1.00    69.50    V-points
        1    69.50   138.50     1.00    69.50    V-points
        2   138.50   207.50     1.00    69.50    V-points
        3   207.50   275.50     1.00    69.50    V-points
        4     0.50    69.50    69.50   138.50    V-points
        5    69.50   138.50    69.50   138.50    V-points
        6   138.50   207.50    69.50   138.50    V-points
        7   207.50   275.50    69.50   138.50    V-points
        8     0.50    69.50   138.50   207.50    V-points
        9    69.50   138.50   138.50   207.50    V-points
       10   138.50   207.50   138.50   207.50    V-points
       11   207.50   275.50   138.50   207.50    V-points
       12     0.50    69.50   207.50   275.00    V-points
       13    69.50   138.50   207.50   275.00    V-points
       14   138.50   207.50   207.50   275.00    V-points
       15   207.50   275.50   207.50   275.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     162
               HaloSizeJ(1) =     162
                TileSide(1) =      75
                TileSize(1) =    5625


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

 Variable               Grid  West Edge    South Edge   East Edge    North Edge
 ---------              ----  ----------   ----------   ----------   ----------

 zeta                     1   Chapman Imp  Chapman Imp  Chapman Imp  Chapman Imp

 ubar                     1   Flather      Flather      Flather      Flather

 vbar                     1   Flather      Flather      Flather      Flather

 u                        1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 v                        1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 temp                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 salt                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 tke                      1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 Activated C-preprocessing Options:

 SSC12               South Sea of China Internal tidal System, 1/12 degree resolution
 ADD_FSOBC           Adding tidal elevation to processed OBC data.
 ADD_M2OBC           Adding tidal currents to processed OBC data.
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_SPONGE          Analytical enhanced viscosity/diffusion sponge.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 AVERAGES            Writing out time-averaged nonlinear model fields.
 BULK_FLUXES         Surface bulk fluxes parameterization.
 CURVGRID            Orthogonal curvilinear grid.
 DIAGNOSTICS_TS      Computing and writing tracer diagnostic terms.
 DIAGNOSTICS_UV      Computing and writing momentum diagnostic terms.
 DIFF_GRID           Horizontal diffusion coefficient scaled by grid size.
 DIURNAL_SRFLUX      Modulate shortwave radiation by the local diurnal cycle.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 EMINUSP             Compute Salt Flux using E-P.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 MIX_GEO_TS          Mixing of tracers along geopotential surfaces.
 MIX_GEO_UV          Mixing of momentum along geopotential surfaces.
 MPI                 MPI distributed-memory configuration.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 QCORRECTION         Surface net heat flux correction.
 K_GSCHEME           Third-order upstream advection of TKE fields.
 RADIATION_2D        Use tangential phase speed in radiation conditions.
 RAMP_TIDES          Ramping tidal forcing for one day.
 RI_SPLINES          Parabolic Spline Reconstruction for Richardson Number.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 SPLINES_VDIFF       Parabolic Spline Reconstruction for Vertical Diffusion.
 SPLINES_VVISC       Parabolic Spline Reconstruction for Vertical Viscosity.
 SSH_TIDES           Add tidal elevation to SSH climatology.
 TS_C4HADVECTION     Fourth-order centered horizontal advection of tracers.
 TS_SVADVECTION      Parabolic splines vertical advection of tracers.
 TS_DIF2             Harmonic mixing of tracers.
 UV_ADV              Advection of momentum.
 UV_COR              Coriolis term.
 UV_U3HADVECTION     Third-order upstream horizontal advection of 3D momentum.
 UV_C4VADVECTION     Fourth-order centered vertical advection of momentum.
 UV_QDRAG            Quadratic bottom stress.
 UV_TIDES            Add tidal currents to 2D momentum climatologies.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 VISC_GRID           Horizontal viscosity coefficient scaled by grid size.

 Process Information:

 Node #  0 (pid=   16064) is active.
 Node # 14 (pid=   16078) is active.
 Node #  1 (pid=   16065) is active.
 Node #  5 (pid=   16069) is active.
 Node #  8 (pid=   16072) is active.
 Node # 10 (pid=   16074) is active.
 Node #  2 (pid=   16066) is active.
 Node #  4 (pid=   16068) is active.
 Node # 15 (pid=   16079) is active.
 Node #  6 (pid=   16070) is active.
 Node #  7 (pid=   16071) is active.
 Node # 13 (pid=   16077) is active.
 Node #  9 (pid=   16073) is active.
 Node # 11 (pid=   16075) is active.
 Node # 12 (pid=   16076) is active.
 Node #  3 (pid=   16067) is active.

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

 Vertical S-coordinate System, Grid 01:

 level   S-coord     Cs-curve   Z   at hmin       at hc    half way     at hmax

    32   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    31  -0.0312500  -0.0001116       -0.188      -0.235      -0.800      -1.136
    30  -0.0625000  -0.0004504       -0.377      -0.472      -2.279      -3.630
    29  -0.0937500  -0.0010281       -0.567      -0.711      -4.473      -7.555
    28  -0.1250000  -0.0018650       -0.757      -0.951      -7.442     -13.031
    27  -0.1562500  -0.0029903       -0.949      -1.194     -11.272     -20.233
    26  -0.1875000  -0.0044433       -1.143      -1.440     -16.083     -29.396
    25  -0.2187500  -0.0062746       -1.338      -1.688     -22.025     -40.823
    24  -0.2500000  -0.0085478       -1.534      -1.939     -29.287     -54.895
    23  -0.2812500  -0.0113419       -1.733      -2.194     -38.107     -72.083
    22  -0.3125000  -0.0147533       -1.934      -2.454     -48.773     -92.966
    21  -0.3437500  -0.0188997       -2.138      -2.720     -61.637    -118.248
    20  -0.3750000  -0.0239231       -2.346      -2.992     -77.122    -148.778
    19  -0.4062500  -0.0299948       -2.557      -3.272     -95.741    -185.581
    18  -0.4375000  -0.0373204       -2.774      -3.561    -118.110    -229.889
    17  -0.4687500  -0.0461459       -2.997      -3.862    -144.963    -283.172
    16  -0.5000000  -0.0567646       -3.227      -4.176    -177.177    -347.187
    15  -0.5312500  -0.0695249       -3.466      -4.506    -215.794    -424.018
    14  -0.5625000  -0.0848395       -3.714      -4.855    -262.047    -516.135
    13  -0.5937500  -0.1031943       -3.975      -5.227    -317.391    -626.446
    12  -0.6250000  -0.1251587       -4.251      -5.626    -383.526    -758.359
    11  -0.6562500  -0.1513952       -4.543      -6.057    -462.434    -915.837
    10  -0.6875000  -0.1826680       -4.856      -6.526    -556.399   -1103.456
     9  -0.7187500  -0.2198500       -5.192      -7.039    -668.032   -1326.437
     8  -0.7500000  -0.2639239       -5.556      -7.604    -800.270   -1590.663
     7  -0.7812500  -0.3159778       -5.951      -8.229    -956.366   -1902.645
     6  -0.8125000  -0.3771881       -6.384      -8.923   -1139.838   -2269.423
     5  -0.8437500  -0.4487876       -6.858      -9.694   -1354.372   -2698.375
     4  -0.8750000  -0.5320097       -7.378     -10.553   -1603.654   -3196.882
     3  -0.9062500  -0.6280041       -7.950     -11.507   -1891.123   -3771.824
     2  -0.9375000  -0.7377147       -8.576     -12.564   -2219.601   -4428.850
     1  -0.9687500  -0.8617150       -9.259     -13.728   -2590.801   -5171.391
     0  -1.0000000  -1.0000000      -10.000     -15.000   -3004.710   -5999.420

 Time Splitting Weights for Grid 01:    ndtfast =  72    nfast =  98
 ==================================

    Primary            Secondary            Accumulated to Current Step

  1-0.0002697221589845 0.0138888888888889-0.0002697221589845 0.0138888888888889
  2-0.0005179738466200 0.0138926350299859-0.0007876960056045 0.0277815239188748
  3-0.0007447551074202 0.0138998291111890-0.0015324511130247 0.0416813530300638
  4-0.0009500661144945 0.0139101729321254-0.0024825172275192 0.0555915259621891
  5-0.0011339073179268 0.0139233682948267-0.0036164245454460 0.0695148942570158
  6-0.0012962796525070 0.0139391170075756-0.0049127041979530 0.0834540112645914
  7-0.0014371848048139 0.0139571208916382-0.0063498890027669 0.0974111321562297
  8-0.0015566255396494 0.0139770817917051-0.0079065145424162 0.1113882139479348
  9-0.0016546060858249 0.0139987015908669-0.0095611206282412 0.1253869155388017
 10-0.0017311325812999 0.0140216822309478-0.0112922532095411 0.1394085977697495
 11-0.0017862135776710 0.0140457257390214-0.0130784667872121 0.1534543235087709
 12-0.0018198606040138 0.0140705342609335-0.0148983273912259 0.1675248577697044
 13-0.0018320887900766 0.0140958101026559-0.0167304161813025 0.1816206678723603
 14-0.0018229175488246 0.0141212557802959-0.0185533337301271 0.1957419236526562
 15-0.0017923713183374 0.0141465740795851-0.0203457050484645 0.2098884977322413
 16-0.0017404803630569 0.0141714681256731-0.0220861854115213 0.2240599658579144
 17-0.0016672816343881 0.0141956414640489-0.0237534670459094 0.2382556073219634
 18-0.0015728196906507 0.0142187981534154-0.0253262867365601 0.2524744054753787
 19-0.0014571476763831 0.0142406428713411-0.0267834344129432 0.2667150483467199
 20-0.0013203283609977 0.0142608810335131-0.0281037627739409 0.2809759293802330
 21-0.0011624352367880 0.0142792189274159-0.0292661980107289 0.2952551483076488
 22-0.0009835536762878 0.0142953638612601-0.0302497516870167 0.3095505121689089
 23-0.0007837821489820 0.0143090243289864-0.0310335338359988 0.3238595364978953
 24-0.0005632334973686 0.0143199101921667-0.0315967673333673 0.3381794466900619
 25-0.0003220362723731 0.0143277328796301-0.0319188036057405 0.3525071795696921
 26-0.0000603361281147 0.0143322056056353-0.0319791397338552 0.3668393851753274
 27 0.0002217027239766 0.0143330436074147-0.0317574370098786 0.3811724287827420
 28 0.0005238960016901 0.0143299644029150-0.0312335410081885 0.3955023931856571
 29 0.0008460377792135 0.0143226880695582-0.0303875032289750 0.4098250812552152
 30 0.0011878988549600 0.0143109375448469-0.0291996043740150 0.4241360188000621
 31 0.0015492250600437 0.0142944389496391-0.0276503793139713 0.4384304577497012
 32 0.0019297355074029 0.0142729219349163-0.0257206438065684 0.4527033796846174
 33 0.0023291207815727 0.0142461200528690-0.0233915230249958 0.4669494997374865
 34 0.0027470410691043 0.0142137711531249-0.0206444819558915 0.4811632708906114
 35 0.0031831242296341 0.0141756178049429-0.0174613577262575 0.4953388886955543
 36 0.0036369638075999 0.0141314077461980-0.0138243939186576 0.5094702964417523
 37 0.0041081169846063 0.0140808943599814-0.0097162769340513 0.5235511908017337
 38 0.0045961024724372 0.0140238371796396-0.0051201744616140 0.5375750279813734
 39 0.0051003983467181 0.0139600024230780-0.0000197761148959 0.5515350304044513
 40 0.0056204398212250 0.0138891635571513 0.0056006637063291 0.5654241939616026
 41 0.0061556169628429 0.0138111018929677 0.0117562806691720 0.5792352958545702
 42 0.0067052723471717 0.0137256072129282 0.0184615530163437 0.5929609030674984
 43 0.0072686986547809 0.0136324784303286 0.0257302516711247 0.6065933814978269
 44 0.0078451362081123 0.0135315242823455 0.0335753878792370 0.6201249057801724
 45 0.0084337704490309 0.0134225640572328 0.0420091583282679 0.6335474698374053
 46 0.0090337293570240 0.0133054283565518 0.0510428876852919 0.6468528981939571
 47 0.0096440808080493 0.0131799598932598 0.0606869684933412 0.6600328580872169
 48 0.0102638298740300 0.0130460143264814 0.0709507983673713 0.6730788724136982
 49 0.0108919160629993 0.0129034611337865 0.0818427144303705 0.6859823335474847
 50 0.0115272104998922 0.0127521845218004 0.0933699249302627 0.6987345180692851
 51 0.0121685130479867 0.0125920843759686 0.1055384379782495 0.7113266024452537
 52 0.0128145493709923 0.0124230772503021 0.1183529873492417 0.7237496796955558
 53 0.0134639679357865 0.0122450973979272 0.1318169552850283 0.7359947770934830
 54 0.0141153369558014 0.0120580978432635 0.1459322922408297 0.7480528749367465
 55 0.0147671412750559 0.0118620514966551 0.1606994335158856 0.7599149264334017
 56 0.0154177791928386 0.0116569523122794 0.1761172127087242 0.7715718787456811
 57 0.0160655592290373 0.0114428164901566 0.1921827719377615 0.7830146952358378
 58 0.0167086968301175 0.0112196837230866 0.2088914687678790 0.7942343789589245
 59 0.0173453110157492 0.0109876184893350 0.2262367797836282 0.8052219974482595
 60 0.0179734209660816 0.0107467113918941 0.2442102007497098 0.8159687088401536
 61 0.0185909425496662 0.0104970805451429 0.2628011432993760 0.8264657893852965
 62 0.0191956847920284 0.0102388730097309 0.2819968280914044 0.8367046623950274
 63 0.0197853462848869 0.0099722662765083 0.3017821743762913 0.8466769286715357
 64 0.0203575115360218 0.0096974698003293 0.3221396859123131 0.8563743984718649
 65 0.0209096472597907 0.0094147265845512 0.3430493331721038 0.8657891250564161
 66 0.0214390986082931 0.0091243148170541 0.3644884317803970 0.8749134398734703
 67 0.0219430853431833 0.0088265495586056 0.3864315171235803 0.8837399894320759
 68 0.0224186979481308 0.0085217844843947 0.4088502150717111 0.8922617739164705
 69 0.0228628936819303 0.0082104136795596 0.4317131087536414 0.9004721875960301
 70 0.0232724925722585 0.0078928734895328 0.4549856013258999 0.9083650610855628
 71 0.0236441733500801 0.0075696444260292 0.4786297746759800 0.9159347055115920
 72 0.0239744693247019 0.0072412531295003 0.5026042440006819 0.9231759586410923
 73 0.0242597641994746 0.0069082743888794 0.5268640082001566 0.9300842330299717
 74 0.0244962878281439 0.0065713332194423 0.5513602960283005 0.9366555662494139
 75 0.0246801119118487 0.0062311069996069 0.5760404079401492 0.9428866732490209
 76 0.0248071456367686 0.0058883276674979 0.6008475535769178 0.9487750009165188
 77 0.0248731312524188 0.0055437839780984 0.6257206848293365 0.9543187848946172
 78 0.0248736395905937 0.0051983238218148 0.6505943244199301 0.9595171087164319
 79 0.0248040655249589 0.0048528566052787 0.6753983899448891 0.9643699653217107
 80 0.0246596233712912 0.0045083556952099 0.7000580133161803 0.9688783210169206
 81 0.0244353422283667 0.0041658609261642 0.7244933555445470 0.9730441819430847
 82 0.0241260612594974 0.0038264811729924 0.7486194168040444 0.9768706631160772
 83 0.0237264249147164 0.0034913969888327 0.7723458417187609 0.9803620601049099
 84 0.0232308780936106 0.0031618633094617 0.7955767198123714 0.9835239234143716
 85 0.0226336612488020 0.0028392122248282 0.8182103810611735 0.9863631356391998
 86 0.0219288054300778 0.0025248558185948 0.8401391864912513 0.9888879914577946
 87 0.0211101272691677 0.0022202890765104 0.8612493137604190 0.9911082805343050
 88 0.0201712239051703 0.0019270928644386 0.8814205376655893 0.9930353733987436
 89 0.0191054678506273 0.0016469369768668 0.9005260055162165 0.9946823103756104
 90 0.0179060017982467 0.0013815832567192 0.9184320073144632 0.9960638936323296
 91 0.0165657333682728 0.0011328887872991 0.9349977406827360 0.9971967824196287
 92 0.0150773297965061 0.0009028091571842 0.9500750704792420 0.9980995915768129
 93 0.0134332125629703 0.0006934017988994 0.9635082830422124 0.9987929933757124
 94 0.0116255519612283 0.0005068294021915 0.9751338350034406 0.9992998227779039
 95 0.0096462616083458 0.0003453634027300 0.9847800966117864 0.9996451861806338
 96 0.0074869928955037 0.0002113875470585 0.9922670895072901 0.9998565737276923
 97 0.0051391293792585 0.0001074015346210 0.9974062188865486 0.9999639752623133
 98 0.0025937811134509 0.0000360247376868 0.9999999999999996 1.0000000000000000

 ndtfast, nfast =   72  98   nfast/ndtfast =  1.36111

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

    1.000000000000 1.025151749473 0.512575874736 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.24456

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

 Minimum X-grid spacing, DXmin =  8.31956155E+00 km
 Maximum X-grid spacing, DXmax =  9.23851059E+00 km
 Minimum Y-grid spacing, DYmin =  8.32752431E+00 km
 Maximum Y-grid spacing, DYmax =  9.23803707E+00 km
 Minimum Z-grid spacing, DZmin =  1.87946450E-01 m
 Maximum Z-grid spacing, DZmax =  8.28028490E+02 m

 Minimum barotropic Courant Number =  7.65400315E-03
 Maximum barotropic Courant Number =  2.03630778E-01
 Maximum Coriolis   Courant Number =  2.30020631E-02

 Horizontal mixing scaled by grid size, GRDMAX =  9.23827383E+00 km

 NLM: GET_STATE - Read state initial conditions,             t =     0 00:00:00
                   (Grid 01, File: ssc12_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -1.48966948E-01 Max =  3.61947780E-01)
                - vertically integrated u-momentum component
                   (Min = -6.00885394E-01 Max =  4.97978379E-01)
                - vertically integrated v-momentum component
                   (Min = -9.50324788E-01 Max =  3.26700159E-01)
                - u-momentum component
                   (Min = -1.05035297E+00 Max =  8.54699225E-01)
                - v-momentum component
                   (Min = -1.33850665E+00 Max =  1.10605943E+00)
                - potential temperature
                   (Min =  1.14820875E+00 Max =  2.94470631E+01)
                - salinity
                   (Min =  3.24921652E+01 Max =  3.49794372E+01)
    GET_NGFLD   - tidal period
                   (Min =  4.30819210E+04 Max =  9.67262421E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Min =  1.05619027E-04 Max =  2.04697708E+01)
    GET_2DFLD   - tidal elevation phase angle
                   (Min =  8.09572363E-06 Max =  6.28317033E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Min =  2.46404577E-06 Max =  6.28318082E+00)
    GET_2DFLD   - tidal current phase angle
                   (Min =  1.82233990E-06 Max =  6.28318114E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Min =  3.85134858E-06 Max =  1.81928232E+01)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Min = -1.04224347E+00 Max =  1.16227576E+00)
forrtl: severe (408): fort: (3): Subscript #1 of the array TVAL has value 0 which is less than the lower bound of 1

Image              PC                Routine            Line        Source             
oceanG             0000000001D7107D  Unknown               Unknown  Unknown
oceanG             0000000001D6FB85  Unknown               Unknown  Unknown
oceanG             0000000001D16409  Unknown               Unknown  Unknown
oceanG             0000000001CD4A8F  Unknown               Unknown  Unknown
oceanG             0000000001CD4E92  Unknown               Unknown  Unknown
oceanG             00000000018CBDBF  get_cycle_                224  get_cycle.f90
oceanG             000000000163732F  inquire_                  277  inquire.f90
oceanG             00000000015CA9C0  get_2dfld_                102  get_2dfld.f90
oceanG             00000000012C728B  get_data_                 100  get_data.f90
oceanG             00000000006930BF  initial_                  282  initial.f90
oceanG             00000000004131ED  ocean_control_mod         130  ocean_control.f90
oceanG             000000000041235F  MAIN__                     95  master.f90
oceanG             00000000004121FC  Unknown               Unknown  Unknown
libc.so.6          00002AE0D133CCDD  Unknown               Unknown  Unknown
oceanG             00000000004120F9  Unknown               Unknown  Unknown

And then I sew the get_cycle.f90

Code: Select all

222 ELSE
223   i=MIN(ntime,Tindex+1)
224   Tend=Tval(i)
225 END IF

that's may mean that the time index of forcing file is wrong.Then I see the forcing file time

Code: Select all

the first forcing file :ssc12_blkm1.nc
lrf_time = -2, -1.75, -1.5, -1.25, -1, -0.75, -0.5, -0.25, 0, 0.25, 0.5, 
    0.75, 1, 1.25, 1.5, 1.75, 2, 2.25, 2.5, 2.75, 3, 3.25, 3.5, 3.75, 4, 
    4.25, 4.5, 4.75, 5, 5.25, 5.5, 5.75, 6, 6.25, 6.5, 6.75, 7, 7.25, 7.5, 
    7.75, 8, 8.25, 8.5, 8.75, 9, 9.25, 9.5, 9.75, 10, 10.25, 10.5, 10.75, 11, 
    11.25, 11.5, 11.75, 12, 12.25, 12.5, 12.75, 13, 13.25, 13.5, 13.75, 14, 
    14.25, 14.5, 14.75, 15, 15.25, 15.5, 15.75, 16, 16.25, 16.5, 16.75, 17, 
    17.25, 17.5, 17.75, 18, 18.25, 18.5, 18.75, 19, 19.25, 19.5, 19.75, 20, 
    20.25, 20.5, 20.75, 21, 21.25, 21.5, 21.75, 22, 22.25, 22.5, 22.75, 23, 
    23.25, 23.5, 23.75, 24, 24.25, 24.5, 24.75, 25, 25.25, 25.5, 25.75, 26, 
    26.25, 26.5, 26.75, 27, 27.25, 27.5, 27.75, 28, 28.25, 28.5, 28.75, 29, 
    29.25, 29.5, 29.75, 30, 30.25, 30.5, 30.75, 31, 31.25, 31.5, 31.75, 32, 
    32.25, 32.5, 32.75 ;

Code: Select all

the last forcing file :ssc12_blkm12.nc
lrf_time = 332, 332.25, 332.5, 332.75, 333, 333.25, 333.5, 333.75, 334, 
    334.25, 334.5, 334.75, 335, 335.25, 335.5, 335.75, 336, 336.25, 336.5, 
    336.75, 337, 337.25, 337.5, 337.75, 338, 338.25, 338.5, 338.75, 339, 
    339.25, 339.5, 339.75, 340, 340.25, 340.5, 340.75, 341, 341.25, 341.5, 
    341.75, 342, 342.25, 342.5, 342.75, 343, 343.25, 343.5, 343.75, 344, 
    344.25, 344.5, 344.75, 345, 345.25, 345.5, 345.75, 346, 346.25, 346.5, 
    346.75, 347, 347.25, 347.5, 347.75, 348, 348.25, 348.5, 348.75, 349, 
    349.25, 349.5, 349.75, 350, 350.25, 350.5, 350.75, 351, 351.25, 351.5, 
    351.75, 352, 352.25, 352.5, 352.75, 353, 353.25, 353.5, 353.75, 354, 
    354.25, 354.5, 354.75, 355, 355.25, 355.5, 355.75, 356, 356.25, 356.5, 
    356.75, 357, 357.25, 357.5, 357.75, 358, 358.25, 358.5, 358.75, 359, 
    359.25, 359.5, 359.75, 360, 360.25, 360.5, 360.75, 361, 361.25, 361.5, 
    361.75, 362, 362.25, 362.5, 362.75, 363, 363.25, 363.5, 363.75, 364, 
    364.25, 364.5, 364.75, 365, 365.25, 365.5, 365.75, 366, 366.25, 366.5, 
    366.75 ;

Code: Select all

the first bry file:ssc12_brym1.nc
temp_time = -16, 14, 44 ;

Code: Select all

the last bry file:ssc12_brym12.nc
temp_time = 318, 348, 378 ;

Code: Select all

the initial file:ssc12_ini.nc
ocean_time = 0 ;

Code: Select all

and the relative parameters:
  NTIMES == 87600
  DT == 360.0d0
  NDTFAST == 72
  DSTART =  0                     ! days
  TIDE_START =  0.0d0                      ! days
  TIME_REF =  20090101.0d0                     ! yyyymmdd.dd

I cant find where the error exist.one point is when I creat the forcing file with the NCEP data set ,I compute the previous month and next month. I don't know why but the scripts of Roms_tools of agrif version does that.
How should I do? Looking forward to your reply

[kate]

I don't think the problem is in the time indexing. What is at line 100 of your get_data.F?

Code: Select all

  get_data_                 100  get_data.f90

Which field is it trying to read there? You're right in that it's trying to read the time associated with that variable.

[lindandan]

I don't think the problem is in the time indexing. What is at line 100 of your get_data.F?

Code: Select all

  get_data_                 100  get_data.f90

Which field is it trying to read there? You're right in that it's trying to read the time associated with that variable.

the position of the get_data.f90 line 100 is try to read the variable Uwind of Surface wind components.

Code: Select all

!
!-------------------------------------------------------------
!Surface wind components.
!-------------------------------------------------------------
!
100     CALL get_2dfld (ng, iNLM, idUair, ncFRCid(idUair,ng),        &
101    &             nFfiles(ng), FRC(1,ng), update(1),           &
102    &             LBi, UBi, LBj, UBj, 2, 1,                    &
103    &             GRID(ng) % rmask,                            &
104    &             FORCES(ng) % UwindG)
105     IF (exit_flag.ne.NoError)  RETURN
106     CALL get_2dfld (ng, iNLM, idVair, ncFRCid(idVair,ng),        &

in the get_2dfld.f90

Code: Select all

!
!  If appropriate, inquire about the contents of input NetCDF file and
!  fill information arrays.
!
      IF (Linquire) THEN
        job=2
 102       CALL inquire (ng, model, job, Iout, Irec, 1, ifield, ncid,      &
     &                Lmulti, nfiles, S)
        IF (exit_flag.ne.NoError) RETURN
      END IF
!

in the inquire.f90

Code: Select all

!
!  Determine initial time record to read and cycling switch.
!
        IF (foundit) THEN
   277       CALL get_cycle (ng, model, ifield, job, Lmulti,               &
     &                    ncfile, Tname(ifield), Nrec, tdays(ng),       &
     &                    Tid, Liocycle, Clength, Trec,                 &
     &                    Tstr, Tend, Tmin, Tmax,  Tscale)
          IF (exit_flag.ne.NoError) RETURN
        END IF
!

I have check the variable Uwind, every time index has the resonable value,and in the forcing file , the time index is all the same.And the cycle attribute cycle_lenth is set to 0.Now I have no idea why it's wrong.Just a question, whether the forcing_time can be negative.

[kate]

You can't have negative time, you can't have a cycle_length of zero, and you can't have this:

in the forcing file , the time index is all the same.

[lindandan]

You can't have negative time, you can't have a cycle_length of zero, and you can't have this:

in the forcing file , the time index is all the same.

I have remove the cycle attribute, and just use bulkm2.nc and brym2.nc. but also the same problem.How to check the error.

[kate]

No need to double post - I answered your other one. But what I'd really do is run this in a debugger and figure out exactly what's going on. This can be solved in serial mode, so the free gdb would work.

[lindandan]

No need to double post - I answered your other one. But what I'd really do is run this in a debugger and figure out exactly what's going on. This can be solved in serial mode, so the free gdb would work.

Sorry, I am just too worry about my research. I used the blkm2.nc and brym2.nc which time index from 29 to 61 also cause this question. And I also change the cycle_length to 28 days of Feb., it can read the wind,but it down with error that Tend>Tstart. At last, I remove the attribute of cycle because I saw a page that you say when split the ncfile with time the model canot use cycle. But it doesnot work. I will use serial mode to try. Thanks very much.