Changeset 7837 for branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

Timestamp:

2017-03-27T10:50:24+02:00 (7 years ago)

Author:

mattmartin

Message:

First version of generic obs oper which works at NEMO3.6 with other GO6 branches.

File:

• branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90 (modified) (6 diffs)

branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

@@ -31 +31 @@
    USE sbcdcy,        ONLY : &    ! For calculation of where it is night-time
       & sbc_dcy, nday_qsr
+   USE obs_grid,      ONLY : & 
+      & obs_level_search     
 
    IMPLICIT NONE
@@ -53 +55 @@
 CONTAINS
 
+
    SUBROUTINE obs_prof_opt( prodatqc, kt, kpi, kpj, kpk,          &
       &                     kit000, kdaystp,                      &
-      &                     pvar1, pvar2, pgdept, pgdepw,
-      &                     pmask1, pmask2, &
+      &                     pvar1, pvar2, pgdept, pgdepw,         &
+      &                     pmask1, pmask2,                       &  
       &                     plam1, plam2, pphi1, pphi2,           &
       &                     k1dint, k2dint, kdailyavtypes )
@@ -154 +157 @@
       INTEGER ::   iend
       INTEGER ::   iobs
+      INTEGER ::   iin, ijn, ikn, ik   ! looping indices over interpolation nodes 
+      INTEGER ::   inum_obs
       INTEGER, DIMENSION(imaxavtypes) :: &
          & idailyavtypes
@@ -161 +166 @@
          & igrdj1, &
          & igrdj2
+      INTEGER, ALLOCATABLE, DIMENSION(:) :: iv_indic
+
       REAL(KIND=wp) :: zlam
       REAL(KIND=wp) :: zphi
@@ -188 +195 @@
          & zgphi2
       REAL(KIND=wp), DIMENSION(1) :: zmsk_1, zmsk_2   
-      REAL(KIND=wp), DIMENSION(:,:,:), ALLOCATABLE :: interp_corner       
+      REAL(KIND=wp), DIMENSION(:,:,:), ALLOCATABLE :: interp_corner
 
       LOGICAL :: ld_dailyav
@@ -627 +634 @@
           
          ENDIF 
+
+      ENDDO
 
       ! Deallocate the data for interpolation