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 8506 – NEMO

Changeset 8506


Ignore:
Timestamp:
2017-09-07T12:28:42+02:00 (7 years ago)
Author:
clem
Message:

changes in style - part5 - continue changing init routines

Location:
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3
Files:
4 edited

Legend:

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

    r8505 r8506  
    176176   INTEGER , PUBLIC ::   jiceprt          !: debug j-point 
    177177 
    178    !                                     !!** ice-init namelist (namiceini) ** 
    179                                           ! -- iceist -- ! 
    180    LOGICAL , PUBLIC ::   ln_limini        ! initialization or not 
    181    LOGICAL , PUBLIC ::   ln_limini_file   ! Ice initialization state from 2D netcdf file 
    182    REAL(wp), PUBLIC ::   rn_thres_sst     ! threshold water temperature for initial sea ice 
    183    REAL(wp), PUBLIC ::   rn_hts_ini_n     ! initial snow thickness in the north 
    184    REAL(wp), PUBLIC ::   rn_hts_ini_s     ! initial snow thickness in the south 
    185    REAL(wp), PUBLIC ::   rn_hti_ini_n     ! initial ice thickness in the north 
    186    REAL(wp), PUBLIC ::   rn_hti_ini_s     ! initial ice thickness in the south 
    187    REAL(wp), PUBLIC ::   rn_ati_ini_n     ! initial leads area in the north 
    188    REAL(wp), PUBLIC ::   rn_ati_ini_s     ! initial leads area in the south 
    189    REAL(wp), PUBLIC ::   rn_smi_ini_n     ! initial salinity  
    190    REAL(wp), PUBLIC ::   rn_smi_ini_s     ! initial salinity 
    191    REAL(wp), PUBLIC ::   rn_tmi_ini_n     ! initial temperature 
    192    REAL(wp), PUBLIC ::   rn_tmi_ini_s     ! initial temperature 
    193     
    194178   !                                     !!** ice-dynamics namelist (namicedyn) ** 
    195179                                          ! -- iceadv -- ! 
     
    252236   REAL(wp), PUBLIC ::   rn_simin         !: minimum ice salinity [PSU] 
    253237 
    254    !                                     !!** ice-mechanical redistribution namelist (namiceitdme) 
    255    REAL(wp), PUBLIC ::   rn_cs            !: fraction of shearing energy contributing to ridging             
    256    INTEGER , PUBLIC ::   nn_partfun       !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007) 
    257    REAL(wp), PUBLIC ::   rn_gstar         !: fractional area of young ice contributing to ridging 
    258    REAL(wp), PUBLIC ::   rn_astar         !: equivalent of G* for an exponential participation function 
    259    LOGICAL , PUBLIC ::   ln_ridging       !: ridging of ice or not                         
    260    REAL(wp), PUBLIC ::   rn_hstar         !: thickness that determines the maximal thickness of ridged ice 
    261    REAL(wp), PUBLIC ::   rn_por_rdg       !: initial porosity of ridges (0.3 regular value) 
    262    REAL(wp), PUBLIC ::   rn_fsnowrdg      !: fractional snow loss to the ocean during ridging 
    263    REAL(wp), PUBLIC ::   rn_fpondrdg      !: fractional melt pond loss to the ocean during ridging 
    264    LOGICAL , PUBLIC ::   ln_rafting       !: rafting of ice or not                         
    265    REAL(wp), PUBLIC ::   rn_hraft         !: threshold thickness (m) for rafting / ridging  
    266    REAL(wp), PUBLIC ::   rn_craft         !: coefficient for smoothness of the hyperbolic tangent in rafting 
    267    REAL(wp), PUBLIC ::   rn_fsnowrft      !: fractional snow loss to the ocean during ridging 
    268    REAL(wp), PUBLIC ::   rn_fpondrft      !: fractional snow loss to the ocean during rafting 
    269  
    270238   ! MV MP 2016 
    271239   !                                     !!** melt pond namelist (namicemp) 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceist.F90

    r8498 r8506  
    4646   INTEGER , PARAMETER ::   jp_smi = 6           ! index of ice sali at T-point 
    4747   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   si  ! structure of input fields (file informations, fields read) 
     48   ! 
     49   ! ** namelist (namiceini) ** 
     50   LOGICAL  ::   ln_limini        ! initialization or not 
     51   LOGICAL  ::   ln_limini_file   ! Ice initialization state from 2D netcdf file 
     52   REAL(wp) ::   rn_thres_sst     ! threshold water temperature for initial sea ice 
     53   REAL(wp) ::   rn_hts_ini_n     ! initial snow thickness in the north 
     54   REAL(wp) ::   rn_hts_ini_s     ! initial snow thickness in the south 
     55   REAL(wp) ::   rn_hti_ini_n     ! initial ice thickness in the north 
     56   REAL(wp) ::   rn_hti_ini_s     ! initial ice thickness in the south 
     57   REAL(wp) ::   rn_ati_ini_n     ! initial leads area in the north 
     58   REAL(wp) ::   rn_ati_ini_s     ! initial leads area in the south 
     59   REAL(wp) ::   rn_smi_ini_n     ! initial salinity  
     60   REAL(wp) ::   rn_smi_ini_s     ! initial salinity 
     61   REAL(wp) ::   rn_tmi_ini_n     ! initial temperature 
     62   REAL(wp) ::   rn_tmi_ini_s     ! initial temperature 
     63    
    4864   !!---------------------------------------------------------------------- 
    4965   !! NEMO/ICE 4.0 , NEMO Consortium (2017) 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerdgrft.F90

    r8500 r8506  
    3636   PUBLIC   ice_rdgrft_alloc         ! called by ice_init  
    3737 
    38    !----------------------------------------------------------------------- 
    3938   ! Variables shared among ridging subroutines 
    40    !----------------------------------------------------------------------- 
    4139   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   asum     ! sum of total ice and open water area 
    4240   REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   aksum    ! ratio of area removed to area ridged 
     
    5149   REAL(wp), PARAMETER ::   krdgmin = 1.1_wp    ! min ridge thickness multiplier 
    5250   REAL(wp), PARAMETER ::   kraft   = 0.5_wp    ! rafting multipliyer 
    53    ! 
    5451   REAL(wp)            ::   zdrho               !  
    5552   ! 
     53   ! ** namelist (namiceitdme) ** 
     54   REAL(wp) ::   rn_cs            !: fraction of shearing energy contributing to ridging             
     55   INTEGER  ::   nn_partfun       !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007) 
     56   REAL(wp) ::   rn_gstar         !: fractional area of young ice contributing to ridging 
     57   REAL(wp) ::   rn_astar         !: equivalent of G* for an exponential participation function 
     58   LOGICAL  ::   ln_ridging       !: ridging of ice or not                         
     59   REAL(wp) ::   rn_hstar         !: thickness that determines the maximal thickness of ridged ice 
     60   REAL(wp) ::   rn_por_rdg       !: initial porosity of ridges (0.3 regular value) 
     61   REAL(wp) ::   rn_fsnowrdg      !: fractional snow loss to the ocean during ridging 
     62   REAL(wp) ::   rn_fpondrdg      !: fractional melt pond loss to the ocean during ridging 
     63   LOGICAL  ::   ln_rafting       !: rafting of ice or not                         
     64   REAL(wp) ::   rn_hraft         !: threshold thickness (m) for rafting / ridging  
     65   REAL(wp) ::   rn_craft         !: coefficient for smoothness of the hyperbolic tangent in rafting 
     66   REAL(wp) ::   rn_fsnowrft      !: fractional snow loss to the ocean during ridging 
     67   REAL(wp) ::   rn_fpondrft      !: fractional snow loss to the ocean during rafting 
    5668   ! 
    5769   !!---------------------------------------------------------------------- 
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90

    r8498 r8506  
    448448         END DO 
    449449         ! 
    450          ! Change thickness to volume 
     450         ! Change thickness to volume (replaces routine ice_var_eqv2glo) 
    451451         v_i_1d(1:nidx)   = ht_i_1d(1:nidx) * a_i_1d(1:nidx) 
    452452         v_s_1d(1:nidx)   = ht_s_1d(1:nidx) * a_i_1d(1:nidx) 
Note: See TracChangeset for help on using the changeset viewer.