88 |
CALL READ_REC_XYZ_RL( 'gsNm1', gSNm1, 1,myIter, myThid) |
CALL READ_REC_XYZ_RL( 'gsNm1', gSNm1, 1,myIter, myThid) |
89 |
CALL READ_REC_XY_RL ('etaN', etaN, 1,myIter, myThid) |
CALL READ_REC_XY_RL ('etaN', etaN, 1,myIter, myThid) |
90 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
91 |
CALL READ_REC_XY_RL |
CALL READ_REC_XY_RL ('etaNm1', etaNm1, 1,myIter, myThid) |
|
& ( 'etaNm1', etaNm1, 1,myIter, myThid) |
|
92 |
CALL READ_REC_XYZ_RL( 'uVelD', uVelD, 1,myIter, myThid) |
CALL READ_REC_XYZ_RL( 'uVelD', uVelD, 1,myIter, myThid) |
93 |
CALL READ_REC_XYZ_RL( 'vVelD', vVelD, 1,myIter, myThid) |
CALL READ_REC_XYZ_RL( 'vVelD', vVelD, 1,myIter, myThid) |
94 |
CALL READ_REC_XYZ_RL( 'uNm1', uNM1, 1,myIter, myThid) |
CALL READ_REC_XYZ_RL( 'uNm1', uNM1, 1,myIter, myThid) |
124 |
CALL MDSREADFIELD(fn,prec,'RL', 1,etaN,12*Nr+1,myThid) |
CALL MDSREADFIELD(fn,prec,'RL', 1,etaN,12*Nr+1,myThid) |
125 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
126 |
IF ( nonlinFreeSurf.GE.0) |
IF ( nonlinFreeSurf.GE.0) |
127 |
& CALL MDSREADFIELD(fn,prec,'RL',1,etaNm1,12*Nr+2,myThid) |
& CALL MDSREADFIELD(fn,prec,'RL',1,etaH,12*Nr+2,myThid) |
128 |
#endif |
#endif |
129 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
130 |
WRITE(fn,'(A,I10.10)') 'pickup_cd.',myIter |
WRITE(fn,'(A,I10.10)') 'pickup_cd.',myIter |
134 |
CALL MDSREADFIELD(fn,prec,'RL',Nr,vNM1, 4,myThid) |
CALL MDSREADFIELD(fn,prec,'RL',Nr,vNM1, 4,myThid) |
135 |
CALL MDSREADFIELD(fn,prec,'RL',Nr,guCD, 5,myThid) |
CALL MDSREADFIELD(fn,prec,'RL',Nr,guCD, 5,myThid) |
136 |
CALL MDSREADFIELD(fn,prec,'RL',Nr,gvCD, 6,myThid) |
CALL MDSREADFIELD(fn,prec,'RL',Nr,gvCD, 6,myThid) |
|
#ifdef NONLIN_FRSURF |
|
|
IF ( nonlinFreeSurf.LT.0) |
|
|
& CALL MDSREADFIELD(fn,prec,'RL', 1,etaNm1,6*Nr+1,myThid) |
|
|
#else |
|
137 |
CALL MDSREADFIELD(fn,prec,'RL', 1,etaNm1,6*Nr+1,myThid) |
CALL MDSREADFIELD(fn,prec,'RL', 1,etaNm1,6*Nr+1,myThid) |
|
#endif |
|
138 |
#endif /* INCLUDE_CD_CODE */ |
#endif /* INCLUDE_CD_CODE */ |
139 |
|
|
140 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
178 |
_EXCH_XYZ_R8(gs , myThid ) |
_EXCH_XYZ_R8(gs , myThid ) |
179 |
_EXCH_XYZ_R8(gsNM1 , myThid ) |
_EXCH_XYZ_R8(gsNM1 , myThid ) |
180 |
_EXCH_XY_R8 (etaN, myThid ) |
_EXCH_XY_R8 (etaN, myThid ) |
181 |
_EXCH_XY_R8( etaNm1, myThid ) |
_EXCH_XY_R8( etaH, myThid ) |
182 |
|
|
183 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
184 |
c**** CALL EXCH_DUV_XYZ_RL(uVelD,vVelD,.TRUE.,myThid) |
c**** CALL EXCH_DUV_XYZ_RL(uVelD,vVelD,.TRUE.,myThid) |
190 |
c _EXCH_XYZ_R8( vNM1, myThid ) |
c _EXCH_XYZ_R8( vNM1, myThid ) |
191 |
_EXCH_XYZ_R8( guCD, myThid ) |
_EXCH_XYZ_R8( guCD, myThid ) |
192 |
_EXCH_XYZ_R8( gvCD, myThid ) |
_EXCH_XYZ_R8( gvCD, myThid ) |
193 |
|
_EXCH_XY_R8( etaNm1, myThid ) |
194 |
#endif |
#endif |
195 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
196 |
IF ( nonHydrostatic ) THEN |
IF ( nonHydrostatic ) THEN |
356 |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaN,12*Nr+1, |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaN,12*Nr+1, |
357 |
& myIter,myThid) |
& myIter,myThid) |
358 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
359 |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaNm1,12*Nr+2, |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaH,12*Nr+2, |
360 |
& myIter,myThid) |
& myIter,myThid) |
361 |
#endif |
#endif |
362 |
#ifdef INCLUDE_CD_CODE |
#ifdef INCLUDE_CD_CODE |
371 |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,vNM1, 4,myIter,myThid) |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,vNM1, 4,myIter,myThid) |
372 |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,guCD, 5,myIter,myThid) |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,guCD, 5,myIter,myThid) |
373 |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,gvCD, 6,myIter,myThid) |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,gvCD, 6,myIter,myThid) |
|
#ifndef NONLIN_FRSURF |
|
374 |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaNm1,6*Nr+1, |
CALL MDSWRITEFIELD(fn,prec,lgf,'RL', 1,etaNm1,6*Nr+1, |
375 |
& myIter,myThid) |
& myIter,myThid) |
|
#endif |
|
376 |
#endif /* INCLUDE_CD_CODE */ |
#endif /* INCLUDE_CD_CODE */ |
377 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
378 |
IF ( nonHydrostatic ) THEN |
IF ( nonHydrostatic ) THEN |