3dVarBio/obs_vec.f90 at master · inogs/3dVarBio · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
subroutine obs_vec


  !---------------------------------------------------------------------------
  !                                                                          !
  !    Copyright 2018 Anna Teruzzi, OGS, Trieste                         !
  !                                                                          !
  !    This file is part of 3DVarBio.
  !    3DVarBio is based on OceanVar (Dobricic, 2006)                                          !
  !                                                                          !
  !    3DVarBio is  free software: you can redistribute it and/or modify.     !
  !    it under the terms of the GNU General Public License as published by  !
  !    the Free Software Foundation, either version 3 of the License, or     !
  !    (at your option) any later version.                                   !
  !                                                                          !
  !    3DVarBio is  distributed in the hope that it will be useful,           !
  !    but WITHOUT ANY WARRANTY; without even the implied warranty of        !
  !    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         !
  !    GNU General Public License for more details.                          !
  !                                                                          !
  !    You should have received a copy of the GNU General Public License     !
  !    along with OceanVar.  If not, see <http://www.gnu.org/licenses/>.       !
  !                                                                          !
  !---------------------------------------------------------------------------

  !-----------------------------------------------------------------------
  !                                                                      !
  ! Create the observational vector                                      !
  !                                                                      !
  ! Version 1: A. Teruzzi 2018                                           !
  !-----------------------------------------------------------------------


  use set_knd
  use drv_str
  use obs_str
  use mpi_str

  implicit none

  INTEGER(i4)    ::  k, i

  ! -------
  ! Define observational vector

  obs%no = sat%nc + arg%no

  if(MyId .eq. 0) &
       write(drv%dia,*) ' Total number of observations: ', obs%no

  ALLOCATE ( obs%inc(obs%no)) ; obs%inc = huge(obs%inc(1))
  ALLOCATE ( obs%amo(obs%no)) ; obs%amo = huge(obs%amo(1))
  ALLOCATE ( obs%res(obs%no)) ; obs%res = huge(obs%res(1))
  ALLOCATE ( obs%err(obs%no)) ; obs%err = huge(obs%err(1))
  ALLOCATE ( obs%gra(obs%no)) ; obs%gra = huge(obs%gra(1))


 k=0

 if (drv%argo_obs .eq. 1) then
    ! ARGO observations
    do i=1,arg%no
       if(arg%flc(i).eq.1)then
          k=k+1
          obs%res(k) = arg%res(i)
          obs%err(k) = arg%err(i)
       endif
    enddo

    DEALLOCATE(arg%res, arg%err)

 endif

 ! Observations of chlorophyll
 if(drv%sat_obs .eq. 1) then
  do i=1,sat%no
    if(sat%flc(i).eq.1)then
       k=k+1
       obs%res(k) = sat%res(i)
       obs%err(k) = sat%err(i)
    endif
  enddo

  DEALLOCATE(sat%res, sat%err)

 endif


end subroutine obs_vec