1. Manu's Change:

Bottom topography deeper than 5500m set tup 5500m: Data there from Levitus is not reliable.
Closed walls along eastern, western, and nothern boundaries.
Radiation/nudging at the southern boundaries (to SODA climatology)

# elif defined NEPD
 
/*
**                     N  E  P  D
*/
#define OUT_DOUBLE
#define AVERAGES
#define AVERAGES_FLUXES
#ifdef MANU_DIAG
#define DIAGNOSTICS_UV
#define DIAGNOSTICS_TS
#endif
#define TS_U3HADVECTION
#define TS_SVADVECTION
#define UV_ADV
#define DJ_GRADPS
#define UV_COR
#define UV_VIS2
#define MIX_S_UV
#define TS_DIF2
#define MIX_S_TS
#define SALINITY
#define NONLIN_EOS
#define CURVGRID
#define MASKING
#define SOLVE3D
#define SPLINES
 
#ifdef NL_MODEL
#  define UV_QDRAG
#  define LMD_MIXING
#  define QCORRECTION
#  define SCORRECTION
#else
#  define UV_LDRAG
#  undef  LMD_MIXING
#endif
 
#undef  MY25_MIXING
#ifdef LMD_MIXING
# define DIURNAL_SRFLUX
# define SOLAR_SOURCE
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_NONLOCAL
# define LMD_SKPP
#endif
 
#undef  M2CLIMATOLOGY
#undef  M3CLIMATOLOGY
#undef  TCLIMATOLOGY
#undef  M2CLM_NUDGING
#undef  M3CLM_NUDGING
#undef  TCLM_NUDGING
! MANU: changed this below
#undef CLOSED_OBC
 
#ifdef CONVOLUTION
# define CLOSED_OBC
#endif
 
#ifdef CLOSED_OBC
# define NORTHERN_WALL
# define SOUTHERN_WALL
# define EASTERN_WALL
# define WESTERN_WALL
#else
# define NORTHERN_WALL
# define WESTERN_WALL
# define EASTERN_WALL
 
#define SPONGE
#define SOUTH_VOLCONS
 
#define SOUTH_FSGRADIENT
#define SOUTH_M2RADIATION
#define SOUTH_M2NUDGING
#define SOUTH_M3RADIATION
#define SOUTH_M3NUDGING
#define SOUTH_TRADIATION
#define SOUTH_TNUDGING
 
#endif
 
#define ANA_BSFLUX
#define ANA_BTFLUX
#if defined ADJOINT || defined TANGENT
#  define FORWARD_READ
#  define FORWARD_MIXING
#else
#  define FORWARD_RHS
#  define FORWARD_WRITE
#  define FORWARD_MIXING
#endif
 
 


2. Model and Domain

3. Preliminary Results:

Benchmark estimate



2. Model


Nx * Ny * Nz
RSM
121 * 66 * 28
ROMS
481 * 261 * 30
  • Extent; The domain covers the whole North Pacific Ocean from eastern Asia to west coast coast of North America. The domain includes the Kuroshio Extension region, East/Japan Sea, Sea of Okhotsk, Bering Sea, Gulf of Alaska, and the California Current System. Please see the domain below.
  1. Lonitudes : 114.240 - 245.0400
  2. Latitudes : 18.489 - 65.001


RSM

  • Domain
  1. dx (zontal resolution) 1.09 degree (fixed)
  2. dy (meridional resolution) : 0.98 degree in the southern boundary (equatorward) ~ 0.442 degree in the northern boundary (poleward).
  3. Note: In IBM machine such as in NCAR's bluevista and blueice, the number of grid points has to be a product of any power of 2,3,5 but the power of 3 has to be either 0 or 1 (not more than 1). This leads to IGRD and JGRD of the RSM domain to be 120 and 65.
  4. Here is the grid.nc file for RSM: rsm-NP6-grid.nc
  5. Optimal dt for RSM integration:
  6. DELTAT_REG_SPRING=360
    DELTAT_REG_SUMMER=450
    DELTAT_REG_FALL=360
    DELTAT_REG_WINTER=270
  7. > [[code]]
    > this fits RSM in BLUEVISTA
    > xdef 121 linear 114.240 1.090
    > ydef 66 levels
    > 18.489 19.471 20.447 21.417 22.380 23.337 24.287 25.230 26.166 27.094
    > 28.015 28.928 29.832 30.729 31.618 32.498 33.369 34.232 35.086 35.931
    > 36.768 37.595 38.413 39.222 40.022 40.813 41.594 42.366 43.129 43.882
    > 44.626 45.360 46.085 46.800 47.507 48.203 48.891 49.569 50.238 50.897
    > 51.548 52.189 52.821 53.444 54.058 54.663 55.259 55.846 56.425 56.995
    > 57.556 58.109 58.654 59.190 59.717 60.237 60.748 61.252 61.747 62.235
    > 62.715 63.187 63.652 64.109 64.559 65.001
    > zdef 12 levels
    > 1000 925 850 700 600 500 400 300 250 200
    > 150 100
    > [[code]]
    > [[code]]
    ===__Details of the chosen options for RSM__=== 
    * Please see the define.h file in RSM [[file:define.h]]
    * A few important things:
    # **Orogrpahy: 4 min global mountain**. #define MTN4MIN /* use global 4 min mountain instead of 8 min */
    # **Convection Parameterization: Relaxed Arakawa-Schubert scheme**. #define RAS /* Relaxed Arakawa Schubert scheme */
    # **No spetram damping by Kanamaru and Kanamitsu (2007).**
    #undef STDAMP /* Spectral tendency damping scheme for rsm */
    #undef STDAMP3 /* area average log ps correction to be used with stdamp */
    ----
     
    ==ROMS== 
    * __**Domain**__
    # The ROMS domain as constructed ( see the modified the make_grid.m [[file:make_grid_roms_NP6.m]]) such that **[[code]]**
    > NX ROMS =(NX RSM - 1) * 4 + 1=  120 * 4 + 1 = 481
    > NY ROMS =(NY RSM - 1) * 4 + 1=  65 * 4 + 1 = 261
    > **[[code]]**
    # dx: **0.2725 degree** (fixed, recall 1.09/4 = 0.2725)
    # dy: **0.2455 degree** (in the southern boundary) ~ **0.1105 degree** (in the northern boundary)
    # Each grid box of RSM has been divided equally in zonal and meridional direction by dx/4 and dy/4 respectively, where dx and dy are RSM's horizontal resolution. And the column at NX ROMS=481 and the row at NY ROMS= 261 were the same values at NX RSM 121 and NY RSM =66. Thus, 16 grid boxes of ROMS correspond to 1 RSM grid box. Here is the grid.nc file for ROMS, [[http://ecpc.ucsd.edu/%7Ehseo/scoar/roms-NP6-grid.nc|roms-NP6-grid.nc]].
     
     
    Details of the chosen options for ROMS (3.0)
    * ===**Climatological spinup solution:**===
    # ===16 yr spinup using COADS climatological surface flux forcing and WOA05 climatoliogical T/S profiles. I have used IRD's Roms_Tool to gerenerate these files.=== 
    # 
    # Geostrophic U/V are computed from T/S with the reference depth of 1000m
    # Topography was smoothed with rfactor of 0.25.
    # hmin=50m, hmax_coast=500m
    # theta_b=5, theta_s=0.4, hc=10
    # ===Explicit tracer diffusion and viscosity coefficients: VISC2=25 m2/s, TNU2=25 m2/s, (plus additional implicit diffusion and viscosity from 3rd order upstream biased advection scheme used in the simulation)=== 
    # ===**I also modified the Funtionals/ana_hmixcoef.h to increase the sponge layer (where the viscosity increases exponentially toward 1200 m2/s) The number of sponge layers are 4 grid points at each open boundary side. See** __**ana_hmixcoef.h**__===
    [[code format="cpp"]]
    ana_hmixcoef.h
    # if defined NP6
    ! sponge layer
    #  if defined UV_VIS2
          DO i=IstrR,IendR
            DO j=JstrR,JendR
              visc2_r(i,j)=visc2(ng)                                   &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(j,r8)**2/16.0_r8)        &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(i,r8)**2/16.0_r8)        &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(Mm(ng)-j,r8)**2/16.0_r8) &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(Lm(ng)-i,r8)**2/16.0_r8)
              visc2_p(i,j)=visc2(ng)                                   &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(j,r8)**2/16.0_r8)        &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(i,r8)**2/16.0_r8)        &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(Mm(ng)-j,r8)**2/16.0_r8) &
         &  +(1200.0_r8-visc2(ng))*exp(-REAL(Lm(ng)-i,r8)**2/16.0_r8)
            END DO
          END DO
    #  endif
    #  if defined TS_DIF2
          DO i=IstrR,IendR
            DO j=JstrR,JendR
              diff2(i,j,itemp)= diff2(i,j,itemp)                     &
         & +(1200.0_r8-diff2(i,j,itemp))*exp(-REAL(j,r8)**2/16.0_r8) &
         & +(1200.0_r8-diff2(i,j,itemp))*exp(-REAL(Mm(ng)-j,r8)**2/16.0_r8) &
         & +(1200.0_r8-diff2(i,j,itemp))*exp(-REAL(Lm(ng)-i,r8)**2/16.0_r8) &
         & +(1200.0_r8-diff2(i,j,itemp))*exp(-REAL(i,r8)**2/16.0_r8)
              diff2(i,j,isalt)= diff2(i,j,isalt)                     &
         & +(1200.0_r8-diff2(i,j,isalt))*exp(-REAL(j,r8)**2/16.0_r8)       &
         & +(1200.0_r8-diff2(i,j,isalt))*exp(-REAL(Mm(ng)-j,r8)**2/16.0_r8) &
         & +(1200.0_r8-diff2(i,j,isalt))*exp(-REAL(Lm(ng)-i,r8)**2/16.0_r8) &
         & +(1200.0_r8-diff2(i,j,isalt))*exp(-REAL(i,r8)**2/16.0_r8)
            END DO
          END DO
    #  endif
    # endif
     
4. Important options I used for spinup forced run. (here is the np6_spinup.h )
/* SOUTH */
#define  SOUTH_M3CLAMPED    /* Southen edge, 3D momentum, clamped condition */
#define  SOUTH_TCLAMPED     /* Southen edge, tracers, clamped condition */
#define  SOUTH_M2FLATHER    /* Southen edge, 2D momentum, Flather condition */
#define  SOUTH_FSCHAPMAN    /* Southen edge, free-surface, Chapman condition */
 
/* NORTH */
#define  NORTH_M3CLAMPED    /* Southen edge, 3D momentum, clamped condition */
#define  NORTH_TCLAMPED     /* Southen edge, tracers, clamped condition */
#define  NORTH_M2FLATHER    /* Southen edge, 2D momentum, Flather condition */
#define  NORTH_FSCHAPMAN    /* Southen edge, free-surface, Chapman condition */
 
/* EAST */
#define  EAST_M3CLAMPED    /* Southen edge, 3D momentum, clamped condition */
#define  EAST_TCLAMPED     /* Southen edge, tracers, clamped condition */
#define  EAST_M2FLATHER    /* Southen edge, 2D momentum, Flather condition */
#define  EAST_FSCHAPMAN    /* Southen edge, free-surface, Chapman condition */
 
/* WEST */
#define  WEST_M3CLAMPED    /* Southen edge, 3D momentum, clamped condition */
#define  WEST_TCLAMPED     /* Southen edge, tracers, clamped condition */
#define  WEST_M2FLATHER    /* Southen edge, 2D momentum, Flather condition */
#define  WEST_FSCHAPMAN    /* Southen edge, free-surface, Chapman condition */
 
#undef BULK_FLUXES     /* turn ON or OFF bulk fluxes computation */
#undef COOL_SKIN       /* turn ON or OFF cool skin correction */
#undef LONGWAVE_OUT
#undef SWFLUX

SCOAR Coupled Simulation:

  • Use IC from the end solution of spinup. Same boundary condition from WOA05, the same OBC clamping.
  • IC (representing December 30) looks like this
  • NP6_IC.zeta.jpg
    SSH aftere 16 year spinup
    NP6_IC.sst.jpg
    NP6_IC.sss.jpg
    SSS aftere 16 year spinup
  • But use BULK_FLUX, COOL_SKIN, LONGWAVE_OUT and SWFLUX in the coupled run. (See the np6_scoar.h



More detailed RSM and ROMS Grid comparison
RSM:
Lonitudes : 114.240 : 1.090 : 245.0400
Latitudes :
18.489 19.471 20.447 21.417 22.380 23.337 24.287 25.230 26.166 27.094
28.015 28.928 29.832 30.729 31.618 32.498 33.369 34.232 35.086 35.931
36.768 37.595 38.413 39.222 40.022 40.813 41.594 42.366 43.129 43.882
44.626 45.360 46.085 46.800 47.507 48.203 48.891 49.569 50.238 50.897
51.548 52.189 52.821 53.444 54.058 54.663 55.259 55.846 56.425 56.995
57.556 58.109 58.654 59.190 59.717 60.237 60.748 61.252 61.747 62.235
62.715 63.187 63.652 64.109 64.559 65.001

ROMS:
Longitudes: 18.4890 : 0.2725 : 65.0010
Latitudes:
18.4890 18.7345 18.9800 19.2255 19.4710 19.7150 19.9590 20.2030 20.4470 20.6895 20.9320
21.1745 21.4170 21.6578 21.8985 22.1393 22.3800 22.6193 22.8585 23.0977 23.3370 23.5745
23.8120 24.0495 24.2870 24.5227 24.7585 24.9943 25.2300 25.4640 25.6980 25.9320 26.1660
26.3980 26.6300 26.8620 27.0940 27.3242 27.5545 27.7848 28.0150 28.2432 28.4715 28.6998
28.9280 29.1540 29.3800 29.6060 29.8320 30.0562 30.2805 30.5048 30.7290 30.9512 31.1735
31.3957 31.6180 31.8380 32.0580 32.2780 32.4980 32.7157 32.9335 33.1512 33.3690 33.5847
33.8005 34.0162 34.2320 34.4455 34.6590 34.8725 35.0860 35.2972 35.5085 35.7197 35.9310
36.1402 36.3495 36.5588 36.7680 36.9748 37.1815 37.3882 37.5950 37.7995 38.0040 38.2085
38.4130 38.6152 38.8175 39.0198 39.2220 39.4220 39.6220 39.8220 40.0220 40.2197 40.4175
40.6153 40.8130 41.0083 41.2035 41.3987 41.5940 41.7870 41.9800 42.1730 42.3660 42.5568
42.7475 42.9382 43.1290 43.3173 43.5055 43.6937 43.8820 44.0680 44.2540 44.4400 44.6260
44.8095 44.9930 45.1765 45.3600 45.5412 45.7225 45.9038 46.0850 46.2638 46.4425 46.6212
46.8000 46.9767 47.1535 47.3302 47.5070 47.6810 47.8550 48.0290 48.2030 48.3750 48.5470
48.7190 48.8910 49.0605 49.2300 49.3995 49.5690 49.7362 49.9035 50.0708 50.2380 50.4027
50.5675 50.7323 50.8970 51.0598 51.2225 51.3852 51.5480 51.7082 51.8685 52.0288 52.1890
52.3470 52.5050 52.6630 52.8210 52.9767 53.1325 53.2883 53.4440 53.5975 53.7510 53.9045
54.0580 54.2092 54.3605 54.5117 54.6630 54.8120 54.9610 55.1100 55.2590 55.4057 55.5525
55.6992 55.8460 55.9907 56.1355 56.2802 56.4250 56.5675 56.7100 56.8525 56.9950 57.1352
57.2755 57.4157 57.5560 57.6942 57.8325 57.9708 58.1090 58.2452 58.3815 58.5178 58.6540
58.7880 58.9220 59.0560 59.1900 59.3217 59.4535 59.5853 59.7170 59.8470 59.9770 60.1070
60.2370 60.3648 60.4925 60.6202 60.7480 60.8740 61.0000 61.1260 61.2520 61.3758 61.4995
61.6232 61.7470 61.8690 61.9910 62.1130 62.2350 62.3550 62.4750 62.5950 62.7150 62.8330
62.9510 63.0690 63.1870 63.3032 63.4195 63.5358 63.6520 63.7662 63.8805 63.9947 64.1090
64.2215 64.3340 64.4465 64.5590 64.6695 64.7800 64.8905 65.0010