- Timestamp:
- 2017-12-14T11:29:10+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/src/domain.f90
r7200 r9033 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 162 & ppa2, ppkth2, ppacr2, rn_wd_ref_depth, ln_wd 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_wd 275 WRITE(numout,*) ' WAD Reference depth) rn_wd_ref_depth = ', rn_wd_ref_depth 274 276 WRITE(numout,*) ' min depth of the ocean (>0) or rn_hmin = ', rn_hmin 275 277 WRITE(numout,*) ' min number of ocean level (<0) ' … … 431 433 IF( ln_isfcav ) THEN ; icav = 1 ; ELSE ; icav = 0 ; ENDIF 432 434 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 433 445 ! 434 446 ! !== horizontal mesh !
Note: See TracChangeset
for help on using the changeset viewer.