/[MITgcm]/MITgcm/pkg/mom_common/mom_calc_hfacz.F
ViewVC logotype

Annotation of /MITgcm/pkg/mom_common/mom_calc_hfacz.F

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

Revision 1.1 - (hide annotations) (download)
Fri May 14 17:43:11 2004 UTC (20 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint53d_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint57s_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint55h_post, checkpoint57g_pre, checkpoint54b_post, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint57y_pre, checkpoint55, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint54f_post, checkpoint57v_post, checkpoint55g_post, checkpoint55f_post, checkpoint57r_post, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint53g_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint57q_post, checkpoint58g_post, checkpoint57z_post, checkpoint57c_post, checkpoint55e_post, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post 
Created pkg/mom_common/
 o changed pkg_groups and pkg_depend
 o moved duplicate files in mom_vecinv and mom_fluxform to mom_common/
This is in preparation for checking in Smagorinsky and Leith type variable
viscosities.
1 adcroft 1.1 C $Header: $
2     C $Name: $
3    
4     #include "MOM_COMMON_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: MOM_CALC_HFACZ
8    
9     C !INTERFACE: ==========================================================
10     SUBROUTINE MOM_CALC_HFACZ(
11     I bi,bj,k,
12     O hFacZ,r_hFacZ,
13     I myThid)
14    
15     C !DESCRIPTION:
16     C Calculates the fractional thickness at vorticity points
17    
18     C !USES: ===============================================================
19     IMPLICIT NONE
20     #include "SIZE.h"
21     #include "GRID.h"
22    
23     C !INPUT PARAMETERS: ===================================================
24     C bi,bj :: tile indices
25     C k :: vertical level
26     C myThid :: thread number
27     INTEGER bi,bj,k
28     INTEGER myThid
29    
30     C !OUTPUT PARAMETERS: ==================================================
31     C hFacZ :: fractional thickness at vorticity points
32     C r_hFacZ :: reciprocal
33     _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
34     _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
35    
36     C !LOCAL VARIABLES: ====================================================
37     C i,j :: loop indices
38     INTEGER I,J
39     _RL hFacZOpen
40     CEOP
41    
42     C-- Calculate open water fraction at vorticity points
43    
44     DO i=1-Olx,sNx+Olx
45     hFacZ(i,1-Oly)=0.
46     r_hFacZ(i,1-Oly)=0.
47     ENDDO
48    
49     DO j=2-Oly,sNy+Oly
50     hFacZ(1-Olx,j)=0.
51     r_hFacZ(1-Olx,j)=0.
52     DO i=2-Olx,sNx+Olx
53     hFacZOpen=min(_hFacW(i,j,k,bi,bj),
54     & _hFacW(i,j-1,k,bi,bj))
55     hFacZOpen=min(_hFacS(i,j,k,bi,bj),hFacZOpen)
56     hFacZOpen=min(_hFacS(i-1,j,k,bi,bj),hFacZOpen)
57     hFacZ(i,j)=hFacZOpen
58     IF (hFacZ(i,j).EQ.0.) THEN
59     r_hFacZ(i,j)=0.
60     ELSE
61     r_hFacZ(i,j)=1./hFacZ(i,j)
62     ENDIF
63     ENDDO
64     ENDDO
65    
66     RETURN
67     END
  ViewVC Help
Powered by ViewVC 1.1.22