Changeset 11951 for NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/ice.F90

Timestamp:

2019-11-22T16:52:49+01:00 (4 years ago)

Author:

stefryn

Message:

fixed bug in stress calculation

File:

• NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/ice.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/ice.F90

@@ -226 +226 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   shear_i     !: Shear of the velocity field                  [s-1]
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   aniso_11, aniso_12   !: structure tensor elements
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   rdg_conv
    !
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   t_bo        !: Sea-Ice bottom temperature [Kelvin]     
@@ -409 +410 @@
          &      stress1_i(jpi,jpj) , stress2_i(jpi,jpj) , stress12_i(jpi,jpj) ,                      &
          &      delta_i  (jpi,jpj) , divu_i   (jpi,jpj) , shear_i   (jpi,jpj) ,                      &
-         &      aniso_11 (jpi,jpj) , aniso_12 (jpi,jpj) , STAT=ierr(ii) )
+         &      aniso_11 (jpi,jpj) , aniso_12 (jpi,jpj) , rdg_conv  (jpi,jpj) , STAT=ierr(ii) )
 
       ii = ii + 1