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

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

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

Revision 1.15 - (show annotations) (download)
Sat Sep 10 18:30:06 2005 UTC (18 years, 8 months ago) by edhill
Branch: MAIN
Changes since 1.14: +2 -1 lines 
 o various changes to mnc including:
   - all files use the new "BASENAME[[.ITER].{t|f}NUM].nc" format
   - output can now be grouped so that all files within a group
       change the ITER portion of their names in lock-step together
   - can now read ("global") PER-FACE (in addition to PER-TILE) files
       and works with both EXCH1 and EXCH2 (but needs more testing)
   - writing works for all verification test cases w/ g77 on Linux
1 C $Header: /u/gcmpack/MITgcm/model/src/write_grid.F,v 1.14 2005/09/06 17:50:56 jmc Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
8 CBOP
9 C !ROUTINE: WRITE_GRID
10 C !INTERFACE:
11 SUBROUTINE WRITE_GRID(
12 I myThid )
13
14 C !DESCRIPTION:
15 C Writes the model geometry and grid arrays to file(s) using which
16 C ever I/O package is available/active. This includes horizontal
17 C and vertical grid and finite volume modulating factors (hFacs).
18 C Essentailly, it dumps the contents of GRID.h to disk.
19
20 C !CALLING SEQUENCE:
21 C INITIALIZE_FIXED
22 C |
23 C |-- WRITE_GRID
24
25 IMPLICIT NONE
26 #include "SIZE.h"
27 #include "EEPARAMS.h"
28 #include "PARAMS.h"
29 #ifdef ALLOW_MNC
30 #include "MNC_PARAMS.h"
31 #endif
32 #include "GRID.h"
33
34 C !INPUT/OUTPUT PARAMETERS:
35 C myThid - Number of this instances
36 INTEGER myThid
37
38 C !LOCAL VARIABLES:
39 character*(4) bfname
40 C tmpfld - Temporary array used to compute & write Total Depth
41 C has to be in common for multi threading
42 C ?aja: not sure why the COMMON block is necessary - should ask JMC?
43 c COMMON / LOCAL_INI_MASKS_ETC / tmpfld
44 _RS tmpfld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45 INTEGER i,j,bi,bj
46 CEOP
47
48 C Calculate the fluid thickness in R coordinates as seen by the
49 C model. This may differ from what the user specified due to
50 C partial cells, etc.
51 DO bj = myByLo(myThid), myByHi(myThid)
52 DO bi = myBxLo(myThid), myBxHi(myThid)
53 DO j=1-Oly,sNy+Oly
54 DO i=1-Olx,sNx+Olx
55 C Total fluid column thickness (r_unit) :
56 tmpfld(i,j,bi,bj) = Ro_surf(i,j,bi,bj) - R_low(i,j,bi,bj)
57 ENDDO
58 ENDDO
59 ENDDO
60 ENDDO
61
62 _BEGIN_MASTER( myThid )
63
64 C Write horizontal grid arrays
65 CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)
66 CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)
67 CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)
68 CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)
69 CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)
70 CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)
71 CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)
72 CALL WRITE_FLD_XY_RS( 'RAZ',' ',rAz,0,myThid)
73 CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)
74 CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)
75 CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)
76 CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)
77 IF ( usingCurvilinearGrid ) THEN
78 CALL WRITE_FLD_XY_RS( 'AngleCS',' ',angleCosC,0,myThid)
79 CALL WRITE_FLD_XY_RS( 'AngleSN',' ',angleSinC,0,myThid)
80 ENDIF
81 C Write 3D geometry arrays
82 CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,0,myThid)
83 CALL WRITE_FLD_XYZ_RS( 'hFacC',' ',hFacC,0,myThid)
84 CALL WRITE_FLD_XYZ_RS( 'hFacW',' ',hFacW,0,myThid)
85 CALL WRITE_FLD_XYZ_RS( 'hFacS',' ',hFacS,0,myThid)
86 IF (buoyancyRelation .EQ. 'ATMOSPHERIC')
87 & CALL WRITE_FLD_XY_RS( 'topo_P',' ',Ro_surf,0,myThid)
88
89 _END_MASTER(myThid)
90
91 #ifdef ALLOW_MNC
92 _BEGIN_MASTER( myThid )
93 IF (useMNC) THEN
94 bfname='grid'
95
96 C Write the GRID.h variables to a file
97 CALL MNC_CW_SET_UDIM(bfname, 0, myThid)
98 CALL MNC_CW_SET_CITER(bfname, 2, -1, -1, -1, myThid)
99 CALL MNC_CW_RL_W('D',bfname,0,0,'RC',rC,myThid)
100 CALL MNC_CW_RL_W('D',bfname,0,0,'RF',rF,myThid)
101 CALL MNC_CW_RL_W('D',bfname,0,0,'RU',rF(2),myThid)
102 CALL MNC_CW_RL_W('D',bfname,0,0,'RL',rF,myThid)
103 CALL MNC_CW_RL_W('D',bfname,0,0,'drC',drC,myThid)
104 CALL MNC_CW_RL_W('D',bfname,0,0,'drF',drF,myThid)
105 CALL MNC_CW_RL_W('D',bfname,0,0,'XC',xC,myThid)
106 CALL MNC_CW_RL_W('D',bfname,0,0,'YC',yC,myThid)
107 c CALL MNC_CW_RL_W('D',bfname,0,0,'XU',xG,myThid)
108 c CALL MNC_CW_RL_W('D',bfname,0,0,'YU',yC,myThid)
109 c CALL MNC_CW_RL_W('D',bfname,0,0,'XV',xC,myThid)
110 c CALL MNC_CW_RL_W('D',bfname,0,0,'YV',yG,myThid)
111 CALL MNC_CW_RL_W('D',bfname,0,0,'XG',xG,myThid)
112 CALL MNC_CW_RL_W('D',bfname,0,0,'YG',yG,myThid)
113 CALL MNC_CW_RL_W('D',bfname,0,0,'dxC',dxC,myThid)
114 CALL MNC_CW_RL_W('D',bfname,0,0,'dyC',dyC,myThid)
115 CALL MNC_CW_RL_W('D',bfname,0,0,'dxF',dxF,myThid)
116 CALL MNC_CW_RL_W('D',bfname,0,0,'dyF',dyF,myThid)
117 CALL MNC_CW_RL_W('D',bfname,0,0,'dxG',dxG,myThid)
118 CALL MNC_CW_RL_W('D',bfname,0,0,'dyG',dyG,myThid)
119 CALL MNC_CW_RL_W('D',bfname,0,0,'dxV',dxV,myThid)
120 CALL MNC_CW_RL_W('D',bfname,0,0,'dyU',dyU,myThid)
121 CALL MNC_CW_RL_W('D',bfname,0,0,'rA',rA,myThid)
122 CALL MNC_CW_RL_W('D',bfname,0,0,'rAw',rAw,myThid)
123 CALL MNC_CW_RL_W('D',bfname,0,0,'rAs',rAs,myThid)
124 CALL MNC_CW_RL_W('D',bfname,0,0,'rAz',rAz,myThid)
125 IF ( usingCurvilinearGrid ) THEN
126 CALL MNC_CW_RL_W('D',bfname,0,0,'AngleCS',angleCosC,myThid)
127 CALL MNC_CW_RL_W('D',bfname,0,0,'AngleSN',angleSinC,myThid)
128 ENDIF
129 CALL MNC_CW_RL_W('D',bfname,0,0,'fCori',fCori,myThid)
130 CALL MNC_CW_RL_W('D',bfname,0,0,'fCoriG',fCoriG,myThid)
131 CALL MNC_CW_RL_W('D',bfname,0,0,'R_low',R_low,myThid)
132 CALL MNC_CW_RL_W('D',bfname,0,0,'Ro_surf',Ro_surf,myThid)
133 CALL MNC_CW_RL_W('D',bfname,0,0,'Depth',tmpfld,myThid)
134 CALL MNC_CW_RL_W('D',bfname,0,0,'HFacC',HFacC,myThid)
135 CALL MNC_CW_RL_W('D',bfname,0,0,'HFacW',HFacW,myThid)
136 CALL MNC_CW_RL_W('D',bfname,0,0,'HFacS',HFacS,myThid)
137
138 ENDIF
139 _END_MASTER(myThid)
140 #endif /* ALLOW_MNC */
141
142 RETURN
143 END
144
145 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
146
  ViewVC Help
Powered by ViewVC 1.1.22