Changeset 8516


Ignore:
Timestamp:
2017-09-08T19:53:20+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part5 - I think I can see the end of the tunnel

Location:
branches/2017/dev_r8183_ICEMODEL/NEMOGCM
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/namelist_ice_cfg

    r8512 r8516  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! LIM3 configuration namelist: Overwrites SHARED/namelist_ice_lim3_ref 
     2!! ESIM configuration namelist: Overwrites SHARED/namelist_ice_lim3_ref 
    33!!              1 - Generic parameters                 (namice_run) 
    44!!              2 - Ice thickness discretization       (namice_itd) 
    5 !!              3 - Ice ridging/rafting                (namice_rdgrft) 
    6 !!              4 - Ice rheology                       (namice_rhg) 
    7 !!              5 - Ice advection                      (namice_adv) 
    8 !!              6 - Ice thermodynamics                 (namice_thd) 
    9 !!              7 - Ice salinity                       (namice_sal) 
    10 !!              8 - Ice melt ponds                     (namice_mp) 
    11 !!              9 - Ice initialization                 (namice_ini) 
    12 !!             10 - Ice/snow albedos                   (namice_alb) 
    13 !!             11 - Ice diagnostics                    (namice_dia) 
     5!!              3 - Ice dynamics                       (namice_dyn) 
     6!!              4 - Ice ridging/rafting                (namice_rdgrft) 
     7!!              5 - Ice rheology                       (namice_rhg) 
     8!!              6 - Ice advection                      (namice_adv) 
     9!!              7 - Ice thermodynamics                 (namice_thd) 
     10!!              8 - Ice salinity                       (namice_sal) 
     11!!              9 - Ice melt ponds                     (namice_mp) 
     12!!             10 - Ice initialization                 (namice_ini) 
     13!!             11 - Ice/snow albedos                   (namice_alb) 
     14!!             12 - Ice diagnostics                    (namice_dia) 
    1415!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1516! 
     
    2021!------------------------------------------------------------------------------ 
    2122&namice_itd     !   Ice discretization 
     23!------------------------------------------------------------------------------ 
     24/ 
     25!------------------------------------------------------------------------------ 
     26&namice_dyn     !   Ice dynamics 
    2227!------------------------------------------------------------------------------ 
    2328/ 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r8515 r8516  
    33!!              1 - Generic parameters                 (namice_run) 
    44!!              2 - Ice thickness discretization       (namice_itd) 
    5 !!              3 - Ice ridging/rafting                (namice_rdgrft) 
    6 !!              4 - Ice rheology                       (namice_rhg) 
    7 !!              5 - Ice advection                      (namice_adv) 
    8 !!              6 - Ice thermodynamics                 (namice_thd) 
    9 !!              7 - Ice salinity                       (namice_sal) 
    10 !!              8 - Ice melt ponds                     (namice_mp) 
    11 !!              9 - Ice initialization                 (namice_ini) 
    12 !!             10 - Ice/snow albedos                   (namice_alb) 
    13 !!             11 - Ice diagnostics                    (namice_dia) 
     5!!              3 - Ice dynamics                       (namice_dyn) 
     6!!              4 - Ice ridging/rafting                (namice_rdgrft) 
     7!!              5 - Ice rheology                       (namice_rhg) 
     8!!              6 - Ice advection                      (namice_adv) 
     9!!              7 - Ice thermodynamics                 (namice_thd) 
     10!!              8 - Ice salinity                       (namice_sal) 
     11!!              9 - Ice melt ponds                     (namice_mp) 
     12!!             10 - Ice initialization                 (namice_ini) 
     13!!             11 - Ice/snow albedos                   (namice_alb) 
     14!!             12 - Ice diagnostics                    (namice_dia) 
    1415!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1516! 
     
    3637   rn_himean        =   2.0           !  expected domain-average ice thickness (m) 
    3738   rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
     39/ 
     40!------------------------------------------------------------------------------ 
     41&namice_dyn     !   Ice dynamics 
     42!------------------------------------------------------------------------------ 
     43   ln_icedyn        = .true.          !  ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
     44      nn_icedyn     =   2             !     switch for ice dynamics    
     45                                      !        2: total 
     46                                      !        1: advection only (no diffusion, no ridging/rafting) 
     47                                      !        0: advection only (as 1 but with prescribed velocity, bypass rheology) 
     48         rn_uice    =   0.00001       !           prescribed ice u-velocity 
     49         rn_vice    =  -0.00001       !           prescribed ice v-velocity 
     50   rn_ishlat        =   2.            !  free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 
     51   rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-) 
     52   ln_landfast      = .false.         !  landfast ice parameterization (T or F)                            
     53      rn_gamma      =   0.15          !     fraction of ocean depth that ice must reach to initiate landfast 
     54                                      !        recommended range: [0.1 ; 0.25] 
     55      rn_icebfr     =  10.            !     maximum bottom stress per unit area of contact (N/m2)                  
     56                                      !        a very large value ensures ice velocity=0 even with a small contact area 
     57                                      !        recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2) 
     58      rn_lfrelax    =   1.e-5         !     relaxation time scale to reach static friction (s-1)                  
    3859/ 
    3960!------------------------------------------------------------------------------ 
     
    6788&namice_rhg     !   Ice rheology 
    6889!------------------------------------------------------------------------------ 
    69    rn_ishlat        =   2.            !  free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 
    70    rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-) 
    71    rn_creepl        =   1.0e-12       !  creep limit (s-1) 
    72    rn_ecc           =   2.0           !  eccentricity of the elliptical yield curve           
    73    nn_nevp          = 120             !  number of EVP subcycles                              
    74    rn_relast        =   0.333         !  ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
    75                                       !     advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
    76    ln_landfast      = .false.         !  landfast ice parameterization (T or F)                            
    77       rn_gamma      =   0.15          !     fraction of ocean depth that ice must reach to initiate landfast 
    78                                       !        recommended range: [0.1 ; 0.25] 
    79       rn_icebfr     =  10.            !     maximum bottom stress per unit area of contact (N/m2)                  
    80                                       !        a very large value ensures ice velocity=0 even with a small contact area 
    81                                       !        recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2) 
    82       rn_lfrelax    =   1.e-5         !     relaxation time scale to reach static friction (s-1)                  
     90   ln_rhg_EVP       = .true.          !  EVP rheology 
     91      rn_creepl     =   1.0e-12       !     creep limit (s-1) 
     92      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve           
     93      nn_nevp       = 120             !     number of EVP subcycles                              
     94      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
     95                                      !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
    8396/ 
    8497!------------------------------------------------------------------------------ 
    8598&namice_adv     !   Ice advection 
    8699!------------------------------------------------------------------------------ 
    87    ln_icedyn        = .true.          !  ice dynamics (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 
    88       nn_icedyn     =   2             !     switch for ice dynamics    
    89                                       !        2: total 
    90                                       !        1: advection only (no diffusion, no ridging/rafting) 
    91                                       !        0: advection only (as 1 but with prescribed velocity, bypass rheology) 
    92          rn_uice    =   0.00001       !           prescribed ice u-velocity 
    93          rn_vice    =  -0.00001       !           prescribed ice v-velocity 
    94                       ! -- ice_adv -- ! 
    95100   ln_adv_Pra       = .false.         !  Advection scheme (Prather) 
    96101   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho) 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r8514 r8516  
    163163   CHARACTER(len=256), PUBLIC ::   cn_icerst_indir !: ice restart input directory 
    164164   CHARACTER(len=256), PUBLIC ::   cn_icerst_outdir!: ice restart output directory 
     165 
     166   !                                     !!** ice-itd namelist (namice_itd) ** 
     167   REAL(wp), PUBLIC ::   rn_himin         !: minimum ice thickness 
    165168    
    166    !                                     !!** ice-diagnostics namelist (namice_dia) ** 
    167    LOGICAL , PUBLIC ::   ln_icediachk     !: flag for ice diag (T) or not (F) 
    168    LOGICAL , PUBLIC ::   ln_icediahsb     !: flag for ice diag (T) or not (F) 
    169    LOGICAL , PUBLIC ::   ln_icectl        !: flag for sea-ice points output (T) or not (F) 
    170    INTEGER , PUBLIC ::   iiceprt          !: debug i-point 
    171    INTEGER , PUBLIC ::   jiceprt          !: debug j-point 
    172  
    173    !                                     !!** ice-dynamics namelist (namicedyn) ** 
     169   !                                     !!** ice-dynamics namelist (namice_dyn) ** 
    174170   LOGICAL , PUBLIC ::   ln_icedyn        !: flag for ice dynamics (T) or not (F) 
    175171   INTEGER , PUBLIC ::   nn_icedyn        !: flag for ice dynamics 
    176172   REAL(wp), PUBLIC ::   rn_uice          !: prescribed u-vel (case nn_icedyn=0) 
    177173   REAL(wp), PUBLIC ::   rn_vice          !: prescribed v-vel (case nn_icedyn=0) 
    178                                           ! -- iceadv -- ! 
    179    LOGICAL , PUBLIC ::   ln_adv_Pra       !: Prather        advection scheme 
    180    LOGICAL , PUBLIC ::   ln_adv_UMx       !: Ultimate-Macho advection scheme 
    181    INTEGER , PUBLIC ::   nn_UMx           !: order of the UMx advection scheme    
    182                                           ! -- icerdgrft -- ! 
     174   REAL(wp), PUBLIC ::   rn_ishlat        !: lateral boundary condition for sea-ice 
     175   REAL(wp), PUBLIC ::   rn_cio           !: drag coefficient for oceanic stress 
     176   LOGICAL , PUBLIC ::   ln_landfast      !: landfast ice parameterization (T or F)  
     177   REAL(wp), PUBLIC ::   rn_gamma         !: fraction of ocean depth that ice must reach to initiate landfast ice 
     178   REAL(wp), PUBLIC ::   rn_icebfr        !: maximum bottom stress per unit area of contact (landfast ice)  
     179   REAL(wp), PUBLIC ::   rn_lfrelax       !: relaxation time scale (s-1) to reach static friction (landfast ice)  
     180   ! 
     181   !                                     !!** ice-rdige/raft namelist (namice_rdgrft) ** 
    183182   LOGICAL , PUBLIC ::   ln_str_H79       !: ice strength parameterization (Hibler79) 
    184183   REAL(wp), PUBLIC ::   rn_pstar         !: determines ice strength, Hibler JPO79 
     
    186185   LOGICAL , PUBLIC ::   ln_str_R75       !: ice strength parameterization (Rothrock75) 
    187186   REAL(wp), PUBLIC ::   rn_perdg         !: ridging work divided by pot. energy change in ridging 
    188                                           ! -- icerhg -- ! 
    189    REAL(wp), PUBLIC ::   rn_ishlat        !: lateral boundary condition for sea-ice 
    190    REAL(wp), PUBLIC ::   rn_cio           !: drag coefficient for oceanic stress 
     187   ! 
     188   !                                     !!** ice-rheology namelist (namice_rhg) ** 
     189   LOGICAL , PUBLIC ::   ln_rhg_EVP       !: EVP rheology 
    191190   REAL(wp), PUBLIC ::   rn_creepl        !: creep limit : has to be under 1.0e-9 
    192191   REAL(wp), PUBLIC ::   rn_ecc           !: eccentricity of the elliptical yield curve 
    193192   INTEGER , PUBLIC ::   nn_nevp          !: number of iterations for subcycling 
    194193   REAL(wp), PUBLIC ::   rn_relast        !: ratio => telast/rdt_ice (1/3 or 1/9 depending on nb of subcycling nevp)  
    195    LOGICAL , PUBLIC ::   ln_landfast      !: landfast ice parameterization (T or F)  
    196    REAL(wp), PUBLIC ::   rn_gamma         !: fraction of ocean depth that ice must reach to initiate landfast ice 
    197    REAL(wp), PUBLIC ::   rn_icebfr        !: maximum bottom stress per unit area of contact (landfast ice)  
    198    REAL(wp), PUBLIC ::   rn_lfrelax       !: relaxation time scale (s-1) to reach static friction (landfast ice)  
    199  
     194   ! 
     195   !                                     !!** ice-advection namelist (namice_adv) ** 
     196   LOGICAL , PUBLIC ::   ln_adv_Pra       !: Prather        advection scheme 
     197   LOGICAL , PUBLIC ::   ln_adv_UMx       !: Ultimate-Macho advection scheme 
     198   INTEGER , PUBLIC ::   nn_UMx           !: order of the UMx advection scheme    
     199   ! 
    200200   !                                     !!** ice-thermodynamics namelist (namice_thd) ** 
    201201   LOGICAL , PUBLIC ::   ln_icethd        !: flag for ice thermo (T) or not (F) 
     
    221221   REAL(wp), PUBLIC ::   rn_vfraz         !: threshold drift speed for collection of bottom frazil ice 
    222222   REAL(wp), PUBLIC ::   rn_Cfraz         !: squeezing coefficient for collection of bottom frazil ice 
    223                                           ! -- iceitd_th -- ! 
    224    REAL(wp), PUBLIC ::   rn_himin         !: minimum ice thickness 
    225                                           ! --  -- ! 
     223                                          ! -- icethd -- ! 
    226224   INTEGER , PUBLIC ::   nn_iceflx        !: Redistribute heat flux over ice categories 
    227225   !                                      !   =-1  Do nothing (needs N(cat) fluxes) 
     
    254252   REAL(wp), PUBLIC ::   rn_hpnd          !: prescribed pond depth    (0<rn_hpnd<1), only if nn_pnd_scheme = 0 
    255253   ! END MV MP 2016 
     254   !                                     !!** ice-diagnostics namelist (namice_dia) ** 
     255   LOGICAL , PUBLIC ::   ln_icediachk     !: flag for ice diag (T) or not (F) 
     256   LOGICAL , PUBLIC ::   ln_icediahsb     !: flag for ice diag (T) or not (F) 
     257   LOGICAL , PUBLIC ::   ln_icectl        !: flag for sea-ice points output (T) or not (F) 
     258   INTEGER , PUBLIC ::   iiceprt          !: debug i-point 
     259   INTEGER , PUBLIC ::   jiceprt          !: debug j-point 
    256260 
    257261   !                                     !!** some other parameters  
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceadv.F90

    r8514 r8516  
    3636   PUBLIC   ice_adv        ! called by icestp 
    3737   PUBLIC   ice_adv_init   ! called by icestp 
    38  
    39    INTEGER, PUBLIC ::              nice_dyn   !: choice of the type of advection scheme 
    40    !                                  ! associated indices: 
    41    INTEGER, PUBLIC, PARAMETER ::   np_dynNO   = 0   !: no ice dynamics and ice advection 
    42    INTEGER, PUBLIC, PARAMETER ::   np_dynFULL = 1   !: full ice dynamics  (rheology + advection + ridging/rafting + correction) 
    43    INTEGER, PUBLIC, PARAMETER ::   np_dyn     = 2   !: no ridging/rafting (rheology + advection                   + correction) 
    44    INTEGER, PUBLIC, PARAMETER ::   np_dynPURE = 3   !: pure dynamics      (rheology + advection)  
    4538 
    4639   !! * Substitution 
     
    250243      INTEGER ::   ios   ! Local integer output status for namelist read 
    251244      !! 
    252       NAMELIST/namice_adv/ ln_icedyn, nn_icedyn, rn_uice, rn_vice, ln_adv_Pra, ln_adv_UMx, nn_UMx 
     245      NAMELIST/namice_adv/ ln_adv_Pra, ln_adv_UMx, nn_UMx 
    253246      !!------------------------------------------------------------------- 
    254247      ! 
     
    267260         WRITE(numout,*) '~~~~~~~~~~~~' 
    268261         WRITE(numout,*) '   Namelist namice_adv' 
    269          WRITE(numout,*) '      Ice dynamics       (T) or not (F)                         ln_icedyn  = ', ln_icedyn 
    270          WRITE(numout,*) '         associated switch                                      nn_icedyn  = ', nn_icedyn 
    271          WRITE(numout,*) '            =2 all processes (default option)' 
    272          WRITE(numout,*) '            =1 advection only (no ridging/rafting)' 
    273          WRITE(numout,*) '            =0 advection only with prescribed velocity given by ' 
    274          WRITE(numout,*) '               a uniform field               (u,v)_ice = (rn_uice,rn_vice) = (', rn_uice,',', rn_vice,')' 
    275262         WRITE(numout,*) '      advection scheme for ice transport (limtrp)' 
    276263         WRITE(numout,*) '         type of advection scheme (Prather)                     ln_adv_Pra = ', ln_adv_Pra  
     
    279266      ENDIF 
    280267      ! 
    281       !                             ! set the choice of ice dynamics 
    282       IF( lk_c1d .OR. .NOT. ln_icedyn ) THEN 
    283          nice_dyn = np_dynNO                    !--- no dynamics 
    284       ELSE 
    285          SELECT CASE( nn_icedyn ) 
    286          CASE( 2 )                     
    287             IF( nn_monocat /= 2 ) THEN          !--- full dynamics (rheology + advection + ridging/rafting + correction) 
    288                nice_dyn = np_dynFULL 
    289             ELSE 
    290                nice_dyn = np_dyn                !--- dynamics without ridging/rafting 
    291             ENDIF 
    292          CASE( 0 , 1 )                          !--- dynamics without ridging/rafting and correction  
    293             nice_dyn = np_dynPURE 
    294          END SELECT 
    295       ENDIF 
    296       !                                         !--- simple conservative piling, comparable with LIM2 
    297       l_piling = nn_icedyn == 1 .OR. ( nn_monocat == 2  .AND.  jpl == 1 ) 
    298       ! 
    299268      IF ( ( ln_adv_Pra .AND. ln_adv_UMx ) .OR. ( .NOT.ln_adv_Pra .AND. .NOT.ln_adv_UMx ) ) THEN 
    300269         CALL ctl_stop( 'ice_adv_init: choose one and only one ice advection scheme (ln_adv_Pra or ln_adv_UMx)' ) 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerhg.F90

    r8514 r8516  
    2020   USE icerhg_evp     ! sea-ice: EVP rheology 
    2121   USE icectl         ! sea-ice: control prints 
    22    USE icevar         ! sea-ice: operations 
    2322   ! 
    2423   USE lbclnk         ! lateral boundary conditions - MPP exchanges 
     
    5554      INTEGER, INTENT(in) ::   kt     ! ice time step 
    5655      !! 
    57       INTEGER  ::   jl   ! dummy loop indices 
    5856      REAL(wp) :: zvi_b, zsmv_b, zei_b, zfs_b, zfw_b, zft_b  
    5957      !!-------------------------------------------------------------------- 
     
    6664         WRITE(numout,*)'~~~~~~~' 
    6765      ENDIF 
     66      !                             ! -- conservation test 
     67      IF( ln_icediachk   )   CALL ice_cons_hsm(0, 'icerhg', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    6868 
    69       CALL ice_var_agg(1)           ! -- aggregate ice categories 
    70       ! 
    71       !                             ! -- conservation test 
    72       IF( ln_icediachk )   CALL ice_cons_hsm(0, 'icerhg', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    73       !                       
    74       IF( ln_landfast ) THEN        ! -- Landfast ice parameterization: define max bottom friction 
    75          tau_icebfr(:,:) = 0._wp 
    76          DO jl = 1, jpl 
    77             WHERE( ht_i(:,:,jl) > ht_n(:,:) * rn_gamma )   tau_icebfr(:,:) = tau_icebfr(:,:) + a_i(:,:,jl) * rn_icebfr 
    78          END DO 
    79       ENDIF 
    80        
    8169      ! ----------------------- 
    8270      ! Rheology (ice dynamics) 
    8371      ! -----------------------    
    84       IF( nn_icedyn /= 0 ) THEN     ! -- Ice dynamics 
    85          ! 
    86          CALL ice_rhg_evp( kt, stress1_i, stress2_i, stress12_i, u_ice, v_ice, shear_i, divu_i, delta_i ) 
    87          ! 
    88       ELSE                          ! -- prescribed uniform velocity 
    89          ! 
    90          u_ice(:,:) = rn_uice * umask(:,:,1) 
    91          v_ice(:,:) = rn_vice * vmask(:,:,1) 
    92          !!CALL RANDOM_NUMBER(u_ice(:,:)) 
    93          !!CALL RANDOM_NUMBER(v_ice(:,:)) 
    94          ! 
    95       ENDIF 
     72      CALL ice_rhg_evp( kt, stress1_i, stress2_i, stress12_i, u_ice, v_ice, shear_i, divu_i, delta_i ) 
    9673      ! 
    97       !                                                   !- conservation test 
     74      !                             !- conservation test 
    9875      IF( ln_icediachk   )   CALL ice_cons_hsm(1, 'icerhg', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    9976      IF( ln_ctl         )   CALL ice_prt3D  ('icerhg')   !- Control prints 
     
    11794      INTEGER ::   ios   ! Local integer output status for namelist read 
    11895      !! 
    119       NAMELIST/namice_rhg/  rn_ishlat  , rn_cio   , rn_creepl, rn_ecc , nn_nevp, rn_relast,   & 
    120          &                ln_landfast, rn_gamma , rn_icebfr, rn_lfrelax 
     96      NAMELIST/namice_rhg/  ln_rhg_EVP, rn_creepl, rn_ecc , nn_nevp, rn_relast 
    12197      !!------------------------------------------------------------------- 
    12298      ! 
     
    135111         WRITE(numout,*) '~~~~~~~~~~~~' 
    136112         WRITE(numout,*) '   Namelist namice_rhg' 
    137         WRITE(numout,*) '      rheology EVP (icerhg_evp)' 
    138          WRITE(numout,*) '         lateral boundary condition for sea ice dynamics          rn_ishlat     = ', rn_ishlat 
    139          WRITE(numout,*) '         drag coefficient for oceanic stress                      rn_cio        = ', rn_cio 
     113         WRITE(numout,*) '      rheology EVP (icerhg_evp)                                   ln_rhg_EVP    = ', ln_rhg_EVP 
    140114         WRITE(numout,*) '         creep limit                                              rn_creepl     = ', rn_creepl 
    141115         WRITE(numout,*) '         eccentricity of the elliptical yield curve               rn_ecc        = ', rn_ecc 
    142116         WRITE(numout,*) '         number of iterations for subcycling                      nn_nevp       = ', nn_nevp 
    143117         WRITE(numout,*) '         ratio of elastic timescale over ice time step            rn_relast     = ', rn_relast 
    144          WRITE(numout,*) '      Landfast: param (T or F)                                    ln_landfast   = ', ln_landfast 
    145          WRITE(numout,*) '         fraction of ocean depth that ice must reach              rn_gamma      = ', rn_gamma 
    146          WRITE(numout,*) '         maximum bottom stress per unit area of contact           rn_icebfr     = ', rn_icebfr 
    147          WRITE(numout,*) '         relax time scale (s-1) to reach static friction          rn_lfrelax    = ', rn_lfrelax 
    148118      ENDIF 
    149       ! 
    150       IF     (      rn_ishlat == 0.                ) THEN   ;   IF(lwp) WRITE(numout,*) '   ===>>>   ice lateral  free-slip' 
    151       ELSEIF (      rn_ishlat == 2.                ) THEN   ;   IF(lwp) WRITE(numout,*) '   ===>>>   ice lateral  no-slip' 
    152       ELSEIF ( 0. < rn_ishlat .AND. rn_ishlat < 2. ) THEN   ;   IF(lwp) WRITE(numout,*) '   ===>>>   ice lateral  partial-slip' 
    153       ELSEIF ( 2. < rn_ishlat                      ) THEN   ;   IF(lwp) WRITE(numout,*) '   ===>>>   ice lateral  strong-slip' 
    154       ENDIF 
    155       ! 
    156       IF( .NOT. ln_landfast )   tau_icebfr(:,:) = 0._wp     ! NO Landfast ice : set to zero one for all 
    157119      ! 
    158120   END SUBROUTINE ice_rhg_init 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

    r8515 r8516  
    3737   USE icerhg         ! sea-ice: rheology 
    3838   USE iceadv         ! sea-ice: advection 
     39   USE icedyn         ! sea-ice: dynamics 
    3940   USE icethd         ! sea-ice: thermodynamics 
    4041   USE icerdgrft      ! sea-ice: ridging/rafting 
     
    144145         CALL ice_rst_opn( kt )     ! Open Ice restart file (if necessary)  
    145146         ! 
    146          SELECT CASE( nice_dyn )    
    147          CASE ( np_dynFULL )          !==  all dynamical processes  ==! 
    148             CALL ice_rhg   ( kt )          ! -- rheology   
    149             CALL ice_adv   ( kt )          ! -- advection of ice 
    150             CALL ice_rdgrft( kt )          ! -- ridging/rafting  
    151             CALL ice_cor   ( kt , 1 )      ! -- Corrections 
    152          CASE ( np_dyn )              !==  pure dynamics only ==!   (no ridging/rafting)   (nono cat. case 2) 
    153             CALL ice_rhg   ( kt )          ! -- rheology   
    154             CALL ice_adv   ( kt )          ! -- advection of ice 
    155             CALL ice_cor   ( kt , 1 )      ! -- Corrections 
    156          CASE ( np_dynPURE )          !==  pure dynamics only ==!   (nn_icedyn= 0 or 1 ) 
    157             CALL ice_rhg   ( kt )          ! -- rheology   
    158             CALL ice_adv   ( kt )          ! -- advection of ice 
    159          END SELECT 
     147         IF( ln_icedyn .AND. .NOT.lk_c1d )   CALL ice_dyn( kt )            ! -- Ice dynamics 
    160148         ! 
    161149         !                          !==  lateral boundary conditions  ==! 
    162150#if defined key_agrif 
    163          IF( .NOT. Agrif_Root()     )   CALL agrif_interp_lim3('T')   ! -- AGRIF  
     151         IF( .NOT. Agrif_Root()     )        CALL agrif_interp_lim3('T')   ! -- AGRIF  
    164152#endif 
    165          IF( ln_icethd .AND. ln_bdy )   CALL bdy_ice( kt )            ! -- bdy ice thermo 
     153         IF( ln_icethd .AND. ln_bdy )        CALL bdy_ice( kt )            ! -- bdy ice thermo 
    166154         ! 
    167155         ! 
     
    186174         !    qprec_ice, qevap_ice 
    187175         !------------------------------------------------------! 
    188          CALL ice_forcing_flx( kt, ksbc ) 
     176                                    CALL ice_forcing_flx( kt, ksbc ) 
    189177 
    190178         !----------------------------! 
     
    269257      ! 
    270258      CALL ice_itd_init                ! ice thickness distribution initialization 
     259      ! 
     260      CALL ice_dyn_init                ! set ice dynamics parameters 
    271261      ! 
    272262      CALL ice_rdgrft_init             ! set ice ridging/rafting parameters 
Note: See TracChangeset for help on using the changeset viewer.