- Timestamp:
- 2017-09-08T17:09:25+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90
r8500 r8514 191 191 o_i(:,:,:) = oa_i(:,:,:) * z1_a_i(:,:,:) !--- ice age 192 192 193 IF( nn_icesal == 2 ) THEN!--- salinity (with a minimum value imposed everywhere)193 IF( nn_icesal == 2 ) THEN !--- salinity (with a minimum value imposed everywhere) 194 194 WHERE( v_i(:,:,:) > epsi20 ) ; sm_i(:,:,:) = MAX( rn_simin , smv_i(:,:,:) / v_i(:,:,:) ) 195 195 ELSEWHERE ; sm_i(:,:,:) = rn_simin … … 289 289 !!------------------------------------------------------------------ 290 290 291 !!gm much much more secure to defined when reading nn_icesal in the namelist integers =1, 2, 3 with explicit names292 !! for example np_Scst_noProfile = 1 ; np_Svar_linProfile = 2 ; np_Scst_fixProfile293 294 291 !!gm Question: Remove the option 3 ? How many years since it last use ? 295 292 296 293 SELECT CASE ( nn_icesal ) 297 294 ! 298 ! !---------------------------------------!299 CASE( 1 ) ! constant salinity in time and space !300 ! !---------------------------------------!295 ! !---------------------------------------! 296 CASE( 1 ) ! constant salinity in time and space ! 297 ! !---------------------------------------! 301 298 s_i (:,:,:,:) = rn_icesal 302 299 sm_i(:,:,:) = rn_icesal 303 300 ! 304 ! !---------------------------------------------!305 CASE( 2 ) ! time varying salinity with linear profile !306 ! !---------------------------------------------!301 ! !---------------------------------------------! 302 CASE( 2 ) ! time varying salinity with linear profile ! 303 ! !---------------------------------------------! 307 304 ! 308 305 ALLOCATE( z_slope_s(jpi,jpj,jpl) , zalpha(jpi,jpj,jpl) ) … … 343 340 DEALLOCATE( z_slope_s , zalpha ) 344 341 ! 345 ! !-------------------------------------------!346 CASE( 3 ) ! constant salinity with a fix profile ! (Schwarzacher (1959) multiyear salinity profile347 ! !-------------------------------------------! (mean = 2.30)342 ! !-------------------------------------------! 343 CASE( 3 ) ! constant salinity with a fix profile ! (Schwarzacher (1959) multiyear salinity profile 344 ! !-------------------------------------------! (mean = 2.30) 348 345 ! 349 346 sm_i(:,:,:) = 2.30_wp … … 415 412 SELECT CASE ( nn_icesal ) 416 413 ! 417 ! !---------------------------------------!418 CASE( 1 ) ! constant salinity in time and space !419 ! !---------------------------------------!414 ! !---------------------------------------! 415 CASE( 1 ) ! constant salinity in time and space ! 416 ! !---------------------------------------! 420 417 s_i_1d(:,:) = rn_icesal 421 418 ! 422 ! !---------------------------------------------!423 CASE( 2 ) ! time varying salinity with linear profile !424 ! !---------------------------------------------!419 ! !---------------------------------------------! 420 CASE( 2 ) ! time varying salinity with linear profile ! 421 ! !---------------------------------------------! 425 422 ! 426 423 ALLOCATE( z_slope_s(jpij) ) … … 445 442 DEALLOCATE( z_slope_s ) 446 443 447 ! !-------------------------------------------!448 CASE( 3 ) ! constant salinity with a fix profile ! (Schwarzacher (1959) multiyear salinity profile449 ! !-------------------------------------------! (mean = 2.30)444 ! !-------------------------------------------! 445 CASE( 3 ) ! constant salinity with a fix profile ! (Schwarzacher (1959) multiyear salinity profile 446 ! !-------------------------------------------! (mean = 2.30) 450 447 ! 451 448 sm_i_1d(:) = 2.30_wp … … 572 569 !! ice thickness distribution follows a gaussian law 573 570 !! around the concentration of the most likely ice thickness 574 !! (similar as iceist .F90)571 !! (similar as iceistate.F90) 575 572 !! 576 573 !! ** Method: Iterative procedure
Note: See TracChangeset
for help on using the changeset viewer.