Changeset 9051
- Timestamp:
- 2017-12-14T15:47:12+01:00 (5 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg
- Files:
-
- 6 edited
-
field_def.xml (modified) (1 diff)
-
namelist_cfg (modified) (1 diff)
-
namelist_ref (modified) (3 diffs)
-
src/dom_oce.f90 (modified) (1 diff)
-
src/domain.f90 (modified) (3 diffs)
-
src/domzgr.f90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/field_def.xml
r9019 r9051 98 98 <field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" /> 99 99 <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> 129 101 130 102 <!-- Tides --> -
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/namelist_cfg
r9033 r9051 58 58 ppkth2 = 999999. ! 59 59 ppacr2 = 999999. ! 60 rn_wd_ref_depth = 0.0 ! Reference Depth for WAD cases (0 all other cases)61 60 / 62 61 !----------------------------------------------------------------------- -
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/namelist_ref
r9033 r9051 170 170 ppkth2 = 48.029893720000 ! 171 171 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 174 181 / 175 182 !----------------------------------------------------------------------- … … 265 272 ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component 266 273 ! =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) 267 279 ! Sea-ice : 268 280 nn_ice = 2 ! =0 no ice boundary condition , 269 281 ! =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") 271 283 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 272 284 ! =1 levitating ice with mass and salt exchange but no presure effect … … 658 670 bn_tem = 'amm12_bdyT_tra', 24 , 'votemper', .true. , .false. , 'daily' , '' , '' , '' 659 671 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' , '' , '' , '' 660 676 ! for lim3 661 677 ! 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 31 31 ! !!* Namelist namdom : time & space domain * 32 32 INTEGER , PUBLIC :: nn_bathy !: = 0/1 ,compute/read the bathymetry file 33 LOGICAL PUBLIC :: ln_wd !: Is it a WAD domain? T/F Default T34 REAL(wp), PUBLIC :: rn_wd_ref_depth !: Reference depth for Wet and Dry35 33 REAL(wp), PUBLIC :: rn_bathy !: depth of flat bottom (active if nn_bathy=0; if =0 depth=jpkm1) 36 34 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 160 160 & ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 161 161 & ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 162 & ppa2, ppkth2, ppacr2 , rn_wd_ref_depth, ln_wd162 & ppa2, ppkth2, ppacr2 163 163 164 164 … … 272 272 WRITE(numout,*) ' flag read/compute bathymetry nn_bathy = ', nn_bathy 273 273 WRITE(numout,*) ' Depth (if =0 bathy=jpkm1) rn_bathy = ', rn_bathy 274 WRITE(numout,*) ' WAD Domain? ln_wd = ', ln_wd275 WRITE(numout,*) ' WAD Reference depth) rn_wd_ref_depth = ', rn_wd_ref_depth276 274 WRITE(numout,*) ' min depth of the ocean (>0) or rn_hmin = ', rn_hmin 277 275 WRITE(numout,*) ' min number of ocean level (<0) ' … … 433 431 IF( ln_isfcav ) THEN ; icav = 1 ; ELSE ; icav = 0 ; ENDIF 434 432 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 dep436 !this is later corrected for with specified min depth bg user for above greoid437 ! WAD points438 !where (z2d (:,:).lte.1e-5) z2d(:,:) = -10.0439 where (tmask (:,:,1).eq.0) z2d(:,:) = 0.0440 IF( ln_wd ) THEN441 CALL iom_rstput( 0, 0, inum, 'rn_wd_ref_depth' , rn_wd_ref_depth ) ! replace this later with variable442 CALL iom_rstput( 0, 0, inum, 'ht_wd', z2d ) ! ht_wd443 ENDIF444 445 433 ! 446 434 ! !== horizontal mesh ! -
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90
r9033 r9051 1957 1957 WRITE(numout,*) ' Bottom cell (Zb) (m) = H*rn_zb_a + rn_zb_b' 1958 1958 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 1987 1960 hift(:,:) = rn_sbot_min ! set the minimum depth for the s-coordinate 1988 1961 hifu(:,:) = rn_sbot_min … … 2105 2078 hbatt(:,:) = zenv(:,:) 2106 2079 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')2109 2080 CALL ctl_warn( ' s-coordinates are tapered in vicinity of the Equator' ) 2110 2081 DO jj = 1, jpj
Note: See TracChangeset
for help on using the changeset viewer.
