Changeset 11308 for NEMO/branches/2019/ENHANCE-02_ISF_domcfg/src/domain.F90

Timestamp:

2019-07-19T15:04:36+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_domcfg: create domisf to isolate isf related work and add new isf geometry computation used by UKESM (ticket #2142)

Location:

NEMO/branches/2019/ENHANCE-02_ISF_domcfg

Files:

• . (copied) (copied from NEMO/branches/2019/ENHANCE-03_domcfg)
• src/domain.F90 (modified) (3 diffs)

NEMO/branches/2019/ENHANCE-02_ISF_domcfg/src/domain.F90

@@ -105 +105 @@
          &             ln_cfmeta, ln_iscpl
       NAMELIST/namdom/ nn_bathy, cn_topo, cn_bath, cn_lon, cn_lat, nn_interp,                        &
-         &             rn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh, rn_hmin, rn_isfhmin,           &
-         &             rn_atfp , rn_rdt   , ln_crs      , jphgr_msh ,                  &
+         &             rn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh, rn_hmin,                       &
+         &             rn_atfp , rn_rdt   , ln_crs      , jphgr_msh ,                                &
          &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m,                         &
          &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh,                  &
@@ -204 +204 @@
          WRITE(numout,*) '      min depth of the ocean    (>0) or    rn_hmin   = ', rn_hmin
          WRITE(numout,*) '      min number of ocean level (<0)       '
-         WRITE(numout,*) '      treshold to open the isf cavity   rn_isfhmin   = ', rn_isfhmin, ' (m)'
          WRITE(numout,*) '      minimum thickness of partial      rn_e3zps_min = ', rn_e3zps_min, ' (m)'
          WRITE(numout,*) '         step level                     rn_e3zps_rat = ', rn_e3zps_rat
@@ -415 +414 @@
          END DO
       END DO
-      CALL iom_rstput( 0, 0, inum, 'bathy_metry'   , z2d , ktype = jp_r4 )
+      CALL iom_rstput( 0, 0, inum, 'bathy_metry_e3'   , z2d , ktype = jp_r4 )
+      DO jj = 1,jpj
+         DO ji = 1,jpi
+            z2d (ji,jj) = SUM ( e3t_0(ji,jj, 1:mikt(ji,jj)-1 ) ) * ssmask(ji,jj) 
+         END DO
+      END DO
+      CALL iom_rstput( 0, 0, inum, 'isf_draft_e3'   , z2d , ktype = jp_r4 )
+      CALL iom_rstput( 0, 0, inum, 'isf_draft'   , risfdep , ktype = jp_r4 )
+      CALL iom_rstput( 0, 0, inum, 'bathy_metry'   , bathy , ktype = jp_r4 )
+      CALL iom_rstput( 0, 0, inum, 'hw',bathy-risfdep, ktype = jp_r4 )
+      CALL iom_rstput( 0, 0, inum, 'mhw',mbkt*ssmask-mikt*ssmask, ktype = jp_i4 )
       !
       !                              !== closed sea ==!