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 9660 – NEMO

Changeset 9660


Ignore:
Timestamp:
2018-05-28T14:09:12+02:00 (6 years ago)
Author:
clem
Message:

remove boundary conditions for ice dynamics. It makes the ice rheology numerically unstable. However the effect of this removal on sea ice boundary conditions is still unclear

Location:
NEMO/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/ICE/icedyn_rhg_evp.F90

    r9656 r9660  
    498498            CALL agrif_interp_ice( 'V' ) 
    499499#endif 
    500             IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
     500!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
     501!            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
    501502            ! 
    502503            DO jj = 2, jpjm1 
    … …  
    546547            CALL agrif_interp_ice( 'U' ) 
    547548#endif 
    548             IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
     549!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
     550!            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
    549551            ! 
    550552         ELSE ! odd iterations 
    … …  
    596598            CALL agrif_interp_ice( 'U' ) 
    597599#endif 
    598             IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
     600!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
     601!            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
    599602            ! 
    600603            DO jj = 2, jpjm1 
    … …  
    644647            CALL agrif_interp_ice( 'V' ) 
    645648#endif 
    646             IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
     649!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
     650!            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
    647651            ! 
    648652         ENDIF 

  • NEMO/trunk/tests/SAS_BIPER/EXP00/1_namelist_ice_cfg

    r9019 r9660  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! ESIM namelist:   
     2!! SI3 namelist:   
    33!!              1 - Generic parameters                 (nampar) 
    44!!              2 - Ice thickness discretization       (namitd) 

  • NEMO/trunk/tests/SAS_BIPER/EXP00/namelist_ice_cfg

    r9019 r9660  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! ESIM namelist:   
     2!! SI3 namelist:   
    33!!              1 - Generic parameters                 (nampar) 
    44!!              2 - Ice thickness discretization       (namitd) 
Note: See TracChangeset for help on using the changeset viewer.