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

Annotation 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 - (hide annotations) (download)
Mon Oct 27 22:32:55 2003 UTC (20 years, 8 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 heimbach 1.4 C $Header: /u/gcmpack/MITgcm/pkg/grdchk/grdchk_get_obcs_mask.F,v 1.3 2003/07/18 21:10:16 heimbach Exp $
2 heimbach 1.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 heimbach 1.4 #if (defined (ALLOW_GRDCHK) && defined (ALLOW_OBCS_CONTROL))
35 heimbach 1.2 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 heimbach 1.3 if ( grdchkvarindex .EQ. 11 ) then
68 heimbach 1.2 #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 heimbach 1.3 call active_read_xz_loc( fname, tmpfldxz, iobcs,
74 heimbach 1.2 & .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 heimbach 1.3 else if ( grdchkvarindex .EQ. 12 ) then
91 heimbach 1.2 #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 heimbach 1.3 call active_read_xz_loc( fname, tmpfldxz, iobcs,
97 heimbach 1.2 & .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 heimbach 1.3 else if ( grdchkvarindex .EQ. 13 ) then
114 heimbach 1.2 #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 heimbach 1.3 call active_read_yz_loc( fname, tmpfldyz, iobcs,
120 heimbach 1.2 & .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 heimbach 1.3 else if ( grdchkvarindex .EQ. 14 ) then
137 heimbach 1.2 #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 heimbach 1.3 call active_read_yz_loc( fname, tmpfldyz, iobcs,
143 heimbach 1.2 & .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 heimbach 1.3 endif
159 heimbach 1.2
160     _END_MASTER( mythid )
161    
162     _BARRIER
163    
164 heimbach 1.4 #endif /* ALLOW_GRDCHK */
165 heimbach 1.2
166     end
167    
  ViewVC Help
Powered by ViewVC 1.1.22