how to keep constant temp when running tide only modal?

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
xiaozhu557
Posts: 61
Joined: Fri Sep 11, 2009 1:48 pm
Location: nmefc

how to keep constant temp when running tide only modal?

#1 Post by xiaozhu557 » Tue Feb 10, 2015 3:34 am

I am running the tides only model and keeping constant temp. But I have got the surface temperature as follow. I also found the temperature of all boundaries in the bottom layers are not constant. What's the reason for this?

Please find the attachment for my model's configure.
Attachments
Temperature on the bottom layer
Temperature on the bottom layer
SST
SST

xiaozhu557
Posts: 61
Joined: Fri Sep 11, 2009 1:48 pm
Location: nmefc

Re: how to keep constant temp when running tide only modal?

#2 Post by xiaozhu557 » Tue Feb 10, 2015 7:44 am

Though I have not defined TS_FIXED, I have set the initial conditions, OBC to T0(ng), and set the surface and bottom fluxes to 0, by ana_initial, ana_tobc.h, ana_stflux.h and ana_btflux.h.

So I really can not understand why the temperature can be changed still during the LOOP?

d.kobashi
Posts: 60
Joined: Tue Sep 28, 2010 11:59 pm
Location: Texas A&M University

Re: how to keep constant temp when running tide only modal?

#3 Post by d.kobashi » Tue Feb 10, 2015 11:58 am

You don't need to define TS_FIXED for tide-only modeling. Just constant tempe/salt as you specified.
It looks to me that your problem may be something to do with open boundaries.
What's your cppdefs and ocean.in?

Also, why did you edit ana_tobc.h? I don't think I did that for tide-only modeling.

Cheers,

DJ@TAMU

User avatar
kate
Posts: 3780
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: how to keep constant temp when running tide only modal?

#4 Post by kate » Tue Feb 10, 2015 4:50 pm

For a really simple tides-only problem you can turn off SOLVE3D. Then you won't have any temperature to be evolving. The one thing you would have to change in ocean.in is to set DT to your fast DT and NDTFAST to 1.

xiaozhu557
Posts: 61
Joined: Fri Sep 11, 2009 1:48 pm
Location: nmefc

Re: how to keep constant temp when running tide only modal?

#5 Post by xiaozhu557 » Thu Feb 12, 2015 1:21 am

d.kobashi wrote:You don't need to define TS_FIXED for tide-only modeling. Just constant tempe/salt as you specified.
It looks to me that your problem may be something to do with open boundaries.
What's your cppdefs and ocean.in?

Also, why did you edit ana_tobc.h? I don't think I did that for tide-only modeling.

Cheers,

DJ@TAMU
My cppdefs are including but not limited as follow,
UV_* ANA_BSFLUX ANA_BTFLUX VISC_GRID MIX_S_UV ANA_FSOBC ANA_M2OBC ANA_M3OBC ANA_TOBC ANA_INITIAL ANA_STFLUX ANA_SMFLUX, RADIATION_2D SPLINES SOLVE3D CURVGRID MASKING DJ_GRADPS.

For other cppdefs, you can check the attached file "batch-41443.out", as follow,

Code: Select all

xcn069

 Model Input Parameters:  ROMS/TOMS version 3.7  
                          Monday - February 9, 2015 -  4:17:44 PM
 -----------------------------------------------------------------------------

 The new SCS Operational Forecasting System, 1/30 degree resolution

 Operating system : Linux
 CPU/hardware     : x86_64
 Compiler system  : ifort
 Compiler command : /export/home/zhuxm/mpich-3.1/bin/mpif90
 Compiler flags   : -heap-arrays -fp-model precise -ip -O3 -free -free -free

 Input Script  : 

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

 Local Root    : /mnt/swgfs/zhuxm/roms/ROMS_trunk_20141023
 Header Dir    : /export/home/zhuxm/swgfs/roms/ofs_scs
 Header file   : scs.h
 Analytical Dir: /export/home/zhuxm/swgfs/roms/ofs_scs

 Resolution, Grid 01: 0985x0793x030,  Parallel Nodes: 128,  Tiling: 016x008


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

   10000000  ntimes            Number of timesteps for 3-D equations.
     90.000  dt                Timestep size (s) for 3-D equations.
         30  ndtfast           Number of timesteps for 2-D equations between
                                 each 3D timestep.
          1  ERstr             Starting ensemble/perturbation run number.
          1  ERend             Ending ensemble/perturbation run number.
          0  nrrec             Number of restart records to read from disk.
          T  LcycleRST         Switch to recycle time-records in restart file.
        480  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).
         40  nHIS              Number of timesteps between the writing fields
                                 into history file.
     172800  ndefHIS           Number of timesteps between creation of new
                                 history files.
          1  ntsAVG            Starting timestep for the accumulation of output
                                 time-averaged data.
          0  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
     172800  ndefAVG           Number of timesteps between creation of new
                                 time-averaged file.
 1.0000E+01  nl_visc2          NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for momentum.
          T  LuvSponge         Turning ON  sponge on horizontal momentum.
          F  LtracerSponge(01) Turning OFF sponge on tracer 01: temp
          F  LtracerSponge(02) Turning OFF sponge on tracer 02: salt
 1.0000E-05  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 1.0000E-05  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-03  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 1.0000E-03  Zob               Bottom roughness (m).
          2  Vtransform        S-coordinate transformation equation.
          4  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.
     30.000  Tcline            S-coordinate surface/bottom layer width (m) used
                                 in vertical coordinate stretching.
   1027.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).
20000101.00  time_ref          Reference time for units attribute (yyyymmdd.dd)
 3.6000E+02  Tnudg(01)         Nudging/relaxation time scale (days)
                                 for tracer 01: temp
 3.6000E+02  Tnudg(02)         Nudging/relaxation time scale (days)
                                 for tracer 02: salt
 3.6000E+02  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 3.6000E+02  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 3.6000E+02  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 1.2000E+02  obcfac            Factor between passive and active
                                 open boundary conditions.
          T  VolCons(1)        NLM western  edge boundary volume conservation.
          T  VolCons(2)        NLM southern edge boundary volume conservation.
          T  VolCons(3)        NLM eastern  edge boundary volume conservation.
          T  VolCons(4)        NLM northern edge boundary volume conservation.
     14.000  T0                Background potential temperature (C) constant.
     35.000  S0                Background salinity (PSU) constant.
   1027.000  R0                Background density (kg/m3) used in linear Equation
                                 of State.
 1.7000E-04  Tcoef             Thermal expansion coefficient (1/Celsius).
 0.0000E+00  Scoef             Saline contraction coefficient (1/PSU).
      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  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(idTvar)      Write out averaged tracer 01: temp
          T  Aout(idTvar)      Write out averaged tracer 02: salt

 Output/Input Files:

             Output Restart File:  scs_rst.nc
        Prefix for History Files:  scs_his0
       Prefix for Averages Files:  scs_avg
                 Input Grid File:  scs_grd.nc
           Input Forcing File 01:  scs_tide_20000101.nc

 Tile partition information for Grid 01:  0985x0793x0030  tiling: 016x008

     tile     Istr     Iend     Jstr     Jend     Npts


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

     tile     Xmin     Xmax     Ymin     Ymax     grid

 RHO-points
   U-points
   V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     148
               HaloSizeJ(1) =     224
                TileSide(1) =     106
                TileSize(1) =    7208


 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:

 SCS                 The new SCS Operational Forecasting System, 1/30 degree resolution
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_FSOBC           Analytical free-surface boundary conditions.
 ANA_INITIAL         Analytical initial conditions.
 ANA_M2OBC           Analytical 2D momentum boundary conditions.
 ANA_M3OBC           Analytical 3D momentum boundary conditions.
 ANA_SMFLUX          Analytical kinematic surface momentum flux.
 ANA_STFLUX          Analytical kinematic surface temperature flux.
 ANA_TOBC            Analytical tracers boundary conditions.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 AVERAGES            Writing out time-averaged nonlinear model fields.
 CURVGRID            Orthogonal curvilinear grid.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 INLINE_2DIO         Processing 3D IO level by level to reduce memory needs.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 NONLINEAR           Nonlinear Model.
 !NONLIN_EOS         Linear Equation of State for seawater.
 NO_WRITE_GRID       Not Writing grid arrays into NetCDF ouput files.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 RADIATION_2D        Use tangential phase speed in radiation conditions.
 RAMP_TIDES          Ramping tidal forcing for one day.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SOLVE3D             Solving 3D Primitive Equations.
 SPLINES             Conservative parabolic spline reconstruction.
 SSH_TIDES           Add tidal elevation to SSH climatology.
 TIDES_ASTRO         Compute astronomical (long-period tide) corrections.
 TS_C4HADVECTION     Fourth-order centered horizontal advection of tracers.
 TS_C4VADVECTION     Fourth-order centered vertical advection 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:


 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

    30   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    29  -0.0333333  -0.0002307       -0.143      -0.503      -1.676      -2.372
    28  -0.0666667  -0.0009291       -0.286      -1.014      -4.742      -7.537
    27  -0.1000000  -0.0021144       -0.430      -1.532      -9.256     -15.608
    26  -0.1333333  -0.0038189       -0.574      -2.057     -15.313     -26.779
    25  -0.1666667  -0.0060892       -0.719      -2.591     -23.052     -41.328
    24  -0.2000000  -0.0089871       -0.864      -3.135     -32.657     -59.624
    23  -0.2333333  -0.0125913       -1.009      -3.689     -44.362     -82.137
    22  -0.2666667  -0.0169996       -1.155      -4.255     -58.457    -109.450
    21  -0.3000000  -0.0223307       -1.302      -4.835     -75.296    -142.273
    20  -0.3333333  -0.0287278       -1.449      -5.431     -95.303    -181.459
    19  -0.3666667  -0.0363614       -1.597      -6.045    -118.986    -228.028
    18  -0.4000000  -0.0454338       -1.747      -6.682    -146.947    -283.187
    17  -0.4333333  -0.0561830       -1.897      -7.343    -179.892    -348.356
    16  -0.4666667  -0.0688882       -2.049      -8.033    -218.652    -425.203
    15  -0.5000000  -0.0838752       -2.203      -8.758    -264.196    -515.673
    14  -0.5333333  -0.1015226       -2.358      -9.523    -317.648    -622.026
    13  -0.5666667  -0.1222681       -2.516     -10.334    -380.310    -746.874
    12  -0.6000000  -0.1466152       -2.676     -11.199    -453.679    -893.225
    11  -0.6333333  -0.1751391       -2.839     -12.127    -539.465   -1064.512
    10  -0.6666667  -0.2084930       -3.006     -13.127    -639.610   -1264.635
     9  -0.7000000  -0.2474113       -3.177     -14.211    -756.296   -1497.978
     8  -0.7333333  -0.2927113       -3.352     -15.391    -891.953   -1769.421
     7  -0.7666667  -0.3452901       -3.532     -16.679   -1049.250   -2084.319
     6  -0.8000000  -0.4061157       -3.719     -18.092   -1231.062   -2448.452
     5  -0.8333333  -0.4762095       -3.912     -19.643   -1440.427   -2867.917
     4  -0.8666667  -0.5566161       -4.112     -21.349   -1680.449   -3348.952
     3  -0.9000000  -0.6483588       -4.320     -23.225   -1954.172   -3897.665
     2  -0.9333333  -0.7523753       -4.537     -25.286   -2264.382   -4519.654
     1  -0.9666667  -0.8694298       -4.764     -27.541   -2613.351   -5219.481
     0  -1.0000000  -1.0000000       -5.000     -30.000   -3002.500   -6000.000

 Time Splitting Weights for Grid 01:    ndtfast =  30    nfast =  42
 ==================================

    Primary            Secondary            Accumulated to Current Step

  1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
  2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
  3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
  4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
  5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
  6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
  7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
  8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
  9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
 10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
 11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
 12 0.0043042707117221 0.0335513677379153-0.0022367614257356 0.4026337758325831
 13 0.0060106451121704 0.0334078920475245 0.0037738836864348 0.4360416678801076
 14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
 15 0.0099910761708920 0.0329439123123590 0.0216737068001213 0.5021931174029188
 16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429148
 17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276978
 18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780105
 19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976138
 21 0.0256737146312911 0.0297153720153352 0.1342125541712342 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524359 0.7189082721405746
 23 0.0316613792205220 0.0279045862015855 0.1945237931729578 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236653 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893981 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
 27 0.0432936737565710 0.0230872677537126 0.3506756411451924 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772280 0.8685254094681656
 29 0.0481720587108284 0.0201154903974257 0.4447073467880565 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256415 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325521 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348098 0.9524525304259084
 34 0.0531891349131379 0.0115726061288397 0.7060109510479478 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213239 0.9738247715198165
 36 0.0503158038019030 0.0080591141492892 0.8085423793232269 0.9818838856691057
 37 0.0473537721847153 0.0063819206892258 0.8558961515079423 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497611 0.9930692679747334
 39 0.0376397765791564 0.0033640675316746 0.9367177506289175 0.9964333355064080
 40 0.0305543017255206 0.0021094083123694 0.9672720523544381 0.9985427438187774
 41 0.0217382098544504 0.0010909315881854 0.9890102622088885 0.9996336754069628
 42 0.0109897377911118 0.0003663245930371 1.0000000000000004 0.9999999999999999

 ndtfast, nfast =   30  42   nfast/ndtfast =  1.40000

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

    1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

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

 Minimum X-grid spacing, DXmin =  3.23632803E+00 km
 Maximum X-grid spacing, DXmax =  9.17057348E+00 km
 Minimum Y-grid spacing, DYmin =  3.59466836E+00 km
 Maximum Y-grid spacing, DYmax =  8.96506546E+00 km
 Minimum Z-grid spacing, DZmin =  1.43021931E-01 m
 Maximum Z-grid spacing, DZmax =  7.80518607E+02 m

 Minimum barotropic Courant Number =  3.73923787E-03
 Maximum barotropic Courant Number =  2.81036039E-01
 Maximum Coriolis   Courant Number =  6.21345144E-03

 Horizontal mixing scaled by grid size, GRDMAX =  9.05374848E+00 km
    GET_NGFLD   - tidal period
                   (Min =  4.30819210E+04 Max =  9.67262421E+04)
    GET_2DFLD   - tidal elevation amplitude
                   (Min =  2.52492159E-05 Max =  1.81496844E+01)
    GET_2DFLD   - tidal elevation phase angle
                   (Min =  7.93344988E-06 Max =  6.28318465E+00)
    GET_2DFLD   - tidal current inclination angle
                   (Min =  3.70293851E-07 Max =  6.28318519E+00)
    GET_2DFLD   - tidal current phase angle
                   (Min =  1.48087661E-06 Max =  6.28318352E+00)
    GET_2DFLD   - maximum tidal current, ellipse major axis
                   (Min =  1.46277260E-06 Max =  3.88200810E+00)
    GET_2DFLD   - minimum tidal current, ellipse minor axis
                   (Min = -4.84188479E-01 Max =  4.11129338E-01)

 Basin information for Grid 01:

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

 Initial basin volumes: TotVolume =  3.7659566873E+16 m3
                        MinVolume =  1.7915665593E+06 m3
                        MaxVolume =  5.3941176727E+10 m3
                          Max/Min =  3.0108385561E+04

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


   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  2.142222E+04  2.142222E+04  3.772388E+16
          (000,000,00)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
      DEF_HIS   - creating history file, Grid 01: scs_his0_0001.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      1     0 00:01:30  1.658685E-30  2.142222E+04  2.142222E+04  3.772388E+16
          (666,679,01)  1.379357E-16  1.050798E-16  0.000000E+00  1.323787E-14
I edited ana_tobc.h file just wanting to set the four(east,west,north,south) temperature and salinity boundary conditions to T0(ng) and S0(ng). I have attached it to this post.

Thanks for your suggestions.
Attachments
ana_tobc.h
(11.45 KiB) Downloaded 62 times

xiaozhu557
Posts: 61
Joined: Fri Sep 11, 2009 1:48 pm
Location: nmefc

Re: how to keep constant temp when running tide only modal?

#6 Post by xiaozhu557 » Thu Feb 12, 2015 1:23 am

kate wrote:For a really simple tides-only problem you can turn off SOLVE3D. Then you won't have any temperature to be evolving. The one thing you would have to change in ocean.in is to set DT to your fast DT and NDTFAST to 1.
I just want to run 3D model, not 2D, though I know 2D is enough for tide only case.

Thanks. Any other suggestions still?

Post Reply