Ignore:
Timestamp:
2020-05-28T16:21:45+02:00 (5 months ago)
Author:
clem
Message:

deal with possible out of bounds ice concentration at initial time step and remove unecessary initialization of before values

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icestp.F90

    r12919 r12988  
    202202         IF( lrst_ice )                 CALL ice_rst_write( kt )      ! -- Ice restart file  
    203203         ! 
    204          IF( ln_icectl )                CALL ice_ctl( kt )            ! -- alerts in case of model crash 
     204         IF( ln_icectl )                CALL ice_ctl( kt )            ! -- Control checks 
    205205         ! 
    206206      ENDIF   ! End sea-ice time step only 
     
    223223      !! ** purpose :   Initialize sea-ice parameters 
    224224      !!---------------------------------------------------------------------- 
    225       INTEGER :: ji, jj, ierr 
     225      INTEGER :: jl, ierr 
    226226      !!---------------------------------------------------------------------- 
    227227      IF(lwp) WRITE(numout,*) 
     
    247247      IF( ierr /= 0 )   CALL ctl_stop('STOP', 'ice_init : unable to allocate ice arrays') 
    248248      ! 
     249      !                                ! set max concentration in both hemispheres 
     250      WHERE( gphit(:,:) > 0._wp )   ;   rn_amax_2d(:,:) = rn_amax_n  ! NH 
     251      ELSEWHERE                     ;   rn_amax_2d(:,:) = rn_amax_s  ! SH 
     252      END WHERE 
     253      ! 
    249254      CALL ice_itd_init                ! ice thickness distribution initialization 
    250255      ! 
     
    252257      ! 
    253258      !                                ! Initial sea-ice state 
    254  
    255259      IF ( ln_rstart .OR. nn_iceini_file == 2 ) THEN 
    256260         CALL ice_rst_read                      ! start from a restart file 
    257261      ELSE 
    258262         CALL ice_istate_init 
    259          CALL ice_istate( nit000 )              ! start from rest: sea-ice deduced from sst 
     263         CALL ice_istate( nit000 )              ! start from rest or read a file 
    260264      ENDIF 
    261265      CALL ice_var_glo2eqv 
     
    275279      tn_ice(:,:,:) = t_su(:,:,:)      ! initialisation of surface temp for coupled simu 
    276280      ! 
    277       !                                ! set max concentration in both hemispheres 
    278       WHERE( gphit(:,:) > 0._wp )   ;   rn_amax_2d(:,:) = rn_amax_n  ! NH 
    279       ELSEWHERE                     ;   rn_amax_2d(:,:) = rn_amax_s  ! SH 
    280       END WHERE 
    281  
    282281      IF( ln_rstart )   CALL iom_close( numrir )  ! close input ice restart file 
    283282      ! 
Note: See TracChangeset for help on using the changeset viewer.