Changeset 6394

Timestamp:

2016-03-17T18:05:01+01:00 (8 years ago)

Author:

jcastill

Message:

Bug fix when coupling sea surface height with ln_apr_dyn - at the first time step the pressure can not be used in the equations as it is not defined yet

File:

• branches/UKMO/dev_r5518_ww3_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_ww3_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -2200 +2200 @@
          !                          ! removed inverse barometer ssh when Patm
          !                          forcing is used (for sea-ice dynamics)
-         IF( ln_apr_dyn ) THEN   ;   ztmp1(:,:) = sshb(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) )
-         ELSE                    ;   ztmp1(:,:) = sshn(:,:)
+         IF( ln_apr_dyn ) THEN
+             IF( kt /= nit000 ) THEN
+                 ztmp1(:,:) = sshb(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) )
+             ELSE
+                 ztmp1(:,:) = sshb(:,:)
+             ENDIF
+         ELSE
+             ztmp1(:,:) = sshn(:,:)
          ENDIF
          CALL cpl_snd( jps_ssh   , isec, RESHAPE ( ztmp1            , (/jpi,jpj,1/) ), info )
-
       ENDIF
       !                                                        ! SSS