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

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

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

Revision 1.2 - (hide 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 gforget 1.2 C $Header: /u/gcmpack/MITgcm/model/src/taueddy_init_varia.F,v 1.1 2006/02/15 03:52:53 heimbach Exp $
2 heimbach 1.1 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 gforget 1.2 #ifdef ALLOW_EDDYPSI
33 heimbach 1.1
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 gforget 1.2 eddyPsiX(i,j,k,bi,bj)=0.
49     eddyPsiY(i,j,k,bi,bj)=0.
50 heimbach 1.1 ENDDO
51     ENDDO
52     ENDDO
53     ENDDO
54     ENDDO
55    
56 gforget 1.2 IF ( eddyPsixFile .NE. ' ' ) THEN
57     CALL MDSREADFIELD(eddyPsiXFile,precFloat32,
58     & 'RS',Nr,eddyPsiX,1,myThid)
59 heimbach 1.1 ENDIF
60 gforget 1.2 IF ( eddyPsiyFile .NE. ' ' ) THEN
61     CALL MDSREADFIELD(eddyPsiYFile,precFloat32,
62     & 'RS',Nr,eddyPsiY,1,myThid)
63 heimbach 1.1 ENDIF
64    
65 gforget 1.2 CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
66 heimbach 1.1
67     #endif /* ALLOW_AUTODIFF_TAMC */
68    
69     END
  ViewVC Help
Powered by ViewVC 1.1.22