pkg/admtlm/admtlm_init_fixed.F

C $Header: /u/gcmpack/MITgcm/pkg/admtlm/admtlm_init_fixed.F,v 1.3 2007/10/08 23:45:51 jmc Exp $
C $Name:  $

#include "ADMTLM_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	C $Header: /u/gcmpack/MITgcm/pkg/admtlm/admtlm_init_fixed.F,v 1.3 2007/10/08 23:45:51 jmc Exp $
2	C $Name: $
3
4	#include "ADMTLM_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	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	DO k=1,Nr
43	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	END DO
47	c
48	DO j=1,sNy
49	DO i=1,sNx
50	wetan(i,j,bi,bj) = 1. / (0.1 _d 0)**2
51	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
61	end