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 4298 for branches – NEMO

Changeset 4298 for branches


Ignore:
Timestamp:
2013-11-20T19:24:53+01:00 (10 years ago)
Author:
clem
Message:
 
Location:
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90

    r4293 r4298  
    7171         &     'use more ocean levels or less ice/snow layers/categories.' ) 
    7272 
    73       !                                ! Open the namelist file  
    74       CALL ctl_opn( numnam_ice, 'namelist_ice', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 
     73                                       ! Open the reference and configuration namelist files and namelist output file  
     74      CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )  
     75      CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 
     76      CALL ctl_opn( numoni,         'output.namelist.ice', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 
    7577      ! 
    7678      CALL ice_run                     ! set some ice run parameters 
     
    129131      !!------------------------------------------------------------------- 
    130132      NAMELIST/namicerun/ cn_icerst_in, cn_icerst_out, ln_limdyn, amax, cai, cao, ln_nicep, ln_limdiahsb, ln_limdiaout 
     133      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    131134      !!------------------------------------------------------------------- 
    132135      !                     
    133       REWIND( numnam_ice )                ! Read Namelist namicerun  
    134       READ  ( numnam_ice , namicerun ) 
    135       ! 
    136       !IF( lk_mpp .AND. ln_nicep ) THEN 
    137       !   ln_nicep = .FALSE. 
    138       !   CALL ctl_warn( 'ice_run : specific control print for LIM3 desactivated with MPI' ) 
    139       !ENDIF 
     136      REWIND( numnam_ice_ref )              ! Namelist namicerun in reference namelist : Parameters for ice 
     137      READ  ( numnam_ice_ref, namicerun, IOSTAT = ios, ERR = 901) 
     138901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in reference namelist', lwp ) 
     139 
     140      REWIND( numnam_ice_cfg )              ! Namelist namicerun in configuration namelist : Parameters for ice 
     141      READ  ( numnam_ice_cfg, namicerun, IOSTAT = ios, ERR = 902 ) 
     142902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in configuration namelist', lwp ) 
     143      WRITE ( numoni, namicerun ) 
     144      ! 
     145      IF( lk_mpp .AND. ln_nicep ) THEN 
     146         ln_nicep = .FALSE. 
     147         CALL ctl_warn( 'ice_run : specific control print for LIM3 desactivated with MPI' ) 
     148      ENDIF 
    140149      ! 
    141150      IF(lwp) THEN                        ! control print 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90

    r4269 r4298  
    5353      !-------------------------------------------------------------------- 
    5454      ijpij = SIZE(zhti,1) 
    55       zht_i(1:ijpij,1:jpl) = 0.d0 
    56       zht_s(1:ijpij,1:jpl) = 0.d0 
    57       za_i (1:ijpij,1:jpl) = 0.d0 
     55      zht_i(1:ijpij,1:jpl) = 0._wp 
     56      zht_s(1:ijpij,1:jpl) = 0._wp 
     57      za_i (1:ijpij,1:jpl) = 0._wp 
    5858 
    5959      !------------------------------------------------------------------------------------ 
     
    101101             
    102102            ! initialisation of ice variables for each try 
    103             zht_i(ji,1:jpl) = 0.d0 
    104             za_i (ji,1:jpl) = 0.d0 
     103            zht_i(ji,1:jpl) = 0._wp 
     104            za_i (ji,1:jpl) = 0._wp 
    105105             
    106106            ! *** case very thin ice: fill only category 1 
     
    114114               ! Fill ice thicknesses except the last one (i_fill) by (hmax-hmin)/2  
    115115               DO jl = 1, i_fill - 1 
    116                   zht_i(ji,jl) = ( zhi_max(jl) + zhi_max(jl-1) ) / 2. 
     116                  zht_i(ji,jl) = ( zhi_max(jl) + zhi_max(jl-1) ) * 0.5_wp 
    117117               END DO 
    118118                
     
    130130               DO jl = 1, i_fill - 1 
    131131                  IF ( jl == jl0 ) CYCLE 
    132                   zarg           = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) / 2. ) 
     132                  zarg           = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) * 0.5_wp ) 
    133133                  za_i(ji,jl) =   za_i (ji,jl0) * EXP(-zarg**2) 
    134134               END DO 
     
    160160            ztest_4 = 1 
    161161            DO jl = 1, i_fill 
    162                IF ( za_i(ji,jl) < 0.0d0 ) ztest_4 = 0 
     162               IF ( za_i(ji,jl) < 0._wp ) ztest_4 = 0 
    163163            END DO 
    164164             
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r4293 r4298  
    14771477      !! ** input   :   Namelist namiceitdme 
    14781478      !!------------------------------------------------------------------- 
     1479      INTEGER :: ios                 ! Local integer output status for namelist read 
    14791480      NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft,&  
    14801481         Gstar, astar,                                & 
     
    14841485      !!------------------------------------------------------------------- 
    14851486      ! 
    1486       REWIND( numnam_ice )                   ! read namiceitdme namelist 
    1487       READ  ( numnam_ice , namiceitdme) 
     1487      REWIND( numnam_ice_ref )              ! Namelist namicetdme in reference namelist : Ice mechanical ice redistribution 
     1488      READ  ( numnam_ice_ref, namiceitdme, IOSTAT = ios, ERR = 901) 
     1489901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitdme in reference namelist', lwp ) 
     1490 
     1491      REWIND( numnam_ice_cfg )              ! Namelist namiceitdme in configuration namelist : Ice mechanical ice redistribution 
     1492      READ  ( numnam_ice_cfg, namiceitdme, IOSTAT = ios, ERR = 902 ) 
     1493902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitdme in configuration namelist', lwp ) 
     1494      WRITE ( numoni, namiceitdme ) 
    14881495      ! 
    14891496      IF (lwp) THEN                          ! control print 
Note: See TracChangeset for help on using the changeset viewer.