/[MITgcm]/MITgcm/pkg/ctrl/ctrl_set_globfld_xyz.F
ViewVC logotype

Contents of /MITgcm/pkg/ctrl/ctrl_set_globfld_xyz.F

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

Revision 1.8 - (show annotations) (download)
Tue Nov 16 05:42:12 2004 UTC (19 years, 5 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint56b_post, checkpoint57y_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint59, checkpoint58, checkpoint57, checkpoint56, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint57l_post, checkpoint57t_post, checkpoint57v_post, checkpoint57f_post, checkpoint57a_post, checkpoint57h_pre, checkpoint58w_post, checkpoint57h_post, checkpoint57y_pre, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, checkpoint57c_pre, checkpoint58r_post, checkpoint58n_post, checkpoint57e_post, checkpoint59a, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint58v_post, checkpoint56a_post, checkpoint58l_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint58g_post, checkpoint58x_post, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint57a_pre, checkpoint57o_post, checkpoint57k_post, checkpoint57w_post, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post
Changes since 1.7: +1 -12 lines 
More on dsvd vs. MITgcm interfacing
o handling of g_, ad, via admtlm_vector (mds...vector)
o use ctrl_pack/unpack for admtlm_vector I/O
o use optimcycle for dsvd iteration
o make sure norm is w.r.t. derived quantities
1 C
2 C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_set_globfld_xyz.F,v 1.7 2003/11/06 22:05:08 heimbach Exp $
3 C $Name: $
4
5 #include "CTRL_CPPOPTIONS.h"
6
7 subroutine ctrl_set_globfld_xyz(
8 I fname, ivartype, mythid )
9
10 c ==================================================================
11 c SUBROUTINE ctrl_set_globfld_xyz
12 c ==================================================================
13 c
14 c o initialise field
15 c
16 c started: heimbach@mit.edu, 16-Aug-2001
17 c
18 c ==================================================================
19
20 implicit none
21
22 c == global variables ==
23
24 #include "EEPARAMS.h"
25 #include "SIZE.h"
26 #include "PARAMS.h"
27 #include "GRID.h"
28
29 #include "ctrl.h"
30 #include "optim.h"
31
32 c == routine arguments ==
33
34 character*( 80) fname
35 integer ivartype
36 integer mythid
37
38 c == local variables ==
39
40 integer bi,bj
41 integer ip,jp
42 integer i,j,k
43 integer itlo,ithi
44 integer jtlo,jthi
45 integer jmin,jmax
46 integer imin,imax
47 integer irec
48
49 _RL globfld3d( snx,nsx,npx,sny,nsy,npy,nr )
50
51 c == external ==
52
53 c == end of interface ==
54
55 jtlo = 1
56 jthi = nsy
57 itlo = 1
58 ithi = nsx
59 jmin = 1
60 jmax = sny
61 imin = 1
62 imax = snx
63
64 c Initialise temporary file
65 do k = 1,nr
66 do jp = 1,nPy
67 do bj = jtlo,jthi
68 do j = jmin,jmax
69 do ip = 1,nPx
70 do bi = itlo,ithi
71 do i = imin,imax
72 globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0
73 enddo
74 enddo
75 enddo
76 enddo
77 enddo
78 enddo
79 enddo
80
81 c-- Only the master thread will do I/O.
82 _BEGIN_MASTER( mythid )
83
84 do irec = 1, ncvarrecs(ivartype)
85 call MDSWRITEFIELD_3D_GL( fname, ctrlprec, 'RL',
86 & Nr, globfld3d,
87 & irec, optimcycle, mythid)
88 enddo
89
90 _END_MASTER( mythid )
91
92 end
  ViewVC Help
Powered by ViewVC 1.1.22