1 |
jmc |
1.1 |
C $Header: /u/gcmpack/MITgcm/pkg/ocn_compon_interf/cpl_ini_vars.F,v 1.1 2003/12/15 02:49:09 jmc Exp $ |
2 |
|
|
C $Name: $ |
3 |
|
|
|
4 |
|
|
#include "CPP_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
CStartOfInterface |
7 |
|
|
SUBROUTINE OCN_CPL_READ_PICKUP( myIter, myThid ) |
8 |
|
|
C *==========================================================* |
9 |
|
|
C | SUBROUTINE OCN_CPL_READ_PICKUP |
10 |
|
|
C | o Reads fields from a pickup file for a restart |
11 |
|
|
C *==========================================================* |
12 |
|
|
C *==========================================================* |
13 |
|
|
IMPLICIT NONE |
14 |
|
|
|
15 |
|
|
C == Global variables == |
16 |
|
|
#include "SIZE.h" |
17 |
|
|
#include "EEPARAMS.h" |
18 |
|
|
#include "PARAMS.h" |
19 |
|
|
#include "CPL_PARAMS.h" |
20 |
|
|
#include "OCNCPL.h" |
21 |
|
|
|
22 |
|
|
C == Routine arguments == |
23 |
|
|
C myIter :: time-step number |
24 |
|
|
C myThid :: Thread number for this instance of the routine |
25 |
|
|
INTEGER myIter |
26 |
|
|
INTEGER myThid |
27 |
|
|
CEndOfInterface |
28 |
|
|
|
29 |
|
|
#ifdef COMPONENT_MODULE |
30 |
|
|
C == Local variables == |
31 |
|
|
INTEGER prec |
32 |
|
|
CHARACTER*(MAX_LEN_FNAM) fn |
33 |
|
|
|
34 |
|
|
WRITE(fn,'(A,I10.10)') 'pickup_cpl.',myIter |
35 |
|
|
prec = precFloat64 |
36 |
|
|
|
37 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,atmSLPr , 1,myThid) |
38 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,HeatFlux , 2,myThid) |
39 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,qShortWave, 3,myThid) |
40 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,tauX , 4,myThid) |
41 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,tauY , 5,myThid) |
42 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,FWFlux , 6,myThid) |
43 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,iceSaltFlx, 7,myThid) |
44 |
|
|
CALL MDSREADFIELD(fn,prec,'RL',1,seaIceMass, 8,myThid) |
45 |
|
|
|
46 |
|
|
c CALL MDSREADFIELD(fn,prec,'RL',1,ocMxlD2cpl, 9,myThid) |
47 |
|
|
c CALL MDSREADFIELD(fn,prec,'RL',1,SSTocn2cpl,10,myThid) |
48 |
|
|
c CALL MDSREADFIELD(fn,prec,'RL',1,SSSocn2cpl,11,myThid) |
49 |
|
|
c CALL MDSREADFIELD(fn,prec,'RL',1,vSqocn2cpl,12,myThid) |
50 |
|
|
|
51 |
|
|
C-- Fill in the halo region: |
52 |
|
|
_EXCH_XY_R8( atmSLPr , myThid ) |
53 |
|
|
_EXCH_XY_R8( HeatFlux , myThid ) |
54 |
|
|
_EXCH_XY_R8( qShortWave, myThid ) |
55 |
|
|
CALL EXCH_UV_XY_RL(tauX, tauY, .TRUE., myThid ) |
56 |
|
|
_EXCH_XY_R8( FWFlux , myThid ) |
57 |
|
|
_EXCH_XY_R8( iceSaltFlx, myThid ) |
58 |
|
|
_EXCH_XY_R8( seaIceMass, myThid ) |
59 |
|
|
|
60 |
|
|
#endif /* COMPONENT_MODULE */ |
61 |
|
|
|
62 |
|
|
RETURN |
63 |
|
|
END |