/[MITgcm]/MITgcm/model/src/taueddy_init_varia.F
ViewVC logotype

Contents of /MITgcm/model/src/taueddy_init_varia.F

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

Revision 1.2 - (show annotations) (download)
Fri May 30 02:45:43 2008 UTC (16 years ago) by gforget
Branch: MAIN
CVS Tags: checkpoint60, checkpoint61, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i
Changes since 1.1: +11 -30 lines 
o bridging the gap between eddy stress and GM.
  -> eddyTau is replaced with eddyPsi (eddyTau = f x rho0 x eddyPsi)
      along with a change in CPP option (now ALLOW_EDDYPSI).
  -> when using GM w/ GM_AdvForm:
      The total eddy streamfunction (Psi = eddyPsi + K x Slope)
      is applied either in the tracer Eq. or in momentum Eq.
      depending on data.gmredi (intro. GM_InMomAsStress).
  -> ALLOW_EDDYPSI_CONTROL for estimation purpose.
  The key modifications are in model/src/taueddy_external_forcing.F
  pkg/gmredi/gmredi_calc_*F pkg/gmredi/gmredi_*transport.F
1 C $Header: /u/gcmpack/MITgcm/model/src/taueddy_init_varia.F,v 1.1 2006/02/15 03:52:53 heimbach Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: TAUEDDY_INIT_VARIA
9 C !INTERFACE:
10 SUBROUTINE TAUEDDY_INIT_VARIA( myThid )
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE TAUEDDY_INIT_VARIA
14 C | o Initialise to zero some active arrays
15 C *==========================================================*
16 C \ev
17
18 C !USES:
19 IMPLICIT NONE
20 C === Global variables ===
21 #include "SIZE.h"
22 #include "EEPARAMS.h"
23 #include "PARAMS.h"
24 #include "GRID.h"
25 #include "FFIELDS.h"
26
27 C !INPUT/OUTPUT PARAMETERS:
28 C == Routine arguments ==
29 C myThid - Number of this instance of INI_UVEL
30 INTEGER myThid
31
32 #ifdef ALLOW_EDDYPSI
33
34 C !LOCAL VARIABLES:
35 C == Local variables ==
36 C bi,bj - Loop counters
37 C I,J,K
38 INTEGER bi, bj
39 INTEGER I, J, K
40 CEOP
41
42 C-- Over all tiles
43 DO bj = myByLo(myThid), myByHi(myThid)
44 DO bi = myBxLo(myThid), myBxHi(myThid)
45 DO K=1,Nr
46 DO J=1-Oly,sNy+Oly
47 DO I=1-Olx,sNx+Olx
48 eddyPsiX(i,j,k,bi,bj)=0.
49 eddyPsiY(i,j,k,bi,bj)=0.
50 ENDDO
51 ENDDO
52 ENDDO
53 ENDDO
54 ENDDO
55
56 IF ( eddyPsixFile .NE. ' ' ) THEN
57 CALL MDSREADFIELD(eddyPsiXFile,precFloat32,
58 & 'RS',Nr,eddyPsiX,1,myThid)
59 ENDIF
60 IF ( eddyPsiyFile .NE. ' ' ) THEN
61 CALL MDSREADFIELD(eddyPsiYFile,precFloat32,
62 & 'RS',Nr,eddyPsiY,1,myThid)
63 ENDIF
64
65 CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
66
67 #endif /* ALLOW_AUTODIFF_TAMC */
68
69 END
  ViewVC Help
Powered by ViewVC 1.1.22