/[MITgcm]/MITgcm/pkg/ebm/ebm_atmosphere.F
ViewVC logotype

Diff of /MITgcm/pkg/ebm/ebm_atmosphere.F

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

revision 1.8 by jmc, Tue Aug 24 14:14:30 2010 UTC revision 1.9 by jmc, Sun Aug 28 21:54:40 2011 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "EBM_OPTIONS.h"  #include "EBM_OPTIONS.h"
5    
6        SUBROUTINE EBM_ATMOSPHERE ( myTime, myiter, myThid )  CBOP 0
7    C !ROUTINE: EBM_ATMOSPHERE
8    
9  C     |==========================================================|  C !INTERFACE:
10  C     | S/R CALCULATE FORCING FROM ENERGY AND MOISTURE           |        SUBROUTINE EBM_ATMOSPHERE ( myTime, myIter, myThid )
11  C     | BALANCE ATMOSPHERE                                       |  
12  C     |==========================================================|  C     !DESCRIPTION:
13    C     *==========================================================*
14    C     | S/R CALCULATE FORCING FROM ENERGY AND MOISTURE
15    C     | BALANCE ATMOSPHERE
16    C     *==========================================================*
17  C      References:  C      References:
18  C      * X. Wang, P. Stone and J. Marotzke, 1999:  C      * X. Wang, P. Stone and J. Marotzke, 1999:
19  C        Global thermohaline circulation. Part I:  C        Global thermohaline circulation. Part I:
# Line 23  C        Destabilization of the thermoha Line 28  C        Destabilization of the thermoha
28  C        by atmospheric eddy transports.  C        by atmospheric eddy transports.
29  C        J. Climate 7(12), 1870-1882  C        J. Climate 7(12), 1870-1882
30    
31    C     !USES:
32        IMPLICIT NONE        IMPLICIT NONE
   
33  C     === Global variables ===  C     === Global variables ===
34  #include "SIZE.h"  #include "SIZE.h"
35  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 37  C     === Global variables === Line 42  C     === Global variables ===
42  # include "tamc_keys.h"  # include "tamc_keys.h"
43  #endif  #endif
44    
45    C     !INPUT PARAMETERS:
46  C     === Routine arguments ===  C     === Routine arguments ===
47  C     myThid - Instance number for this innvocation of CALC_FORCING  C     myThid  :: my Thread Id number
       INTEGER myThid  
       INTEGER myIter  
48        _RL myTime        _RL myTime
49  CEndOfInterface        INTEGER myIter
50          INTEGER myThid
51    CEOP
52    
53  #ifdef ALLOW_EBM  #ifdef ALLOW_EBM
54    C     !LOCAL VARIABLES:
 C     == Local variables ==  
55        _RL ReCountX(1-OLy:sNy+OLy,nSy)        _RL ReCountX(1-OLy:sNy+OLy,nSy)
56        INTEGER bi, bj        INTEGER bi, bj
57        INTEGER i, j        INTEGER i, j
58        INTEGER no_so        INTEGER no_so
59    #ifdef ALLOW_AUTODIFF_TAMC
60        INTEGER iebmkey        INTEGER iebmkey
61        LOGICAL TOP_LAYER  #endif /* ALLOW_AUTODIFF_TAMC */
   
 C--   Top layer only  
 cph      TOP_LAYER = k .EQ. 1  
   
 cph      IF ( TOP_LAYER ) THEN  
62    
63        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
64         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 336  c--   Pacific Line 337  c--   Pacific
337             ENDIF             ENDIF
338            ENDIF            ENDIF
339  #endif  #endif
340              EmPmR(i,j,bi,bj) = EmPmR(i,j,bi,bj)
341         &                     - Run(i,j,bi,bj)*scale_runoff
342            EmPmR(i,j,bi,bj) = EmPmR(i,j,bi,bj)*rhoConstFresh            EmPmR(i,j,bi,bj) = EmPmR(i,j,bi,bj)*rhoConstFresh
343           ENDDO           ENDDO
344          ENDDO          ENDDO
# Line 345  c--   Pacific Line 348  c--   Pacific
348        _EXCH_XY_RS(Qnet , myThid )        _EXCH_XY_RS(Qnet , myThid )
349        _EXCH_XY_RS(EmPmR , myThid )        _EXCH_XY_RS(EmPmR , myThid )
350    
   
351  C      CALL PLOT_FIELD_XYRS( Qnet, 'Qnet' , 1, myThid )  C      CALL PLOT_FIELD_XYRS( Qnet, 'Qnet' , 1, myThid )
352  C      CALL PLOT_FIELD_XYRS( EmPmR, 'EmPmR' , 1, myThid )  C      CALL PLOT_FIELD_XYRS( EmPmR, 'EmPmR' , 1, myThid )
353    
 cph  end of IF TOP_LAYER  
 cph      ENDIF  
   
354  #endif /* ALLOW_EBM */  #endif /* ALLOW_EBM */
355    
356        RETURN        RETURN
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22