- Timestamp:
- 2015-02-06T19:12:57+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90
r5064 r5067 154 154 155 155 ! Default new ice thickness 156 hicol(:,:) = hiccrit157 158 IF( fraz_swi== 1 ) THEN156 hicol(:,:) = rn_hnewice 157 158 IF( nn_frazil == 1 ) THEN 159 159 160 160 !-------------------- … … 175 175 !------------- 176 176 ! C-grid wind stress components 177 ztaux = ( utau_ice(ji-1,jj ) * tmu(ji-1,jj) &178 & + utau_ice(ji ,jj ) * tmu(ji ,jj) ) * 0.5_wp179 ztauy = ( vtau_ice(ji ,jj-1) * tmv(ji ,jj-1) &180 & + vtau_ice(ji ,jj ) * tmv(ji ,jj) ) * 0.5_wp177 ztaux = ( utau_ice(ji-1,jj ) * umask(ji-1,jj ,1) & 178 & + utau_ice(ji ,jj ) * umask(ji ,jj ,1) ) * 0.5_wp 179 ztauy = ( vtau_ice(ji ,jj-1) * vmask(ji ,jj-1,1) & 180 & + vtau_ice(ji ,jj ) * vmask(ji ,jj ,1) ) * 0.5_wp 181 181 ! Square root of wind stress 182 182 ztenagm = SQRT( SQRT( ztaux**2 + ztauy**2 ) ) … … 194 194 ! C-grid ice velocity 195 195 rswitch = MAX( 0._wp, SIGN( 1._wp , at_i(ji,jj) ) ) 196 zvgx = rswitch * ( u_ice(ji-1,jj ) * tmu(ji-1,jj ) + u_ice(ji,jj) * tmu(ji,jj) ) * 0.5_wp197 zvgy = rswitch * ( v_ice(ji ,jj-1) * tmv(ji ,jj-1) + v_ice(ji,jj) * tmv(ji,jj) ) * 0.5_wp196 zvgx = rswitch * ( u_ice(ji-1,jj ) * umask(ji-1,jj ,1) + u_ice(ji,jj) * umask(ji,jj,1) ) * 0.5_wp 197 zvgy = rswitch * ( v_ice(ji ,jj-1) * vmask(ji ,jj-1,1) + v_ice(ji,jj) * vmask(ji,jj,1) ) * 0.5_wp 198 198 199 199 !----------------------------------- … … 319 319 !---------------------- 320 320 DO ji = 1, nbpac 321 zh_newice(ji) = hiccrit322 END DO 323 IF( fraz_swi== 1 ) zh_newice(1:nbpac) = hicol_1d(1:nbpac)321 zh_newice(ji) = rn_hnewice 322 END DO 323 IF( nn_frazil == 1 ) zh_newice(1:nbpac) = hicol_1d(1:nbpac) 324 324 325 325 !---------------------- 326 326 ! Salinity of new ice 327 327 !---------------------- 328 SELECT CASE ( n um_sal )328 SELECT CASE ( nn_icesal ) 329 329 CASE ( 1 ) ! Sice = constant 330 zs_newice(1:nbpac) = bulk_sal330 zs_newice(1:nbpac) = rn_icesal 331 331 CASE ( 2 ) ! Sice = F(z,t) [Vancoppenolle et al (2005)] 332 332 DO ji = 1, nbpac 333 333 ii = MOD( npac(ji) - 1 , jpi ) + 1 334 334 ij = ( npac(ji) - 1 ) / jpi + 1 335 zs_newice(ji) = MIN( 4.606 + 0.91 / zh_newice(ji) , s_i_max , 0.5 * sss_m(ii,ij) )335 zs_newice(ji) = MIN( 4.606 + 0.91 / zh_newice(ji) , rn_simax , 0.5 * sss_m(ii,ij) ) 336 336 END DO 337 337 CASE ( 3 ) ! Sice = F(z) [multiyear ice] … … 386 386 ! A fraction zfrazb of frazil ice is accreted at the ice bottom 387 387 rswitch = 1._wp - MAX( 0._wp, SIGN( 1._wp , - zat_i_1d(ji) ) ) 388 zfrazb = rswitch * ( TANH ( Cfrazb * ( zvrel_1d(ji) - vfrazb ) ) + 1.0 ) * 0.5 *maxfrazb388 zfrazb = rswitch * ( TANH ( rn_Cfrazb * ( zvrel_1d(ji) - rn_vfrazb ) ) + 1.0 ) * 0.5 * rn_maxfrazb 389 389 zv_frazb(ji) = zfrazb * zv_newice(ji) 390 390 zv_newice(ji) = ( 1.0 - zfrazb ) * zv_newice(ji) … … 408 408 ! we keep the excessive volume in memory and attribute it later to bottom accretion 409 409 DO ji = 1, nbpac 410 IF ( za_newice(ji) > ( amax - zat_i_1d(ji) ) ) THEN411 zda_res(ji) = za_newice(ji) - ( amax - zat_i_1d(ji) )410 IF ( za_newice(ji) > ( rn_amax - zat_i_1d(ji) ) ) THEN 411 zda_res(ji) = za_newice(ji) - ( rn_amax - zat_i_1d(ji) ) 412 412 zdv_res(ji) = zda_res (ji) * zh_newice(ji) 413 413 za_newice(ji) = za_newice(ji) - zda_res (ji)
Note: See TracChangeset
for help on using the changeset viewer.