Changeset 9722

Timestamp:

2018-06-01T11:26:52+02:00 (6 years ago)

Author:

flavoni

Message:

add bathy_meter field in domain_cfg.nc input files

File:

• utils/tools/DOMAINcfg/src/domain.f90 (modified) (1 diff)

utils/tools/DOMAINcfg/src/domain.f90

@@ -481 +481 @@
       CALL iom_rstput( 0, 0, inum, 'bottom_level' , REAL( mbkt, wp )*ssmask , ktype = jp_i4 )   ! nb of ocean T-points
       CALL iom_rstput( 0, 0, inum, 'top_level'    , REAL( mikt, wp )*ssmask , ktype = jp_i4 )   ! nb of ocean T-points (ISF)
+      DO jj = 1,jpj
+         DO ji = 1,jpi
+            z2d (ji,jj) = SUM ( e3t_0(ji,jj, 1:mbkt(ji,jj) ) ) * ssmask(ji,jj) 
+         END DO
+      END DO
+      CALL iom_rstput( 0, 0, inum, 'bathy_metry'   , z2d , ktype = jp_r4 )
+
       !
       IF( ln_sco ) THEN             ! s-coordinate: store grid stiffness ratio  (Not required anyway)