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

Changeset 15367


Ignore:
Timestamp:
2021-10-14T10:00:58+02:00 (2 years ago)
Author:
davestorkey
Message:

r4.0-HEAD : fix for wind forcing of icebergs #2728

Location:
NEMO/releases/r4.0/r4.0-HEAD/src/OCE
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/ICB/icbini.F90

    r13350 r15367  
    7272      IF( .NOT. ln_icebergs )   RETURN 
    7373 
     74      ALLOCATE( utau_icb(jpi,jpj) , vtau_icb(jpi,jpj) ) 
    7475      !                          ! allocate gridded fields 
    7576      IF( icb_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'icb_alloc : unable to allocate arrays' ) 

  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/ICB/icbutl.F90

    r14372 r15367  
    7676      ss_e(1:jpi,1:jpj) = sss_m(:,:) 
    7777      fr_e(1:jpi,1:jpj) = fr_i (:,:) 
    78       ua_e(1:jpi,1:jpj) = utau (:,:) * umask(:,:,1) ! maybe mask useless because mask applied in sbcblk 
    79       va_e(1:jpi,1:jpj) = vtau (:,:) * vmask(:,:,1) ! maybe mask useless because mask applied in sbcblk 
     78      ua_e(1:jpi,1:jpj) = utau_icb (:,:) * umask(:,:,1) ! maybe mask useless because mask applied in sbcblk 
     79      va_e(1:jpi,1:jpj) = vtau_icb (:,:) * vmask(:,:,1) ! maybe mask useless because mask applied in sbcblk 
    8080      ! 
    8181      CALL lbc_lnk_icb( 'icbutl', uo_e, 'U', -1._wp, 1, 1 ) 

  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/SBC/sbc_oce.F90

    r13284 r15367  
    111111   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau   , utau_b   !: sea surface i-stress (ocean referential)     [N/m2] 
    112112   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vtau   , vtau_b   !: sea surface j-stress (ocean referential)     [N/m2] 
     113   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau_icb, vtau_icb !: sea surface (i,j)-stress used by icebergs   [N/m2] 
    113114   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   taum              !: module of sea surface stress (at T-point)    [N/m2]  
    114115   !! wndm is used onmpute surface gases exchanges in ice-free ocean or leads 

  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/SBC/sbcmod.F90

    r13481 r15367  
    448448      IF( ln_mixcpl )          CALL sbc_cpl_rcv   ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing 
    449449      ! 
    450       IF ( ln_wave .AND. (ln_tauwoc .OR. ln_tauw) ) CALL sbc_wstress( )      ! Wind stress provided by waves  
     450      IF ( ln_wave .AND. (ln_tauwoc .OR. ln_tauw) ) CALL sbc_wstress( )       ! Wind stress provided by waves  
     451      ! 
     452      IF( ln_icebergs ) utau_icb(:,:) = utau(:,:) ; vtau_icb(:,:) = vtau(:,:) ! save pure wind stresses (with no ice-ocean stress) to be used by icebergs 
    451453      ! 
    452454      !                                            !==  Misc. Options  ==! 
Note: See TracChangeset for help on using the changeset viewer.