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

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

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


Revision 1.1 - (hide annotations) (download)
Tue Nov 8 02:13:17 2005 UTC (18 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58e_post, checkpoint57y_post, checkpoint58h_post, checkpoint57y_pre, checkpoint58j_post, checkpoint58, checkpoint58f_post, checkpoint57x_post, checkpoint58d_post, checkpoint58c_post, checkpoint58a_post, checkpoint58i_post, checkpoint58g_post, checkpoint57z_post, checkpoint58b_post
initialise NH-variables in specific S/R: ini_nh_vars.F

1 jmc 1.1 C $Header: $
2     C $Name: $
3    
4     #include "CPP_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: INI_NH_VARS
8     C !INTERFACE:
9     SUBROUTINE INI_NH_VARS ( myThid )
10     C !DESCRIPTION: \bv
11     C *==========================================================*
12     C | SUBROUTINE INI_NH_VARS
13     C | o Initialise to zero all NH_VARS.h arrays
14     C *==========================================================*
15     C | Sets all the NH State variables & tendencies to zero.
16     C *==========================================================*
17     C \ev
18    
19     C !USES:
20     IMPLICIT NONE
21     C === Global variables ===
22     #include "SIZE.h"
23     #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "NH_VARS.h"
26    
27     C !INPUT/OUTPUT PARAMETERS:
28     C == Routine arguments ==
29     C myThid :: My Thread Id number
30     INTEGER myThid
31     CEOP
32    
33     #ifdef ALLOW_NONHYDROSTATIC
34     C !LOCAL VARIABLES:
35     C == Local variables ==
36     C bi,bj :: tile indices
37     C i,j,k :: loop counters
38     INTEGER bi, bj
39     INTEGER i, j, k
40    
41     c IF ( nonHydrostatic ) THEN
42     C- Note: comment out this IF so that we are sure that all variables in
43     C common blocs (even if never used later) are always initialised
44     DO bj = myByLo(myThid), myByHi(myThid)
45     DO bi = myBxLo(myThid), myBxHi(myThid)
46     DO k=1,Nr
47     DO j=1-Oly,sNy+Oly
48     DO i=1-Olx,sNx+Olx
49     phi_nh(i,j,k,bi,bj) = 0. _d 0
50     gW (i,j,k,bi,bj) = 0. _d 0
51     gWnm1(i,j,k,bi,bj) = 0. _d 0
52     ENDDO
53     ENDDO
54     ENDDO
55     ENDDO
56     ENDDO
57     c ENDIF
58    
59     #endif /* ALLOW_NONHYDROSTATIC */
60    
61     RETURN
62     END

  ViewVC Help
Powered by ViewVC 1.1.22