New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
lib_cray.f90 in branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC – NEMO

source: branches/UKMO/dev_r5785_SSS_obsoper/NEMOGCM/NEMO/OPA_SRC/lib_cray.f90 @ 7773

Last change on this file since 7773 was 7773, checked in by mattmartin, 7 years ago

Committing updates after doing the following:

  • merging the branch dev_r4650_general_vert_coord_obsoper@7763 into this branch
  • updating it so that the following OBS changes were implemented correctly on top of the simplification changes:
    • generalised vertical coordinate for profile obs. This was done so that is now the default option.
    • sst bias correction implemented with the new simplified obs code.
    • included the biogeochemical obs types int he new simplified obs code.
    • included the changes to exclude obs in the boundary for limited area models
    • included other changes for the efficiency of the obs operator to remove global arrays.
File size: 1.1 KB
Line 
1!  Cray subroutines or functions used by OPA model and possibly
2!  not found on other platforms.
3!
4!  check their existence
5! 
6!  wheneq
7!!----------------------------------------------------------------------
8!!  OPA 9.0 , LOCEAN-IPSL (2005)
9!! $Id: lib_cray.f90 7740 2017-02-27 13:18:43Z mattmartin $
10!! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
11!!----------------------------------------------------------------------
12SUBROUTINE lib_cray
13      WRITE(*,*) 'lib_cray: You should not have seen this print! error?'
14END SUBROUTINE lib_cray
15
16SUBROUTINE wheneq ( i, x, j, t, ind, nn )
17        IMPLICIT NONE
18
19        INTEGER , INTENT (  in ) :: i, j
20        INTEGER , INTENT ( out ) :: nn
21        REAL    , INTENT (  in ), DIMENSION (1+(i-1)*j) :: x
22        REAL    , INTENT (  in ) :: t
23        INTEGER , INTENT ( out ), DIMENSION (1+(i-1)*j) :: ind
24        INTEGER :: n, k
25        nn = 0
26        DO n = 1, i
27          k = 1 + (n-1) * j
28          IF ( x ( k) == t ) THEN
29              nn = nn + 1
30              ind (nn) = k
31          ENDIF
32        END DO
33
34END SUBROUTINE wheneq
Note: See TracBrowser for help on using the repository browser.