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

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

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

Revision 1.5 - (show annotations) (download)
Wed Jan 30 04:12:12 2002 UTC (22 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint51k_post, checkpoint47e_post, checkpoint44e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, checkpoint50c_post, checkpoint46f_post, checkpoint48e_post, checkpoint50c_pre, checkpoint44f_post, checkpoint46b_post, checkpoint51o_pre, checkpoint51l_post, checkpoint48i_post, checkpoint46l_pre, chkpt44d_post, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint50b_pre, checkpoint44e_pre, checkpoint51f_post, checkpoint48b_post, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint51n_post, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint51n_pre, checkpoint47d_post, checkpoint46d_pre, checkpoint48d_post, checkpoint48f_post, checkpoint45d_post, checkpoint46j_pre, checkpoint51l_pre, chkpt44a_post, checkpoint44h_pre, checkpoint48h_post, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, chkpt44c_pre, checkpoint48a_post, checkpoint45a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint47j_post, branch-exfmods-tag, checkpoint44g_post, branchpoint-genmake2, checkpoint46e_pre, checkpoint48c_post, checkpoint45b_post, checkpoint46b_pre, checkpoint51i_post, release1_final_v1, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint46, checkpoint47b_post, checkpoint44b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint50g_post, checkpoint45c_post, checkpoint44h_post, checkpoint46g_post, checkpoint50h_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint47f_post, checkpoint50e_post, chkpt44a_pre, checkpoint46i_post, checkpoint46c_post, checkpoint50d_pre, checkpoint46e_post, checkpoint51e_post, checkpoint44b_pre, checkpoint47, checkpoint45, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint51o_post, checkpoint51f_pre, chkpt44c_post, checkpoint48g_post, checkpoint47h_post, checkpoint44f_pre, checkpoint51g_post, checkpoint46d_post, checkpoint50b_post, checkpoint51m_post, checkpoint51a_post
Branch point for: branch-exfmods-curt, release1_final, branch-genmake2, tg2-branch, checkpoint51n_branch
Changes since 1.4: +14 -1 lines 
small changes associated with NonLin_FreeSurf option :
* initialization (ini_psurf.F);
* dump hFac fields (write_state.F);
* avoid unnecessary re-computation (forward_step.F initialise_varia.F calc_surf_dr.F);
1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_psurf.F,v 1.4 2001/09/26 18:09:15 cnh Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: INI_PSURF
8 C !INTERFACE:
9 SUBROUTINE INI_PSURF( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE INI_PSURF |
14 C | o Set model initial free-surface height/pressure. |
15 C *==========================================================*
16 C | There are several options for setting the initial |
17 C | surface displacement (r unit) field. |
18 C | 1. Inline code |
19 C | 2. Two-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 "DYNVARS.h"
31
32 C !INPUT/OUTPUT PARAMETERS:
33 C == Routine arguments ==
34 C myThid - Number of this instance of INI_PSURF
35 INTEGER myThid
36
37 C !LOCAL VARIABLES:
38 C == Local variables ==
39 C bi,bj - Loop counters
40 C I,J
41 INTEGER bi, bj
42 INTEGER I, J
43 CEOP
44
45 C-- Initialise surface position anomaly to zero
46 DO bj = myByLo(myThid), myByHi(myThid)
47 DO bi = myBxLo(myThid), myBxHi(myThid)
48 DO J=1-Oly,sNy+Oly
49 DO I=1-Olx,sNx+Olx
50 etaN(I,J,bi,bj) = 0. _d 0
51 ENDDO
52 ENDDO
53 ENDDO
54 ENDDO
55 C Read an initial state
56 IF (pSurfInitFile .NE. ' ') THEN
57 _BEGIN_MASTER( myThid )
58 CALL READ_FLD_XY_RL( pSurfInitFile, ' ', etaN, 0, myThid )
59 _END_MASTER(myThid)
60 ENDIF
61 C
62 _EXCH_XY_R8(etaN, myThid)
63
64 #ifdef INCLUDE_CD_CODE
65 C-- By default, initialize etaNm1 with etaN :
66 DO bj=myByLo(myThid),myByHi(myThid)
67 DO bi=myBxLo(myThid),myBxHi(myThid)
68 DO J=1-Oly,sNy+Oly
69 DO I=1-Olx,sNx+Olx
70 etaNm1(I,J,bi,bj) = etaN(I,J,bi,bj)
71 ENDDO
72 ENDDO
73 ENDDO
74 ENDDO
75
76 C _EXCH_XY_R8(etaNm1, myThid)
77 #endif
78
79 #ifdef NONLIN_FRSURF
80 C-- By default, initialize etaH with etaN :
81 DO bj=myByLo(myThid),myByHi(myThid)
82 DO bi=myBxLo(myThid),myBxHi(myThid)
83 DO J=1-Oly,sNy+Oly
84 DO I=1-Olx,sNx+Olx
85 etaH(I,J,bi,bj) = etaN(I,J,bi,bj)
86 ENDDO
87 ENDDO
88 ENDDO
89 ENDDO
90 #endif
91
92 RETURN
93 END
  ViewVC Help
Powered by ViewVC 1.1.22