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 9051 – NEMO

Changeset 9051


Ignore:
Timestamp:
2017-12-14T15:47:12+01:00 (6 years ago)
Author:
flavoni
Message:

put back correct routines for TOOLS/DOMAINcfg, do not change these routines in merge development or merge phase

Location:
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/field_def.xml

    r9019 r9051  
    9898         <field id="tpt_dep"      long_name="T-point depth"                  standard_name="depth_below_geoid"   unit="m"   grid_ref="grid_T_3D" /> 
    9999         <field id="e3tdef"       long_name="T-cell thickness deformation"                                       unit="%"   grid_ref="grid_T_3D" /> 
    100        </field_group> 
    101  
    102        <!-- OSMOSIS diagnostics --> 
    103        <field_group id="OSMOSIS_T" grid_ref="grid_T_2D" operation="instant" > 
    104          <field id="ghamt"          long_name="non-local temperature flux"    grid_ref="grid_T_3D"  unit="deg m/s"        /> 
    105          <field id="ghams"          long_name="non-local salinity flux"        grid_ref="grid_T_3D"   unit="psu m/s"        /> 
    106          <field id="zwth0"          long_name="surface non-local temperature flux"                             unit="deg m/s"        /> 
    107          <field id="zws0"          long_name="surface non-local salinity flux"                                     unit="psu m/s"        /> 
    108          <field id="hbl"          long_name="boundary layer depth"                                                unit="m"        /> 
    109          <field id="hbli"          long_name="initial boundary layer depth"                                     unit="m"        /> 
    110          <field id="dstokes"          long_name="stokes drift  depth scale"                                     unit="m"        /> 
    111          <field id="zustke"          long_name="magnitude of stokes drift  at T-points"                     unit="m/s"        /> 
    112          <field id="zwstrc"          long_name="convective velocity scale"                                  unit="m/s"        /> 
    113          <field id="zwstrl"          long_name="langmuir velocity scale"                                   unit="m/s"        /> 
    114          <field id="zustar"          long_name="friction velocity"                                           unit="m/s"        /> 
    115          <field id="zhbl"          long_name="boundary layer depth"                                           unit="m"        /> 
    116          <field id="zhml"          long_name="mixed layer depth"                                           unit="m"        /> 
    117        </field_group> 
    118      
    119        <field_group id="OSMOSIS_U" grid_ref="grid_U_2D" > 
    120          <field id="ghamu"          long_name="non-local u-momentum flux"        grid_ref="grid_U_3D"   unit="m^2/s^2"        /> 
    121          <field id="us_x"          long_name="i component of Stokes drift"                             unit="m/s"        /> 
    122        </field_group> 
    123     
    124        <field_group id="OSMOSIS_V" grid_ref="grid_V_2D" > 
    125          <field id="ghamv"          long_name="non-local v-momentum flux"        grid_ref="grid_V_3D"   unit="m^2/s^2"        /> 
    126          <field id="us_y"          long_name="j component of Stokes drift"                             unit="m/s"        /> 
    127        </field_group> 
    128  
     100      </field_group> 
    129101 
    130102      <!-- Tides --> 
  • branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/namelist_cfg

    r9033 r9051  
    5858   ppkth2      =  999999.              ! 
    5959   ppacr2      =  999999.              ! 
    60    rn_wd_ref_depth   =    0.0          !  Reference Depth for WAD cases (0 all other cases) 
    6160/ 
    6261!----------------------------------------------------------------------- 
  • branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/namelist_ref

    r9033 r9051  
    170170   ppkth2      =       48.029893720000 ! 
    171171   ppacr2      =       13.000000000000 ! 
    172    ln_wd       =             .false.   !  T/F activation of wetting and drying 
    173    rn_wd_ref_depth   =    0.0          !  Reference Depth for WAD cases (0 all other cases) 
     172/ 
     173!----------------------------------------------------------------------- 
     174&namwad        !   Wetting and drying                                   (default F) 
     175!----------------------------------------------------------------------- 
     176   ln_wd       = .false.   !  T/F activation of wetting and drying 
     177   rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells 
     178   rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells 
     179   rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed 
     180   nn_wdit     =  10       !  Max iterations for W/D limiter 
    174181/ 
    175182!----------------------------------------------------------------------- 
     
    265272                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component 
    266273                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component  
     274   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
     275                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
     276                           !  = 0  Average per-category fluxes (forced and coupled mode) 
     277                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
     278                           !  = 2  Redistribute a single flux over categories (coupled mode only) 
    267279                     ! Sea-ice : 
    268280   nn_ice      = 2         !  =0 no ice boundary condition   , 
    269281                           !  =1 use observed ice-cover      , 
    270                            !  =2-3 ice-model used                       ("key_lim3", "key_cice") 
     282                           !  =2 ice-model used                         ("key_lim3", "key_lim2", "key_cice") 
    271283   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    272284                           !  =1 levitating ice with mass and salt exchange but no presure effect 
     
    658670   bn_tem      = 'amm12_bdyT_tra',         24        , 'votemper',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    659671   bn_sal      = 'amm12_bdyT_tra',         24        , 'vosaline',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
     672! for lim2 
     673!   bn_frld    = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
     674!   bn_hicif   = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
     675!   bn_hsnif   = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    660676! for lim3 
    661677!   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
  • branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/src/dom_oce.f90

    r9033 r9051  
    3131   !                                    !!* Namelist namdom : time & space domain * 
    3232   INTEGER , PUBLIC ::   nn_bathy        !: = 0/1 ,compute/read the bathymetry file 
    33    LOGICAL   PUBLIC ::   ln_wd           !: Is it a WAD domain? T/F Default T 
    34    REAL(wp), PUBLIC ::   rn_wd_ref_depth !: Reference depth for Wet and Dry 
    3533   REAL(wp), PUBLIC ::   rn_bathy        !: depth of flat bottom (active if nn_bathy=0; if =0 depth=jpkm1) 
    3634   REAL(wp), PUBLIC ::   rn_hmin         !: minimum ocean depth (>0) or minimum number of ocean levels (<0) 
  • branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/src/domain.f90

    r9033 r9051  
    160160         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m,                         & 
    161161         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh,                  & 
    162          &             ppa2, ppkth2, ppacr2, rn_wd_ref_depth, ln_wd 
     162         &             ppa2, ppkth2, ppacr2 
    163163 
    164164 
     
    272272         WRITE(numout,*) '      flag read/compute bathymetry      nn_bathy     = ', nn_bathy 
    273273         WRITE(numout,*) '      Depth (if =0 bathy=jpkm1)         rn_bathy     = ', rn_bathy 
    274          WRITE(numout,*) '      WAD Domain?                       ln_wd        = ', ln_wd 
    275          WRITE(numout,*) '      WAD Reference depth)              rn_wd_ref_depth     = ', rn_wd_ref_depth 
    276274         WRITE(numout,*) '      min depth of the ocean    (>0) or    rn_hmin   = ', rn_hmin 
    277275         WRITE(numout,*) '      min number of ocean level (<0)       ' 
     
    433431      IF( ln_isfcav ) THEN   ;   icav = 1   ;   ELSE   ;   icav = 0   ;   ENDIF 
    434432      CALL iom_rstput( 0, 0, inum, 'ln_isfcav', REAL( icav, wp), ktype = jp_i4 ) 
    435       z2d(:,:) = hbatt(:,:) ! add back on reference height to get appox dep 
    436                                  !this is later corrected for with specified min depth bg user for above greoid 
    437                                  ! WAD points 
    438       !where (z2d   (:,:).lte.1e-5)  z2d(:,:) = -10.0 
    439       where (tmask  (:,:,1).eq.0)  z2d(:,:) = 0.0 
    440       IF( ln_wd ) THEN 
    441           CALL iom_rstput( 0, 0, inum, 'rn_wd_ref_depth'   , rn_wd_ref_depth  ) ! replace this later with variable 
    442           CALL iom_rstput( 0, 0, inum, 'ht_wd', z2d )        !    ht_wd 
    443       ENDIF 
    444  
    445433      ! 
    446434      !                             !==  horizontal mesh  ! 
  • branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90

    r9033 r9051  
    19571957         WRITE(numout,*) '        Bottom cell (Zb) (m) = H*rn_zb_a + rn_zb_b' 
    19581958      ENDIF 
    1959       ! CEOD 
    1960       ! Add on some reference level to make nagative bathy positive to create 
    1961       ! sensible e3t_0's 
    1962 ! Can used commented code here to force bathy near open bdy not to have any WAD 
    1963 ! points 
    1964  
    1965 !CEOD MAKE points on the perimeter assuming a bondary is here on shallower than 
    1966 !10m so that tide can be specified here. 
    1967 !      DO jj = 1, jpj 
    1968 !         DO ji = 1, jpi 
    1969 !            !CEODIF( mig(ji) <= 2 .or. mjg(jj) <= 2  .or. mig(ji) >= jpiglo-2 .or.  mjg(jj) >=jpjglo-2) then 
    1970 !            IF( mig(ji) <= 10   .or. mjg(jj) <=10  .or.  mjg(ji) >=jpiglo-10+1 .or. mjg(jj) >=jpjglo-10+1   ) then 
    1971 !               IF( bathy(ji,jj) + rn_wd_ref_depth  > 0.0)   then 
    1972 !                  bathy(ji,jj) = max(5.0,bathy(ji,jj)) ! 
    1973 !               ENDIF 
    1974 !               IF( bathy(ji,jj)  + rn_wd_ref_depth   <= 0.0)   then 
    1975 !                  bathy(ji,jj) =  -rn_wd_ref_depth*2._wp ! Make sure it is land 
    1976 !               ENDIF 
    1977 !            ENDIF 
    1978 !         END DO 
    1979 !      END DO 
    1980  
    1981       IF( ln_wd) THEN 
    1982          bathy(:,:) = bathy(:,:) + rn_wd_ref_depth !  Arbitrary reference level for WAD, This needs 
    1983       ELSE 
    1984          bathy(:,:) = bathy(:,:)  !  Arbitrary reference level for WAD, This needs 
    1985       ENDIF 
    1986                                                 !  to be accounted for in istate in the actual model and bdys etc. 
     1959 
    19871960      hift(:,:) = rn_sbot_min                     ! set the minimum depth for the s-coordinate 
    19881961      hifu(:,:) = rn_sbot_min 
     
    21052078      hbatt(:,:) = zenv(:,:)  
    21062079      IF( MINVAL( gphit(:,:) ) * MAXVAL( gphit(:,:) ) <= 0._wp ) THEN 
    2107           IF ( jphgr_msh == 2 .OR. jphgr_msh == 3) CALL ctl_stop( 'dom:zgr_sco:  if jphgr_msh = 2 or 3 and  &  
    2108                                       & s-coordinates stop, if not correction at Equator is applied, but it is wrong') 
    21092080         CALL ctl_warn( ' s-coordinates are tapered in vicinity of the Equator' ) 
    21102081         DO jj = 1, jpj 
Note: See TracChangeset for help on using the changeset viewer.