pkg/admtlm/admtlm_init_fixed.F

b35bd3101a Jean*0001 #include "ADMTLM_OPTIONS.h"
949e225614 Patr*0002
                0003       subroutine admtlm_init_fixed( myThid )
                0004 C     /==========================================================\
                0005 C     | subroutine admtlm_init_fixed                             |
                0006 C     | o This routine assigns final T,S to cost function        |
                0007 C     \==========================================================/
                0008        implicit none
                0009
                0010 C     == Global variables ===
                0011 #include "SIZE.h"
                0012 #include "EEPARAMS.h"
                0013 #include "PARAMS.h"
                0014 #ifdef ALLOW_ADMTLM
4d72283393 Mart*0015 # include "CTRL.h"
949e225614 Patr*0016 # include "ctrl_weights.h"
                0017 #endif
                0018
                0019 C     ======== Routine arguments ======================
                0020 C     myThid - Thread number for this instance of the routine.
                0021       integer myThid
                0022
                0023 #ifdef ALLOW_ADMTLM
                0024 C     ========= Local variables =========================
                0025       integer i, j, k
                0026       integer bi, bj
                0027       _RL dummy
                0028
                0029       dummy = 0.
                0030
                0031       DO bj=myByLo(myThid),myByHi(myThid)
                0032        DO bi=myBxLo(myThid),myBxHi(myThid)
                0033 c
8209108cec Patr*0034 c N.B. !!!
                0035 c Values here should be 1/(std.dev.)^2
                0036 c to be consistent with ctrl_set_*pack_xy* routines
                0037 c N.B. !!!
                0038 c
949e225614 Patr*0039         DO k=1,Nr
8209108cec Patr*0040          wtheta(k,bi,bj) = 1. / (1. _d 0)**2
                0041          wsalt (k,bi,bj) = 1. / (0.1 _d 0)**2
                0042          wuvvel(k,bi,bj) = 1. / (0.01 _d 0)**2
949e225614 Patr*0043         END DO
                0044 c
                0045         DO j=1,sNy
                0046          DO i=1,sNx
8209108cec Patr*0047           wetan(i,j,bi,bj) = 1. / (0.1 _d 0)**2
949e225614 Patr*0048          END DO
                0049         END DO
                0050 c
                0051        END DO
                0052       END DO
                0053
                0054       call active_write_xy('wetan',wetan,1,0,mythid,dummy)
                0055
                0056 #endif
305d315bb8 Jean*0057
949e225614 Patr*0058       end