/[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.48 - (hide 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 dimitri 1.48 C $Header: /u/gcmpack/MITgcm/model/src/ini_forcing.F,v 1.47 2007/10/15 15:28:24 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 jmc 1.42 C | SUBROUTINE INI_FORCING
15     C | o Set model initial forcing fields.
16 cnh 1.22 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 jmc 1.42 C i, j
37 cnh 1.1 INTEGER bi, bj
38 jmc 1.42 INTEGER i, j
39 cnh 1.22 CEOP
40 cnh 1.1
41 jmc 1.42 C- Initialise all arrays in common blocks
42 adcroft 1.10 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 heimbach 1.18 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 jmc 1.27 saltFlux (i,j,bi,bj) = 0. _d 0
51 heimbach 1.18 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 jmc 1.44 pLoad (i,j,bi,bj) = 0. _d 0
55 jmc 1.33 sIceLoad (i,j,bi,bj) = 0. _d 0
56 jmc 1.35 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 jmc 1.44 #ifndef EXCLUDE_FFIELDS_LOAD
62 heimbach 1.24 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 jmc 1.37 saltFlux0 (i,j,bi,bj) = 0. _d 0
71     saltFlux1 (i,j,bi,bj) = 0. _d 0
72 heimbach 1.24 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 jmc 1.42 #ifdef SHORTWAVE_HEATING
77 heimbach 1.24 Qsw0 (i,j,bi,bj) = 0. _d 0
78     Qsw1 (i,j,bi,bj) = 0. _d 0
79     #endif
80     #ifdef ATMOSPHERIC_LOADING
81 jmc 1.44 pLoad0 (i,j,bi,bj) = 0. _d 0
82     pLoad1 (i,j,bi,bj) = 0. _d 0
83 heimbach 1.24 #endif
84 jmc 1.44 #endif /* EXCLUDE_FFIELDS_LOAD */
85 adcroft 1.10 ENDDO
86     ENDDO
87     ENDDO
88     ENDDO
89 jmc 1.42
90 heimbach 1.38 DO bj = myByLo(myThid), myByHi(myThid)
91     DO bi = myBxLo(myThid), myBxHi(myThid)
92 jmc 1.42 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 heimbach 1.38 ELSE
98 jmc 1.42 lambdaThetaClimRelax(i,j,bi,bj) = 0. _d 0
99 heimbach 1.38 ENDIF
100     IF ( doSaltClimRelax .AND.
101 jmc 1.42 & ABS(yC(i,j,bi,bj)).LE.latBandClimRelax ) THEN
102     lambdaSaltClimRelax(i,j,bi,bj) = 1. _d 0/tauSaltClimRelax
103 heimbach 1.38 ELSE
104 jmc 1.42 lambdaSaltClimRelax(i,j,bi,bj) = 0. _d 0
105 heimbach 1.38 ENDIF
106     ENDDO
107     ENDDO
108     ENDDO
109     ENDDO
110 jmc 1.42
111     C- every-one waits before master thread loads from file
112 cnh 1.39 _BARRIER
113 jmc 1.42
114 adcroft 1.16 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 jmc 1.34 ELSEIF ( surfQnetFile .NE. ' ' ) THEN
123     CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
124 adcroft 1.16 ENDIF
125     IF ( EmPmRfile .NE. ' ' ) THEN
126     CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
127 jmc 1.47 c IF ( convertEmP2rUnit.EQ.mass2rUnit ) THEN
128 jmc 1.46 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 jmc 1.47 c ENDIF
140 adcroft 1.16 ENDIF
141 jmc 1.37 IF ( saltFluxFile .NE. ' ' ) THEN
142     CALL READ_FLD_XY_RS( saltFluxFile, ' ', saltFlux, 0, myThid )
143     ENDIF
144 adcroft 1.16 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 heimbach 1.38 IF ( lambdaThetaFile .NE. ' ' ) THEN
151 jmc 1.42 CALL READ_FLD_XY_RS( lambdaThetaFile, ' ',
152 heimbach 1.38 & lambdaThetaClimRelax, 0, myThid )
153     ENDIF
154     IF ( lambdaSaltFile .NE. ' ' ) THEN
155 jmc 1.42 CALL READ_FLD_XY_RS( lambdaSaltFile, ' ',
156 heimbach 1.38 & lambdaSaltClimRelax, 0, myThid )
157     ENDIF
158 heimbach 1.18 #ifdef SHORTWAVE_HEATING
159     IF ( surfQswFile .NE. ' ' ) THEN
160     CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
161 jmc 1.34 IF ( surfQFile .NE. ' ' ) THEN
162     C- Qnet is now (after c54) the net Heat Flux (including SW)
163 jmc 1.42 _BARRIER
164     DO bj = myByLo(myThid), myByHi(myThid)
165     DO bi = myBxLo(myThid), myBxHi(myThid)
166 jmc 1.36 DO j=1-OLy,sNy+OLy
167     DO i=1-OLx,sNx+OLx
168 jmc 1.34 Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
169 jmc 1.36 ENDDO
170 jmc 1.34 ENDDO
171     ENDDO
172 jmc 1.36 ENDDO
173 jmc 1.34 ENDIF
174 heimbach 1.18 ENDIF
175     #endif
176 mlosch 1.23 #ifdef ATMOSPHERIC_LOADING
177     IF ( pLoadFile .NE. ' ' ) THEN
178 jmc 1.44 CALL READ_FLD_XY_RS( pLoadFile, ' ', pLoad, 0, myThid )
179 mlosch 1.23 ENDIF
180     #endif
181 heimbach 1.40
182 jmc 1.43 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 mlosch 1.23 #ifdef SHORTWAVE_HEATING
191 jmc 1.43 CALL EXCH_XY_RS(Qsw , myThid )
192 mlosch 1.23 #endif
193     #ifdef ATMOSPHERIC_LOADING
194 jmc 1.44 CALL EXCH_XY_RS(pLoad , myThid )
195     C CALL PLOT_FIELD_XYRS( pLoad, 'S/R INI_FORCING pLoad',1,myThid)
196 mlosch 1.23 #endif
197 adcroft 1.15 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 cnh 1.1
200     RETURN
201     END
  ViewVC Help
Powered by ViewVC 1.1.22