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

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

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

Revision 1.6 - (show annotations) (download)
Fri Aug 4 00:40:43 2006 UTC (17 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint59, checkpoint58y_post, checkpoint58t_post, checkpoint60, checkpoint61, checkpoint58w_post, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58r_post, checkpoint58n_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, mitgcm_mapl_00, checkpoint58v_post, checkpoint61f, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint61q, checkpoint61z, checkpoint61e, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61x, checkpoint61y
Changes since 1.5: +10 -13 lines 
change commented lines (in case we read initial Phi_nh from a file)
1 C $Header: /u/gcmpack/MITgcm/model/src/ini_pnh.F,v 1.5 2005/11/08 02:14:10 jmc Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: INI_PNH
8 C !INTERFACE:
9 SUBROUTINE INI_PNH( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE INI_PNH
14 C | o Set model initial non-hydrostatic potential.
15 C *==========================================================*
16 C | There are several options for setting the initial
17 C | NH-potential (=Pressure/rhoConst) anomaly:
18 C | 1. Inline code
19 C | 2. Three-dimensional data from a file.
20 C *==========================================================*
21 C \ev
22
23 C !USES:
24 IMPLICIT NONE
25 C === Global variables ===
26 #include "SIZE.h"
27 #include "EEPARAMS.h"
28 #include "PARAMS.h"
29 #include "GRID.h"
30 #include "NH_VARS.h"
31
32 C !INPUT/OUTPUT PARAMETERS:
33 C == Routine arguments ==
34 C myThid :: My Thread Id number
35 INTEGER myThid
36
37 #ifdef ALLOW_NONHYDROSTATIC
38 C !LOCAL VARIABLES:
39 C == Local variables ==
40 INTEGER bi,bj,I,J,K
41 CEOP
42
43
44 C-- Initialise surface position anomaly to zero
45 DO bj = myByLo(myThid), myByHi(myThid)
46 DO bi = myBxLo(myThid), myBxHi(myThid)
47 DO K = 1, Nr
48 DO J=1-Oly,sNy+Oly
49 DO I=1-Olx,sNx+Olx
50 phi_nh(I,J,K,bi,bj) = 0. _d 0
51 ENDDO
52 ENDDO
53 ENDDO
54 ENDDO
55 ENDDO
56
57 C-- Read an initial state
58 c IF (phiNHinitFile .NE. ' ') THEN
59 c _BARRIER
60 c CALL READ_FLD_XY_RL( phiNHinitFile, ' ', phi_nh, 0, myThid )
61 c _EXCH_XYZ_RL(phi_nh, myThid)
62 c ENDIF
63
64 #endif /* ALLOW_NONHYDROSTATIC */
65 RETURN
66 END
  ViewVC Help
Powered by ViewVC 1.1.22