Changeset 5828 for branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90

Timestamp:

2015-10-23T19:03:26+02:00 (9 years ago)

Author:

mathiot

Message:

Fixes for initialisation issues + bug highlited on the MetO machine

File:

• branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90 (modified) (2 diffs)

branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90

@@ -129 +129 @@
          !
          ! initialisation
-         qisf(:,:)        = 0._wp  ; fwfisf(:,:) = 0._wp
-         risf_tsc(:,:,:)  = 0._wp
+         qisf(:,:)        = 0._wp  ; fwfisf  (:,:) = 0._wp
+         risf_tsc(:,:,:)  = 0._wp  ; fwfisf_b(:,:) = 0._wp
          !
          ! define isf tbl tickness, top and bottom indice
@@ -483 +483 @@
                   ! Presumably zdis can never be negative because gammas is very small compared to gammat
                   ! compute s freeze
-                  IF (zsfrz .GE. 0.0_wp) THEN ; zsfrz=(-zbqe-SQRT(zdis))*zaqer
-                  ELSE                        ; zsfrz=(-zbqe+SQRT(zdis))*zaqer
-                  ENDIF
+                  zsfrz=(-zbqe-SQRT(zdis))*zaqer
+                  IF ( zsfrz .LT. 0.0_wp ) zsfrz=(-zbqe+SQRT(zdis))*zaqer
 
                   ! compute t freeze