nudging_sst problem

[yjtack1234]

dear all:
Hello, I have a problem with ROMS. I want to nudge SST satellite data to my model.

So, I define NUDGING_SST and input APARNAM. But, Model blows up as soon as running.

Here is my log file. give me some advice.

Or, is there any way to nudge sst data easier than NUDGING_SST?


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

175200 ntimes Number of timesteps for 3-D equations.
180.000 dt Timestep size (s) for 3-D equations.
10 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.
480 nRST Number of timesteps between the writing of data
into restart fields.
480 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
480 nHIS Number of timesteps between the writing fields
into history file.
480 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
480 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
480 ndefAVG Number of timesteps between creation of new
time-averaged file.
2.0000E+01 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
2.0000E+01 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
2.0000E+01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E+00 blk_ZQ Height (m) of surface air humidity measurement.
2.0000E+00 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
2 lmd_Jwt Jerlov water type.
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
5.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
5.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).
20100101.00 time_ref Reference time for units attribute (yyyymmdd.dd)
1.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
1.0000E+01 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
5.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
5.0000E+00 obcfac Factor between passive and active
open boundary conditions.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T LtracerSrc(01) Processing point sources/Sink on tracer 01: temp
T LtracerSrc(02) Processing point sources/Sink on 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(idTsur) Write out surface net heat 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(idDano) Write out density anomaly.
T Hout(idVvis) Write out vertical viscosity coefficient.
T Hout(idTdif) Write out vertical T-diffusion coefficient.
T Hout(idSdif) Write out vertical S-diffusion coefficient.
T Hout(idHsbl) Write out depth of surface boundary layer.
T Hout(idHbbl) Write out depth of bottom boundary layer.

Output/Input Files:

Output Restart File: ../output/test2/ocean_rst.nc
Prefix for History Files: ../output/test2/ocean_his
Prefix for Averages Files: ../output/test2/ocean_avg
Input Grid File: ../../input2_2/08_ep/1993/roms_grid2_ADD_08_2_ep.nc
Input Nonlinear Initial File: ../../input2_2/08_ep/1993/roms_add2_08_2ep_ini.nc
Input Forcing File 01: ../../input2_2/08_ep/1993/frc_ecmwf_1993Pair_365.nc
Input Forcing File 02: ../../input2_2/08_ep/1993/frc_ecmwf_1993Qair_365.nc
Input Forcing File 03: ../../input2_2/08_ep/1993/frc_ecmwf_1993swradd_365.nc
Input Forcing File 04: ../../input2_2/08_ep/1993/frc_ecmwf_1993Tair_365.nc
Input Forcing File 05: ../../input2_2/08_ep/1993/frc_ecmwf_1993Uwind_365.nc
Input Forcing File 06: ../../input2_2/08_ep/1993/frc_ecmwf_1993Vwind_365.nc
Input Forcing File 07: ../../input2_2/08_ep/1993/roms_river_add2_2_1993.nc
Input Climatology File: ../../input2_2/08_ep/1993/roms_clm.nc
Input Boundary File: ../../input2_2/08_ep/1993/roms_add2_2_bry_SODA-Y1993.nc
Assimilation Parameters File: ./assimilation_sst.in

Tile partition information for Grid 01: 0332x0384x0040 tiling: 008x008

tile Istr Iend Jstr Jend Npts

0 1 40 1 48 76800
1 41 82 1 48 80640
2 83 124 1 48 80640
3 125 166 1 48 80640
4 167 208 1 48 80640
5 209 250 1 48 80640
6 251 292 1 48 80640
7 293 332 1 48 76800
8 1 40 49 96 76800
9 41 82 49 96 80640
10 83 124 49 96 80640
11 125 166 49 96 80640
12 167 208 49 96 80640
13 209 250 49 96 80640
14 251 292 49 96 80640
15 293 332 49 96 76800
16 1 40 97 144 76800
17 41 82 97 144 80640
18 83 124 97 144 80640
19 125 166 97 144 80640
20 167 208 97 144 80640
21 209 250 97 144 80640
22 251 292 97 144 80640
23 293 332 97 144 76800
24 1 40 145 192 76800
25 41 82 145 192 80640
26 83 124 145 192 80640
27 125 166 145 192 80640
28 167 208 145 192 80640
29 209 250 145 192 80640
30 251 292 145 192 80640
31 293 332 145 192 76800
32 1 40 193 240 76800
33 41 82 193 240 80640
34 83 124 193 240 80640
35 125 166 193 240 80640
36 167 208 193 240 80640
37 209 250 193 240 80640
38 251 292 193 240 80640
39 293 332 193 240 76800
40 1 40 241 288 76800
41 41 82 241 288 80640
42 83 124 241 288 80640
43 125 166 241 288 80640
44 167 208 241 288 80640
45 209 250 241 288 80640
46 251 292 241 288 80640
47 293 332 241 288 76800
48 1 40 289 336 76800
49 41 82 289 336 80640
50 83 124 289 336 80640
51 125 166 289 336 80640
52 167 208 289 336 80640
53 209 250 289 336 80640
54 251 292 289 336 80640
55 293 332 289 336 76800
56 1 40 337 384 76800
57 41 82 337 384 80640
58 83 124 337 384 80640
59 125 166 337 384 80640
60 167 208 337 384 80640
61 209 250 337 384 80640
62 251 292 337 384 80640
63 293 332 337 384 76800



Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 108
HaloSizeJ(1) = 120
TileSide(1) = 54
TileSize(1) = 2592


Activated C-preprocessing Options:

EASTASIANSEAS Northwest Pacific ADD
ADD_FSOBC Adding tidal elevation to proccesed OBC data.
ADD_M2OBC Adding tidal currents to proccesed OBC data.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_CLOUD Analytical cloud fraction.
ANA_RAIN Analytical rain fall rate.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged fields.
AVERAGES_AKS Writing out time-averaged vertical S-diffusion.
AVERAGES_AKT Writing out time-averaged vertical T-diffusion.
AVERAGES_FLUXES Writing out time-averaged surface fluxes.
BULK_FLUXES Surface bulk fluxes parametererization.
CURVGRID Orthogonal curvilinear grid.
DIURNAL_SRFLUX Modulate shortwave radiation by the local diurnal cycle.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EAST_FSCHAPMAN Eastern edge, free-surface, Chapman condition.
EAST_M2FLATHER Eastern edge, 2D momentum, Flather condition.
EAST_M3CLAMPED Eastern edge, 3D momentum, Clamped condition.
EAST_TCLAMPED Eastern edge, tracers, Clamped condition.
EAST_VOLCONS Eastern edge, enforce mass conservation.
LMD_BKPP KPP bottom boundary layer mixing.
LMD_CONVEC LMD convective mixing due to shear instability.
LMD_MIXING Large/McWilliams/Doney interior mixing.
LMD_NONLOCAL LMD convective nonlocal transport.
LMD_RIMIX LMD diffusivity due to shear instability.
LMD_SKPP KPP surface boundary layer mixing.
LONGWAVE Compute net longwave radiation internally.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NORTH_VOLCONS Northern edge, enforce mass conservation.
NORTH_FSCHAPMAN Northern edge, free-surface, Chapman condition.
NORTH_M2FLATHER Northern edge, 2D momentum, Flather condition.
NORTH_M3CLAMPED Northern edge, 3D momentum, Clamped condition.
NORTH_TCLAMPED Northern edge, tracers, Clamped condition.
NUDGING_SST Nudging SST observations.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
RADIATION_2D Use tangential phase speed in radiation conditions.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SOUTH_FSCHAPMAN Southern edge, free-surface, Chapman condition.
SOUTH_M2FLATHER Southern edge, 2D momentum, Flather condition.
SOUTH_M3CLAMPED Southern edge, 3D momentum, Clamped condition.
SOUTH_TCLAMPED Southern edge, tracers, Clamped condition.
SOUTH_VOLCONS Southern edge, enforce mass conservation.
SPLINES Conservative parabolic spline reconstruction.
SSH_TIDES Add tidal elevation to SSH climatology.
TCLIMATOLOGY Processing tracer climatology data.
TCLM_NUDGING Nudging toward tracer climatology.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
TS_PSOURCE Tracers point sources and sinks.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_SADVECTION Parabolic splines vertical advection of momentum.
UV_QDRAG Quadratic bottom stress.
UV_PSOURCE Mass point sources and sinks.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WEST_FSCHAPMAN Western edge, free-surface, Chapman condition.
WEST_M2FLATHER Western edge, 2D momentum, Flather condition.
WEST_M3CLAMPED Western edge, 3D momentum, Clamped condition.
WEST_TCLAMPED Western edge, tracers, Clamped condition.
WEST_VOLCONS Western edge, enforce mass conservation.

INITIAL: Configuring and initializing forward nonlinear model ...


Vertical S-coordinate System:

level S-coord Cs-curve at_hmin over_slope at_hmax

40 0.0000000 0.0000000 0.000 0.000 0.000
39 -0.0250000 -0.0017775 -0.125 -4.564 -9.003
38 -0.0500000 -0.0037835 -0.250 -9.699 -19.149
37 -0.0750000 -0.0060921 -0.375 -15.590 -30.805
36 -0.1000000 -0.0087922 -0.500 -22.458 -44.417
35 -0.1250000 -0.0119914 -0.625 -30.573 -60.522
34 -0.1500000 -0.0158197 -0.750 -40.260 -79.769
33 -0.1750000 -0.0204330 -0.875 -51.907 -102.938
32 -0.2000000 -0.0260168 -1.000 -65.977 -130.954
31 -0.2250000 -0.0327879 -1.125 -83.013 -164.901
30 -0.2500000 -0.0409944 -1.250 -103.634 -206.017
29 -0.2750000 -0.0509117 -1.375 -128.527 -255.679
28 -0.3000000 -0.0628318 -1.500 -158.422 -315.345
27 -0.3250000 -0.0770446 -1.625 -194.044 -386.463
26 -0.3500000 -0.0938098 -1.750 -236.040 -470.330
25 -0.3750000 -0.1133185 -1.875 -284.888 -567.901
24 -0.4000000 -0.1356491 -2.000 -340.784 -679.567
23 -0.4250000 -0.1607245 -2.125 -403.534 -804.944
22 -0.4500000 -0.1882840 -2.250 -472.489 -942.729
21 -0.4750000 -0.2178817 -2.375 -546.535 -1090.694
20 -0.5000000 -0.2489214 -2.500 -624.181 -1245.862
19 -0.5250000 -0.2807264 -2.625 -703.739 -1404.853
18 -0.5500000 -0.3126323 -2.750 -783.549 -1564.348
17 -0.5750000 -0.3440788 -2.875 -862.212 -1721.549
16 -0.6000000 -0.3746809 -3.000 -938.766 -1874.531
15 -0.6250000 -0.4042644 -3.125 -1012.775 -2022.425
14 -0.6500000 -0.4328656 -3.250 -1084.332 -2165.414
13 -0.6750000 -0.4607052 -3.375 -1153.986 -2304.597
12 -0.7000000 -0.4881475 -3.500 -1222.648 -2441.797
11 -0.7250000 -0.5156591 -3.625 -1291.484 -2579.342
10 -0.7500000 -0.5437741 -3.750 -1361.826 -2719.902
9 -0.7750000 -0.5730691 -3.875 -1435.115 -2866.355
8 -0.8000000 -0.6041494 -4.000 -1512.863 -3021.726
7 -0.8250000 -0.6376434 -4.125 -1596.639 -3189.154
6 -0.8500000 -0.6742046 -4.250 -1688.076 -3371.902
5 -0.8750000 -0.7145184 -4.375 -1788.885 -3573.394
4 -0.9000000 -0.7593114 -4.500 -1900.880 -3797.260
3 -0.9250000 -0.8093643 -4.625 -2026.012 -4047.400
2 -0.9500000 -0.8655258 -4.750 -2166.401 -4328.051
1 -0.9750000 -0.9287272 -4.875 -2324.371 -4643.868
0 -1.0000000 -1.0000000 -5.000 -2502.500 -5000.000

Time Splitting Weights: ndtfast = 10 nfast = 15

Primary Secondary Accumulated to Current Step

1 0.0014850648765542 0.1000000000000000 0.0014850648765542 0.1000000000000000
2 0.0059386031399402 0.0998514935123446 0.0074236680164944 0.1998514935123446
3 0.0133457074936657 0.0992576331983505 0.0207693755101601 0.2991091267106951
4 0.0236484051180389 0.0979230624489840 0.0444177806281990 0.3970321891596791
5 0.0366959666818608 0.0955582219371801 0.0811137473100598 0.4925904110968592
6 0.0521753389587948 0.0918886252689940 0.1332890862688545 0.5844790363658532
7 0.0695217010484137 0.0866710913731145 0.2028107873172682 0.6711501277389678
8 0.0878091442019233 0.0797189212682732 0.2906199315191915 0.7508690490072409
9 0.1056214752525633 0.0709380068480809 0.3962414067717548 0.8218070558553218
10 0.1209031436506851 0.0603758593228245 0.5171445504224399 0.8821829151781463
11 0.1307902921035067 0.0482855449577560 0.6479348425259466 0.9304684601359023
12 0.1314219308195445 0.0352065157474054 0.7793567733454911 0.9656749758833076
13 0.1177312353577227 0.0220643226654509 0.8970880087032138 0.9877392985487585
14 0.0832169680811588 0.0102911991296786 0.9803049767843727 0.9980304976784371
15 0.0196950232156274 0.0019695023215627 1.0000000000000000 0.9999999999999999

ndtfast, nfast = 10 15 nfast/ndtfast = 1.50000

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

1.000000000000 1.083803072626 0.541901536313 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.00000

Minimum X-grid spacing, DXmin = 7.60025906E+00 km
Maximum X-grid spacing, DXmax = 1.17437346E+01 km
Minimum Y-grid spacing, DYmin = 7.61770607E+00 km
Maximum Y-grid spacing, DYmax = 1.17402060E+01 km
Minimum Z-grid spacing, DZmin = 1.25000000E-01 m
Maximum Z-grid spacing, DZmax = 3.56132476E+02 m

Minimum barotropic Courant Number = 1.66925614E-02
Maximum barotropic Courant Number = 6.46860676E-01
Maximum Coriolis Courant Number = 2.06333617E-02


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(File: roms_add2_08_2ep_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = -2.43266248E+00 Max = 2.69932123E+01)
- salinity
(Min = 2.90106609E+01 Max = 3.52245790E+01)
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 4.85051667E-05 Max = 1.73504087E+01)
GET_2DFLD - tidal elevation phase angle
(Min = 8.99863721E-07 Max = 6.28316973E+00)
GET_2DFLD - tidal current inclination angle
(Min = 4.97278241E-06 Max = 6.28318402E+00)
GET_2DFLD - tidal current phase angle
(Min = 8.71432860E-07 Max = 6.28317602E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.32779225E-06 Max = 2.21848579E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -6.09852072E-01 Max = 4.17190933E-01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 2.50000000E+01 Max = 1.26000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 1.42000000E+02 Max = 2.39000000E+02)
GET_NGFLD - river runoff direction
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 2.50000004E-02 Max = 2.50000004E-02)
GET_NGFLD - river runoff mass transport, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_river_add2_2_1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.85000000E+02 Max = 1.82000000E+04)
GET_NGFLD - river runoff potential temperature, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_river_add2_2_1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 4.08249998E+00 Max = 1.15649996E+01)
GET_NGFLD - river runoff salinity, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_river_add2_2_1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 364 12:00:00
(Rec=0365, Index=1, File: frc_ecmwf_1993Uwind_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -1.00494452E+01 Max = 9.66639423E+00)
GET_2DFLD - surface v-wind component, t = 364 12:00:00
(Rec=0365, Index=1, File: frc_ecmwf_1993Vwind_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -1.30453968E+01 Max = 1.74682617E+01)
GET_2DFLD - surface air pressure, t = 364 12:00:00
(Rec=0365, Index=1, File: frc_ecmwf_1993Pair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 9.91656677E+02 Max = 1.03156396E+03)
GET_2DFLD - solar shortwave radiation flux, t = 364 12:00:00
(Rec=0365, Index=1, File: frc_ecmwf_1993swradd_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 5.76239458E-07 Max = 3.86578192E-05)
GET_2DFLD - surface air temperature, t = 364 12:00:00
(Rec=0365, Index=1, File: frc_ecmwf_1993Tair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -3.20782394E+01 Max = 2.68132095E+01)
GET_2DFLD - surface air relative humidity, t = 364 12:00:00
(Rec=0365, Index=1, File: frc_ecmwf_1993Qair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 2.91812077E-01 Max = 9.57439880E-01)
GET_NGFLD - free-surface western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.30773351E-01 Max = 7.71431506E-01)
GET_NGFLD - free-surface eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.23135282E-01 Max = 1.03932746E+00)
GET_NGFLD - free-surface southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.06237587E-01 Max = 9.48591274E-01)
GET_NGFLD - free-surface northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.87908625E-01 Max = 4.87641655E-02)
GET_NGFLD - 2D u-momentum western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.15328573E-01 Max = 1.49283795E-04)
GET_NGFLD - 2D v-momentum western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.04844360E-02 Max = 1.49208666E-04)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.94660487E-02 Max = 7.35969206E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.78450633E-02 Max = 5.80423266E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.22142262E-02 Max = 4.66335363E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.16974921E-02 Max = 6.95386598E-02)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.43963767E-03 Max = 1.79442171E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.27926676E-02 Max = 3.94947392E-02)
GET_NGFLD - 3D u-momentum western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.39760687E-01 Max = 4.08816604E-02)
GET_NGFLD - 3D v-momentum western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.07088160E-01 Max = 6.54563796E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.74079328E-01 Max = 4.83699049E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.33970777E-01 Max = 3.42267690E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.07057950E-01 Max = 1.15894191E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.03669770E-01 Max = 6.30104647E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.67717470E-02 Max = 5.58414561E-02)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.18416991E-01 Max = 9.26704737E-02)
GET_NGFLD - temperature western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.80255201E+00 Max = 2.60598279E+01)
GET_NGFLD - salinity western boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.02756577E+01 Max = 3.47565802E+01)
GET_NGFLD - temperature eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.34928980E-01 Max = 2.79747379E+01)
GET_NGFLD - salinity eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.24857270E+01 Max = 3.52342948E+01)
GET_NGFLD - temperature southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.04792269E+00 Max = 2.81506726E+01)
GET_NGFLD - salinity southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.37780916E+01 Max = 3.51951563E+01)
GET_NGFLD - temperature northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.53488623E-01 Max = 7.72190237E+00)
GET_NGFLD - salinity northern boundary condition, t = 345 00:00:00
(Rec=0012, Index=2, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.12212846E+01 Max = 3.44721993E+01)
GET_3DFLD - potential temperature climatology, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.84504335E+00 Max = 2.80750065E+01)
GET_3DFLD - salinity climatology, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.82957318E+01 Max = 3.52148349E+01)
GET_2DFLD - sea surface temperature observations, t = 364 12:00:00
(Rec=0073, Index=1, File: roms_5days_sate1993.nc)
(Tmin= 4.5000 Tmax= 364.5000)
(Min = -1.07333331E+00 Max = 2.74470496E+01)
GET_2DFLD - sea surface temperature error variance, t = 364 12:00:00
(Rec=0073, Index=1, File: roms_5days_sate1993.nc)
(Tmin= 4.5000 Tmax= 364.5000)
(Min = 5.00000000E-01 Max = 5.00000000E-01)

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


Initial basin volumes: TotVolume = 2.6588166380E+16 m3
MinVolume = 7.4300392740E+06 m3
MaxVolume = 4.9071852037E+10 m3
Max/Min = 6.6045212182E+03

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00175200)
GET_NGFLD - river runoff mass transport, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_river_add2_2_1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.56000000E+02 Max = 1.05000000E+04)
GET_NGFLD - river runoff potential temperature, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_river_add2_2_1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.47559571E+00 Max = 7.49845791E+00)
GET_NGFLD - river runoff salinity, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_river_add2_2_1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 0 12:00:00
(Rec=0001, Index=2, File: frc_ecmwf_1993Uwind_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -8.87222862E+00 Max = 9.14228916E+00)
GET_2DFLD - surface v-wind component, t = 0 12:00:00
(Rec=0001, Index=2, File: frc_ecmwf_1993Vwind_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -1.48801384E+01 Max = 6.07381105E+00)
GET_2DFLD - surface air pressure, t = 0 12:00:00
(Rec=0001, Index=2, File: frc_ecmwf_1993Pair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 1.01194104E+03 Max = 1.02739099E+03)
GET_2DFLD - solar shortwave radiation flux, t = 0 12:00:00
(Rec=0001, Index=2, File: frc_ecmwf_1993swradd_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 2.09976735E-06 Max = 3.78373290E-05)
GET_2DFLD - surface air temperature, t = 0 12:00:00
(Rec=0001, Index=2, File: frc_ecmwf_1993Tair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -2.41697655E+01 Max = 2.54901142E+01)
GET_2DFLD - surface air relative humidity, t = 0 12:00:00
(Rec=0001, Index=2, File: frc_ecmwf_1993Qair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 4.54055710E-01 Max = 9.27759018E-01)
GET_NGFLD - free-surface western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.52844417E-02 Max = 6.53009295E-01)
GET_NGFLD - free-surface eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.83937883E-01 Max = 1.03107993E+00)
GET_NGFLD - free-surface southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 4.05198028E-01 Max = 8.66083701E-01)
GET_NGFLD - free-surface northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.76816204E-01 Max = -2.41349568E-02)
GET_NGFLD - 2D u-momentum western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.70838303E-02 Max = 2.29165856E-03)
GET_NGFLD - 2D v-momentum western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.33806707E-02 Max = 4.80133258E-04)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.35672237E-02 Max = 8.32791361E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.58092570E-02 Max = 4.67241032E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.68170841E-02 Max = 2.38523732E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.93366955E-02 Max = 2.59837296E-02)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.58531301E-02 Max = 1.68521071E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.57071684E-02 Max = 1.56803027E-02)
GET_NGFLD - 3D u-momentum western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.25821645E-01 Max = 4.49971464E-02)
GET_NGFLD - 3D v-momentum western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.13449953E-01 Max = 4.10675805E-02)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.58205333E-01 Max = 5.45059981E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.81034366E-01 Max = 1.29429299E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.91163975E-01 Max = 1.44907935E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.62255872E-01 Max = 5.42261881E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.85694235E-02 Max = 5.25595157E-02)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.25952294E-01 Max = 3.11934491E-02)
GET_NGFLD - temperature western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.81182228E+00 Max = 2.44203925E+01)
GET_NGFLD - salinity western boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.03083383E+01 Max = 3.46355019E+01)
GET_NGFLD - temperature eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.11165111E+00 Max = 2.62109885E+01)
GET_NGFLD - salinity eastern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.27699372E+01 Max = 3.51476733E+01)
GET_NGFLD - temperature southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.03954967E+00 Max = 2.64998003E+01)
GET_NGFLD - salinity southern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.30762672E+01 Max = 3.51479991E+01)
GET_NGFLD - temperature northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.31919429E-01 Max = 4.76351833E+00)
GET_NGFLD - salinity northern boundary condition, t = 15 00:00:00
(Rec=0001, Index=1, File: roms_add2_2_bry_SODA-Y1993.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.19473957E+01 Max = 3.44690172E+01)
GET_3DFLD - potential temperature climatology, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.43266248E+00 Max = 2.69932123E+01)
GET_3DFLD - salinity climatology, t = 15 00:00:00
(Rec=0001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.90106609E+01 Max = 3.52245790E+01)
GET_2DFLD - sea surface temperature observations, t = 4 12:00:00
(Rec=0001, Index=2, File: roms_5days_sate1993.nc)
(Tmin= 4.5000 Tmax= 364.5000)
(Min = -1.42442767E+00 Max = 2.68864955E+01)
GET_2DFLD - sea surface temperature error variance, t = 4 12:00:00
(Rec=0001, Index=2, File: roms_5days_sate1993.nc)
(Tmin= 4.5000 Tmax= 364.5000)
(Min = 5.00000000E-01 Max = 5.00000000E-01)

STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME

0 0 00:00:00 3.109088E-06 2.134717E+04 2.134717E+04 2.662330E+16
DEF_HIS - creating history file: ../output/test2/ocean_his_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
GET_2DFLD - surface u-wind component, t = 1 12:00:00
(Rec=0002, Index=1, File: frc_ecmwf_1993Uwind_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -9.01124763E+00 Max = 1.14758062E+01)
GET_2DFLD - surface v-wind component, t = 1 12:00:00
(Rec=0002, Index=1, File: frc_ecmwf_1993Vwind_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -7.66629934E+00 Max = 1.14586620E+01)
GET_2DFLD - surface air pressure, t = 1 12:00:00
(Rec=0002, Index=1, File: frc_ecmwf_1993Pair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 1.00684985E+03 Max = 1.03398645E+03)
GET_2DFLD - solar shortwave radiation flux, t = 1 12:00:00
(Rec=0002, Index=1, File: frc_ecmwf_1993swradd_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 1.33823110E-06 Max = 3.99642596E-05)
GET_2DFLD - surface air temperature, t = 1 12:00:00
(Rec=0002, Index=1, File: frc_ecmwf_1993Tair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = -2.28215714E+01 Max = 2.56358318E+01)
GET_2DFLD - surface air relative humidity, t = 1 12:00:00
(Rec=0002, Index=1, File: frc_ecmwf_1993Qair_365.nc)
(Tmin= 0.5000 Tmax= 364.5000)
(Min = 3.65235634E-01 Max = 9.48288193E-01)
480 1 00:00:00 NaN NaN NaN NaN

Blowing-up: Saving latest model state into RESTART file

[kate]

Does it run without nudging? Blowing up in the first timestep is kind of suspicious.

NUDGING_SST was taken out completely at some point for being a bad idea.

[yjtack1234]

Thank you for reply.

Model run without nudging_sst. I found another way that is QCORRECTION. But, I dont understand completely dQdSST variable. So, I input dQdSST=-40 without reason. Model is running well for now. Please give me some tips about dQdSST.

[kate]

Sorry, it was NUDGING_SSH that used to be there. The code I have now has neither NUDGING_SST nor NUDGING_SSH. Which code are you using?

dQdSST is the change in heat flux due to changes in the surface temperature. Surface heat flux is strongly related to surface temperature, so if you have an externally computed heat flux, it could do the wrong thing where your SST doesn't match that used to compute the heat flux.