Changeset 12101 for utils/tools_UKMO_MERGE_2019/DOMAINcfg/src/domwri.F90

Timestamp:

2019-12-06T18:45:39+01:00 (4 years ago)

Author:

mathiot

Message:

merge ENHANCE-03_domcfg and ENHANCE-02_ISF_nemo

File:

• utils/tools_UKMO_MERGE_2019/DOMAINcfg/src/domwri.F90 (copied) (copied from utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/src/domwri.F90) (2 diffs)

utils/tools_UKMO_MERGE_2019/DOMAINcfg/src/domwri.F90

@@ -154 +154 @@
       CALL iom_rstput( 0, 0, inum, 'ff_t', ff_t, ktype = jp_r8 )
       
-      ! note that mbkt is set to 1 over land ==> use surface tmask
+      ! note that mbkt and mikt is set to 1 over land ==> use surface tmask
       zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) , wp )
-      CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 )     !    ! nb of ocean T-points
+      CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 )                              !    ! nb of ocean T-points
+      CALL iom_rstput( 0, 0, inum, 'bathy_metry', bathy(:,:)   * ssmask(:,:), ktype = jp_r8 )   !    ! bathymetry
       zprt(:,:) = ssmask(:,:) * REAL( mikt(:,:) , wp )
-      CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 )       !    ! nb of ocean T-points
-      zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp )
-      CALL iom_rstput( 0, 0, inum, 'isfdraft', zprt, ktype = jp_r8 )   !    ! nb of ocean T-points
-      !                                                         ! vertical mesh
+      CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 )                                !    ! first wet level
+      CALL iom_rstput( 0, 0, inum, 'isfdraft'   , risfdep(:,:) * ssmask(:,:), ktype = jp_r8 )   !    ! ice shelf draft
+      zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) - mikt(:,:) + 1, wp )
+      CALL iom_rstput( 0, 0, inum, 'mhw',zprt, ktype = jp_i4 )
+      CALL iom_rstput( 0, 0, inum, 'hw' ,(bathy-risfdep)*ssmask, ktype = jp_r8 )
+
+      !                  ! vertical mesh
       CALL iom_rstput( 0, 0, inum, 'e3t_0', e3t_0, ktype = jp_r8  )    !    ! scale factors
       CALL iom_rstput( 0, 0, inum, 'e3u_0', e3u_0, ktype = jp_r8  )
@@ -177 +181 @@
       ENDIF
       !
-   !   IF( ll_wd ) CALL iom_rstput( 0, 0, inum, 'ht_0'   , ht_0   , ktype = jp_r8 )
-
       !                                     ! ============================
       CALL iom_close( inum )                !        close the files