- Timestamp:
- 2017-12-14T15:47:12+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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 !
Note: See TracChangeset
for help on using the changeset viewer.