Ignore:
Timestamp:
2017-09-15T20:07:33+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part6 - more clarity (still not finished)

File:
1 edited

Legend:

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

    r8518 r8531  
    2525   USE lib_mpp        ! MPP library 
    2626   USE in_out_manager ! I/O manager 
     27   USE iom            ! I/O manager 
    2728   USE lib_fortran    ! glob_sum 
    2829   USE timing         ! Timing 
     
    209210      !!      dynamics 
    210211      !! 
    211       !! ** Method  :  Read the namice_dyn namelist and check the ice-dynamic 
     212      !! ** Method  :  Read the namdyn namelist and check the ice-dynamic 
    212213      !!       parameter values called at the first timestep (nit000) 
    213214      !! 
    214       !! ** input   :   Namelist namice_dyn 
     215      !! ** input   :   Namelist namdyn 
    215216      !!------------------------------------------------------------------- 
    216217      INTEGER ::   ios, ioptio   ! Local integer output status for namelist read 
    217218      !! 
    218       NAMELIST/namice_dyn/ ln_dynFULL, ln_dynRHGADV, ln_dynADV, rn_uice, rn_vice,  & 
    219          &                 rn_ishlat  , rn_cio   ,                         & 
    220          &                 ln_landfast, rn_gamma , rn_icebfr, rn_lfrelax 
    221       !!------------------------------------------------------------------- 
    222       ! 
    223       REWIND( numnam_ice_ref )         ! Namelist namice_dyn in reference namelist : Ice dynamics 
    224       READ  ( numnam_ice_ref, namice_dyn, IOSTAT = ios, ERR = 901) 
    225 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_dyn in reference namelist', lwp ) 
    226       ! 
    227       REWIND( numnam_ice_cfg )         ! Namelist namice_dyn in configuration namelist : Ice dynamics 
    228       READ  ( numnam_ice_cfg, namice_dyn, IOSTAT = ios, ERR = 902 ) 
    229 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_dyn in configuration namelist', lwp ) 
    230       IF(lwm) WRITE ( numoni, namice_dyn ) 
     219      NAMELIST/namdyn/ ln_dynFULL, ln_dynRHGADV, ln_dynADV, rn_uice, rn_vice,  & 
     220         &             rn_ishlat  ,                                            & 
     221         &             ln_landfast, rn_gamma , rn_icebfr, rn_lfrelax 
     222      !!------------------------------------------------------------------- 
     223      ! 
     224      REWIND( numnam_ice_ref )         ! Namelist namdyn in reference namelist : Ice dynamics 
     225      READ  ( numnam_ice_ref, namdyn, IOSTAT = ios, ERR = 901) 
     226901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn in reference namelist', lwp ) 
     227      ! 
     228      REWIND( numnam_ice_cfg )         ! Namelist namdyn in configuration namelist : Ice dynamics 
     229      READ  ( numnam_ice_cfg, namdyn, IOSTAT = ios, ERR = 902 ) 
     230902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn in configuration namelist', lwp ) 
     231      IF(lwm) WRITE ( numoni, namdyn ) 
    231232      ! 
    232233      IF(lwp) THEN                     ! control print 
     
    234235         WRITE(numout,*) 'ice_dyn_init: ice parameters for ice dynamics ' 
    235236         WRITE(numout,*) '~~~~~~~~~~~~' 
    236          WRITE(numout,*) '   Namelist namice_dyn' 
     237         WRITE(numout,*) '   Namelist namdyn:' 
    237238         WRITE(numout,*) '      Full ice dynamics      (rhg + adv + ridge/raft + corr)  ln_dynFULL   = ', ln_dynFULL 
    238239         WRITE(numout,*) '      No ridge/raft & No cor (rhg + adv)                      ln_dynRHGADV = ', ln_dynRHGADV 
     
    240241         WRITE(numout,*) '           with prescribed velocity given by ' 
    241242         WRITE(numout,*) '               a uniform field               (u,v)_ice = (rn_uice,rn_vice) = (', rn_uice,',', rn_vice,')' 
    242          WRITE(numout,*) '      lateral boundary condition for sea ice dynamics        rn_ishlat     = ', rn_ishlat 
    243          WRITE(numout,*) '      drag coefficient for oceanic stress                    rn_cio        = ', rn_cio 
    244          WRITE(numout,*) '      Landfast: param (T or F)                               ln_landfast   = ', ln_landfast 
    245          WRITE(numout,*) '         fraction of ocean depth that ice must reach         rn_gamma      = ', rn_gamma 
    246          WRITE(numout,*) '         maximum bottom stress per unit area of contact      rn_icebfr     = ', rn_icebfr 
    247          WRITE(numout,*) '         relax time scale (s-1) to reach static friction     rn_lfrelax    = ', rn_lfrelax 
     243         WRITE(numout,*) '      lateral boundary condition for sea ice dynamics         rn_ishlat    = ', rn_ishlat 
     244         WRITE(numout,*) '      Landfast: param (T or F)                                ln_landfast  = ', ln_landfast 
     245         WRITE(numout,*) '         fraction of ocean depth that ice must reach          rn_gamma     = ', rn_gamma 
     246         WRITE(numout,*) '         maximum bottom stress per unit area of contact       rn_icebfr    = ', rn_icebfr 
     247         WRITE(numout,*) '         relax time scale (s-1) to reach static friction      rn_lfrelax   = ', rn_lfrelax 
    248248      ENDIF 
    249249      !                             !== set the choice of ice dynamics ==! 
     
    267267      ENDIF 
    268268      !                                      !--- NO Landfast ice : set to zero once for all 
    269       IF( .NOT. ln_landfast )   tau_icebfr(:,:) = 0._wp  
     269      IF( .NOT. ln_landfast )   tau_icebfr(:,:) = 0._wp 
     270      ! 
     271      CALL ice_rdgrft_init          ! set ice ridging/rafting parameters 
     272      CALL ice_rhg_init             ! set ice rheology parameters 
     273      CALL ice_adv_init             ! set ice advection parameters 
    270274      ! 
    271275   END SUBROUTINE ice_dyn_init 
Note: See TracChangeset for help on using the changeset viewer.