/[MITgcm]/MITgcm/pkg/longstep/longstep_average_3d_fac.F
ViewVC logotype

Contents of /MITgcm/pkg/longstep/longstep_average_3d_fac.F

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

Revision 1.2 - (show annotations) (download)
Mon Jan 18 01:02:56 2010 UTC (14 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint63, checkpoint65p, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62b, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x
Changes since 1.1: +3 -3 lines 
use the right type (RS) for argument "fac'
1 C $Header: /u/gcmpack/MITgcm/pkg/longstep/longstep_average_3d_fac.F,v 1.1 2010/01/12 23:55:48 jahn Exp $
2 C $Name: $
3
4 #include "LONGSTEP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: LONGSTEP_AVERAGE_3D_FAC
8 C !INTERFACE:
9 SUBROUTINE LONGSTEP_AVERAGE_3D_FAC(
10 I lsCnt,
11 U lsFld,
12 I fac,
13 I myNr, myThid)
14
15 C !DESCRIPTION: \bv
16 C *==========================================================*
17 C | SUBROUTINE LONGSTEP_AVERAGE_3D_FAC
18 C | o compute longstep average for one field (incl. overlap)
19 C *==========================================================*
20 C \ev
21
22 C !USES:
23 IMPLICIT NONE
24 C == Global variables ===
25 #include "SIZE.h"
26 #include "EEPARAMS.h"
27 #include "PARAMS.h"
28 #include "GRID.h"
29
30 C !INPUT/OUTPUT PARAMETERS:
31 C == Routine arguments ==
32 C lsCnt :: Counter to update
33 C lsFld :: average to update
34 C myNr :: number of levels to update
35 C myThid :: my Thread Id number
36 INTEGER myNr
37 _RL lsFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
38 _RS fac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
39 INTEGER lsCnt(nSx,nSy)
40 INTEGER myThid
41
42 #ifdef ALLOW_LONGSTEP
43 C !LOCAL VARIABLES:
44 C == Local variables ==
45 INTEGER i,j,k,bi,bj
46 _RL factor
47 CEOP
48
49 DO bj=myByLo(myThid), myByHi(myThid)
50 DO bi=myBxLo(myThid), myBxHi(myThid)
51 DO k = 1,myNr
52 DO j = 1-OLy,sNy+OLy
53 DO i = 1-OLx,sNx+OLx
54 factor = FLOAT(lsCnt(bi,bj))*fac(i,j,k,bi,bj)
55 IF (factor.NE.0. _d 0) factor = 1. _d 0 / factor
56 lsFld(i,j,k,bi,bj) = lsFld(i,j,k,bi,bj)*factor
57 ENDDO
58 ENDDO
59 ENDDO
60 ENDDO
61 ENDDO
62 #endif /* ALLOW_LONGSTEP */
63
64 RETURN
65 END
66
  ViewVC Help
Powered by ViewVC 1.1.22