New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8403 – NEMO

Changeset 8403


Ignore:
Timestamp:
2017-08-03T10:05:37+02:00 (7 years ago)
Author:
deazer
Message:

Add in ROMS WAD option ln_rwd+changes for implicit Bed Friction for ln_wd option
Note no ramp placed on ROMS bed friction yet
CS15mini case added as a Test CASE
at this revision AMM15 with Pure sigma coords barotorpic runs for 4 days without failure
in with ROMS option with 20cm min deoth and 50 vertical levels
Both run for CS15mini
In real domains nothing done on reference level yet so real domains
must have not negative depth points yet.
But a basic test has been done in WAD channel test cases (WAD7)

No changes in Main line source yet. See the MY_SRC sub dir of CS15 and TEST_CASES/WAD
for actual code changes.

Location:
branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG
Files:
28 added
5 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/TEST_CASES/WAD/EXP00/namelist_cfg

    r7609 r8403  
    88   rn_dx = 1000.0 
    99   rn_dz = 1.0 
    10    nn_wad_test = 1 
     10   nn_wad_test = 7 
    1111/ 
    1212!----------------------------------------------------------------------- 
    … …  
    142142&nambdy        !  unstructured open boundaries 
    143143!----------------------------------------------------------------------- 
    144     ln_bdy         = .false.             
    145     nb_bdy         = 0                    !  number of open boundary sets 
     144    ln_bdy         = .true.             
     145    nb_bdy         = 1                    !  number of open boundary sets 
    146146    ln_coords_file = .false.              !  =T : read bdy coordinates from file 
    147147    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files 
    … …  
    394394!----------------------------------------------------------------------- 
    395395   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
    396    ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
     396   ln_dyn_trd  = .true.   ! (T) 3D momentum trend output 
    397397   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
    398398   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet) 
    399399   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
    400400   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
    401    ln_tra_trd  = .false.    ! (T) 3D tracer trend output 
     401   ln_tra_trd  = .true.    ! (T) 3D tracer trend output 
    402402   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
    403403   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
    … …  
    444444&namwad  !   Wetting and drying 
    445445!----------------------------------------------------------------------- 
    446    ln_wd             = .true.   ! T/F activation of wetting and drying 
    447    rn_wdmin1         =  0.3     ! Minimum wet depth on dried cells 
    448    rn_wdmin2         =  0.0001  ! Tolerance of min wet depth on dried cells 
     446   ln_wd             = .false.  ! T/F activation of NOC  wetting and drying scheme 
     447   ln_rwd            = .true.   ! T/F activation of ROMS wetting and drying scheme 
     448   ln_rwd_rmp        = .true.   ! Turn on the limiter 
     449   ln_rwd_bc         = .true.   ! ROMS Baroclinic option 
     450   ln_wd_diag        = .false.   ! T/F activation of diagnostics for ROMS wd scheme 
     451   rn_wdmin0         =  0.50     ! Rmp value for NOCL option 
     452   rn_wdmin1         =  0.150    ! Minimum wet depth on dried cells 
     453   rn_wdmin2         =  0.001    ! Tolerance of min wet depth on dried cells 
     454   rn_ssh_ref        =  3.0     ! reference level 
    449455   rn_wdld           =  2.5     ! Land elevation below which wetting/drying is allowed 
    450456   nn_wdit           =   20     ! Max iterations for W/D limiter 
    451 / 
     457   jn_wd_i           =  22      ! i index of diagnostics 
     458   jn_wd_j           =   3      ! j index of diagnostics 
     459   jn_wd_k           =   3      ! k index of diagnostics 
     460/ 

  • branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/TEST_CASES/WAD/MY_SRC/domain.F90

    r7616 r8403  
    674674      ENDIF 
    675675      ! 
    676       IF( ln_wd ) THEN              ! wetting and drying domain 
     676      IF( ln_wd  .or. ln_rwd ) THEN              ! wetting and drying domain 
    677677         CALL iom_rstput( 0, 0, inum, 'ht_0'   , ht_0   , ktype = jp_r8 ) 
    678678         CALL iom_rstput( 0, 0, inum, 'ht_wd'  , ht_wd  , ktype = jp_r8 ) 

  • branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/TEST_CASES/WAD/MY_SRC/usrdef_istate.F90

    r7616 r8403  
    169169            ! 
    170170      END SELECT 
     171 
     172! subtract the height of z=0 above the geoid (this allows z = 0 to be higher than all points that may become wet)    
     173      pssh(:,:) = pssh(:,:) - rn_ssh_ref 
     174 
    171175      ! 
    172176      ! Apply minimum wetdepth criterion 

  • branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/TEST_CASES/WAD/MY_SRC/usrdef_zgr.F90

    r7616 r8403  
    1919                      & mj0, mj1, glamt, gphit   ! ocean space and time domain 
    2020   USE usrdef_nam     ! User defined : namelist variables 
    21    USE wet_dry ,  ONLY: rn_wdmin1, rn_wdmin2, rn_wdld, ht_wd 
     21   USE wet_dry ,  ONLY: rn_wdmin1, rn_wdmin2, rn_wdld, ht_wd, rn_ht_0,rn_ssh_ref 
    2222   ! 
    2323   USE in_out_manager ! I/O manager 
    … …  
    232232         END SELECT 
    233233      END IF 
    234                ! 
    235       ! 
     234               
     235! increase the depth of the bathymetry by rn_ssh_ref and rn_ht_0    
     236      zht(:,:) = zht(:,:) + rn_ssh_ref + rn_ht_0   
     237 
    236238      ! at u-point: averaging zht 
    237239      DO ji = 1, jpim1 

  • branches/UKMO/ROMS_WAD_7832/NEMOGCM/CONFIG/cfg.txt

    r7715 r8403  
    77ORCA2_OFF_TRC OPA_SRC OFF_SRC TOP_SRC 
    88ORCA2_LIM3_PISCES OPA_SRC LIM_SRC_3 TOP_SRC NST_SRC 
     9CS15mini OPA_SRC 
Note: See TracChangeset for help on using the changeset viewer.