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 6020 for branches/UKMO/icebergs_restart_single_file/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90 – NEMO

Ignore:
Timestamp:
2015-12-08T12:39:53+01:00 (8 years ago)
Author:
timgraham
Message:

Merged with head of trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/icebergs_restart_single_file/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

    r6019 r6020  
    99   !!            3.3  !  2010-04  (M. Leclair, G. Madec)  modified LF-RA 
    1010   !!            - -  !  2010-10  (C. Ethe, G. Madec) TRC-TRA merge (T-S in 4D) 
     11   !!            3.7  !  2014-01  (G. Madec) suppression of curl and hdiv from the restart 
     12   !!             -   !  2014-12  (G. Madec) remove KPP scheme 
    1113   !!---------------------------------------------------------------------- 
    1214 
     
    1820   USE oce             ! ocean dynamics and tracers  
    1921   USE dom_oce         ! ocean space and time domain 
     22   USE sbc_ice         ! only lk_lim3  
    2023   USE phycst          ! physical constants 
     24   USE eosbn2          ! equation of state            (eos bn2 routine) 
     25   USE trdmxl_oce      ! ocean active mixed layer tracers trends variables 
     26   ! 
    2127   USE in_out_manager  ! I/O manager 
    2228   USE iom             ! I/O module 
    23    USE eosbn2          ! equation of state            (eos bn2 routine) 
    24    USE trdmxl_oce      ! ocean active mixed layer tracers trends variables 
    25    USE divcur          ! hor. divergence and curl      (div & cur routines) 
    2629 
    2730   IMPLICIT NONE 
     
    130133                     CALL iom_rstput( kt, nitrst, numrow, 'tb'     , tsb(:,:,:,jp_tem) ) 
    131134                     CALL iom_rstput( kt, nitrst, numrow, 'sb'     , tsb(:,:,:,jp_sal) ) 
    132                      CALL iom_rstput( kt, nitrst, numrow, 'rotb'   , rotb      ) 
    133                      CALL iom_rstput( kt, nitrst, numrow, 'hdivb'  , hdivb     ) 
    134135                     CALL iom_rstput( kt, nitrst, numrow, 'sshb'   , sshb      ) 
    135136                     ! 
     
    138139                     CALL iom_rstput( kt, nitrst, numrow, 'tn'     , tsn(:,:,:,jp_tem) ) 
    139140                     CALL iom_rstput( kt, nitrst, numrow, 'sn'     , tsn(:,:,:,jp_sal) ) 
    140                      CALL iom_rstput( kt, nitrst, numrow, 'rotn'   , rotn      ) 
    141                      CALL iom_rstput( kt, nitrst, numrow, 'hdivn'  , hdivn     ) 
    142141                     CALL iom_rstput( kt, nitrst, numrow, 'sshn'   , sshn      ) 
    143142                     CALL iom_rstput( kt, nitrst, numrow, 'rhop'   , rhop      ) 
    144 #if defined key_zdfkpp 
    145                      CALL iom_rstput( kt, nitrst, numrow, 'rhd'    , rhd       ) 
    146 #endif 
    147143      IF( kt == nitrst ) THEN 
    148144         CALL iom_close( numrow )     ! close the restart file (only at last time step) 
     
    197193   END SUBROUTINE rst_read_open 
    198194 
     195 
    199196   SUBROUTINE rst_read 
    200197      !!----------------------------------------------------------------------  
     
    227224         CALL iom_get( numror, jpdom_autoglo, 'tb'     , tsb(:,:,:,jp_tem) ) 
    228225         CALL iom_get( numror, jpdom_autoglo, 'sb'     , tsb(:,:,:,jp_sal) ) 
    229          CALL iom_get( numror, jpdom_autoglo, 'rotb'   , rotb    ) 
    230          CALL iom_get( numror, jpdom_autoglo, 'hdivb'  , hdivb   ) 
    231226         CALL iom_get( numror, jpdom_autoglo, 'sshb'   , sshb    ) 
    232227      ELSE 
     
    239234      CALL iom_get( numror, jpdom_autoglo, 'sn'     , tsn(:,:,:,jp_sal) ) 
    240235      CALL iom_get( numror, jpdom_autoglo, 'sshn'   , sshn    ) 
    241       IF( iom_varid( numror, 'rotn', ldstop = .FALSE. ) > 0 ) THEN 
    242          CALL iom_get( numror, jpdom_autoglo, 'rotn'   , rotn    ) 
    243          CALL iom_get( numror, jpdom_autoglo, 'hdivn'  , hdivn   ) 
    244       ELSE 
    245          CALL div_cur( 0 )                              ! Horizontal divergence & Relative vorticity 
    246       ENDIF 
    247236      IF( iom_varid( numror, 'rhop', ldstop = .FALSE. ) > 0 ) THEN 
    248237         CALL iom_get( numror, jpdom_autoglo, 'rhop'   , rhop    )   ! now    potential density 
     
    250239         CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) )    
    251240      ENDIF 
    252 #if defined key_zdfkpp 
    253       IF( iom_varid( numror, 'rhd', ldstop = .FALSE. ) > 0 ) THEN 
    254          CALL iom_get( numror, jpdom_autoglo, 'rhd'    , rhd     )   ! now    in situ density anomaly 
    255       ELSE 
    256          CALL eos( tsn, rhd, fsdept_n(:,:,:) )   ! compute rhd 
    257       ENDIF 
    258 #endif 
    259241      ! 
    260242      IF( neuler == 0 ) THEN                                  ! Euler restart (neuler=0) 
     
    262244         ub   (:,:,:)   = un   (:,:,:) 
    263245         vb   (:,:,:)   = vn   (:,:,:) 
    264          rotb (:,:,:)   = rotn (:,:,:) 
    265          hdivb(:,:,:)   = hdivn(:,:,:) 
    266246         sshb (:,:)     = sshn (:,:) 
    267  
     247         ! 
    268248         IF( lk_vvl ) THEN 
    269249            DO jk = 1, jpk 
     
    271251            END DO 
    272252         ENDIF 
    273  
     253         ! 
    274254      ENDIF 
    275255      ! 
Note: See TracChangeset for help on using the changeset viewer.