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

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

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


Revision 1.5 - (hide annotations) (download)
Wed Dec 22 23:31:51 2010 UTC (13 years, 4 months ago) by jmc
Branch: MAIN
Changes since 1.4: +1 -2 lines
clean up previous modif (forgot to remove LOAD_GRID_SPACING_OLD)

1 jmc 1.5 C $Header: /u/gcmpack/MITgcm/model/src/load_grid_spacing.F,v 1.4 2010/12/22 00:05:31 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4 jmc 1.2 c #include "PACKAGES_CONFIG.h"
5 jmc 1.1 #include "CPP_OPTIONS.h"
6    
7     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
8     CBOP
9     C !ROUTINE: LOAD_GRID_SPACING
10     C !INTERFACE:
11     SUBROUTINE LOAD_GRID_SPACING( myThid )
12    
13     C !DESCRIPTION:
14 jmc 1.2 C load grid-spacing (vector array) delX, delY, delR or delRc from file.
15    
16 jmc 1.1 C !USES:
17     IMPLICIT NONE
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PARAMS.h"
21 jmc 1.3 #include "GRID.h"
22 jmc 1.1
23     C !INPUT/OUTPUT PARAMETERS:
24 jmc 1.2 C myThid :: my Thread Id. number
25 jmc 1.1 INTEGER myThid
26     CEOP
27    
28 jmc 1.4 C !FUNCTIONS:
29     INTEGER ILNBLNK
30     EXTERNAL ILNBLNK
31    
32 jmc 1.1 C !LOCAL VARIABLES:
33 jmc 1.3 C msgBuf :: Informational/error message buffer
34 jmc 1.4 INTEGER iLen
35 jmc 1.1 CHARACTER*(MAX_LEN_MBUF) msgBuf
36    
37 jmc 1.2 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
38    
39     _BEGIN_MASTER( myThid )
40    
41 jmc 1.3 C-- X coordinate
42 jmc 1.1 IF ( delXFile .NE. ' ' ) THEN
43     iLen = ILNBLNK(delXFile)
44 jmc 1.4 CALL READ_GLVEC_RL( delXFile, ' ', delX, Nx, 1, myThid )
45 jmc 1.1 WRITE(msgBuf,'(3A)') 'S/R LOAD_GRID_SPACING:',
46     & ' delX loaded from file: ', delXFile(1:iLen)
47     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
48 jmc 1.2 & SQUEEZE_RIGHT , myThid )
49 jmc 1.1 ENDIF
50    
51 jmc 1.3 C-- Y coordinate
52 jmc 1.1 IF ( delYFile .NE. ' ' ) THEN
53     iLen = ILNBLNK(delYFile)
54 jmc 1.4 CALL READ_GLVEC_RL( delYFile, ' ', delY, Ny, 1, myThid )
55 jmc 1.1 WRITE(msgBuf,'(3A)') 'S/R LOAD_GRID_SPACING:',
56     & ' delY loaded from file: ', delYFile(1:iLen)
57     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
58 jmc 1.2 & SQUEEZE_RIGHT , myThid )
59     ENDIF
60    
61 jmc 1.3 C-- vertical coordinate "R"
62 jmc 1.2 IF ( delRFile .NE. ' ' ) THEN
63     iLen = ILNBLNK(delRFile)
64 jmc 1.4 CALL READ_GLVEC_RL( delRFile, ' ', delR, Nr, 1, myThid )
65 jmc 1.2 WRITE(msgBuf,'(3A)') 'S/R LOAD_GRID_SPACING:',
66     & ' delR loaded from file: ', delRFile(1:iLen)
67     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
68     & SQUEEZE_RIGHT , myThid )
69     ENDIF
70    
71     IF ( delRcFile .NE. ' ' ) THEN
72     iLen = ILNBLNK(delRcFile)
73 jmc 1.4 CALL READ_GLVEC_RL( delRcFile, ' ', delRc, Nr+1, 1, myThid )
74 jmc 1.2 WRITE(msgBuf,'(3A)') 'S/R LOAD_GRID_SPACING:',
75     & ' delRc loaded from file: ', delRcFile(1:iLen)
76     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
77     & SQUEEZE_RIGHT , myThid )
78 jmc 1.1 ENDIF
79    
80 jmc 1.3 C-- hybrid sigma vertical coordinate coefficient
81     IF ( hybSigmFile .NE. ' ' ) THEN
82     iLen = ILNBLNK(hybSigmFile)
83 jmc 1.4 CALL READ_GLVEC_RL( hybSigmFile,' ',aHybSigmF,Nr+1, 1,myThid )
84     CALL READ_GLVEC_RL( hybSigmFile,' ',bHybSigmF,Nr+1, 2,myThid )
85 jmc 1.3 WRITE(msgBuf,'(3A)') 'S/R LOAD_GRID_SPACING:',
86     & ' a&b_HybSigmF loaded from file: ', hybSigmFile(1:iLen)
87     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
88     & SQUEEZE_RIGHT , myThid )
89     ENDIF
90    
91 jmc 1.2 _END_MASTER(myThid)
92 jmc 1.1 C-- Everyone else must wait for the parameters to be loaded
93     _BARRIER
94    
95     RETURN
96     END

  ViewVC Help
Powered by ViewVC 1.1.22