New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11986 for NEMO/branches/2019/ENHANCE-02_ISF_domcfg/src/domwri.F90 – NEMO

Ignore:
Timestamp:
2019-11-27T15:43:33+01:00 (5 years ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_domcfg: changes needed after Dave's review

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/ENHANCE-02_ISF_domcfg/src/domwri.F90

    r11201 r11986  
    154154      CALL iom_rstput( 0, 0, inum, 'ff_t', ff_t, ktype = jp_r8 ) 
    155155       
    156       ! note that mbkt is set to 1 over land ==> use surface tmask 
     156      ! note that mbkt and mikt is set to 1 over land ==> use surface tmask 
    157157      zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) , wp ) 
    158       CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 )     !    ! nb of ocean T-points 
     158      CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 )                              !    ! nb of ocean T-points 
     159      CALL iom_rstput( 0, 0, inum, 'bathy_metry', bathy(:,:)   * ssmask(:,:), ktype = jp_r8 )   !    ! bathymetry 
    159160      zprt(:,:) = ssmask(:,:) * REAL( mikt(:,:) , wp ) 
    160       CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 )       !    ! nb of ocean T-points 
    161       zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp ) 
    162       CALL iom_rstput( 0, 0, inum, 'isfdraft', zprt, ktype = jp_r8 )   !    ! nb of ocean T-points 
    163       !                                                         ! vertical mesh 
     161      CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 )                                !    ! first wet level 
     162      CALL iom_rstput( 0, 0, inum, 'isfdraft'   , risfdep(:,:) * ssmask(:,:), ktype = jp_r8 )   !    ! ice shelf draft 
     163      zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) - mikt(:,:) + 1, wp ) 
     164      CALL iom_rstput( 0, 0, inum, 'mhw',zprt, ktype = jp_i4 ) 
     165      CALL iom_rstput( 0, 0, inum, 'hw' ,(bathy-risfdep)*ssmask, ktype = jp_r8 ) 
     166 
     167      !                  ! vertical mesh 
    164168      CALL iom_rstput( 0, 0, inum, 'e3t_0', e3t_0, ktype = jp_r8  )    !    ! scale factors 
    165169      CALL iom_rstput( 0, 0, inum, 'e3u_0', e3u_0, ktype = jp_r8  ) 
Note: See TracChangeset for help on using the changeset viewer.