76 |
C read the raw grid data from files |
C read the raw grid data from files |
77 |
|
|
78 |
C- Cell centered quantities |
C- Cell centered quantities |
79 |
CALL MDSREADFIELD('LONC.bin',readBinaryPrec,'RL',1,XC, 1,myThid) |
CALL MDSREADFIELD('LONC.bin',readBinaryPrec,'RS',1,XC, 1,myThid) |
80 |
CALL MDSREADFIELD('LATC.bin',readBinaryPrec,'RL',1,YC, 1,myThid) |
CALL MDSREADFIELD('LATC.bin',readBinaryPrec,'RS',1,YC, 1,myThid) |
81 |
_EXCH_XY_R4(XC,myThid) |
_EXCH_XY_R4(XC,myThid) |
82 |
_EXCH_XY_R4(YC,myThid) |
_EXCH_XY_R4(YC,myThid) |
83 |
|
|
84 |
CALL MDSREADFIELD('DXF.bin',readBinaryPrec,'RL',1,DXF, 1,myThid) |
CALL MDSREADFIELD('DXF.bin',readBinaryPrec,'RS',1,DXF, 1,myThid) |
85 |
CALL MDSREADFIELD('DYF.bin',readBinaryPrec,'RL',1,DYF, 1,myThid) |
CALL MDSREADFIELD('DYF.bin',readBinaryPrec,'RS',1,DYF, 1,myThid) |
86 |
C !!! _EXCH_OUV_XY_R4(DXF, DYF, unSigned, myThid ) |
C !!! _EXCH_OUV_XY_R4(DXF, DYF, unSigned, myThid ) |
87 |
cs! this is not correct! <= need paired exchange for DXF,DYF |
cs! this is not correct! <= need paired exchange for DXF,DYF |
88 |
_EXCH_XY_R4(DXF,myThid) |
_EXCH_XY_R4(DXF,myThid) |
89 |
_EXCH_XY_R4(DYF,myThid) |
_EXCH_XY_R4(DYF,myThid) |
90 |
|
|
91 |
CALL MDSREADFIELD('RA.bin',readBinaryPrec,'RL',1,RA, 1,myThid) |
CALL MDSREADFIELD('RA.bin',readBinaryPrec,'RS',1,RA, 1,myThid) |
92 |
_EXCH_XY_R4(RA,myThid ) |
_EXCH_XY_R4(RA,myThid ) |
93 |
|
|
94 |
C- Corner quantities |
C- Corner quantities |
95 |
C *********** this are not degbugged ************ |
C *********** this are not degbugged ************ |
96 |
CALL MDSREADFIELD('LONG.bin',readBinaryPrec,'RL',1,XG, 1,myThid) |
CALL MDSREADFIELD('LONG.bin',readBinaryPrec,'RS',1,XG, 1,myThid) |
97 |
CALL MDSREADFIELD('LATG.bin',readBinaryPrec,'RL',1,YG, 1,myThid) |
CALL MDSREADFIELD('LATG.bin',readBinaryPrec,'RS',1,YG, 1,myThid) |
98 |
cs- this block needed by cubed sphere until we write more useful I/O routines |
cs- this block needed by cubed sphere until we write more useful I/O routines |
99 |
bi=3 |
bi=3 |
100 |
bj=1 |
bj=1 |
114 |
CALL EXCH_Z_XY_RS(XG,myThid) |
CALL EXCH_Z_XY_RS(XG,myThid) |
115 |
CALL EXCH_Z_XY_RS(YG,myThid) |
CALL EXCH_Z_XY_RS(YG,myThid) |
116 |
|
|
117 |
CALL MDSREADFIELD('DXV.bin',readBinaryPrec,'RL',1,DXV, 1,myThid) |
CALL MDSREADFIELD('DXV.bin',readBinaryPrec,'RS',1,DXV, 1,myThid) |
118 |
CALL MDSREADFIELD('DYU.bin',readBinaryPrec,'RL',1,DYU, 1,myThid) |
CALL MDSREADFIELD('DYU.bin',readBinaryPrec,'RS',1,DYU, 1,myThid) |
119 |
cs- this block needed by cubed sphere until we write more useful I/O routines |
cs- this block needed by cubed sphere until we write more useful I/O routines |
120 |
C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid) |
C !!! _EXCH_ZUV_XY_R4(DXV, DYU, unSigned, myThid) |
121 |
cs! this is not correct <= need paired exchange for dxv,dyu |
cs! this is not correct <= need paired exchange for dxv,dyu |
135 |
CALL EXCH_Z_XY_RS(DXV,myThid) |
CALL EXCH_Z_XY_RS(DXV,myThid) |
136 |
CALL EXCH_Z_XY_RS(DYU,myThid) |
CALL EXCH_Z_XY_RS(DYU,myThid) |
137 |
|
|
138 |
CALL MDSREADFIELD('RAZ.bin',readBinaryPrec,'RL',1,RAZ, 1,myThid) |
CALL MDSREADFIELD('RAZ.bin',readBinaryPrec,'RS',1,RAZ, 1,myThid) |
139 |
cs- this block needed by cubed sphere until we write more useful I/O routines |
cs- this block needed by cubed sphere until we write more useful I/O routines |
140 |
CALL EXCH_Z_XY_RS(RAZ , myThid ) |
CALL EXCH_Z_XY_RS(RAZ , myThid ) |
141 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
148 |
CALL EXCH_Z_XY_RS(RAZ,myThid) |
CALL EXCH_Z_XY_RS(RAZ,myThid) |
149 |
|
|
150 |
C- Staggered (u,v pairs) quantities |
C- Staggered (u,v pairs) quantities |
151 |
CALL MDSREADFIELD('DXC.bin',readBinaryPrec,'RL',1,DXC, 1,myThid) |
CALL MDSREADFIELD('DXC.bin',readBinaryPrec,'RS',1,DXC, 1,myThid) |
152 |
CALL MDSREADFIELD('DYC.bin',readBinaryPrec,'RL',1,DYC, 1,myThid) |
CALL MDSREADFIELD('DYC.bin',readBinaryPrec,'RS',1,DYC, 1,myThid) |
153 |
CALL EXCH_UV_XY_RS(DXC,DYC,.FALSE.,myThid) |
CALL EXCH_UV_XY_RS(DXC,DYC,.FALSE.,myThid) |
154 |
|
|
155 |
CALL MDSREADFIELD('RAW.bin',readBinaryPrec,'RL',1,RAW, 1,myThid) |
CALL MDSREADFIELD('RAW.bin',readBinaryPrec,'RS',1,RAW, 1,myThid) |
156 |
CALL MDSREADFIELD('RAS.bin',readBinaryPrec,'RL',1,RAS, 1,myThid) |
CALL MDSREADFIELD('RAS.bin',readBinaryPrec,'RS',1,RAS, 1,myThid) |
157 |
cs- this block needed by cubed sphere until we write more useful I/O routines |
cs- this block needed by cubed sphere until we write more useful I/O routines |
158 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
159 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
166 |
cs- end block |
cs- end block |
167 |
CALL EXCH_UV_XY_RS(RAW,RAS,.FALSE.,myThid) |
CALL EXCH_UV_XY_RS(RAW,RAS,.FALSE.,myThid) |
168 |
|
|
169 |
CALL MDSREADFIELD('DXG.bin',readBinaryPrec,'RL',1,DXG, 1,myThid) |
CALL MDSREADFIELD('DXG.bin',readBinaryPrec,'RS',1,DXG, 1,myThid) |
170 |
CALL MDSREADFIELD('DYG.bin',readBinaryPrec,'RL',1,DYG, 1,myThid) |
CALL MDSREADFIELD('DYG.bin',readBinaryPrec,'RS',1,DYG, 1,myThid) |
171 |
cs- this block needed by cubed sphere until we write more useful I/O routines |
cs- this block needed by cubed sphere until we write more useful I/O routines |
172 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
173 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |