- Timestamp:
- 2017-09-07T20:08:11+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceist.F90
r8506 r8512 47 47 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: si ! structure of input fields (file informations, fields read) 48 48 ! 49 ! ** namelist (namice ini) **50 LOGICAL :: ln_ limini ! initialization or not51 LOGICAL :: ln_ limini_file ! Ice initialization state from 2D netcdf file49 ! ** namelist (namice_ini) ** 50 LOGICAL :: ln_iceini ! initialization or not 51 LOGICAL :: ln_iceini_file ! Ice initialization state from 2D netcdf file 52 52 REAL(wp) :: rn_thres_sst ! threshold water temperature for initial sea ice 53 53 REAL(wp) :: rn_hts_ini_n ! initial snow thickness in the north … … 126 126 ! 2) Initialization of sea ice state variables 127 127 !-------------------------------------------------------------------- 128 IF( ln_ limini ) THEN128 IF( ln_iceini ) THEN 129 129 ! 130 IF( ln_ limini_file )THEN130 IF( ln_iceini_file )THEN 131 131 ! 132 132 zht_i_ini(:,:) = si(jp_hti)%fnow(:,:,1) … … 141 141 END WHERE 142 142 ! 143 ELSE ! ln_ limini_file = F143 ELSE ! ln_iceini_file = F 144 144 145 145 !-------------------------------------------------------------------- … … 175 175 END DO 176 176 ! 177 ENDIF ! ln_ limini_file177 ENDIF ! ln_iceini_file 178 178 179 179 zvt_i_ini(:,:) = zht_i_ini(:,:) * zat_i_ini(:,:) ! ice volume … … 403 403 ! END MV MP 2016 404 404 405 ELSE ! if ln_ limini=false405 ELSE ! if ln_iceini=false 406 406 a_i (:,:,:) = 0._wp 407 407 v_i (:,:,:) = 0._wp … … 431 431 h_ip (:,:,:) = 0._wp 432 432 433 ENDIF ! ln_ limini433 ENDIF ! ln_iceini 434 434 435 435 at_i (:,:) = 0.0_wp … … 518 518 !! ** Purpose : Definition of initial state of the ice 519 519 !! 520 !! ** Method : Read the namice ini namelist and check the parameter520 !! ** Method : Read the namice_ini namelist and check the parameter 521 521 !! values called at the first timestep (nit000) 522 522 !! 523 523 !! ** input : 524 !! Namelist namice ini524 !! Namelist namice_ini 525 525 !! 526 526 !! history : … … 537 537 TYPE(FLD_N), DIMENSION(jpfldi) :: slf_i ! array of namelist informations on the fields to read 538 538 ! 539 NAMELIST/namice ini/ ln_limini, ln_limini_file, rn_thres_sst, rn_hts_ini_n, rn_hts_ini_s, &539 NAMELIST/namice_ini/ ln_iceini, ln_iceini_file, rn_thres_sst, rn_hts_ini_n, rn_hts_ini_s, & 540 540 & rn_hti_ini_n, rn_hti_ini_s, rn_ati_ini_n, rn_ati_ini_s, rn_smi_ini_n, & 541 541 & rn_smi_ini_s, rn_tmi_ini_n, rn_tmi_ini_s, & … … 543 543 !!----------------------------------------------------------------------------- 544 544 ! 545 REWIND( numnam_ice_ref ) ! Namelist namice ini in reference namelist : Ice initial state546 READ ( numnam_ice_ref, namice ini, IOSTAT = ios, ERR = 901)547 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice ini in reference namelist', lwp )548 549 REWIND( numnam_ice_cfg ) ! Namelist namice ini in configuration namelist : Ice initial state550 READ ( numnam_ice_cfg, namice ini, IOSTAT = ios, ERR = 902 )551 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice ini in configuration namelist', lwp )552 IF(lwm) WRITE ( numoni, namice ini )545 REWIND( numnam_ice_ref ) ! Namelist namice_ini in reference namelist : Ice initial state 546 READ ( numnam_ice_ref, namice_ini, IOSTAT = ios, ERR = 901) 547 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_ini in reference namelist', lwp ) 548 549 REWIND( numnam_ice_cfg ) ! Namelist namice_ini in configuration namelist : Ice initial state 550 READ ( numnam_ice_cfg, namice_ini, IOSTAT = ios, ERR = 902 ) 551 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_ini in configuration namelist', lwp ) 552 IF(lwm) WRITE ( numoni, namice_ini ) 553 553 554 554 slf_i(jp_hti) = sn_hti ; slf_i(jp_hts) = sn_hts … … 563 563 WRITE(numout,*) 'ice_ist_init : ice parameters inititialisation ' 564 564 WRITE(numout,*) '~~~~~~~~~~~~' 565 WRITE(numout,*) ' Namelist namice ini'566 WRITE(numout,*) ' initialization with ice (T) or not (F) ln_ limini = ', ln_limini567 WRITE(numout,*) ' ice initialization from a netcdf file ln_ limini_file = ', ln_limini_file565 WRITE(numout,*) ' Namelist namice_ini' 566 WRITE(numout,*) ' initialization with ice (T) or not (F) ln_iceini = ', ln_iceini 567 WRITE(numout,*) ' ice initialization from a netcdf file ln_iceini_file = ', ln_iceini_file 568 568 WRITE(numout,*) ' threshold water temp. for initial sea-ice rn_thres_sst = ', rn_thres_sst 569 569 WRITE(numout,*) ' initial snow thickness in the north rn_hts_ini_n = ', rn_hts_ini_n … … 579 579 ENDIF 580 580 581 IF( ln_ limini_file ) THEN ! Ice initialization using input file581 IF( ln_iceini_file ) THEN ! Ice initialization using input file 582 582 ! 583 583 ! set si structure
Note: See TracChangeset
for help on using the changeset viewer.