New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10555 for NEMO/trunk/src – NEMO

Changeset 10555 for NEMO/trunk/src


Ignore:
Timestamp:
2019-01-22T12:36:01+01:00 (5 years ago)
Author:
clem
Message:

choose to update ice velocity every sub-iteration in rheology or not (temporary logical has been added in icedyn_rhg_evp.F90). For now it has been set to false, so there is no update at each sub-iteration

Location:
NEMO/trunk/src/ICE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icedyn_rhg_evp.F90

    r10425 r10555  
    112112      REAL(wp), DIMENSION(:,:), INTENT(  out) ::   pshear_i  , pdivu_i   , pdelta_i      ! 
    113113      !! 
     114      LOGICAL, PARAMETER ::   ll_bdy_substep = .FALSE. ! temporary option to call bdy at each sub-time step (T) 
     115      !                                                                              or only at the main time step (F) 
    114116      INTEGER ::   ji, jj       ! dummy loop indices 
    115117      INTEGER ::   jter         ! local integers 
     
    535537            CALL agrif_interp_ice( 'V' ) 
    536538#endif 
    537             IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
     539            IF( ln_bdy .AND. ll_bdy_substep ) CALL bdy_ice_dyn( 'V' ) 
    538540            ! 
    539541            DO jj = 2, jpjm1 
     
    583585            CALL agrif_interp_ice( 'U' ) 
    584586#endif 
    585             IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
     587            IF( ln_bdy .AND. ll_bdy_substep ) CALL bdy_ice_dyn( 'U' ) 
    586588            ! 
    587589         ELSE ! odd iterations 
     
    633635            CALL agrif_interp_ice( 'U' ) 
    634636#endif 
    635             IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
     637            IF( ln_bdy .AND. ll_bdy_substep ) CALL bdy_ice_dyn( 'U' ) 
    636638            ! 
    637639            DO jj = 2, jpjm1 
     
    681683            CALL agrif_interp_ice( 'V' ) 
    682684#endif 
    683             IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
     685            IF( ln_bdy .AND. ll_bdy_substep ) CALL bdy_ice_dyn( 'V' ) 
    684686            ! 
    685687         ENDIF 
     
    696698      END DO                                              !  end loop over jter  ! 
    697699      !                                                   ! ==================== ! 
     700      ! 
     701      IF( ln_bdy .AND. .NOT.ll_bdy_substep ) THEN 
     702         CALL bdy_ice_dyn( 'U' ) 
     703         CALL bdy_ice_dyn( 'V' ) 
     704      ENDIF 
    698705      ! 
    699706      !------------------------------------------------------------------------------! 
  • NEMO/trunk/src/ICE/icevar.F90

    r10415 r10555  
    567567               DO ji = 1 , jpi 
    568568                  IF( pe_i(ji,jj,jk,jl) < 0._wp .OR. pa_i(ji,jj,jl) < 0._wp ) THEN 
    569                      hfx_res(ji,jj)   = hfx_res(ji,jj) - pe_i(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0 
     569                     hfx_res(ji,jj)   = hfx_res(ji,jj) - pe_i(ji,jj,jk,jl) * r1_rdtice ! W.m-2 >0 
    570570                     pe_i(ji,jj,jk,jl) = 0._wp 
    571571                  ENDIF 
Note: See TracChangeset for help on using the changeset viewer.