/[MITgcm]/MITgcm_contrib/ESMF/global_ocean.128x64x15/code/driver_run.F
ViewVC logotype

Diff of /MITgcm_contrib/ESMF/global_ocean.128x64x15/code/driver_run.F

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

revision 1.2 by cnh, Sun Feb 22 02:56:19 2004 UTC revision 1.4 by cnh, Mon Mar 29 16:19:43 2004 UTC
# Line 45  C     == Global variables == Line 45  C     == Global variables ==
45  #include "OCNCPL.h"  #include "OCNCPL.h"
46    
47  C     == Routine arguments ==  C     == Routine arguments ==
48        _RL atm_HeatFlux(   1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_HeatFlux(   1-OLx:sNx+OLx,1-OLy:sNy+OLy)
49        _RL atm_TauX(       1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_TauX(       1-OLx:sNx+OLx,1-OLy:sNy+OLy)
50        _RL atm_TauY(       1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_TauY(       1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51        _RL atm_Qlatent(    1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_Qlatent(    1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52        _RL atm_Qsensible(  1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_Qsensible(  1-OLx:sNx+OLx,1-OLy:sNy+OLy)
53        _RL atm_Qlongwave(  1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_Qlongwave(  1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54        _RL atm_Qshortwave( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_Qshortwave( 1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55        _RL atm_uVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_uVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy)
56        _RL atm_vVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_vVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy)
57        _RL atm_FWFlux(     1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_FWFlux(     1-OLx:sNx+OLx,1-OLy:sNy+OLy)
58        _RL atm_Hatm(       1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL atm_Hatm(       1-OLx:sNx+OLx,1-OLy:sNy+OLy)
59        _RL ocn_SSTocn(     1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL ocn_SSTocn(     1-OLx:sNx+OLx,1-OLy:sNy+OLy)
60        _RL ocn_Hocn(       1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL ocn_Hocn(       1-OLx:sNx+OLx,1-OLy:sNy+OLy)
61  C      myThid :: Thread number for this instance of the routine  C      myThid :: Thread number for this instance of the routine
62        INTEGER myThid        INTEGER myThid
63        INTEGER myCurrentIter        INTEGER myCurrentIter
# Line 74  CEOP Line 74  CEOP
74        CALL COMP_OCN_PUSH_CONTEXT(1)        CALL COMP_OCN_PUSH_CONTEXT(1)
75    
76        myThid = 1        myThid = 1
77        CLOSE(errorMessageUnit)  !cnhdbg      CLOSE(errorMessageUnit)
78        CLOSE(standardMessageUnit)  !cnhdbg      CLOSE(standardMessageUnit)
79  !     myProcessStr  = '0000'  !     myProcessStr  = '0000'
80        WRITE(fNam,'(A,A)') 'STDOUT.', myProcessStr(1:4)  !cnhdbg      WRITE(fNam,'(A,A)') 'STDOUT.', myProcessStr(1:4)
81        OPEN(standardMessageUnit,FILE=fNam,STATUS='old',POSITION='append')  !cnhdbg      OPEN(standardMessageUnit,FILE=fNam,STATUS='old',POSITION='append')
82        WRITE(fNam,'(A,A)') 'STDERR.', myProcessStr(1:4)  !cnhdbg      WRITE(fNam,'(A,A)') 'STDERR.', myProcessStr(1:4)
83        OPEN(errorMessageUnit,FILE=fNam,STATUS='old',POSITION='append')  !cnhdbg      OPEN(errorMessageUnit,FILE=fNam,STATUS='old',POSITION='append')
84        WRITE(msgBuf,*) 'Entering OCN driver_run'  !cnhdbg      WRITE(msgBuf,*) 'Entering OCN driver_run'
85        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  !cnhdbg      CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
86       &                    SQUEEZE_RIGHT , 1)  !cnhdbg     &                    SQUEEZE_RIGHT , 1)
87    
88  C     Extract the import fields  C     Extract the import fields
89        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
90         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
91          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
92           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
93            HeatFlux(   i,j,bi,bj)=atm_HeatFlux(  i,j,bi,bj)            HeatFlux(   i,j,bi,bj)=atm_HeatFlux(  i,j)
94            TauX(       i,j,bi,bj)=atm_tauX(      i,j,bi,bj)            TauX(       i,j,bi,bj)=atm_tauX(      i,j)
95            TauY(       i,j,bi,bj)=atm_tauY(      i,j,bi,bj)            TauY(       i,j,bi,bj)=atm_tauY(      i,j)
96            Qlatent(    i,j,bi,bj)=atm_qLatent(   i,j,bi,bj)            Qlatent(    i,j,bi,bj)=atm_qLatent(   i,j)
97            Qsensible(  i,j,bi,bj)=atm_qSensible( i,j,bi,bj)            Qsensible(  i,j,bi,bj)=atm_qSensible( i,j)
98            Qlongwave(  i,j,bi,bj)=atm_qLongwave( i,j,bi,bj)            Qlongwave(  i,j,bi,bj)=atm_qLongwave( i,j)
99            Qshortwave( i,j,bi,bj)=atm_qShortwave(i,j,bi,bj)            Qshortwave( i,j,bi,bj)=atm_qShortwave(i,j)
100            uVelGround( i,j,bi,bj)=atm_uVelGround(i,j,bi,bj)            uVelGround( i,j,bi,bj)=atm_uVelGround(i,j)
101            vVelGround( i,j,bi,bj)=atm_vVelGround(i,j,bi,bj)            vVelGround( i,j,bi,bj)=atm_vVelGround(i,j)
102            FWFlux(     i,j,bi,bj)=atm_FWFlux(    i,j,bi,bj)            FWFlux(     i,j,bi,bj)=atm_FWFlux(    i,j)
103            Hatm(       i,j,bi,bj)=atm_Hatm(      i,j,bi,bj)            Hatm(       i,j,bi,bj)=atm_Hatm(      i,j)
104           ENDDO           ENDDO
105          ENDDO          ENDDO
106         ENDDO         ENDDO
# Line 108  C     Extract the import fields Line 108  C     Extract the import fields
108    
109    
110  C     myThid = 1  C     myThid = 1
111  C     myCurrentTime = startTime+startStep*deltaTClock        myCurrentTime = startTime
112  C     myCurrentIter = nIter0+startStep        myCurrentIter = nIter0
113  C     myThid        = 1  C     myThid        = 1
114  C     nTimesteps    = stopStep-startStep  C     nTimesteps    = 1
115    
116        CALL MONITOR( myCurrentIter, myCurrentTime, myThid )        CALL MONITOR( myCurrentIter, myCurrentTime, myThid )
117    
# Line 129  C     nTimesteps    = stopStep-startStep Line 129  C     nTimesteps    = stopStep-startStep
129        CALL FORWARD_STEP_SETUP(    iLoop, myCurrentTime,        CALL FORWARD_STEP_SETUP(    iLoop, myCurrentTime,
130       &                            myCurrentIter, myThid )       &                            myCurrentIter, myThid )
131    
132          WRITE(41,*) 'MAXVAL SSTocn2cpl = ', MAXVAL(SSTocn2cpl(1:snx,1:sny,:,:))
133          WRITE(41,*) 'MINVAL SSTocn2cpl = ', MINVAL(SSTocn2cpl(1:snx,1:sny,:,:))
134    
135  C     Fill the export fields  C     Fill the export fields
136        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
137         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
138          DO j=1-OLy,sNy+OLy          DO j=1,sNy
139           DO i=1-OLx,sNx+OLx           DO i=1,sNx
140            ocn_SSTocn(i,j,bi,bj) = SSTocn2cpl(i,j,bi,bj)            ocn_SSTocn(i,j) = SSTocn2cpl(i,j,bi,bj)
141            ocn_Hocn(  i,j,bi,bj) = R_low(     i,j,bi,bj)            ocn_Hocn(  i,j) = R_low(     i,j,bi,bj)
142           ENDDO           ENDDO
143          ENDDO          ENDDO
144         ENDDO         ENDDO
145        ENDDO        ENDDO
146    
147        CLOSE(errorMessageUnit)        WRITE(41,*) 'D_RUN MAXVAL ocn_SSTocn = ', MAXVAL(ocn_SSTocn(1:snx,1:sny))
148        CLOSE(standardMessageUnit)        WRITE(41,*) 'D_RUN MINVAL ocn_SSTocn = ', MINVAL(ocn_SSTocn(1:snx,1:sny))
149    
150    !cnhdbg      CLOSE(errorMessageUnit)
151    !cnhdbg      CLOSE(standardMessageUnit)
152        CALL COMP_OCN_POP_CONTEXT(1)        CALL COMP_OCN_POP_CONTEXT(1)
153    
154        RETURN        RETURN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22