Changeset 10009 for NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/SAO/sao_read.F90

Timestamp:

2018-07-29T11:23:51+02:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-04): RK3 branch - step II.1 time-level dimension on ssh

File:

• NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/SAO/sao_read.F90 (modified) (5 diffs)

NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/SAO/sao_read.F90

@@ -9 +9 @@
    USE par_kind, ONLY: lc
    USE netcdf
-   USE oce,     ONLY: tsn, sshn
+   USE oce,     ONLY: tsn, ssh
    USE dom_oce, ONLY: nlci, nlcj, nimpp, njmpp, tmask
    USE par_oce, ONLY: jpi, jpj, jpk
@@ -77 +77 @@
 
       IF (TRIM(filename) == 'nofile') THEN
-         tsn (:,:,:,:) = fbrmdi
-         sshn(:,:)     = fbrmdi
+         tsn(:,:,:,:) = fbrmdi
+         ssh(:,:,Nnn) = fbrmdi
       ELSE
          WRITE(numout,*) "Opening :", TRIM(filename)
@@ -133 +133 @@
             WHERE(temp_sshn(:,:) == fill_val) temp_sshn(:,:) = fbrmdi
 
-            ! Initialise tsn, sshn to fbrmdi
+            ! Initialise tsn, ssh to fbrmdi
             tsn(:,:,:,:) = fbrmdi
-            sshn(:,:) = fbrmdi
+            ssh(:,:,Nnn) = fbrmdi
 
             ! Mask out missing data index
             tsn(1:nlci,1:nlcj,1:jpk,1) = temp_tn(:,:,:) * tmask(1:nlci,1:nlcj,1:jpk)
             tsn(1:nlci,1:nlcj,1:jpk,2) = temp_sn(:,:,:) * tmask(1:nlci,1:nlcj,1:jpk)
-            sshn(1:nlci,1:nlcj)        = temp_sshn(:,:) * tmask(1:nlci,1:nlcj,1)
+            ssh(1:nlci,1:nlcj,Nnn)     = temp_sshn(:,:) * tmask(1:nlci,1:nlcj,1)
 
-            ! Remove halo from tmask, tsn, sshn to prevent double obs counting
+            ! Remove halo from tmask, tsn, ssh to prevent double obs counting
             IF (jpi > nlci) THEN
                 tmask(nlci+1:,:,:) = 0
                 tsn(nlci+1:,:,:,1) = 0
                 tsn(nlci+1:,:,:,2) = 0
-                sshn(nlci+1:,:) = 0
+                ssh(nlci+1:,:,Nnn) = 0
             END IF
             IF (jpj > nlcj) THEN
@@ -153 +153 @@
                 tsn(:,nlcj+1:,:,1) = 0
                 tsn(:,nlcj+1:,:,2) = 0
-                sshn(:,nlcj+1:) = 0
+                ssh(:,nlcj+1:,Nnn) = 0
             END IF
 
@@ -161 +161 @@
             ! Mark all as missing data
             tsn(:,:,:,:) = fbrmdi
-            sshn(:,:) = fbrmdi
+            ssh(:,:,Nnn) = fbrmdi
          ENDIF
          ! Close netcdf file