/[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.36 - (hide annotations) (download)
Thu Sep 2 17:39:05 2004 UTC (19 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57b_post, checkpoint56b_post, checkpoint55, checkpoint57, checkpoint56, checkpoint54f_post, checkpoint55i_post, checkpoint55c_post, checkpoint57a_post, checkpoint55g_post, checkpoint55d_post, checkpoint55d_pre, checkpoint55j_post, checkpoint55h_post, checkpoint55b_post, checkpoint55f_post, checkpoint56a_post, checkpoint56c_post, checkpoint57a_pre, checkpoint55a_post, checkpoint55e_post
Changes since 1.35: +7 -3 lines 
fix a bug: bi,bj do loops missing.
1 jmc 1.36 C $Header: /u/gcmpack/MITgcm/model/src/ini_forcing.F,v 1.35 2004/07/18 01:04:23 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     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 adcroft 1.10 ENDDO
87     ENDDO
88     ENDDO
89     ENDDO
90 adcroft 1.15 C
91 adcroft 1.16 _BEGIN_MASTER(myThid)
92     IF ( zonalWindFile .NE. ' ' ) THEN
93     CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
94     ENDIF
95     IF ( meridWindFile .NE. ' ' ) THEN
96     CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
97     ENDIF
98     IF ( surfQFile .NE. ' ' ) THEN
99     CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
100 jmc 1.34 ELSEIF ( surfQnetFile .NE. ' ' ) THEN
101     CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
102 adcroft 1.16 ENDIF
103     IF ( EmPmRfile .NE. ' ' ) THEN
104     CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
105     ENDIF
106     IF ( thetaClimFile .NE. ' ' ) THEN
107     CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
108     ENDIF
109     IF ( saltClimFile .NE. ' ' ) THEN
110     CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
111     ENDIF
112 heimbach 1.18 #ifdef SHORTWAVE_HEATING
113     IF ( surfQswFile .NE. ' ' ) THEN
114     CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
115 jmc 1.34 IF ( surfQFile .NE. ' ' ) THEN
116     C- Qnet is now (after c54) the net Heat Flux (including SW)
117 jmc 1.36 DO bj = 1,nSy
118     DO bi = 1,nSx
119     DO j=1-OLy,sNy+OLy
120     DO i=1-OLx,sNx+OLx
121 jmc 1.34 Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
122 jmc 1.36 ENDDO
123 jmc 1.34 ENDDO
124     ENDDO
125 jmc 1.36 ENDDO
126 jmc 1.34 ENDIF
127 heimbach 1.18 ENDIF
128     #endif
129 mlosch 1.23 #ifdef ATMOSPHERIC_LOADING
130     IF ( pLoadFile .NE. ' ' ) THEN
131     CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )
132     ENDIF
133     #endif
134 adcroft 1.16 _END_MASTER(myThid)
135     C
136 adcroft 1.15 _EXCH_XY_R4(fu , myThid )
137     _EXCH_XY_R4(fv , myThid )
138     _EXCH_XY_R4(Qnet , myThid )
139     _EXCH_XY_R4(EmPmR, myThid )
140     _EXCH_XY_R4(SST , myThid )
141     _EXCH_XY_R4(SSS , myThid )
142 mlosch 1.23 #ifdef SHORTWAVE_HEATING
143 adcroft 1.16 _EXCH_XY_R4(Qsw , myThid )
144 mlosch 1.23 #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 adcroft 1.10
150 adcroft 1.15 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)
152 cnh 1.1
153     RETURN
154     END
  ViewVC Help
Powered by ViewVC 1.1.22