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

Diff of /MITgcm/pkg/aim_compon_interf/atm_store_evmpr.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 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$    C $Name$  
3    
4  #include "AIM_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE ATM_STORE_EVMPR( bi,bj,        SUBROUTINE ATM_STORE_EVMPR( bi,bj,
8       I                             myCurrentTime,       I                             myTime,
9       I                             myCurrentIter,       I                             myIter,
10       I                             myThid )       I                             myThid )
11  C     /==========================================================\  C     /==========================================================\
12  C     | SUBROUTINE ATM_STORE_EVMPR                               |  C     | SUBROUTINE ATM_STORE_EVMPR                               |
# Line 22  C     == Global variables == Line 22  C     == Global variables ==
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #include "CPL_PARAMS.h"  #include "CPL_PARAMS.h"
25    C     == Global variables (main model)
26    #include "FFIELDS.h"
27  C     == Global variables for coupling interface ==  C     == Global variables for coupling interface ==
28  #include "ATMCPL.h"  #include "ATMCPL.h"
 C     == AIMPHYS specific global data ==  
 #include "com_physvar.h"  
29    
30    
31  C     == Routine arguments ==  C     == Routine arguments ==
32  C     bi,bj         - Tile index  C     bi,bj  - Tile index
33  C     myThid        - Thread number for this instance of the routine  C     myThid - Thread number for this instance of the routine
34  C     myCurrentIter - Current timestep number  C     myIter - Current timestep number
35  C     myCurrentTime - Current model time  C     myTime - Current model time
36        INTEGER bi, bj        INTEGER bi, bj
37        _RL     myCurrentTime        _RL     myTime
38        INTEGER myCurrentIter        INTEGER myIter
39        INTEGER myThid        INTEGER myThid
40  CEndOfInterface  CEndOfInterface
41    
42    #ifdef COMPONENT_MODULE
43  C     == Local variables ==  C     == Local variables ==
44  C     I,J,K,II - Loop counters  C     i,j - Loop counters
45  C     conv_precip  :: conversion factor for precip: from g.m-2.s-1 to m/s        INTEGER i,j
       INTEGER I,J,K,II  
       _RL conv_precip  
46        _RL cplTimeFraction        _RL cplTimeFraction
47    
 C     from g/(m^2.s) to m/s :  
       conv_Precip = 1. _d -3 / rhoConstFresh  
   
48  C     o Accumulate fresh water flux that will be exported to the  C     o Accumulate fresh water flux that will be exported to the
49  C       coupling layer. Fresh-water flux is in units of g/m^2/s.  C       coupling layer. Fresh-water flux is in units of g/m^2/s.
50  C       To convert to meters per second divide by density in  C       To convert to meters per second divide by density in
51  C       g/m^3. Positive flux is into the atmosphere (E-P).  C       g/m^3. Positive flux is into the atmosphere (E-P).
52        cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)  
53  c     DO bj=myByLo(myThid),myByHi(myThid)         cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
54  c      DO bi=myBxLo(myThid),myBxHi(myThid)         EvMPrTime(bi,bj) = EvMPrTime(bi,bj) + cplTimeFraction
55          EvMPrTime(bi,bj) = EvMPrTime(bi,bj) + cplTimeFraction         DO j=1,sNy
56          DO J=1,sNy           DO i=1,sNx
57           DO I=1,sNx            EvMPrFlux(i,j,bi,bj) = EvMPrFlux(i,j,bi,bj)
58            II = (sNx)*(J-1)+I       &                         + EmPmR(i,j,bi,bj)*cplTimeFraction
           EvMPrFlux(I,J,bi,bj) = EvMPrFlux(I,J,bi,bj)  
      &                         +( EVAP(II,2,myThid)  
      &                           -PRECNV(II,myThid)  
      &                           -PRECLS(II,myThid)  
      &                           )*conv_precip*cplTimeFraction  
59           ENDDO           ENDDO
60          ENDDO         ENDDO
61  c      ENDDO  
62  c     ENDDO  #endif /* COMPONENT_MODULE */
63    
64        RETURN        RETURN
65        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22