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

Contents 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 - (show annotations) (download)
Fri May 14 17:43:11 2004 UTC (21 years, 9 months 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 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