model/src/ini_forcing.F

C $Header: /u/gcmpack/MITgcm/model/src/ini_forcing.F,v 1.33 2004/07/12 20:19:24 jmc Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: INI_FORCING
C     !INTERFACE:
      SUBROUTINE INI_FORCING( myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE INI_FORCING                                    
C     | o Set model initial forcing fields.                       
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "FFIELDS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid -  Number of this instance of INI_FORCING
      INTEGER myThid

C     !LOCAL VARIABLES:
C     == Local variables ==
C     bi,bj  - Loop counters
C     I,J
      INTEGER bi, bj
      INTEGER  I,  J
CEOP

      _BARRIER

      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        DO j=1-OLy,sNy+OLy
         DO i=1-OLx,sNx+OLx
          fu              (i,j,bi,bj) = 0. _d 0
          fv              (i,j,bi,bj) = 0. _d 0
          Qnet            (i,j,bi,bj) = 0. _d 0
          EmPmR           (i,j,bi,bj) = 0. _d 0
          saltFlux        (i,j,bi,bj) = 0. _d 0
          SST             (i,j,bi,bj) = 0. _d 0
          SSS             (i,j,bi,bj) = 0. _d 0
          Qsw             (i,j,bi,bj) = 0. _d 0
#ifdef ATMOSPHERIC_LOADING
          pload           (i,j,bi,bj) = 0. _d 0
          sIceLoad        (i,j,bi,bj) = 0. _d 0
#endif
          surfaceTendencyU(i,j,bi,bj) = 0. _d 0
          surfaceTendencyV(i,j,bi,bj) = 0. _d 0
          surfaceTendencyT(i,j,bi,bj) = 0. _d 0
          surfaceTendencyS(i,j,bi,bj) = 0. _d 0
          surfaceTendencyTice(i,j,bi,bj) = 0. _d 0
#ifndef ALLOW_EXF
          taux0           (i,j,bi,bj) = 0. _d 0
          taux1           (i,j,bi,bj) = 0. _d 0
          tauy0           (i,j,bi,bj) = 0. _d 0
          tauy1           (i,j,bi,bj) = 0. _d 0
          Qnet0           (i,j,bi,bj) = 0. _d 0
          Qnet1           (i,j,bi,bj) = 0. _d 0
          EmPmR0          (i,j,bi,bj) = 0. _d 0
          EmPmR1          (i,j,bi,bj) = 0. _d 0
          SST0            (i,j,bi,bj) = 0. _d 0
          SST1            (i,j,bi,bj) = 0. _d 0
          SSS0            (i,j,bi,bj) = 0. _d 0
          SSS1            (i,j,bi,bj) = 0. _d 0
#ifdef SHORTWAVE_HEATING          
          Qsw0            (i,j,bi,bj) = 0. _d 0
          Qsw1            (i,j,bi,bj) = 0. _d 0
#endif
#ifdef ATMOSPHERIC_LOADING
          pload0          (i,j,bi,bj) = 0. _d 0
          pload1          (i,j,bi,bj) = 0. _d 0
#endif
#endif
         ENDDO
        ENDDO
       ENDDO
      ENDDO
C
      _BEGIN_MASTER(myThid)
      IF ( zonalWindFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
      ENDIF
      IF ( meridWindFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
      ENDIF
      IF ( surfQFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid )
      ELSEIF ( surfQnetFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid )
      ENDIF
      IF ( EmPmRfile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid )
      ENDIF
      IF ( thetaClimFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( thetaClimFile, ' ', SST, 0, myThid )
      ENDIF
      IF ( saltClimFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( saltClimFile, ' ', SSS, 0, myThid )
      ENDIF
#ifdef SHORTWAVE_HEATING
      IF ( surfQswFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid )
       IF ( surfQFile .NE. ' '  ) THEN
C-     Qnet is now (after c54) the net Heat Flux (including SW)
         DO j=1-OLy,sNy+OLy
          DO i=1-OLx,sNx+OLx
            Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
          ENDDO
         ENDDO
       ENDIF
      ENDIF
#endif
#ifdef ATMOSPHERIC_LOADING
      IF ( pLoadFile .NE. ' '  ) THEN
       CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )
      ENDIF
#endif
      _END_MASTER(myThid)
C
      _EXCH_XY_R4(fu   , myThid )
      _EXCH_XY_R4(fv   , myThid )
      _EXCH_XY_R4(Qnet , myThid )
      _EXCH_XY_R4(EmPmR, myThid )
      _EXCH_XY_R4(SST  , myThid )
      _EXCH_XY_R4(SSS  , myThid )
#ifdef SHORTWAVE_HEATING
      _EXCH_XY_R4(Qsw  , myThid )
#endif
#ifdef ATMOSPHERIC_LOADING
      _EXCH_XY_R4(pload  , myThid )
C     CALL PLOT_FIELD_XYRS( pload, 'S/R INI_FORCING pload',1,myThid)
#endif

C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)

      RETURN
      END
1	jmc	1.34	C $Header: /u/gcmpack/MITgcm/model/src/ini_forcing.F,v 1.33 2004/07/12 20:19: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			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	heimbach	1.18	surfaceTendencyU(i,j,bi,bj) = 0. _d 0
60			surfaceTendencyV(i,j,bi,bj) = 0. _d 0
61			surfaceTendencyT(i,j,bi,bj) = 0. _d 0
62			surfaceTendencyS(i,j,bi,bj) = 0. _d 0
63	dimitri	1.31	surfaceTendencyTice(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			DO j=1-OLy,sNy+OLy
118			DO i=1-OLx,sNx+OLx
119			Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj)
120			ENDDO
121			ENDDO
122			ENDIF
123	heimbach	1.18	ENDIF
124			#endif
125	mlosch	1.23	#ifdef ATMOSPHERIC_LOADING
126			IF ( pLoadFile .NE. ' ' ) THEN
127			CALL READ_FLD_XY_RS( pLoadFile, ' ', pload, 0, myThid )
128			ENDIF
129			#endif
130	adcroft	1.16	_END_MASTER(myThid)
131			C
132	adcroft	1.15	_EXCH_XY_R4(fu , myThid )
133			_EXCH_XY_R4(fv , myThid )
134			_EXCH_XY_R4(Qnet , myThid )
135			_EXCH_XY_R4(EmPmR, myThid )
136			_EXCH_XY_R4(SST , myThid )
137			_EXCH_XY_R4(SSS , myThid )
138	mlosch	1.23	#ifdef SHORTWAVE_HEATING
139	adcroft	1.16	_EXCH_XY_R4(Qsw , myThid )
140	mlosch	1.23	#endif
141			#ifdef ATMOSPHERIC_LOADING
142			_EXCH_XY_R4(pload , myThid )
143			C CALL PLOT_FIELD_XYRS( pload, 'S/R INI_FORCING pload',1,myThid)
144			#endif
145	adcroft	1.10
146	adcroft	1.15	C CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
147			C CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
148	cnh	1.1
149			RETURN
150			END