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

Annotation 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 - (hide annotations) (download)
Tue Feb 26 16:05:07 2002 UTC (22 years, 2 months 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 adcroft 1.1.4.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 jmc 1.1
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