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

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

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


Revision 1.4 - (show annotations) (download)
Wed Dec 22 00:05:31 2010 UTC (13 years, 6 months ago) by jmc
Branch: MAIN
Changes since 1.3: +13 -157 lines
read grid-spacing and ref.profile using S/R READ_GLVEC_RL.

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

  ViewVC Help
Powered by ViewVC 1.1.22