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 10009 for NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/DIA/diawri.F90 – NEMO

Ignore:
Timestamp:
2018-07-29T11:23:51+02:00 (6 years ago)
Author:
gm
Message:

#1911 (ENHANCE-04): RK3 branch - step II.1 time-level dimension on ssh

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/DIA/diawri.F90

    r10001 r10009  
    136136 
    137137      IF( ll_wd ) THEN 
    138          CALL iom_put( "ssh" , (sshn+ssh_ref)*tmask(:,:,1) )   ! sea surface height (brought back to the reference used for wetting and drying) 
     138         CALL iom_put( "ssh" , ( ssh(:,:,Nnn)+ssh_ref)*ssmask(:,:) )   ! sea surface height (brought back to the reference used for wetting and drying) 
    139139      ELSE 
    140          CALL iom_put( "ssh" , sshn )              ! sea surface height 
     140         CALL iom_put( "ssh" ,   ssh(:,:,Nnn) )              ! sea surface height 
    141141      ENDIF 
    142142 
    143143      IF( iom_use("wetdep") )   &                  ! wet depth 
    144          CALL iom_put( "wetdep" , ht_0(:,:) + sshn(:,:) ) 
     144         CALL iom_put( "wetdep" , ht_0(:,:) + ssh(:,:,Nnn) ) 
    145145       
    146146      CALL iom_put( "toce", tsn(:,:,:,jp_tem) )    ! 3D temperature 
    … …  
    771771         CALL histwrite( nid_T, "vovvldef", it, zw3d             , ndim_T , ndex_T  )   ! level thickness deformation 
    772772      ENDIF 
    773       CALL histwrite( nid_T, "sossheig", it, sshn          , ndim_hT, ndex_hT )   ! sea surface height 
     773      CALL histwrite( nid_T, "sossheig", it, ssh(:,:,Nnn)  , ndim_hT, ndex_hT )   ! sea surface height 
    774774      CALL histwrite( nid_T, "sowaflup", it, ( emp-rnf )   , ndim_hT, ndex_hT )   ! upward water flux 
    775775      CALL histwrite( nid_T, "sorunoff", it, rnf           , ndim_hT, ndex_hT )   ! river runoffs 
    … …  
    993993      CALL histwrite( id_i, "votemper", kt, tsn(:,:,:,jp_tem), jpi*jpj*jpk, idex )    ! now temperature 
    994994      CALL histwrite( id_i, "vosaline", kt, tsn(:,:,:,jp_sal), jpi*jpj*jpk, idex )    ! now salinity 
    995       CALL histwrite( id_i, "sossheig", kt, sshn             , jpi*jpj    , idex )    ! sea surface height 
     995      CALL histwrite( id_i, "sossheig", kt, ssh(:,:,Nnn)     , jpi*jpj    , idex )    ! sea surface height 
    996996      CALL histwrite( id_i, "vozocrtx", kt, un               , jpi*jpj*jpk, idex )    ! now i-velocity 
    997997      CALL histwrite( id_i, "vomecrty", kt, vn               , jpi*jpj*jpk, idex )    ! now j-velocity 
Note: See TracChangeset for help on using the changeset viewer.