pkg/admtlm/admtlm_init_fixed.F


#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
        DO k=1,Nr
         wtheta(k,bi,bj) = 1. _d 0
         wsalt (k,bi,bj) = 0.1 _d 0
         wuvvel(k,bi,bj) = 0.01 _d 0
        END DO
c
        DO j=1,sNy
         DO i=1,sNx
          wetan(i,j,bi,bj) = 0.1 _d 0
         END DO
        END DO
c
       END DO
      END DO

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

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