- Timestamp:
- 2016-10-05T09:43:42+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90
r6853 r6989 114 114 CHARACTER (len=50) :: charout 115 115 116 REAL(wp) :: zrhoco ! rau0 * rn_cio 116 117 REAL(wp) :: zdtevp, z1_dtevp ! time step for subcycling 117 118 REAL(wp) :: ecc2, z1_ecc2 ! square of yield ellipse eccenticity … … 203 204 ! 1) define some variables and initialize arrays 204 205 !------------------------------------------------------------------------------! 206 zrhoco = rau0 * rn_cio 207 205 208 ! ecc2: square of yield ellipse eccenticrity 206 209 ecc2 = rn_ecc * rn_ecc … … 357 360 358 361 ! delta at T points 359 zdelta = SQRT( zdiv2 + ( zdt2 + zds2 ) * usecc2 )362 zdelta = SQRT( zdiv2 + ( zdt2 + zds2 ) * z1_ecc2 ) 360 363 361 364 ! P/delta at T points … … 424 427 425 428 ! tau_io/(v_oce - v_ice) 426 zTauO = zaV(ji,jj) * rhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) ) &427 & + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) )429 zTauO = zaV(ji,jj) * zrhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) ) & 430 & + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) ) 428 431 429 432 ! tau_bottom/v_ice … … 469 472 470 473 ! tau_io/(u_oce - u_ice) 471 zTauO = zaU(ji,jj) * rhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) ) &472 & + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) )474 zTauO = zaU(ji,jj) * zrhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) ) & 475 & + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) ) 473 476 474 477 ! tau_bottom/u_ice … … 515 518 516 519 ! tau_io/(u_oce - u_ice) 517 zTauO = zaU(ji,jj) * rhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) ) &518 & + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) )520 zTauO = zaU(ji,jj) * zrhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) ) & 521 & + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) ) 519 522 520 523 ! tau_bottom/u_ice … … 559 562 560 563 ! tau_io/(v_oce - v_ice) 561 zTauO = zaV(ji,jj) * rhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) ) &562 & + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) )564 zTauO = zaV(ji,jj) * zrhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) ) & 565 & + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) ) 563 566 564 567 ! tau_bottom/v_ice … … 651 654 652 655 ! delta at T points 653 zdelta = SQRT( divu_i(ji,jj) * divu_i(ji,jj) + ( zdt2 + zds2 ) * usecc2 )656 zdelta = SQRT( divu_i(ji,jj) * divu_i(ji,jj) + ( zdt2 + zds2 ) * z1_ecc2 ) 654 657 rswitch = 1._wp - MAX( 0._wp, SIGN( 1._wp, -zdelta ) ) ! 0 if delta=0 655 658 delta_i(ji,jj) = zdelta + rn_creepl * rswitch
Note: See TracChangeset
for help on using the changeset viewer.