Changeset 5313 for branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90
- Timestamp:
- 2015-05-29T11:46:03+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90
r5312 r5313 18 18 USE ice ! sea-ice variables 19 19 USE oce , ONLY : snwice_mass, snwice_mass_b 20 USE par_ice ! sea-ice parameters21 20 USE dom_oce ! ocean domain 22 21 USE sbc_oce ! Surface boundary condition: ocean fields … … 27 26 USE wrk_nemo ! work arrays 28 27 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 28 USE limctl 29 29 30 30 IMPLICIT NONE … … 33 33 PUBLIC lim_rst_opn ! routine called by icestep.F90 34 34 PUBLIC lim_rst_write ! routine called by icestep.F90 35 PUBLIC lim_rst_read ! routine called by iceini.F9035 PUBLIC lim_rst_read ! routine called by sbc_lim_init 36 36 37 37 LOGICAL, PUBLIC :: lrst_ice !: logical to control the ice restart write … … 93 93 ENDIF 94 94 ! 95 IF( ln_icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - Beginning the time step - ' ) ! control print 95 96 END SUBROUTINE lim_rst_opn 96 97 … … 172 173 CALL iom_rstput( iter, nitrst, numriw, 'stress2_i' , stress2_i ) 173 174 CALL iom_rstput( iter, nitrst, numriw, 'stress12_i' , stress12_i ) 174 CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass' , snwice_mass ) !clem modif175 CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass_b', snwice_mass_b ) !clem modif175 CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass' , snwice_mass ) 176 CALL iom_rstput( iter, nitrst, numriw, 'snwice_mass_b', snwice_mass_b ) 176 177 177 178 DO jl = 1, jpl … … 313 314 !! ** purpose : read of sea-ice variable restart in a netcdf file 314 315 !!---------------------------------------------------------------------- 315 INTEGER :: ji, jj, jk, jl , indx316 INTEGER :: ji, jj, jk, jl 316 317 REAL(wp) :: zfice, ziter 317 REAL(wp) :: zs_inf, z_slope_s, zsmax, zsmin, zalpha ! local scalars used for the salinity profile318 REAL(wp), POINTER, DIMENSION(:) :: zs_zero319 318 REAL(wp), POINTER, DIMENSION(:,:) :: z2d 320 319 CHARACTER(len=15) :: znam … … 324 323 !!---------------------------------------------------------------------- 325 324 326 CALL wrk_alloc( nlay_i, zs_zero )327 325 CALL wrk_alloc( jpi, jpj, z2d ) 328 326 … … 402 400 CALL iom_get( numrir, jpdom_autoglo, 'stress2_i' , stress2_i ) 403 401 CALL iom_get( numrir, jpdom_autoglo, 'stress12_i', stress12_i ) 404 CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass' , snwice_mass ) !clem modif405 CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass_b', snwice_mass_b ) !clem modif402 CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass' , snwice_mass ) 403 CALL iom_get( numrir, jpdom_autoglo, 'snwice_mass_b', snwice_mass_b ) 406 404 407 405 DO jl = 1, jpl … … 528 526 ! 529 527 ! clem: I do not understand why the following IF is needed 530 ! I suspect something inconsistent in the main code with option n um_sal=1531 IF( n um_sal == 1 ) THEN528 ! I suspect something inconsistent in the main code with option nn_icesal=1 529 IF( nn_icesal == 1 ) THEN 532 530 DO jl = 1, jpl 533 sm_i(:,:,jl) = bulk_sal531 sm_i(:,:,jl) = rn_icesal 534 532 DO jk = 1, nlay_i 535 s_i(:,:,jk,jl) = bulk_sal533 s_i(:,:,jk,jl) = rn_icesal 536 534 END DO 537 535 END DO … … 540 538 !CALL iom_close( numrir ) !clem: closed in sbcice_lim.F90 541 539 ! 542 CALL wrk_dealloc( nlay_i, zs_zero )543 540 CALL wrk_dealloc( jpi, jpj, z2d ) 544 541 !
Note: See TracChangeset
for help on using the changeset viewer.