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

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

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


Revision 1.1.4.1 - (show annotations) (download)
Tue Feb 26 16:05:07 2002 UTC (22 years, 1 month ago) by adcroft
Branch: release1
CVS Tags: release1_p12, release1_p13, release1_p10, release1_p11, release1_p16, release1_p17, release1_p15, release1_p14, release1_p8, release1_p9, release1_p2, release1_p3, release1_p4, release1_p6, release1_p1, release1_p5, release1_p7, release1_p13_pre, release1_p12_pre, release1_chkpt44d_post
Branch point for: release1_50yr
Changes since 1.1: +2 -2 lines
Merging changes on MAIN between checkpoint43 and checkpoint43a-release1mods
Command: cvs -q update -jcheckpoint43 -jcheckpoint43a-release1mods -d -P

These changes are most of the changes between c43 and c44 except those
that occured after "12:45 11 Jan 2002". As far as I can tell it is
checkpoint43 with the following mods:

  o fix bug in mom_vi_del2uv
  o select when filters are applied ; add options to zonal_filter (data.zonfilt)  o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection)
  o update AIM experiments (NCEP input files)
  o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf)
  o added some stuff for AD
  o Jamar wet-points

This update does not contain the following mods that are in checkpoint44

  o bug fix in pkg/generic_advdiff/
    - thread related bug, bi,bj arguments in vertical advection routines
  o some changes to pkg/autodiff, pkg/cost, pkg/exf, pkg/ecco,
    verification/carbon and model/src/ related to adjoint
  o some new Matlab scripts for diagnosing model density
    - utils/matlab/dens_poly3.m and ini_poly3.m

The list of exclusions is accurate based on a "cvs diff". The list of
inclusions is based on the record in doc/tag-index which may not be complete.

1 C $Header: /u/gcmpack/MITgcm/verification/aim.5l_cs/code/shap_filt_relvort3.F,v 1.1 2002/01/09 00:28:56 jmc Exp $
2 C $Name: checkpoint43a-release1mods $
3
4 #include "CPP_OPTIONS.h"
5
6 SUBROUTINE SHAP_FILT_RELVORT3(
7 I bi,bj,k,
8 I uFld, vFld, hFacZ,
9 O vort3,
10 I myThid)
11 IMPLICIT NONE
12 C /==========================================================\
13 C | S/R SHAP_FILT_RELVORT3 |
14 C |==========================================================|
15 C \==========================================================/
16
17 C == Global variables ==
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "GRID.h"
22 C == Routine arguments ==
23 C myThid - Instance number for this innvocation of CALC_MOM_RHS
24 INTEGER bi,bj,k
25 _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
26 _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
27 _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28 _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29 INTEGER myThid
30
31 #ifdef ALLOW_SHAP_FILT
32 C == Local variables ==
33 INTEGER i,j
34 _RL AZcorner
35
36 DO J=2-Oly,sNy+Oly
37 DO I=2-Olx,sNx+Olx
38
39 C Horizontal curl of flow field - ignoring lopping factors
40 vort3(I,J)=
41 & recip_rAz(I,J,bi,bj)*(
42 & vFld(I,J)*dyc(I,J,bi,bj)
43 & -vFld(I-1,J)*dyc(I-1,J,bi,bj)
44 & -uFld(I,J)*dxc(I,J,bi,bj)
45 & +uFld(I,J-1)*dxc(I,J-1,bi,bj)
46 & )
47
48 C Horizontal curl of flow field - including lopping factors
49 c IF (hFacZ(i,j).NE.0.) THEN
50 c vort3(I,J)=
51 c & recip_rAz(I,J,bi,bj)*(
52 c & vFld(I,J)*dyc(I,J,bi,bj)*_hFacW(i,j,k,bi,bj)
53 c & -vFld(I-1,J)*dyc(I-1,J,bi,bj)*_hFacW(i-1,j,k,bi,bj)
54 c & -uFld(I,J)*dxc(I,J,bi,bj)*_hFacS(i,j,k,bi,bj)
55 c & +uFld(I,J-1)*dxc(I,J-1,bi,bj)*_hFacS(i,j-1,k,bi,bj)
56 c & )
57 c & /hFacZ(i,j)
58 c ELSE
59 c vort3(I,J)=0.
60 c ENDIF
61
62 ENDDO
63 ENDDO
64
65 C Special stuff for Cubed Sphere
66 IF (useCubedSphereExchange) THEN
67 AZcorner = 0.75 _d 0
68 I=1
69 J=1
70 vort3(I,J)=
71 c & +recip_rAz(I,J,bi,bj)*(
72 & +recip_rA(I,J,bi,bj)/AZcorner*(
73 & vFld(I,J)*dyc(I,J,bi,bj)
74 & -uFld(I,J)*dxc(I,J,bi,bj)
75 & +uFld(I,J-1)*dxc(I,J-1,bi,bj)
76 & )
77 cph & -vFld(I-1,J)*dyc(I-1,J,bi,bj)
78 I=sNx+1
79 J=1
80 vort3(I,J)=
81 c & +recip_rAz(I,J,bi,bj)*(
82 & +recip_rA(I-1,J,bi,bj)/AZcorner*(
83 & -vFld(I-1,J)*dyc(I-1,J,bi,bj)
84 & -uFld(I,J)*dxc(I,J,bi,bj)
85 & +uFld(I,J-1)*dxc(I,J-1,bi,bj)
86 & )
87 cph & vFld(I,J)*dyc(I,J,bi,bj)
88 I=1
89 J=sNy+1
90 vort3(I,J)=
91 c & +recip_rAz(I,J,bi,bj)*(
92 & +recip_rA(I,J-1,bi,bj)/AZcorner*(
93 & vFld(I,J)*dyc(I,J,bi,bj)
94 & -uFld(I,J)*dxc(I,J,bi,bj)
95 & +uFld(I,J-1)*dxc(I,J-1,bi,bj)
96 & )
97 cph & -vFld(I-1,J)*dyc(I-1,J,bi,bj)
98 I=sNx+1
99 J=sNy+1
100 vort3(I,J)=
101 c & +recip_rAz(I,J,bi,bj)*(
102 & +recip_rA(I-1,J-1,bi,bj)/AZcorner*(
103 & -vFld(I-1,J)*dyc(I-1,J,bi,bj)
104 & -uFld(I,J)*dxc(I,J,bi,bj)
105 & +uFld(I,J-1)*dxc(I,J-1,bi,bj)
106 & )
107 cph & vFld(I,J)*dyc(I,J,bi,bj)
108 ENDIF
109
110 #endif /* ALLOW_SHAP_FILT */
111
112 RETURN
113 END

  ViewVC Help
Powered by ViewVC 1.1.22