- Timestamp:
- 2017-09-08T18:19:17+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90
r8514 r8515 36 36 PRIVATE 37 37 38 PUBLIC ice_istate ! called by icestp.F90 38 PUBLIC ice_istate ! called by icestp.F90 39 PUBLIC ice_istate_init ! called by icestp.F90 39 40 40 41 INTEGER , PARAMETER :: jpfldi = 6 ! maximum number of files to read … … 102 103 103 104 IF(lwp) WRITE(numout,*) 104 IF(lwp) WRITE(numout,*) 'ice_istate 105 IF(lwp) WRITE(numout,*) '~~~~~~~~~~ 105 IF(lwp) WRITE(numout,*) 'ice_istate: sea-ice initialization ' 106 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 106 107 107 108 !-------------------------------------------------------------------- … … 109 110 !-------------------------------------------------------------------- 110 111 ! 111 CALL ice_istate_init112 113 112 ! init surface temperature 114 113 DO jl = 1, jpl … … 154 153 !----------------------------- 155 154 ! assign initial thickness, concentration, snow depth and salinity to an hemisphere-dependent array 156 DO jj = 1, jpj 157 DO ji = 1, jpi 158 IF( ff_t(ji,jj) >= 0._wp ) THEN 159 zht_i_ini(ji,jj) = rn_hti_ini_n * zswitch(ji,jj) 160 zht_s_ini(ji,jj) = rn_hts_ini_n * zswitch(ji,jj) 161 zat_i_ini(ji,jj) = rn_ati_ini_n * zswitch(ji,jj) 162 zts_u_ini(ji,jj) = rn_tmi_ini_n * zswitch(ji,jj) 163 zsm_i_ini(ji,jj) = rn_smi_ini_n * zswitch(ji,jj) 164 ztm_i_ini(ji,jj) = rn_tmi_ini_n * zswitch(ji,jj) 165 ELSE 166 zht_i_ini(ji,jj) = rn_hti_ini_s * zswitch(ji,jj) 167 zht_s_ini(ji,jj) = rn_hts_ini_s * zswitch(ji,jj) 168 zat_i_ini(ji,jj) = rn_ati_ini_s * zswitch(ji,jj) 169 zts_u_ini(ji,jj) = rn_tmi_ini_s * zswitch(ji,jj) 170 zsm_i_ini(ji,jj) = rn_smi_ini_s * zswitch(ji,jj) 171 ztm_i_ini(ji,jj) = rn_tmi_ini_s * zswitch(ji,jj) 172 ENDIF 173 END DO 174 END DO 155 WHERE( ff_t(:,:) >= 0._wp ) 156 zht_i_ini(:,:) = rn_hti_ini_n * zswitch(:,:) 157 zht_s_ini(:,:) = rn_hts_ini_n * zswitch(:,:) 158 zat_i_ini(:,:) = rn_ati_ini_n * zswitch(:,:) 159 zts_u_ini(:,:) = rn_tmi_ini_n * zswitch(:,:) 160 zsm_i_ini(:,:) = rn_smi_ini_n * zswitch(:,:) 161 ztm_i_ini(:,:) = rn_tmi_ini_n * zswitch(:,:) 162 ELSEWHERE 163 zht_i_ini(:,:) = rn_hti_ini_s * zswitch(:,:) 164 zht_s_ini(:,:) = rn_hts_ini_s * zswitch(:,:) 165 zat_i_ini(:,:) = rn_ati_ini_s * zswitch(:,:) 166 zts_u_ini(:,:) = rn_tmi_ini_s * zswitch(:,:) 167 zsm_i_ini(:,:) = rn_smi_ini_s * zswitch(:,:) 168 ztm_i_ini(:,:) = rn_tmi_ini_s * zswitch(:,:) 169 END WHERE 175 170 ! 176 171 ENDIF ! ln_iceini_file … … 554 549 slf_i(jp_ati) = sn_ati ; slf_i(jp_tsu) = sn_tsu 555 550 slf_i(jp_tmi) = sn_tmi ; slf_i(jp_smi) = sn_smi 556 557 ! Define the initial parameters 558 ! ------------------------- 559 560 IF(lwp) THEN 551 ! 552 ! 553 IF(lwp) THEN ! control print 561 554 WRITE(numout,*) 562 555 WRITE(numout,*) 'ice_istate_init: ice parameters inititialisation ' 563 556 WRITE(numout,*) '~~~~~~~~~~~~~~~' 564 557 WRITE(numout,*) ' Namelist namice_ini' 565 WRITE(numout,*) ' initialization with ice (T) or not (F) ln_iceini = ', ln_iceini566 WRITE(numout,*) ' ice initialization from a netcdf file ln_iceini_file = ', ln_iceini_file567 WRITE(numout,*) ' threshold water temp. for initial sea-icern_thres_sst = ', rn_thres_sst568 WRITE(numout,*) ' initial snow thickness in the north rn_hts_ini_n = ', rn_hts_ini_n569 WRITE(numout,*) ' initial snow thickness in the south rn_hts_ini_s = ', rn_hts_ini_s570 WRITE(numout,*) ' initial ice thickness in the north rn_hti_ini_n = ', rn_hti_ini_n571 WRITE(numout,*) ' initial ice thickness in the south rn_hti_ini_s = ', rn_hti_ini_s572 WRITE(numout,*) ' initial ice concentr. in the north rn_ati_ini_n = ', rn_ati_ini_n573 WRITE(numout,*) ' initial ice concentr. in the north rn_ati_ini_s = ', rn_ati_ini_s574 WRITE(numout,*) ' initial ice salinity in the north rn_smi_ini_n = ', rn_smi_ini_n575 WRITE(numout,*) ' initial ice salinity in the south rn_smi_ini_s = ', rn_smi_ini_s576 WRITE(numout,*) ' initial ice/snw temp in the north rn_tmi_ini_n = ', rn_tmi_ini_n577 WRITE(numout,*) ' initial ice/snw temp in the south rn_tmi_ini_s = ', rn_tmi_ini_s558 WRITE(numout,*) ' initialization with ice (T) or not (F) ln_iceini = ', ln_iceini 559 WRITE(numout,*) ' ice initialization from a netcdf file ln_iceini_file = ', ln_iceini_file 560 WRITE(numout,*) ' max delta ocean temp. above Tfreeze with initial ice rn_thres_sst = ', rn_thres_sst 561 WRITE(numout,*) ' initial snow thickness in the north rn_hts_ini_n = ', rn_hts_ini_n 562 WRITE(numout,*) ' initial snow thickness in the south rn_hts_ini_s = ', rn_hts_ini_s 563 WRITE(numout,*) ' initial ice thickness in the north rn_hti_ini_n = ', rn_hti_ini_n 564 WRITE(numout,*) ' initial ice thickness in the south rn_hti_ini_s = ', rn_hti_ini_s 565 WRITE(numout,*) ' initial ice concentr. in the north rn_ati_ini_n = ', rn_ati_ini_n 566 WRITE(numout,*) ' initial ice concentr. in the north rn_ati_ini_s = ', rn_ati_ini_s 567 WRITE(numout,*) ' initial ice salinity in the north rn_smi_ini_n = ', rn_smi_ini_n 568 WRITE(numout,*) ' initial ice salinity in the south rn_smi_ini_s = ', rn_smi_ini_s 569 WRITE(numout,*) ' initial ice/snw temp in the north rn_tmi_ini_n = ', rn_tmi_ini_n 570 WRITE(numout,*) ' initial ice/snw temp in the south rn_tmi_ini_s = ', rn_tmi_ini_s 578 571 ENDIF 579 572
Note: See TracChangeset
for help on using the changeset viewer.