/[MITgcm]/MITgcm/model/src/ini_forcing.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_forcing.F

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

revision 1.21 by cnh, Sun Feb 4 14:38:47 2001 UTC revision 1.36 by jmc, Thu Sep 2 17:39:05 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CStartOfInterface  CBOP
8    C     !ROUTINE: INI_FORCING
9    C     !INTERFACE:
10        SUBROUTINE INI_FORCING( myThid )        SUBROUTINE INI_FORCING( myThid )
 C     /==========================================================\  
 C     | SUBROUTINE INI_FORCING                                   |  
 C     | o Set model initial forcing fields.                      |  
 C     \==========================================================/  
       IMPLICIT NONE  
11    
12    C     !DESCRIPTION: \bv
13    C     *==========================================================*
14    C     | SUBROUTINE INI_FORCING                                    
15    C     | o Set model initial forcing fields.                      
16    C     *==========================================================*
17    C     \ev
18    
19    C     !USES:
20          IMPLICIT NONE
21  C     === Global variables ===  C     === Global variables ===
22  #include "SIZE.h"  #include "SIZE.h"
23  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 18  C     === Global variables === Line 25  C     === Global variables ===
25  #include "GRID.h"  #include "GRID.h"
26  #include "FFIELDS.h"  #include "FFIELDS.h"
27    
28    C     !INPUT/OUTPUT PARAMETERS:
29  C     == Routine arguments ==  C     == Routine arguments ==
30  C     myThid -  Number of this instance of INI_FORCING  C     myThid -  Number of this instance of INI_FORCING
31        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
32    
33    C     !LOCAL VARIABLES:
34  C     == Local variables ==  C     == Local variables ==
35  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
36  C     I,J  C     I,J
37        INTEGER bi, bj        INTEGER bi, bj
38        INTEGER  I,  J        INTEGER  I,  J
39    CEOP
40    
41        _BARRIER        _BARRIER
42    
# Line 39  C     I,J Line 48  C     I,J
48            fv              (i,j,bi,bj) = 0. _d 0            fv              (i,j,bi,bj) = 0. _d 0
49            Qnet            (i,j,bi,bj) = 0. _d 0            Qnet            (i,j,bi,bj) = 0. _d 0
50            EmPmR           (i,j,bi,bj) = 0. _d 0            EmPmR           (i,j,bi,bj) = 0. _d 0
51              saltFlux        (i,j,bi,bj) = 0. _d 0
52            SST             (i,j,bi,bj) = 0. _d 0            SST             (i,j,bi,bj) = 0. _d 0
53            SSS             (i,j,bi,bj) = 0. _d 0            SSS             (i,j,bi,bj) = 0. _d 0
54            Qsw             (i,j,bi,bj) = 0. _d 0            Qsw             (i,j,bi,bj) = 0. _d 0
55            surfaceTendencyU(i,j,bi,bj) = 0. _d 0  #ifdef ATMOSPHERIC_LOADING
56            surfaceTendencyV(i,j,bi,bj) = 0. _d 0            pload           (i,j,bi,bj) = 0. _d 0
57            surfaceTendencyT(i,j,bi,bj) = 0. _d 0            sIceLoad        (i,j,bi,bj) = 0. _d 0
58            surfaceTendencyS(i,j,bi,bj) = 0. _d 0  #endif
59              surfaceForcingU(i,j,bi,bj) = 0. _d 0
60              surfaceForcingV(i,j,bi,bj) = 0. _d 0
61              surfaceForcingT(i,j,bi,bj) = 0. _d 0
62              surfaceForcingS(i,j,bi,bj) = 0. _d 0
63              surfaceForcingTice(i,j,bi,bj) = 0. _d 0
64    #ifndef ALLOW_EXF
65              taux0           (i,j,bi,bj) = 0. _d 0
66              taux1           (i,j,bi,bj) = 0. _d 0
67              tauy0           (i,j,bi,bj) = 0. _d 0
68              tauy1           (i,j,bi,bj) = 0. _d 0
69              Qnet0           (i,j,bi,bj) = 0. _d 0
70              Qnet1           (i,j,bi,bj) = 0. _d 0
71              EmPmR0          (i,j,bi,bj) = 0. _d 0
72              EmPmR1          (i,j,bi,bj) = 0. _d 0
73              SST0            (i,j,bi,bj) = 0. _d 0
74              SST1            (i,j,bi,bj) = 0. _d 0
75              SSS0            (i,j,bi,bj) = 0. _d 0
76              SSS1            (i,j,bi,bj) = 0. _d 0
77    #ifdef SHORTWAVE_HEATING          
78              Qsw0            (i,j,bi,bj) = 0. _d 0
79              Qsw1            (i,j,bi,bj) = 0. _d 0
80    #endif
81    #ifdef ATMOSPHERIC_LOADING
82              pload0          (i,j,bi,bj) = 0. _d 0
83              pload1          (i,j,bi,bj) = 0. _d 0
84    #endif
85    #endif
86           ENDDO           ENDDO
87          ENDDO          ENDDO
88         ENDDO         ENDDO
# Line 60  C Line 97  C
97        ENDIF        ENDIF
98        IF ( surfQFile .NE. ' '  ) THEN        IF ( surfQFile .NE. ' '  ) THEN
99         CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )         CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
100          ELSEIF ( surfQnetFile .NE. ' '  ) THEN
101           CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
102        ENDIF        ENDIF
103        IF ( EmPmRfile .NE. ' '  ) THEN        IF ( EmPmRfile .NE. ' '  ) THEN
104         CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )         CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
# Line 73  C Line 112  C
112  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
113        IF ( surfQswFile .NE. ' '  ) THEN        IF ( surfQswFile .NE. ' '  ) THEN
114         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )         CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
115           IF ( surfQFile .NE. ' '  ) THEN
116    C-     Qnet is now (after c54) the net Heat Flux (including SW)
117            DO bj = 1,nSy
118             DO bi = 1,nSx
119              DO j=1-OLy,sNy+OLy
120               DO i=1-OLx,sNx+OLx
121                Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
122               ENDDO
123              ENDDO
124             ENDDO
125            ENDDO
126           ENDIF
127          ENDIF
128    #endif
129    #ifdef ATMOSPHERIC_LOADING
130          IF ( pLoadFile .NE. ' '  ) THEN
131           CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )
132        ENDIF        ENDIF
133  #endif  #endif
134        _END_MASTER(myThid)        _END_MASTER(myThid)
# Line 83  C Line 139  C
139        _EXCH_XY_R4(EmPmR, myThid )        _EXCH_XY_R4(EmPmR, myThid )
140        _EXCH_XY_R4(SST  , myThid )        _EXCH_XY_R4(SST  , myThid )
141        _EXCH_XY_R4(SSS  , myThid )        _EXCH_XY_R4(SSS  , myThid )
142    #ifdef SHORTWAVE_HEATING
143        _EXCH_XY_R4(Qsw  , myThid )        _EXCH_XY_R4(Qsw  , myThid )
144    #endif
145    #ifdef ATMOSPHERIC_LOADING
146          _EXCH_XY_R4(pload  , myThid )
147    C     CALL PLOT_FIELD_XYRS( pload, 'S/R INI_FORCING pload',1,myThid)
148    #endif
149    
150  C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)  C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
151  C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)  C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22