Changeset 717 for trunk/NEMO/OPA_SRC/restart.F90
- Timestamp:
- 2007-10-16T13:03:55+02:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/restart.F90
r709 r717 19 19 USE phycst ! physical constants 20 20 USE daymod ! calendar 21 USE ice_oce ! ice variables22 21 USE cpl_oce, ONLY : lk_cpl ! 23 22 USE in_out_manager ! I/O manager … … 138 137 CALL iom_rstput( kt, nitrst, numrow, 'rotn' , rotn ) 139 138 CALL iom_rstput( kt, nitrst, numrow, 'hdivn' , hdivn ) 140 141 #if defined key_ice_lim142 CALL iom_rstput( kt, nitrst, numrow, 'nfice' , REAL( nfice, wp) ) ! ice computation frequency143 CALL iom_rstput( kt, nitrst, numrow, 'sst_io' , sst_io )144 CALL iom_rstput( kt, nitrst, numrow, 'sss_io' , sss_io )145 CALL iom_rstput( kt, nitrst, numrow, 'u_io' , u_io )146 CALL iom_rstput( kt, nitrst, numrow, 'v_io' , v_io )147 # if defined key_coupled148 CALL iom_rstput( kt, nitrst, numrow, 'alb_ice', alb_ice )149 # endif150 #endif151 #if defined key_flx_bulk_monthly || defined key_flx_bulk_daily || defined key_flx_core152 CALL iom_rstput( kt, nitrst, numrow, 'nfbulk' , REAL( nfbulk, wp) ) ! bulk computation frequency153 CALL iom_rstput( kt, nitrst, numrow, 'gsst' , gsst )154 #endif155 139 156 140 IF( nn_dynhpg_rst == 1 .OR. lk_vvl ) THEN … … 204 188 !! has been stored in the restart file. 205 189 !!---------------------------------------------------------------------- 206 REAL(wp) :: zcoef, zkt, zrdt, zrdttra1, zndastp , znfice, znfbulk190 REAL(wp) :: zcoef, zkt, zrdt, zrdttra1, zndastp 207 191 #if defined key_ice_lim 208 192 INTEGER :: ji, jj … … 299 283 ENDIF 300 284 301 !!sm: TO BE MOVED IN NEW SURFACE MODULE...302 303 #if defined key_ice_lim304 ! Louvain La Neuve Sea Ice Model305 IF( iom_varid( numror, 'nfice' ) > 0 ) then306 CALL iom_get( numror , 'nfice' , znfice ) ! ice computation frequency307 CALL iom_get( numror, jpdom_autoglo, 'sst_io' , sst_io )308 CALL iom_get( numror, jpdom_autoglo, 'sss_io' , sss_io )309 CALL iom_get( numror, jpdom_autoglo, 'u_io' , u_io )310 CALL iom_get( numror, jpdom_autoglo, 'v_io' , v_io )311 # if defined key_coupled312 CALL iom_get( numror, jpdom_autoglo, 'alb_ice', alb_ice )313 # endif314 IF( znfice /= REAL( nfice, wp ) ) THEN ! if nfice changed between 2 runs315 zcoef = REAL( nfice-1, wp ) / znfice316 sst_io(:,:) = zcoef * sst_io(:,:)317 sss_io(:,:) = zcoef * sss_io(:,:)318 u_io (:,:) = zcoef * u_io (:,:)319 v_io (:,:) = zcoef * v_io (:,:)320 ENDIF321 ELSE322 IF(lwp) WRITE(numout,*)323 IF(lwp) WRITE(numout,*) 'rst_read : LLN sea Ice Model => Ice initialization'324 IF(lwp) WRITE(numout,*)325 zcoef = REAL( nfice-1, wp )326 sst_io(:,:) = zcoef *( tn(:,:,1) + rt0 ) !!bug a explanation is needed here!327 sss_io(:,:) = zcoef * sn(:,:,1)328 zcoef = 0.5 * REAL( nfice-1, wp )329 DO jj = 2, jpj330 DO ji = fs_2, jpi ! vector opt.331 u_io(ji,jj) = zcoef * ( un(ji-1,jj ,1) + un(ji-1,jj-1,1) )332 v_io(ji,jj) = zcoef * ( vn(ji ,jj-1,1) + vn(ji-1,jj-1,1) )333 END DO334 END DO335 # if defined key_coupled336 alb_ice(:,:) = 0.8 * tmask(:,:,1)337 # endif338 ENDIF339 #endif340 #if defined key_flx_bulk_monthly || defined key_flx_bulk_daily || defined key_flx_core341 ! Louvain La Neuve Sea Ice Model342 IF( iom_varid( numror, 'nfbulk' ) > 0 ) THEN343 CALL iom_get( numror , 'nfbulk', znfbulk ) ! bulk computation frequency344 CALL iom_get( numror, jpdom_autoglo, 'gsst' , gsst )345 IF( znfbulk /= REAL(nfbulk, wp) ) THEN ! if you change nfbulk between 2 runs346 zcoef = REAL( nfbulk-1, wp ) / znfbulk347 gsst(:,:) = zcoef * gsst(:,:)348 ENDIF349 ELSE350 IF(lwp) WRITE(numout,*)351 IF(lwp) WRITE(numout,*) 'rst_read : LLN sea Ice Model => Ice initialization'352 IF(lwp) WRITE(numout,*)353 gsst(:,:) = REAL( nfbulk - 1, wp )*( tn(:,:,1) + rt0 )354 ENDIF355 #endif356 357 !!sm: end of TO BE MOVED IN NEW SURFACE MODULE...358 359 285 IF( iom_varid( numror, 'rhd' ) > 0 ) THEN 360 286 CALL iom_get( numror, jpdom_autoglo, 'rhd' , rhd )
Note: See TracChangeset
for help on using the changeset viewer.