/[MITgcm]/MITgcm/verification/aim.5l_cs/code/mom_vi_hfacz_diss.F
ViewVC logotype

Contents of /MITgcm/verification/aim.5l_cs/code/mom_vi_hfacz_diss.F

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


Revision 1.2 - (show annotations) (download)
Sun Jul 13 19:26:05 2003 UTC (16 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint57m_post, checkpoint52l_pre, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint56b_post, checkpoint57y_post, checkpoint52j_pre, checkpoint51o_pre, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint54d_post, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint54e_post, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint51l_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint52k_post, checkpoint59, checkpoint58, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint58f_post, checkpoint52f_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint51t_post, checkpoint58t_post, checkpoint51n_post, checkpoint55i_post, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, checkpoint51s_post, checkpoint57t_post, checkpoint55c_post, checkpoint51j_post, checkpoint52e_pre, checkpoint57v_post, checkpoint57f_post, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint64, checkpoint65, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint57a_post, checkpoint57h_pre, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint52b_pre, checkpoint54b_post, checkpoint58w_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint57y_pre, checkpoint55g_post, checkpoint51q_post, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint52f_pre, checkpoint55d_post, checkpoint58e_post, checkpoint54a_pre, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint53c_post, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, branchpoint-genmake2, checkpoint54a_post, checkpoint55h_post, checkpoint58n_post, checkpoint51r_post, checkpoint51i_post, checkpoint57e_post, checkpoint55b_post, checkpoint53a_post, checkpoint65o, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint52d_post, checkpoint53g_post, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint52a_pre, checkpoint62b, checkpoint58v_post, checkpoint52i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint56a_post, checkpoint64y, checkpoint64x, checkpoint58l_post, checkpoint64z, checkpoint53f_post, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint57h_done, checkpoint52j_post, checkpoint57j_post, checkpoint57f_pre, checkpoint61f, checkpoint58g_post, branch-netcdf, checkpoint52l_post, checkpoint58x_post, checkpoint61n, checkpoint52n_post, checkpoint53b_pre, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint51e_post, checkpoint57a_pre, checkpoint55a_post, checkpoint57o_post, checkpoint51o_post, checkpoint61q, checkpoint57k_post, checkpoint51f_pre, checkpoint61z, checkpoint53b_post, checkpoint52a_post, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint51g_post, ecco_c52_e35, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint51m_post, checkpoint53d_pre, checkpoint58s_post, checkpoint55e_post, checkpoint61g, checkpoint61d, checkpoint54c_post, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint51p_post, checkpoint61x, checkpoint61y, checkpoint51u_post, HEAD
Branch point for: branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.1: +24 -1 lines
return to the standard def. of hFacZ at the corner of the C-S. grid.
(safer if topography step at the corner, but not the case here in AIM-cs)

1 C $Header: /u/gcmpack/MITgcm/verification/aim.5l_cs/code/mom_vi_hfacz_diss.F,v 1.1 2002/01/09 00:28:56 jmc Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 SUBROUTINE MOM_VI_HFACZ_DISS(
7 I bi,bj,k,
8 O hFacZ,r_hFacZ,
9 I myThid)
10 IMPLICIT NONE
11 C
12 C Compute hFactor (and reciprol) at the corner (Z-point)
13 C used for vorticity,divergence form of viscous term
14 C (also in Shap_S2 formulation)
15 C
16
17 C == Global variables ==
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "GRID.h"
22
23 C == Routine arguments ==
24 INTEGER bi,bj,k
25 _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26 _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27 _RS hFacZOpen
28 INTEGER myThid
29
30 C == Local variables ==
31 INTEGER I,J
32
33 C-- Calculate open water fraction at vorticity points
34
35 DO i=1-Olx,sNx+Olx
36 hFacZ(i,1-Oly)=0.
37 r_hFacZ(i,1-Oly)=0.
38 ENDDO
39
40 DO j=2-Oly,sNy+Oly
41 hFacZ(1-Olx,j)=0.
42 r_hFacZ(1-Olx,j)=0.
43 DO i=2-Olx,sNx+Olx
44 hFacZOpen=MIN(_hFacW(i,j,k,bi,bj)
45 & + _hFacW(i,j-1,k,bi,bj)
46 & , _hFacS(i,j,k,bi,bj)
47 & + _hFacS(i-1,j,k,bi,bj)
48 & )
49 hFacZ(i,j)=0.5*hFacZOpen
50 IF (hFacZ(i,j).EQ.0.) THEN
51 r_hFacZ(i,j)=0.
52 ELSE
53 r_hFacZ(i,j)=1./hFacZ(i,j)
54 ENDIF
55 ENDDO
56 ENDDO
57
58 C- Special stuff for Cubed Sphere
59 C above formula is ambiguous when only 3 edges instead of 4,
60 C => return to min-3 definition at the Cubed-Sphere corners
61 IF (useCubedSphereExchange) THEN
62 DO j=1,sNy+1,sNy
63 DO i=1,sNx+1,sNx
64 hFacZOpen=MIN( _hFacW(i,j,k,bi,bj)
65 & , _hFacW(i,j-1,k,bi,bj)
66 & , _hFacS(i,j,k,bi,bj)
67 & , _hFacS(i-1,j,k,bi,bj)
68 & )
69 hFacZ(i,j)=hFacZOpen
70 IF (hFacZOpen.EQ.0.) THEN
71 r_hFacZ(i,j)=0.
72 ELSE
73 r_hFacZ(i,j)=1./hFacZOpen
74 ENDIF
75 ENDDO
76 ENDDO
77 ENDIF
78
79 RETURN
80 END

  ViewVC Help
Powered by ViewVC 1.1.22