/[MITgcm]/MITgcm/pkg/grdchk/grdchk_get_obcs_mask.F
ViewVC logotype

Contents of /MITgcm/pkg/grdchk/grdchk_get_obcs_mask.F

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

Revision 1.4 - (show annotations) (download)
Mon Oct 27 22:32:55 2003 UTC (20 years, 7 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint58e_post, checkpoint57v_post, checkpoint52n_post, checkpoint52j_post, checkpoint53d_post, checkpoint58u_post, checkpoint58w_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint52e_post, checkpoint57s_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint51q_post, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, hrcube_1, checkpoint57e_post, branch-netcdf, checkpoint52d_pre, checkpoint52l_post, checkpoint55h_post, checkpoint58n_post, checkpoint51r_post, checkpoint58x_post, checkpoint52k_post, checkpoint52b_pre, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint58t_post, checkpoint58h_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint58q_post, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint51o_post, checkpoint51p_post, checkpoint58j_post, checkpoint52a_pre, checkpoint59a, checkpoint55f_post, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint52d_post, eckpoint57e_pre, checkpoint52a_post, checkpoint57h_done, checkpoint58f_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint57x_post, checkpoint57n_post, checkpoint52c_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, ecco_c52_e35, hrcube5, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint55e_post, checkpoint58k_post, checkpoint52i_post, checkpoint52j_pre, checkpoint58v_post, checkpoint53f_post, checkpoint55a_post, checkpoint51t_post, checkpoint53d_pre, checkpoint54c_post, checkpoint58s_post, checkpoint58p_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint57h_post, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: branch-nonh, netcdf-sm0
Changes since 1.3: +3 -3 lines 
o cleaning ALLOW_GRADIENT_CHECK -> ALLOW_GRDCHK
o cleaning some ALLOW_TANGENTLINEAR_RUN -> ALLOW_AUTODIFF
o bug fix in find_alpha.F for MDJWF:
  - modif. to alpha = 1/D*( dN/dT - rho*dD/Dt) to account for
    change rho -> rho-rhoConst
  - replace call find_rho to find_rhonum
1 C $Header: /u/gcmpack/MITgcm/pkg/grdchk/grdchk_get_obcs_mask.F,v 1.3 2003/07/18 21:10:16 heimbach Exp $
2
3 #include "CTRL_CPPOPTIONS.h"
4
5
6 subroutine grdchk_get_obcs_mask( mythid )
7
8 c ==================================================================
9 c SUBROUTINE grdchk_get_obcs_mask
10 c ==================================================================
11 c
12 c o Get obcs masks from file
13 c
14 c started: heimbach@mit.edu: 22-Apr-2003
15 c
16 c ==================================================================
17 c SUBROUTINE grdchk_get_obcs_mask
18 c ==================================================================
19
20 implicit none
21
22 c == global variables ==
23
24 #include "EEPARAMS.h"
25 #include "SIZE.h"
26 #include "GRID.h"
27 #include "ctrl.h"
28 #include "grdchk.h"
29
30 c == routine arguments ==
31
32 integer mythid
33
34 #if (defined (ALLOW_GRDCHK) && defined (ALLOW_OBCS_CONTROL))
35 c == local variables ==
36
37 integer bi,bj
38 integer i,j,k
39 integer irec,iobcs
40 integer itlo,ithi
41 integer jtlo,jthi
42 integer jmin,jmax
43 integer imin,imax
44
45 _RL dummy
46
47 character*( 80) fname
48
49 c == external ==
50
51 integer ilnblnk
52 external ilnblnk
53
54 c == end of interface ==
55
56 jtlo = 1
57 jthi = nsy
58 itlo = 1
59 ithi = nsx
60 jmin = 1
61 jmax = sny
62 imin = 1
63 imax = snx
64
65 _BEGIN_MASTER( mythid )
66
67 if ( grdchkvarindex .EQ. 11 ) then
68 #ifdef ALLOW_OBCSN_CONTROL
69 write(fname(1:80),'(80a)') ' '
70 write(fname(1:80),'(a)') 'maskobcsn'
71 c
72 do iobcs = 1,nobcs
73 call active_read_xz_loc( fname, tmpfldxz, iobcs,
74 & .false., .false., 0, mythid, dummy)
75 c
76 do bj = jtlo,jthi
77 do bi = itlo,ithi
78 do k = 1,nr
79 do i = imin,imax
80 grdchk_maskxz(i,k,bi,bj,iobcs) =
81 & tmpfldxz(i,k,bi,bj)
82 enddo
83 enddo
84 enddo
85 enddo
86 c
87 enddo
88 #endif
89
90 else if ( grdchkvarindex .EQ. 12 ) then
91 #ifdef ALLOW_OBCSS_CONTROL
92 write(fname(1:80),'(80a)') ' '
93 write(fname(1:80),'(a)') 'maskobcss'
94 c
95 do iobcs = 1,nobcs
96 call active_read_xz_loc( fname, tmpfldxz, iobcs,
97 & .false., .false., 0, mythid, dummy)
98 c
99 do bj = jtlo,jthi
100 do bi = itlo,ithi
101 do k = 1,nr
102 do i = imin,imax
103 grdchk_maskxz(i,k,bi,bj,iobcs) =
104 & tmpfldxz(i,k,bi,bj)
105 enddo
106 enddo
107 enddo
108 enddo
109 c
110 enddo
111 #endif
112
113 else if ( grdchkvarindex .EQ. 13 ) then
114 #ifdef ALLOW_OBCSW_CONTROL
115 write(fname(1:80),'(80a)') ' '
116 write(fname(1:80),'(a)') 'maskobcsw'
117 c
118 do iobcs = 1,nobcs
119 call active_read_yz_loc( fname, tmpfldyz, iobcs,
120 & .false., .false., 0, mythid, dummy)
121 c
122 do bj = jtlo,jthi
123 do bi = itlo,ithi
124 do k = 1,nr
125 do j = jmin,jmax
126 grdchk_maskyz(j,k,bi,bj,iobcs) =
127 & tmpfldyz(j,k,bi,bj)
128 enddo
129 enddo
130 enddo
131 enddo
132 c
133 enddo
134 #endif
135
136 else if ( grdchkvarindex .EQ. 14 ) then
137 #ifdef ALLOW_OBCSE_CONTROL
138 write(fname(1:80),'(80a)') ' '
139 write(fname(1:80),'(a)') 'maskobcse'
140 c
141 do iobcs = 1,nobcs
142 call active_read_yz_loc( fname, tmpfldyz, iobcs,
143 & .false., .false., 0, mythid, dummy)
144 c
145 do bj = jtlo,jthi
146 do bi = itlo,ithi
147 do k = 1,nr
148 do j = jmin,jmax
149 grdchk_maskyz(j,k,bi,bj,iobcs) =
150 & tmpfldyz(j,k,bi,bj)
151 enddo
152 enddo
153 enddo
154 enddo
155 c
156 enddo
157 #endif
158 endif
159
160 _END_MASTER( mythid )
161
162 _BARRIER
163
164 #endif /* ALLOW_GRDCHK */
165
166 end
167
  ViewVC Help
Powered by ViewVC 1.1.22