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 921 for trunk/NEMO/LIM_SRC_3/iceini.F90 – NEMO

Ignore:
Timestamp:
2008-05-13T10:28:52+02:00 (16 years ago)
Author:
rblod
Message:

Correct indentation and print for debug in LIM3, see ticket #134, step I

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/iceini.F90

    r888 r921  
    3232 
    3333   !! * Share Module variables 
    34    LOGICAL , PUBLIC  ::     & !!! ** init namelist (namicerun) ** 
    35       ln_limdyn   = .TRUE., & !: flag for ice dynamics (T) or not (F) 
    36       ln_nicep    = .TRUE.    !: flag for sea-ice points output (T) or not (F) 
    3734   INTEGER , PUBLIC  ::   &  !: 
    3835      nstart ,            &  !: iteration number of the begining of the run  
     
    4138      numit                  !: iteration number 
    4239   REAL(wp), PUBLIC  ::   &  !: 
    43       hsndif = 0.e0    ,  &  !: computation of temp. in snow (0) or not (9999) 
    44       hicdif = 0.e0    ,  &  !: computation of temp. in ice (0) or not (9999) 
    45       tpstot           ,  &  !: time of the run in seconds 
    46       cai    = 1.40e-3 ,  &  !: atmospheric drag over sea ice 
    47       cao    = 1.00e-3       !: atmospheric drag over ocean 
    48    REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !: 
    49       acrit  = (/ 1.e-06 , 1.e-06 /)    !: minimum fraction for leads in  
    50       !                                   !  north and south hemisphere 
     40      tpstot                 !: time of the run in seconds 
    5141   !!---------------------------------------------------------------------- 
    5242   !!   LIM 3.0,  UCL-ASTR-LOCEAN-IPSL (2008)  
     
    7262 
    7363      CALL ice_run                    !  read in namelist some run parameters 
    74                   
     64 
    7565      ! Louvain la Neuve Ice model 
    7666      IF( nacc == 1 ) THEN 
    77           dtsd2   = nn_fsbc * rdtmin * 0.5 
    78           rdt_ice = nn_fsbc * rdtmin 
     67         dtsd2   = nn_fsbc * rdtmin * 0.5 
     68         rdt_ice = nn_fsbc * rdtmin 
    7969      ELSE 
    80           dtsd2   = nn_fsbc * rdt * 0.5 
    81           rdt_ice = nn_fsbc * rdt 
     70         dtsd2   = nn_fsbc * rdt * 0.5 
     71         rdt_ice = nn_fsbc * rdt 
    8272      ENDIF 
    8373 
    8474      CALL lim_msh                    ! ice mesh initialization 
    85       
     75 
    8676      CALL lim_itd_ini                ! initialize the ice thickness 
    87                                       ! distribution 
     77      ! distribution 
    8878      ! Initial sea-ice state 
    8979      IF( .NOT.ln_rstart ) THEN 
     
    9282         CALL lim_istate              ! start from rest: sea-ice deduced from sst 
    9383         CALL lim_var_agg(1)          ! aggregate category variables in 
    94                                       ! bulk variables 
     84         ! bulk variables 
    9585         CALL lim_var_glo2eqv         ! convert global variables in equivalent 
    96                                       ! variables 
     86         ! variables 
    9787      ELSE 
    9888         CALL lim_rst_read            ! start from a restart file 
     
    10898      alb_ice(:,:,:) = albege(:,:)      ! sea-ice albedo 
    10999# endif 
    110        
     100 
    111101      nstart = numit  + nn_fsbc       
    112102      nitrun = nitend - nit000 + 1  
     
    138128      REWIND ( numnam_ice ) 
    139129      READ   ( numnam_ice , namicerun ) 
     130      ln_nicep = ln_nicep .AND. lwp 
    140131      IF(lwp) THEN 
    141132         WRITE(numout,*) 
     
    150141         WRITE(numout,*) '   Several ice points in the ice or not in ocean.output = ', ln_nicep 
    151142      ENDIF 
    152        
     143 
    153144   END SUBROUTINE ice_run 
    154145 
    155146   SUBROUTINE lim_itd_ini 
    156         !!------------------------------------------------------------------ 
    157         !!                ***  ROUTINE lim_itd_ini *** 
    158         !! ** Purpose : 
    159         !!            Initializes the ice thickness distribution 
    160         !! ** Method  : 
    161         !!            Very simple. Currently there are no ice types in the 
    162         !!            model... 
    163         !! 
    164         !! ** Arguments : 
    165         !!           kideb , kiut : Starting and ending points on which the 
    166         !!                         the computation is applied 
    167         !! 
    168         !! ** Inputs / Ouputs : (global commons) 
    169         !! 
    170         !! ** External : 
    171         !! 
    172         !! ** References : 
    173         !! 
    174         !! ** History : 
    175         !!           (12-2005) Martin Vancoppenolle 
    176         !! 
    177         !!------------------------------------------------------------------ 
    178         !! * Arguments 
    179  
    180        !! * Local variables 
    181        INTEGER ::   jl,       &   ! ice category dummy loop index 
    182                     jm            ! ice types    dummy loop index 
    183  
    184        REAL(wp)  ::           &  ! constant values 
    185           zeps      =  1.0e-10,   & ! 
    186           zc1                 ,   & ! 
    187           zc2                 ,   & ! 
    188           zc3                 ,   & ! 
    189           zx1 
    190  
    191        WRITE(numout,*) 'lim_itd_ini : Initialization of ice thickness distribution ' 
    192        WRITE(numout,*) '~~~~~~~~~~~~' 
    193  
    194 !!-- End of declarations 
    195 !!------------------------------------------------------------------------------ 
    196  
    197 !------------------------------------------------------------------------------! 
    198 ! 1) Ice thickness distribution parameters initialization     
    199 !------------------------------------------------------------------------------! 
     147      !!------------------------------------------------------------------ 
     148      !!                ***  ROUTINE lim_itd_ini *** 
     149      !! ** Purpose : 
     150      !!            Initializes the ice thickness distribution 
     151      !! ** Method  : 
     152      !!            Very simple. Currently there are no ice types in the 
     153      !!            model... 
     154      !! 
     155      !! ** Arguments : 
     156      !!           kideb , kiut : Starting and ending points on which the 
     157      !!                         the computation is applied 
     158      !! 
     159      !! ** Inputs / Ouputs : (global commons) 
     160      !! 
     161      !! ** External : 
     162      !! 
     163      !! ** References : 
     164      !! 
     165      !! ** History : 
     166      !!           (12-2005) Martin Vancoppenolle 
     167      !! 
     168      !!------------------------------------------------------------------ 
     169      !! * Arguments 
     170 
     171      !! * Local variables 
     172      INTEGER ::   jl,       &   ! ice category dummy loop index 
     173         jm            ! ice types    dummy loop index 
     174 
     175      REAL(wp)  ::           &  ! constant values 
     176         zeps      =  1.0e-10,   & ! 
     177         zc1                 ,   & ! 
     178         zc2                 ,   & ! 
     179         zc3                 ,   & ! 
     180         zx1 
     181 
     182      WRITE(numout,*) 'lim_itd_ini : Initialization of ice thickness distribution ' 
     183      WRITE(numout,*) '~~~~~~~~~~~~' 
     184 
     185      !!-- End of declarations 
     186      !!------------------------------------------------------------------------------ 
     187 
     188      !------------------------------------------------------------------------------! 
     189      ! 1) Ice thickness distribution parameters initialization     
     190      !------------------------------------------------------------------------------! 
    200191 
    201192      !- Types boundaries (integer) 
     
    266257      tn_ice(:,:,:) = t_su(:,:,:) 
    267258 
    268     END SUBROUTINE lim_itd_ini 
     259   END SUBROUTINE lim_itd_ini 
    269260 
    270261#else 
Note: See TracChangeset for help on using the changeset viewer.