/[MITgcm]/MITgcm/pkg/admtlm/admtlm_dsvd2model.F
ViewVC logotype

Annotation of /MITgcm/pkg/admtlm/admtlm_dsvd2model.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (hide annotations) (download)
Fri Apr 29 10:36:45 2005 UTC (19 years, 1 month ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57o_post, checkpoint57v_post, checkpoint57m_post, checkpoint57s_post, checkpoint57k_post, checkpoint57i_post, checkpoint57r_post, checkpoint57h_done, checkpoint57n_post, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, checkpoint57j_post, checkpoint57h_pre, checkpoint57l_post, checkpoint57h_post
o Extend svd to full state.
o Pass control variables as arguments.

1 heimbach 1.1
2     #include "CPP_OPTIONS.h"
3    
4     subroutine admtlm_dsvd2model( tmpstate, mythid )
5     C /==========================================================\
6     C | subroutine admtlm_map |
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 "adcost.h"
17     # include "g_cost.h"
18     #endif
19    
20     C ======== Routine arguments ======================
21     C myThid - Thread number for this instance of the routine.
22     integer myThid
23     _RL tmpstate (snx,sny,nsx,nsy,4*Nr+1)
24    
25     #ifdef ALLOW_ADMTLM
26     C ========= Local variables =========================
27     integer i, j, k
28     integer bi, bj
29    
30     DO bj=myByLo(myThid),myByHi(myThid)
31     DO bi=myBxLo(myThid),myBxHi(myThid)
32     DO j=1,sNy
33     DO i=1,sNx
34     DO k=1,4*Nr
35     g_objf_state_final(i,j,bi,bj,k) =
36     & tmpstate(i,j,bi,bj,k)
37     END DO
38     g_objf_state_final(i,j,bi,bj,4*Nr+1) =
39     & tmpstate(i,j,bi,bj,4*Nr+1)
40     END DO
41     END DO
42     END DO
43     END DO
44     #endif
45    
46     end

  ViewVC Help
Powered by ViewVC 1.1.22