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/icestp.F90

    r8518 r8531  
    142142         ! --- ice dynamics and advection  --- ! 
    143143         !-------------------------------------! 
    144          CALL ice_diag0             ! set diag of mass, heat and salt fluxes to 0 
     144         CALL diag_set0             ! set diag of mass, heat and salt fluxes to 0 
    145145         CALL ice_rst_opn( kt )     ! Open Ice restart file (if necessary)  
    146146         ! 
     
    243243      IF(lwm) CALL ctl_opn( numoni, 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 ) 
    244244      ! 
    245       CALL ice_run_init                ! set some ice run parameters 
     245      CALL par_init                ! set some ice run parameters 
    246246      ! 
    247247      !                                ! Allocate the ice arrays (sbc_ice already allocated in sbc_init) 
     
    268268      CALL ice_var_glo2eqv 
    269269      ! 
     270      CALL ice_forcing_init            ! set ice-ocean and ice-atm. coupling parameters 
     271      ! 
    270272      IF( ln_icedyn ) THEN 
    271273         CALL ice_dyn_init             ! set ice dynamics parameters 
    272          CALL ice_rdgrft_init          ! set ice ridging/rafting parameters 
    273          CALL ice_rhg_init             ! set ice rheology parameters 
    274          CALL ice_adv_init             ! set ice advection parameters 
    275274      ENDIF 
    276  
     275      ! 
    277276      IF( ln_icethd ) THEN 
    278277         CALL ice_thd_init             ! set ice thermodynics parameters 
    279          CALL ice_thd_sal_init         ! set ice salinity parameters 
    280278      ENDIF    
    281279      ! 
     
    299297 
    300298 
    301    SUBROUTINE ice_run_init 
     299   SUBROUTINE par_init 
    302300      !!------------------------------------------------------------------- 
    303       !!                  ***  ROUTINE ice_run_init *** 
     301      !!                  ***  ROUTINE par_init *** 
    304302      !! 
    305303      !! ** Purpose :   Definition some run parameter for ice model 
    306304      !! 
    307       !! ** Method  :   Read the namice_run namelist and check the parameter 
     305      !! ** Method  :   Read the nampar namelist and check the parameter 
    308306      !!                values called at the first timestep (nit000) 
    309307      !! 
    310       !! ** input   :   Namelist namice_run 
     308      !! ** input   :   Namelist nampar 
    311309      !!------------------------------------------------------------------- 
    312310      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    313       NAMELIST/namice_run/ jpl, nlay_i, nlay_s, nn_monocat, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s,  & 
    314          &                 cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir 
     311      NAMELIST/nampar/ jpl, nlay_i, nlay_s, nn_monocat, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s,  & 
     312         &             cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir 
    315313      !!------------------------------------------------------------------- 
    316314      ! 
    317       REWIND( numnam_ice_ref )      ! Namelist namice_run in reference namelist : Parameters for ice 
    318       READ  ( numnam_ice_ref, namice_run, IOSTAT = ios, ERR = 901) 
    319 901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_run in reference namelist', lwp ) 
    320  
    321       REWIND( numnam_ice_cfg )      ! Namelist namice_run in configuration namelist : Parameters for ice 
    322       READ  ( numnam_ice_cfg, namice_run, IOSTAT = ios, ERR = 902 ) 
    323 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_run in configuration namelist', lwp ) 
    324       IF(lwm) WRITE ( numoni, namice_run ) 
     315      REWIND( numnam_ice_ref )      ! Namelist nampar in reference namelist : Parameters for ice 
     316      READ  ( numnam_ice_ref, nampar, IOSTAT = ios, ERR = 901) 
     317901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in reference namelist', lwp ) 
     318 
     319      REWIND( numnam_ice_cfg )      ! Namelist nampar in configuration namelist : Parameters for ice 
     320      READ  ( numnam_ice_cfg, nampar, IOSTAT = ios, ERR = 902 ) 
     321902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampar in configuration namelist', lwp ) 
     322      IF(lwm) WRITE ( numoni, nampar ) 
    325323      ! 
    326324      IF(lwp) THEN                  ! control print 
    327325         WRITE(numout,*) 
    328          WRITE(numout,*) 'ice_run_init : ice share parameters for dynamics/advection/thermo of sea-ice' 
    329          WRITE(numout,*) ' ~~~~~~' 
    330          WRITE(numout,*) '   Namelist namice_run : ' 
     326         WRITE(numout,*) 'par_init: ice parameters shared among all the routines' 
     327         WRITE(numout,*) ' ~~~~~~~' 
     328         WRITE(numout,*) '   Namelist nampar: ' 
    331329         WRITE(numout,*) '      number of ice  categories                              jpl    = ', jpl 
    332330         WRITE(numout,*) '      number of ice  layers                                  nlay_i = ', nlay_i 
    333331         WRITE(numout,*) '      number of snow layers                                  nlay_s = ', nlay_s 
    334332         WRITE(numout,*) '      virtual ITD mono-category param (1-4) or not (0)   nn_monocat = ', nn_monocat 
    335          WRITE(numout,*) '      Ice dynamics       (T) or not (F)                  ln_icedyn = ', ln_icedyn 
    336          WRITE(numout,*) '      Ice thermodynamics (T) or not (F)                  ln_icethd = ', ln_icethd 
     333         WRITE(numout,*) '      Ice dynamics       (T) or not (F)                   ln_icedyn = ', ln_icedyn 
     334         WRITE(numout,*) '      Ice thermodynamics (T) or not (F)                   ln_icethd = ', ln_icethd 
    337335         WRITE(numout,*) '      maximum ice concentration for NH                              = ', rn_amax_n  
    338336         WRITE(numout,*) '      maximum ice concentration for SH                              = ', rn_amax_s 
     
    346344      ENDIF 
    347345      IF ( jpl == 1 .AND. nn_monocat == 0 ) THEN 
    348          CALL ctl_stop( 'STOP', 'ice_run_init : if jpl=1 then nn_monocat should be between 1 and 4' ) 
     346         CALL ctl_stop( 'STOP', 'par_init : if jpl=1 then nn_monocat should be between 1 and 4' ) 
    349347      ENDIF 
    350348      ! 
    351       IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('ice_run_init: online conservation check does not work with BDY') 
     349      IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('par_init: online conservation check does not work with BDY') 
    352350      ! 
    353351      rdt_ice   = REAL(nn_fsbc) * rdt          !--- sea-ice timestep and inverse 
     
    358356      r1_nlay_s = 1._wp / REAL( nlay_s, wp ) 
    359357      ! 
    360    END SUBROUTINE ice_run_init 
     358   END SUBROUTINE par_init 
    361359 
    362360 
     
    393391 
    394392 
    395    SUBROUTINE ice_diag0 
    396       !!---------------------------------------------------------------------- 
    397       !!                  ***  ROUTINE ice_diag0  *** 
     393   SUBROUTINE diag_set0 
     394      !!---------------------------------------------------------------------- 
     395      !!                  ***  ROUTINE diag_set0  *** 
    398396      !! 
    399397      !! ** purpose :  set ice-ocean and ice-atm. fluxes to zeros at the beggining 
     
    442440      tau_icebfr(:,:) = 0._wp; ! landfast ice param only (clem: important to keep the init here) 
    443441       
    444    END SUBROUTINE ice_diag0 
     442   END SUBROUTINE diag_set0 
    445443 
    446444#else 
Note: See TracChangeset for help on using the changeset viewer.