Changeset 14982 for NEMO/trunk/src/OCE/OBS/diaobs.F90

Timestamp:

2021-06-11T16:52:03+02:00 (3 years ago)

Author:

hadcv

Message:

#2665: Various fixes for code enabled with key_qco/key_linssh

File:

• NEMO/trunk/src/OCE/OBS/diaobs.F90 (modified) (6 diffs)

NEMO/trunk/src/OCE/OBS/diaobs.F90

@@ -99 +99 @@
    CHARACTER(len=8), PUBLIC, DIMENSION(:), ALLOCATABLE ::   cobstypesprof, cobstypessurf   !: Profile & surface obs types
 
+#  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
    !! NEMO/OCE 4.0 , NEMO Consortium (2018)
@@ -623 +624 @@
       INTEGER :: jtype             ! Data loop variable
       INTEGER :: jvar              ! Variable number
-      INTEGER :: ji, jj            ! Loop counters
+      INTEGER :: ji, jj, jk        ! Loop counters
       REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: &
          & zprofvar                ! Model values for variables in a prof ob
@@ -634 +635 @@
          & zglam,    &             ! Model longitudes for prof variables
          & zgphi                   ! Model latitudes for prof variables
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zdept, zdepw
 
       !-----------------------------------------------------------------------
@@ -650 +652 @@
 
       IF ( nproftypes > 0 ) THEN
+
+         ALLOCATE( zdept(jpi,jpj,jpk), zdepw(jpi,jpj,jpk) )
+         DO jk = 1, jpk
+            zdept(:,:,jk) = gdept(:,:,jk,Kmm)
+            zdepw(:,:,jk) = gdepw(:,:,jk,Kmm)
+         END DO
 
          DO jtype = 1, nproftypes
@@ -687 +695 @@
                   &               nit000, idaystp, jvar,                   &
                   &               zprofvar(:,:,:,jvar),                    &
-                  &               gdept(:,:,:,Kmm), gdepw(:,:,:,Kmm),      &
+                  &               zdept(:,:,:), zdepw(:,:,:),      &
                   &               zprofmask(:,:,:,jvar),                   &
                   &               zglam(:,:,jvar), zgphi(:,:,jvar),        &
@@ -697 +705 @@
 
          END DO
+
+         DEALLOCATE( zdept, zdepw )
 
       ENDIF