pkg/admtlm/admtlm_init_fixed.F

C $Header:  $
C $Name:  $

#include "CPP_OPTIONS.h"

      subroutine admtlm_init_fixed( myThid )
C     /==========================================================\
C     | subroutine admtlm_init_fixed                             |
C     | o This routine assigns final T,S to cost function        |
C     \==========================================================/
       implicit none

C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#ifdef ALLOW_ADMTLM
# include "ctrl.h"
# include "ctrl_weights.h"
#endif

C     ======== Routine arguments ======================
C     myThid - Thread number for this instance of the routine.
      integer myThid

#ifdef ALLOW_ADMTLM
C     ========= Local variables =========================
      integer i, j, k
      integer bi, bj
      _RL dummy

      dummy = 0.

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
c
c N.B. !!!
c Values here should be 1/(std.dev.)^2
c to be consistent with ctrl_set_*pack_xy* routines
c N.B. !!!
c
        DO k=1,Nr
         wtheta(k,bi,bj) = 1. / (1. _d 0)**2
         wsalt (k,bi,bj) = 1. / (0.1 _d 0)**2
         wuvvel(k,bi,bj) = 1. / (0.01 _d 0)**2
        END DO
c
        DO j=1,sNy
         DO i=1,sNx
          wetan(i,j,bi,bj) = 1. / (0.1 _d 0)**2
         END DO
        END DO
c
       END DO
      END DO

      call active_write_xy('wetan',wetan,1,0,mythid,dummy)

#endif

      end
1	jmc	1.3	C $Header: $
2			C $Name: $
3	heimbach	1.1
4			#include "CPP_OPTIONS.h"
5
6			subroutine admtlm_init_fixed( myThid )
7			C /==========================================================\
8			C \| subroutine admtlm_init_fixed \|
9			C \| o This routine assigns final T,S to cost function \|
10			C \==========================================================/
11			implicit none
12
13			C == Global variables ===
14			#include "SIZE.h"
15			#include "EEPARAMS.h"
16			#include "PARAMS.h"
17			#ifdef ALLOW_ADMTLM
18			# include "ctrl.h"
19			# include "ctrl_weights.h"
20			#endif
21
22			C ======== Routine arguments ======================
23			C myThid - Thread number for this instance of the routine.
24			integer myThid
25
26			#ifdef ALLOW_ADMTLM
27			C ========= Local variables =========================
28			integer i, j, k
29			integer bi, bj
30			_RL dummy
31
32			dummy = 0.
33
34			DO bj=myByLo(myThid),myByHi(myThid)
35			DO bi=myBxLo(myThid),myBxHi(myThid)
36			c
37	heimbach	1.2	c N.B. !!!
38			c Values here should be 1/(std.dev.)^2
39			c to be consistent with ctrl_set_pack_xy routines
40			c N.B. !!!
41			c
42	heimbach	1.1	DO k=1,Nr
43	heimbach	1.2	wtheta(k,bi,bj) = 1. / (1. _d 0)**2
44			wsalt (k,bi,bj) = 1. / (0.1 _d 0)**2
45			wuvvel(k,bi,bj) = 1. / (0.01 _d 0)**2
46	heimbach	1.1	END DO
47			c
48			DO j=1,sNy
49			DO i=1,sNx
50	heimbach	1.2	wetan(i,j,bi,bj) = 1. / (0.1 _d 0)**2
51	heimbach	1.1	END DO
52			END DO
53			c
54			END DO
55			END DO
56
57			call active_write_xy('wetan',wetan,1,0,mythid,dummy)
58
59			#endif
60	jmc	1.3
61	heimbach	1.1	end