Ignore:
Timestamp:
2013-07-16T15:56:48+02:00 (7 years ago)
Author:
clem
Message:

change min thickness by min volume to set up u_ice=u_oce

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r3938 r3978  
    447447         CALL lbc_lnk( zs12(:,:), 'F', 1. ) 
    448448 
     449!#if defined key_bdy 
     450!         ! clem: change zs1, zs2, zs12 at the boundary for each iteration 
     451!         CALL bdy_ice_lim_dyn( 2, zs1, zs2, zs12 ) 
     452!         CALL lbc_lnk( zs1 (:,:), 'T', 1. ) 
     453!         CALL lbc_lnk( zs2 (:,:), 'T', 1. ) 
     454!         CALL lbc_lnk( zs12(:,:), 'F', 1. ) 
     455!#endif          
     456 
    449457         ! Ice internal stresses (Appendix C of Hunke and Dukowicz, 2002) 
    450458         DO jj = k_j1+1, k_jpj-1 
     
    578586!#if defined key_bdy 
    579587!         ! clem: change u_ice and v_ice at the boundary for each iteration 
    580 !         CALL bdy_ice_lim_dyn() 
     588!         CALL bdy_ice_lim_dyn( 1 ) 
    581589!#endif          
    582590 
     
    607615!CDIR NOVERRCHK 
    608616         DO ji = fs_2, fs_jpim1 
    609             zindb  = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )  
    610             zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 
     617            !zindb  = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )  
     618            !zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 
     619            zdummy = vt_i(ji,jj) 
    611620            IF ( zdummy .LE. hminrhg ) THEN 
    612621               u_ice(ji,jj) = u_oce(ji,jj) 
     
    619628      CALL lbc_lnk( v_ice(:,:), 'V', -1. )  
    620629 
     630#if defined key_bdy 
    621631      ! clem: change u_ice and v_ice at the boundary 
    622 #if defined key_bdy 
    623       CALL bdy_ice_lim_dyn() 
     632      CALL bdy_ice_lim_dyn( 1 ) 
    624633#endif          
    625634 
    626635      DO jj = k_j1+1, k_jpj-1  
    627636         DO ji = fs_2, fs_jpim1 
    628             zindb  = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )  
    629             zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 
     637            !zindb  = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )  
     638            !zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 
     639            zdummy = vt_i(ji,jj) 
    630640            IF ( zdummy .LE. hminrhg ) THEN 
    631641               v_ice1(ji,jj)  = 0.5*( (v_ice(ji,jj)+v_ice(ji,jj-1))*e1t(ji+1,jj)   & 
     
    650660            !- zdd(:,:), zdt(:,:): divergence and tension at centre  
    651661            !- zds(:,:): shear on northeast corner of grid cells 
    652             zindb  = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )  
    653             zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 
    654  
     662            !zindb  = MAX( 0.0, SIGN( 1.0, at_i(ji,jj) - 1.0e-6 ) )  
     663            !zdummy = zindb * vt_i(ji,jj) / MAX(at_i(ji,jj) , 1.0e-06 ) 
     664            zdummy = vt_i(ji,jj) 
    655665            IF ( zdummy .LE. hminrhg ) THEN 
    656666 
Note: See TracChangeset for help on using the changeset viewer.