Ignore:
Timestamp:
2017-09-12T20:46:13+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part6 - one more round

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

    r8516 r8517  
    258258      CALL ice_itd_init                ! ice thickness distribution initialization 
    259259      ! 
    260       CALL ice_dyn_init                ! set ice dynamics parameters 
    261       ! 
    262       CALL ice_rdgrft_init             ! set ice ridging/rafting parameters 
    263       ! 
    264       CALL ice_rhg_init                ! set ice rheology parameters 
    265       ! 
    266       CALL ice_adv_init                ! set ice advection parameters 
    267       ! 
    268       CALL ice_thd_init                ! set ice thermodynics parameters 
    269       ! 
    270       CALL ice_thd_sal_init            ! set ice salinity parameters 
    271         
     260      IF( ln_icedyn )   THEN 
     261         CALL ice_dyn_init             ! set ice dynamics parameters 
     262         CALL ice_rdgrft_init          ! set ice ridging/rafting parameters 
     263         CALL ice_rhg_init             ! set ice rheology parameters 
     264         CALL ice_adv_init             ! set ice advection parameters 
     265      ENDIF 
     266 
     267      IF( ln_icethd ) THEN 
     268         CALL ice_thd_init             ! set ice thermodynics parameters 
     269         CALL ice_thd_sal_init         ! set ice salinity parameters 
     270      ENDIF 
     271    
    272272      ! MV MP 2016 
    273273      CALL lim_mp_init                 ! set melt ponds parameters 
     
    283283      CALL ice_var_glo2eqv 
    284284      ! 
    285       CALL ice_update_init                 ! ice surface boundary condition 
    286       ! 
    287       CALL ice_alb_init                    ! ice surface albedo 
    288       ! 
    289       CALL ice_dia_init                    ! initialization for diags 
    290       ! 
    291       fr_i  (:,:)   = at_i(:,:)         ! initialisation of sea-ice fraction 
    292       tn_ice(:,:,:) = t_su(:,:,:)       ! initialisation of surface temp for coupled simu 
    293       ! 
     285      CALL ice_update_init             ! ice surface boundary condition 
     286      ! 
     287      CALL ice_alb_init                ! ice surface albedo 
     288      ! 
     289      CALL ice_dia_init                ! initialization for diags 
     290      ! 
     291      fr_i  (:,:)   = at_i(:,:)        ! initialisation of sea-ice fraction 
     292      tn_ice(:,:,:) = t_su(:,:,:)      ! initialisation of surface temp for coupled simu 
     293      ! 
     294      !                                ! set max concentration in both hemispheres 
    294295      WHERE( gphit(:,:) > 0._wp )   ;   rn_amax_2d(:,:) = rn_amax_n  ! NH 
    295296      ELSEWHERE                     ;   rn_amax_2d(:,:) = rn_amax_s  ! SH 
     
    306307      !! 
    307308      !! ** Method  :   Read the namice_run namelist and check the parameter 
    308       !!              values called at the first timestep (nit000) 
     309      !!                values called at the first timestep (nit000) 
    309310      !! 
    310311      !! ** input   :   Namelist namice_run 
    311312      !!------------------------------------------------------------------- 
    312313      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    313       NAMELIST/namice_run/ jpl, nlay_i, nlay_s, nn_monocat, rn_amax_n, rn_amax_s, cn_icerst_in, cn_icerst_indir,   & 
    314          &                cn_icerst_out, cn_icerst_outdir 
     314      NAMELIST/namice_run/ jpl, nlay_i, nlay_s, nn_monocat, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s,  & 
     315         &                 cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir 
    315316      !!------------------------------------------------------------------- 
    316317      ! 
     
    333334         WRITE(numout,*) '      number of snow layers                                  nlay_s = ', nlay_s 
    334335         WRITE(numout,*) '      virtual ITD mono-category param (1-4) or not (0)   nn_monocat = ', nn_monocat 
     336         WRITE(numout,*) '      Ice dynamics       (T) or not (F)                  ln_icedyn  = ', ln_icedyn 
     337         WRITE(numout,*) '      Ice thermodynamics (T) or not (F)                  ln_icethd  = ', ln_icethd 
    335338         WRITE(numout,*) '      maximum ice concentration for NH                              = ', rn_amax_n  
    336339         WRITE(numout,*) '      maximum ice concentration for SH                              = ', rn_amax_s 
    337340      ENDIF 
    338341      ! 
    339       !                                         !--- check consistency 
     342      !                                        !--- check consistency 
    340343      IF ( jpl > 1 .AND. nn_monocat == 1 ) THEN 
    341344         nn_monocat = 0 
     
    349352      IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('ice_run_init: online conservation check does not work with BDY') 
    350353      ! 
    351       rdt_ice   = REAL(nn_fsbc) * rdt           !--- sea-ice timestep and inverse 
     354      rdt_ice   = REAL(nn_fsbc) * rdt          !--- sea-ice timestep and inverse 
    352355      r1_rdtice = 1._wp / rdt_ice 
    353       IF( lwp ) WRITE(numout,*) '   ice timestep rdt_ice  = ', rdt_ice 
    354       ! 
    355       r1_nlay_i = 1._wp / REAL( nlay_i, wp )    !--- inverse of nlay_i and nlay_s 
     356      IF( lwp ) WRITE(numout,*) '   ice timestep rdt_ice = ', rdt_ice 
     357      ! 
     358      r1_nlay_i = 1._wp / REAL( nlay_i, wp )   !--- inverse of nlay_i and nlay_s 
    356359      r1_nlay_s = 1._wp / REAL( nlay_s, wp ) 
    357360      ! 
Note: See TracChangeset for help on using the changeset viewer.