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

Changeset 4067


Ignore:
Timestamp:
2013-10-16T16:44:36+02:00 (10 years ago)
Author:
smasson
Message:

dev_v3_4_STABLE_2012: bugfix for mixed oce-ice coupling, see #1110 and #1137

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r3915 r4067  
    382382      ! 
    383383      IF( TRIM( sn_rcv_tau%cldes ) /= 'oce and ice' ) THEN        ! 'oce and ice' case ocean stress on ocean mesh used 
    384          srcv(jpr_itz1:jpr_itz2)%laction = .FALSE.    ! ice components not received (itx1 and ity1 used later) 
     384         srcv(jpr_itx1:jpr_itz2)%laction = .FALSE.    ! ice components not received 
    385385         srcv(jpr_itx1)%clgrid = 'U'                  ! ocean stress used after its transformation 
    386386         srcv(jpr_ity1)%clgrid = 'V'                  ! i.e. it is always at U- & V-points for i- & j-comp. resp. 
     
    502502      ! Allocate taum part of frcv which is used even when not received as coupling field 
    503503      IF ( .NOT. srcv(jpr_taum)%laction ) ALLOCATE( frcv(jpr_taum)%z3(jpi,jpj,srcv(jn)%nct) ) 
     504      ! Allocate itx1 and ity1 as they are used in sbc_cpl_ice_tau even if srcv(jpr_itx1)%laction = .FALSE. 
     505      IF( k_ice /= 0 ) THEN 
     506         IF ( .NOT. srcv(jpr_itx1)%laction ) ALLOCATE( frcv(jpr_itx1)%z3(jpi,jpj,srcv(jn)%nct) ) 
     507         IF ( .NOT. srcv(jpr_ity1)%laction ) ALLOCATE( frcv(jpr_ity1)%z3(jpi,jpj,srcv(jn)%nct) ) 
     508      END IF 
    504509 
    505510      ! ================================ ! 
Note: See TracChangeset for help on using the changeset viewer.