Changeset 1244 for trunk/NEMO/OPA_SRC/SBC/sbcssm.F90

Timestamp:

2009-01-06T12:13:44+01:00 (15 years ago)

Author:

ctlod

Message:

improve the LIM 3.0 stability in using sea surface height, see ticket: #286

File:

• trunk/NEMO/OPA_SRC/SBC/sbcssm.F90 (modified) (8 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcssm.F90

@@ -64 +64 @@
          sst_m(:,:) = tn(:,:,1)
          sss_m(:,:) = sn(:,:,1)
+         ssh_m(:,:) = sshn(:,:)
          !
       ELSE
@@ -78 +79 @@
                CALL iom_get( numror, jpdom_autoglo, 'sst_m'  , sst_m  )   !   "         "    temperature (T-point)
                CALL iom_get( numror, jpdom_autoglo, 'sss_m'  , sss_m  )   !   "         "    salinity    (T-point)
+               CALL iom_get( numror, jpdom_autoglo, 'ssh_m'  , ssh_m  )   !   "         "    height      (T-point)
                !
                IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN      ! nn_fsbc has changed between 2 runs
@@ -87 +89 @@
                   sst_m(:,:) = zcoef * sst_m(:,:)
                   sss_m(:,:) = zcoef * sss_m(:,:)
+                  ssh_m(:,:) = zcoef * ssh_m(:,:)
                ELSE
                   IF(lwp) WRITE(numout,*) '~~~~~~~   mean fields read in the ocean restart file'
@@ -97 +100 @@
                sst_m(:,:) = zcoef * tn(:,:,1)
                sss_m(:,:) = zcoef * sn(:,:,1)
+               ssh_m(:,:) = zcoef * sshn(:,:)
             ENDIF
             !                                             ! ---------------------------------------- !
@@ -105 +109 @@
             sst_m(:,:) = 0.e0
             sss_m(:,:) = 0.e0
+            ssh_m(:,:) = 0.e0
          ENDIF
          !                                                ! ---------------------------------------- !
@@ -113 +118 @@
          sst_m(:,:) = sst_m(:,:) + tn(:,:,1)
          sss_m(:,:) = sss_m(:,:) + sn(:,:,1)
+         ssh_m(:,:) = ssh_m(:,:) + sshn(:,:)
 
          !                                                ! ---------------------------------------- !
@@ -122 +128 @@
             ssu_m(:,:) = ssu_m(:,:) * zcoef           ! mean suface current  [m/s]
             ssv_m(:,:) = ssv_m(:,:) * zcoef           !
+            ssh_m(:,:) = ssh_m(:,:) * zcoef           ! mean SSH             [m]
             !
          ENDIF
@@ -137 +144 @@
             CALL iom_rstput( kt, nitrst, numrow, 'sst_m'  , sst_m  )
             CALL iom_rstput( kt, nitrst, numrow, 'sss_m'  , sss_m  )
+            CALL iom_rstput( kt, nitrst, numrow, 'ssh_m'  , ssh_m  )
             !
          ENDIF