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

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

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

Revision 1.48 - (show annotations) (download)
Wed Nov 28 00:18:17 2007 UTC (16 years, 6 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint59q, checkpoint59p, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59k
Changes since 1.47: +1 -7 lines 
added pkg/salt_plume by gathering code, which were previously
spread around various files in model/inc and model/src
results remain numerically identical to before, as a first step
towards adding more options, etc., to this package
1 C $Header: /u/gcmpack/MITgcm/model/src/ini_forcing.F,v 1.47 2007/10/15 15:28:24 jmc Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: INI_FORCING
9 C !INTERFACE:
10 SUBROUTINE INI_FORCING( myThid )
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 ===
22 #include "SIZE.h"
23 #include "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "GRID.h"
26 #include "FFIELDS.h"
27
28 C !INPUT/OUTPUT PARAMETERS:
29 C == Routine arguments ==
30 C myThid - Number of this instance of INI_FORCING
31 INTEGER myThid
32
33 C !LOCAL VARIABLES:
34 C == Local variables ==
35 C bi,bj - Loop counters
36 C i, j
37 INTEGER bi, bj
38 INTEGER i, j
39 CEOP
40
41 C- Initialise all arrays in common blocks
42 DO bj = myByLo(myThid), myByHi(myThid)
43 DO bi = myBxLo(myThid), myBxHi(myThid)
44 DO j=1-OLy,sNy+OLy
45 DO i=1-OLx,sNx+OLx
46 fu (i,j,bi,bj) = 0. _d 0
47 fv (i,j,bi,bj) = 0. _d 0
48 Qnet (i,j,bi,bj) = 0. _d 0
49 EmPmR (i,j,bi,bj) = 0. _d 0
50 saltFlux (i,j,bi,bj) = 0. _d 0
51 SST (i,j,bi,bj) = 0. _d 0
52 SSS (i,j,bi,bj) = 0. _d 0
53 Qsw (i,j,bi,bj) = 0. _d 0
54 pLoad (i,j,bi,bj) = 0. _d 0
55 sIceLoad (i,j,bi,bj) = 0. _d 0
56 surfaceForcingU(i,j,bi,bj) = 0. _d 0
57 surfaceForcingV(i,j,bi,bj) = 0. _d 0
58 surfaceForcingT(i,j,bi,bj) = 0. _d 0
59 surfaceForcingS(i,j,bi,bj) = 0. _d 0
60 surfaceForcingTice(i,j,bi,bj) = 0. _d 0
61 #ifndef EXCLUDE_FFIELDS_LOAD
62 taux0 (i,j,bi,bj) = 0. _d 0
63 taux1 (i,j,bi,bj) = 0. _d 0
64 tauy0 (i,j,bi,bj) = 0. _d 0
65 tauy1 (i,j,bi,bj) = 0. _d 0
66 Qnet0 (i,j,bi,bj) = 0. _d 0
67 Qnet1 (i,j,bi,bj) = 0. _d 0
68 EmPmR0 (i,j,bi,bj) = 0. _d 0
69 EmPmR1 (i,j,bi,bj) = 0. _d 0
70 saltFlux0 (i,j,bi,bj) = 0. _d 0
71 saltFlux1 (i,j,bi,bj) = 0. _d 0
72 SST0 (i,j,bi,bj) = 0. _d 0
73 SST1 (i,j,bi,bj) = 0. _d 0
74 SSS0 (i,j,bi,bj) = 0. _d 0
75 SSS1 (i,j,bi,bj) = 0. _d 0
76 #ifdef SHORTWAVE_HEATING
77 Qsw0 (i,j,bi,bj) = 0. _d 0
78 Qsw1 (i,j,bi,bj) = 0. _d 0
79 #endif
80 #ifdef ATMOSPHERIC_LOADING
81 pLoad0 (i,j,bi,bj) = 0. _d 0
82 pLoad1 (i,j,bi,bj) = 0. _d 0
83 #endif
84 #endif /* EXCLUDE_FFIELDS_LOAD */
85 ENDDO
86 ENDDO
87 ENDDO
88 ENDDO
89
90 DO bj = myByLo(myThid), myByHi(myThid)
91 DO bi = myBxLo(myThid), myBxHi(myThid)
92 DO j=1-Oly,sNy+Oly
93 DO i=1-Olx,sNx+Olx
94 IF ( doThetaClimRelax .AND.
95 & ABS(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN
96 lambdaThetaClimRelax(i,j,bi,bj) = 1. _d 0/tauThetaClimRelax
97 ELSE
98 lambdaThetaClimRelax(i,j,bi,bj) = 0. _d 0
99 ENDIF
100 IF ( doSaltClimRelax .AND.
101 & ABS(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN
102 lambdaSaltClimRelax(i,j,bi,bj) = 1. _d 0/tauSaltClimRelax
103 ELSE
104 lambdaSaltClimRelax(i,j,bi,bj) = 0. _d 0
105 ENDIF
106 ENDDO
107 ENDDO
108 ENDDO
109 ENDDO
110
111 C- every-one waits before master thread loads from file
112 _BARRIER
113
114 IF ( zonalWindFile .NE. ' ' ) THEN
115 CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
116 ENDIF
117 IF ( meridWindFile .NE. ' ' ) THEN
118 CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
119 ENDIF
120 IF ( surfQFile .NE. ' ' ) THEN
121 CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
122 ELSEIF ( surfQnetFile .NE. ' ' ) THEN
123 CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
124 ENDIF
125 IF ( EmPmRfile .NE. ' ' ) THEN
126 CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
127 c IF ( convertEmP2rUnit.EQ.mass2rUnit ) THEN
128 C- EmPmR is now (after c59h) expressed in kg/m2/s (fresh water mass flux)
129 _BARRIER
130 DO bj = myByLo(myThid), myByHi(myThid)
131 DO bi = myBxLo(myThid), myBxHi(myThid)
132 DO j=1-Oly,sNy+Oly
133 DO i=1-Olx,sNx+Olx
134 EmPmR(i,j,bi,bj) = EmPmR(i,j,bi,bj)*rhoConstFresh
135 ENDDO
136 ENDDO
137 ENDDO
138 ENDDO
139 c ENDIF
140 ENDIF
141 IF ( saltFluxFile .NE. ' ' ) THEN
142 CALL READ_FLD_XY_RS( saltFluxFile, ' ', saltFlux, 0, myThid )
143 ENDIF
144 IF ( thetaClimFile .NE. ' ' ) THEN
145 CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
146 ENDIF
147 IF ( saltClimFile .NE. ' ' ) THEN
148 CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
149 ENDIF
150 IF ( lambdaThetaFile .NE. ' ' ) THEN
151 CALL READ_FLD_XY_RS( lambdaThetaFile, ' ',
152 & lambdaThetaClimRelax, 0, myThid )
153 ENDIF
154 IF ( lambdaSaltFile .NE. ' ' ) THEN
155 CALL READ_FLD_XY_RS( lambdaSaltFile, ' ',
156 & lambdaSaltClimRelax, 0, myThid )
157 ENDIF
158 #ifdef SHORTWAVE_HEATING
159 IF ( surfQswFile .NE. ' ' ) THEN
160 CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
161 IF ( surfQFile .NE. ' ' ) THEN
162 C- Qnet is now (after c54) the net Heat Flux (including SW)
163 _BARRIER
164 DO bj = myByLo(myThid), myByHi(myThid)
165 DO bi = myBxLo(myThid), myBxHi(myThid)
166 DO j=1-OLy,sNy+OLy
167 DO i=1-OLx,sNx+OLx
168 Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
169 ENDDO
170 ENDDO
171 ENDDO
172 ENDDO
173 ENDIF
174 ENDIF
175 #endif
176 #ifdef ATMOSPHERIC_LOADING
177 IF ( pLoadFile .NE. ' ' ) THEN
178 CALL READ_FLD_XY_RS( pLoadFile, ' ', pLoad, 0, myThid )
179 ENDIF
180 #endif
181
182 CALL EXCH_UV_XY_RS( fu,fv, .TRUE., myThid )
183 CALL EXCH_XY_RS( Qnet , myThid )
184 CALL EXCH_XY_RS( EmPmR, myThid )
185 CALL EXCH_XY_RS( saltFlux, myThid )
186 CALL EXCH_XY_RS( SST , myThid )
187 CALL EXCH_XY_RS( SSS , myThid )
188 CALL EXCH_XY_RS( lambdaThetaClimRelax, myThid )
189 CALL EXCH_XY_RS( lambdaSaltClimRelax , myThid )
190 #ifdef SHORTWAVE_HEATING
191 CALL EXCH_XY_RS(Qsw , myThid )
192 #endif
193 #ifdef ATMOSPHERIC_LOADING
194 CALL EXCH_XY_RS(pLoad , myThid )
195 C CALL PLOT_FIELD_XYRS( pLoad, 'S/R INI_FORCING pLoad',1,myThid)
196 #endif
197 C CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
198 C CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
199
200 RETURN
201 END
  ViewVC Help
Powered by ViewVC 1.1.22