/[MITgcm]/MITgcm/pkg/aim_compon_interf/cpl_ini_vars.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_compon_interf/cpl_ini_vars.F

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

revision 1.1 by jmc, Mon Dec 15 02:44:48 2003 UTC revision 1.2 by jmc, Fri May 21 20:00:48 2004 UTC
# Line 23  C     == Global variables == Line 23  C     == Global variables ==
23  #include "EEPARAMS.h"  #include "EEPARAMS.h"
24  #include "PARAMS.h"  #include "PARAMS.h"
25  #include "GRID.h"  #include "GRID.h"
26    #include "CPL_PARAMS.h"
27  #include "ATMCPL.h"  #include "ATMCPL.h"
28    
29  C     == Routine arguments ==  C     == Routine arguments ==
# Line 30  C     myThid - Thread number for this in Line 31  C     myThid - Thread number for this in
31        INTEGER myThid        INTEGER myThid
32  CEndOfInterface  CEndOfInterface
33    
34    #ifdef COMPONENT_MODULE
35  C     == Local variables ==  C     == Local variables ==
36  C     I,J,K,bi,bj - Loop counters  C     I,J,bi,bj - Loop counters
37        INTEGER prec        INTEGER I,J,bi,bj
       CHARACTER*(MAX_LEN_FNAM) fn  
       INTEGER I,J,K,bi,bj  
38    
39        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
40         DO bi=myBxLo(myThid),myByLo(myThid)         DO bi=myBxLo(myThid),myByLo(myThid)
41          DO J=1-Oly,sNy+Oly          DO J=1-Oly,sNy+Oly
42           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
43  C         Export field  C-        Export field
44            HeatFlux(I,J,bi,bj)   = 0.            atmSLPr(   I,J,bi,bj) = 0.
45              HeatFlux(  I,J,bi,bj) = 0.
46              qShortWave(I,J,bi,bj) = 0.
47    c         qSensible( I,J,bi,bj) = 0.
48    c         qLatent(   I,J,bi,bj) = 0.
49    c         qLongWave( I,J,bi,bj) = 0.
50    c         uVelGround(I,J,bi,bj) = 0.
51    c         vVelGround(I,J,bi,bj) = 0.
52              tauX(      I,J,bi,bj) = 0.
53              tauY(      I,J,bi,bj) = 0.
54            EvMPrFlux( I,J,bi,bj) = 0.            EvMPrFlux( I,J,bi,bj) = 0.
55            RunOffFlux(I,J,bi,bj) = 0.            RunOffFlux(I,J,bi,bj) = 0.
56            tauX(    I,J,bi,bj)   = 0.            RunOffEnFx(I,J,bi,bj) = 0.
57            tauY(    I,J,bi,bj)   = 0.            iceSaltFlx(I,J,bi,bj) = 0.
58            uVelGround(I,J,bi,bj) = 0.            seaIceMass(I,J,bi,bj) = 0.
59            vVelGround(I,J,bi,bj) = 0.  C-        Import field
60            Qsensible(I,J,bi,bj)  = 0.  c         ocMxlD (I,J,bi,bj)    = 0.
           Qshortwave(I,J,bi,bj) = 0.  
           Qlongwave(I,J,bi,bj)  = 0.  
           Qlatent(I,J,bi,bj)    = 0.  
 C         Import field  
61            SSTocn (I,J,bi,bj)    = 0.            SSTocn (I,J,bi,bj)    = 0.
62              SSSocn (I,J,bi,bj)    = 0.
63              vSqocn (I,J,bi,bj)    = 0.
64           ENDDO           ENDDO
65          ENDDO          ENDDO
66            atmSLPrTime(bi,bj)    = 0.
67          HeatFluxTime(bi,bj)   = 0.          HeatFluxTime(bi,bj)   = 0.
68          EvMPrTime(bi,bj)      = 0.          qShortwaveTime(bi,bj) = 0.
69          RunOffTime(bi,bj)     = 0.  c       qSensibleTime(bi,bj)  = 0.
70    c       qLatentTime(bi,bj)    = 0.
71    c       qLongwaveTime(bi,bj)  = 0.
72    c       uVelGroundTime(bi,bj) = 0.
73    c       vVelGroundTime(bi,bj) = 0.
74          tauXTime(bi,bj)       = 0.          tauXTime(bi,bj)       = 0.
75          tauYTime(bi,bj)       = 0.          tauYTime(bi,bj)       = 0.
76          uVelGroundTime(bi,bj) = 0.          EvMPrTime(bi,bj)      = 0.
77          vVelGroundTime(bi,bj) = 0.          RunOffTime(bi,bj)     = 0.
78          qSensibleTime(bi,bj)  = 0.          ROEnFxTime(bi,bj)     = 0.
79          qLatentTime(bi,bj)    = 0.          saltFxTime(bi,bj)     = 0.
80          qLongwaveTime(bi,bj)  = 0.          seaIceTime(bi,bj)     = 0.
         qShortwaveTime(bi,bj) = 0.  
81         ENDDO         ENDDO
82        ENDDO        ENDDO
83    
84        IF ( startTime .NE. 0. .OR.  nIter0 .NE. 0 ) THEN        IF ( startTime .NE. 0. .OR. nIter0 .NE. 0 ) THEN
85         prec=precFloat64         CALL ATM_CPL_READ_PICKUP( nIter0, myThid )
        WRITE(fn,'(A,I10.10)') 'pickup_cpl.',nIter0  
        CALL MDSREADFIELD(fn,prec,'RL',1,SSTocn    , 1,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,HeatFlux  , 2,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,tauX      , 3,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,tauY      , 4,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,EvMPrFlux , 5,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,RunOffFlux, 6,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,uVelGround, 7,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,vVelGround, 8,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,qLatent   , 9,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,qSensible ,10,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,qLongwave ,11,myThid)  
        CALL MDSREADFIELD(fn,prec,'RL',1,qShortwave,12,myThid)  
86        ENDIF        ENDIF
87    
88    #endif /* COMPONENT_MODULE */
89    
90        RETURN        RETURN
91        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22