/[MITgcm]/MITgcm/pkg/shap_filt/shap_filt_tracer_s2.F
ViewVC logotype

Diff of /MITgcm/pkg/shap_filt/shap_filt_tracer_s2.F

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

revision 1.4 by jmc, Mon Mar 4 01:32:55 2002 UTC revision 1.5 by jmc, Sat Jun 26 01:15:30 2004 UTC
# Line 8  C     !ROUTINE: SHAP_FILT_TRACER_S2 Line 8  C     !ROUTINE: SHAP_FILT_TRACER_S2
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE SHAP_FILT_TRACER_S2(        SUBROUTINE SHAP_FILT_TRACER_S2(
10       U           field, tmpFld,       U           field, tmpFld,
11       I           kSize, myTime, myThid )       I           nShapTr, kSize, myTime, myThid )
12  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
13  C     *==========================================================*  C     *==========================================================*
14  C     | S/R SHAP_FILT_TRACER_S2  C     | S/R SHAP_FILT_TRACER_S2
# Line 33  C     !INPUT/OUTPUT PARAMETERS: Line 33  C     !INPUT/OUTPUT PARAMETERS:
33  C     == Routine arguments  C     == Routine arguments
34  C     field :: cell-centered 2D field on which filter applies  C     field :: cell-centered 2D field on which filter applies
35  C     tmpFld :: working temporary array  C     tmpFld :: working temporary array
36    C     nShapTr :: (total) power of the filter for this tracer
37  C     kSize :: length of 3rd Dim : either =1 (2D field) or =Nr (3D field)  C     kSize :: length of 3rd Dim : either =1 (2D field) or =Nr (3D field)
38  C     myTime :: Current time in simulation  C     myTime :: Current time in simulation
39  C     myThid :: Thread number for this instance of SHAP_FILT_TRACER_S2  C     myThid :: Thread number for this instance of SHAP_FILT_TRACER_S2
40        INTEGER kSize        INTEGER nShapTr, kSize
41        _RL field(1-OLx:sNx+OLx,1-OLy:sNy+OLy,kSize,nSx,nSy)        _RL field(1-OLx:sNx+OLx,1-OLy:sNy+OLy,kSize,nSx,nSy)
42        _RL tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,kSize,nSx,nSy)        _RL tmpFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,kSize,nSx,nSy)
43        _RL     myTime        _RL     myTime
# Line 51  C     == Local variables == Line 52  C     == Local variables ==
52        _RL tmpGrd(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL tmpGrd(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
53  CEOP  CEOP
54    
55        IF (nShapT.gt.0) THEN        IF (nShapTr.gt.0) THEN
56  C-------  C-------
57  C  Apply computational filter ^(nShap-nShapPhys) without grid factor  C  Apply computational filter ^(nShap-nShapPhys) without grid factor
58  C  then apply Physical filter ^nShapPhys  with grid factors  C  then apply Physical filter ^nShapPhys  with grid factors
59  C-------  C-------
60          nShapComput = nShapT - nShapTrPhys          nShapComput = nShapTr - nShapTrPhys
61    
62          DO bj=myByLo(myThid),myByHi(myThid)          DO bj=myByLo(myThid),myByHi(myThid)
63           DO bi=myBxLo(myThid),myBxHi(myThid)           DO bi=myBxLo(myThid),myBxHi(myThid)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22