168 |
mnc_iter = myIter |
mnc_iter = myIter |
169 |
|
|
170 |
C Write dynvars using the MNC package |
C Write dynvars using the MNC package |
171 |
CALL MNC_CW_RL_W_D('pickup',0,0,'iter',-1,mnc_iter, myThid) |
CALL MNC_CW_SET_UDIM('pickup', -1, myThid) |
172 |
CALL MNC_CW_RL_W_D('pickup',0,0,'U', 0, uVel, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'iter',mnc_iter, myThid) |
173 |
CALL MNC_CW_RL_W_D('pickup',0,0,'V', 0, vVel, myThid) |
CALL MNC_CW_SET_UDIM('pickup', 0, myThid) |
174 |
CALL MNC_CW_RL_W_D('pickup',0,0,'T', 0, theta, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'U', uVel, myThid) |
175 |
CALL MNC_CW_RL_W_D('pickup',0,0,'S', 0, salt, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'V', vVel, myThid) |
176 |
CALL MNC_CW_RL_W_D('pickup',0,0,'Eta', 0, etaN, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'T', theta, myThid) |
177 |
|
CALL MNC_CW_RL_W_D('pickup',0,0,'S', salt, myThid) |
178 |
|
CALL MNC_CW_RL_W_D('pickup',0,0,'Eta', etaN, myThid) |
179 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
180 |
IF ( nonlinFreeSurf.GE.0) |
IF ( nonlinFreeSurf.GE.0) |
181 |
& CALL MNC_CW_RL_W_D('pickup',0,0,'EtaNH', 0, etaH, myThid) |
& CALL MNC_CW_RL_W_D('pickup',0,0,'EtaNH', etaH, myThid) |
182 |
#endif |
#endif |
183 |
CALL MNC_CW_RL_W_D('pickup',0,0,'Unm1', 0, gUnm1, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'Unm1', gUnm1, myThid) |
184 |
CALL MNC_CW_RL_W_D('pickup',0,0,'Vnm1', 0, gVnm1, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'Vnm1', gVnm1, myThid) |
185 |
CALL MNC_CW_RL_W_D('pickup',0,0,'Tnm1', 0, gTnm1, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'Tnm1', gTnm1, myThid) |
186 |
CALL MNC_CW_RL_W_D('pickup',0,0,'Snm1', 0, gSnm1, myThid) |
CALL MNC_CW_RL_W_D('pickup',0,0,'Snm1', gSnm1, myThid) |
187 |
|
|
188 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
189 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
220 |
|
|
221 |
C Read variables from the pickup file |
C Read variables from the pickup file |
222 |
CALL MNC_FILE_CLOSE_ALL_MATCHING('pickup', myThid) |
CALL MNC_FILE_CLOSE_ALL_MATCHING('pickup', myThid) |
223 |
CALL MNC_CW_RL_R_D('pickup',0,0,'iter',1,mnc_iter, myThid) |
CALL MNC_CW_SET_UDIM('pickup', 1, myThid) |
224 |
CALL MNC_CW_RL_R_D('pickup',0,0,'U',1,uVel, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'iter',mnc_iter, myThid) |
225 |
CALL MNC_CW_RL_R_D('pickup',0,0,'V',1,vVel, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'U',uVel, myThid) |
226 |
CALL MNC_CW_RL_R_D('pickup',0,0,'T',1,theta, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'V',vVel, myThid) |
227 |
CALL MNC_CW_RL_R_D('pickup',0,0,'S',1,salt, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'T',theta, myThid) |
228 |
CALL MNC_CW_RL_R_D('pickup',0,0,'Eta',1,etaN, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'S',salt, myThid) |
229 |
|
CALL MNC_CW_RL_R_D('pickup',0,0,'Eta',etaN, myThid) |
230 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
231 |
IF ( nonlinFreeSurf.GE.0) |
IF ( nonlinFreeSurf.GE.0) |
232 |
& CALL MNC_CW_RL_R_D('pickup',0,0,'EtaNH', 1, etaH, myThid) |
& CALL MNC_CW_RL_R_D('pickup',0,0,'EtaNH', etaH, myThid) |
233 |
#endif |
#endif |
234 |
CALL MNC_CW_RL_R_D('pickup',0,0,'Unm1',1,gUnm1, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'Unm1',gUnm1, myThid) |
235 |
CALL MNC_CW_RL_R_D('pickup',0,0,'Vnm1',1,gVnm1, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'Vnm1',gVnm1, myThid) |
236 |
CALL MNC_CW_RL_R_D('pickup',0,0,'Tnm1',1,gTnm1, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'Tnm1',gTnm1, myThid) |
237 |
CALL MNC_CW_RL_R_D('pickup',0,0,'Snm1',1,gSnm1, myThid) |
CALL MNC_CW_RL_R_D('pickup',0,0,'Snm1',gSnm1, myThid) |
238 |
ENDIF |
ENDIF |
239 |
#endif /* ALLOW_MNC */ |
#endif /* ALLOW_MNC */ |
240 |
|
|