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

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

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

Revision 1.37 - (hide annotations) (download)
Tue Jan 4 15:57:22 2005 UTC (19 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57d_post, checkpoint57f_post, checkpoint57c_post, checkpoint57c_pre, checkpoint57e_post, eckpoint57e_pre, checkpoint57f_pre
Changes since 1.36: +7 -1 lines 
add saltFlux(+File) forcing capability with periodicExternalForcing
1 jmc 1.37 C $Header: /u/gcmpack/MITgcm/model/src/ini_forcing.F,v 1.36 2004/09/02 17:39:05 jmc Exp $
2 cnh 1.22 C $Name: $
3 cnh 1.1
4 edhill 1.26 #include "PACKAGES_CONFIG.h"
5 cnh 1.12 #include "CPP_OPTIONS.h"
6 cnh 1.1
7 cnh 1.22 CBOP
8     C !ROUTINE: INI_FORCING
9     C !INTERFACE:
10 cnh 1.1 SUBROUTINE INI_FORCING( myThid )
11 cnh 1.22
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 adcroft 1.14 IMPLICIT NONE
21 cnh 1.1 C === Global variables ===
22     #include "SIZE.h"
23     #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "GRID.h"
26     #include "FFIELDS.h"
27    
28 cnh 1.22 C !INPUT/OUTPUT PARAMETERS:
29 cnh 1.1 C == Routine arguments ==
30     C myThid - Number of this instance of INI_FORCING
31     INTEGER myThid
32    
33 cnh 1.22 C !LOCAL VARIABLES:
34 cnh 1.1 C == Local variables ==
35     C bi,bj - Loop counters
36 adcroft 1.17 C I,J
37 cnh 1.1 INTEGER bi, bj
38 adcroft 1.17 INTEGER I, J
39 cnh 1.22 CEOP
40 cnh 1.1
41 cnh 1.6 _BARRIER
42    
43 adcroft 1.10 DO bj = myByLo(myThid), myByHi(myThid)
44     DO bi = myBxLo(myThid), myBxHi(myThid)
45     DO j=1-OLy,sNy+OLy
46     DO i=1-OLx,sNx+OLx
47 heimbach 1.18 fu (i,j,bi,bj) = 0. _d 0
48     fv (i,j,bi,bj) = 0. _d 0
49     Qnet (i,j,bi,bj) = 0. _d 0
50     EmPmR (i,j,bi,bj) = 0. _d 0
51 jmc 1.27 saltFlux (i,j,bi,bj) = 0. _d 0
52 heimbach 1.18 SST (i,j,bi,bj) = 0. _d 0
53     SSS (i,j,bi,bj) = 0. _d 0
54     Qsw (i,j,bi,bj) = 0. _d 0
55 mlosch 1.23 #ifdef ATMOSPHERIC_LOADING
56     pload (i,j,bi,bj) = 0. _d 0
57 jmc 1.33 sIceLoad (i,j,bi,bj) = 0. _d 0
58 mlosch 1.23 #endif
59 jmc 1.35 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 edhill 1.30 #ifndef ALLOW_EXF
65 heimbach 1.24 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 jmc 1.37 saltFlux0 (i,j,bi,bj) = 0. _d 0
74     saltFlux1 (i,j,bi,bj) = 0. _d 0
75 heimbach 1.24 SST0 (i,j,bi,bj) = 0. _d 0
76     SST1 (i,j,bi,bj) = 0. _d 0
77     SSS0 (i,j,bi,bj) = 0. _d 0
78     SSS1 (i,j,bi,bj) = 0. _d 0
79     #ifdef SHORTWAVE_HEATING
80     Qsw0 (i,j,bi,bj) = 0. _d 0
81     Qsw1 (i,j,bi,bj) = 0. _d 0
82     #endif
83     #ifdef ATMOSPHERIC_LOADING
84     pload0 (i,j,bi,bj) = 0. _d 0
85     pload1 (i,j,bi,bj) = 0. _d 0
86     #endif
87     #endif
88 adcroft 1.10 ENDDO
89     ENDDO
90     ENDDO
91     ENDDO
92 adcroft 1.15 C
93 adcroft 1.16 _BEGIN_MASTER(myThid)
94     IF ( zonalWindFile .NE. ' ' ) THEN
95     CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
96     ENDIF
97     IF ( meridWindFile .NE. ' ' ) THEN
98     CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
99     ENDIF
100     IF ( surfQFile .NE. ' ' ) THEN
101     CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
102 jmc 1.34 ELSEIF ( surfQnetFile .NE. ' ' ) THEN
103     CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
104 adcroft 1.16 ENDIF
105     IF ( EmPmRfile .NE. ' ' ) THEN
106     CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
107     ENDIF
108 jmc 1.37 IF ( saltFluxFile .NE. ' ' ) THEN
109     CALL READ_FLD_XY_RS( saltFluxFile, ' ', saltFlux, 0, myThid )
110     ENDIF
111 adcroft 1.16 IF ( thetaClimFile .NE. ' ' ) THEN
112     CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
113     ENDIF
114     IF ( saltClimFile .NE. ' ' ) THEN
115     CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
116     ENDIF
117 heimbach 1.18 #ifdef SHORTWAVE_HEATING
118     IF ( surfQswFile .NE. ' ' ) THEN
119     CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
120 jmc 1.34 IF ( surfQFile .NE. ' ' ) THEN
121     C- Qnet is now (after c54) the net Heat Flux (including SW)
122 jmc 1.36 DO bj = 1,nSy
123     DO bi = 1,nSx
124     DO j=1-OLy,sNy+OLy
125     DO i=1-OLx,sNx+OLx
126 jmc 1.34 Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
127 jmc 1.36 ENDDO
128 jmc 1.34 ENDDO
129     ENDDO
130 jmc 1.36 ENDDO
131 jmc 1.34 ENDIF
132 heimbach 1.18 ENDIF
133     #endif
134 mlosch 1.23 #ifdef ATMOSPHERIC_LOADING
135     IF ( pLoadFile .NE. ' ' ) THEN
136     CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )
137     ENDIF
138     #endif
139 adcroft 1.16 _END_MASTER(myThid)
140     C
141 adcroft 1.15 _EXCH_XY_R4(fu , myThid )
142     _EXCH_XY_R4(fv , myThid )
143     _EXCH_XY_R4(Qnet , myThid )
144     _EXCH_XY_R4(EmPmR, myThid )
145 jmc 1.37 _EXCH_XY_R4( saltFlux, myThid )
146 adcroft 1.15 _EXCH_XY_R4(SST , myThid )
147     _EXCH_XY_R4(SSS , myThid )
148 mlosch 1.23 #ifdef SHORTWAVE_HEATING
149 adcroft 1.16 _EXCH_XY_R4(Qsw , myThid )
150 mlosch 1.23 #endif
151     #ifdef ATMOSPHERIC_LOADING
152     _EXCH_XY_R4(pload , myThid )
153     C CALL PLOT_FIELD_XYRS( pload, 'S/R INI_FORCING pload',1,myThid)
154     #endif
155 adcroft 1.10
156 adcroft 1.15 C CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
157     C CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
158 cnh 1.1
159     RETURN
160     END
  ViewVC Help
Powered by ViewVC 1.1.22