pkg/dic/dic_ini_forcing.F

C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_ini_forcing.F,v 1.23 2014/08/30 00:34:02 jmc Exp $
C $Name:  $

#include "DIC_OPTIONS.h"

CBOP
C !ROUTINE: DIC_INI_FORCING

C !INTERFACE: ==========================================================
      SUBROUTINE DIC_INI_FORCING( myThid )

C !DESCRIPTION:
c initializes forcing fields to zero
c (or to reasonable values)

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DIC_VARS.h"
#include "DIC_LOAD.h"

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
      INTEGER myThid
CEOP

#ifdef ALLOW_DIC

c !LOCAL VARIABLES: ===================================================
       INTEGER bi,bj,i,j
#if (defined (READ_PAR) && defined (USE_QSW))
       CHARACTER*(MAX_LEN_MBUF) msgBuf
#endif

C First call requires that we initialize everything to zero for safety
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
         DIC_ldRec(bi,bj) = 0
       ENDDO
      ENDDO
      CALL LEF_ZERO( dicwind0,myThid )
      CALL LEF_ZERO( dicwind1,myThid )
      CALL LEF_ZERO( atmosp0,myThid )
      CALL LEF_ZERO( atmosp1,myThid )
      CALL LEF_ZERO( silica0,myThid )
      CALL LEF_ZERO( silica1,myThid )
      CALL LEF_ZERO( ice0,myThid )
      CALL LEF_ZERO( ice1,myThid )
#ifdef READ_PAR
      CALL LEF_ZERO( par0,myThid )
      CALL LEF_ZERO( par1,myThid )
#endif
#ifdef ALLOW_FE
      CALL LEF_ZERO( feinput0,myThid )
      CALL LEF_ZERO( feinput1,myThid )
#endif
#ifdef LIGHT_CHL
      CALL LEF_ZERO( chlinput,myThid )
#endif

#ifdef READ_PAR
#ifdef USE_QSW
        WRITE(msgBuf,'(2A)')
     &       ' DIC_INI_FORCING: You can not use READ_PAR ',
     &       ' and USE_QSW'
        CALL PRINT_ERROR( msgBuf, myThid )
        STOP 'ABNORMAL END: S/R DIC_INI_FORCING'
#endif
#endif

c set reasonable values to those that need at least something
       DO bj = myByLo(myThid), myByHi(myThid)
        DO bi = myBxLo(myThid), myBxHi(myThid)
         DO j=1-OLy,sNy+OLy
          DO i=1-OLx,sNx+OLx
             WIND(i,j,bi,bj)   = 5. _d 0*maskC(i,j,1,bi,bj)
             AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
             SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
             FIce(i,j,bi,bj)   = 0. _d 0
             FluxCO2(i,j,bi,bj)= 0. _d 0
#ifdef READ_PAR
             PAR(i,j,bi,bj)    = 100. _d 0*maskC(i,j,1,bi,bj)
#endif
#ifdef LIGHT_CHL
C If the chlorophyll climatology is not provided, use this default value.
             CHL(i,j,bi,bj)    = 1. _d -2*maskC(i,j,1,bi,bj)
#endif
#ifdef ALLOW_FE
             InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
#endif
          ENDDO
         ENDDO
        ENDDO
       ENDDO

#endif /* ALLOW_DIC */
      RETURN
      END
1	jmc	1.24	C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_ini_forcing.F,v 1.23 2014/08/30 00:34:02 jmc Exp $
2	jmc	1.3	C $Name: $
3
4	stephd	1.1	#include "DIC_OPTIONS.h"
5
6	stephd	1.2	CBOP
7	dfer	1.16	C !ROUTINE: DIC_INI_FORCING
8	stephd	1.2
9			C !INTERFACE: ==========================================================
10	stephd	1.1	SUBROUTINE DIC_INI_FORCING( myThid )
11
12	stephd	1.2	C !DESCRIPTION:
13			c initializes forcing fields to zero
14	stephd	1.5	c (or to reasonable values)
15	stephd	1.2
16			C !USES: ===============================================================
17	stephd	1.1	IMPLICIT NONE
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20			#include "PARAMS.h"
21			#include "GRID.h"
22	dfer	1.14	#include "DIC_VARS.h"
23	stephd	1.1	#include "DIC_LOAD.h"
24
25	stephd	1.2	C !INPUT PARAMETERS: ===================================================
26			C myThid :: thread number
27	stephd	1.1	INTEGER myThid
28	stephd	1.2	CEOP
29	stephd	1.1
30	jmc	1.22	#ifdef ALLOW_DIC
31	stephd	1.5
32			c !LOCAL VARIABLES: ===================================================
33			INTEGER bi,bj,i,j
34	jmc	1.21	#if (defined (READ_PAR) && defined (USE_QSW))
35	stephd	1.6	CHARACTER*(MAX_LEN_MBUF) msgBuf
36	jmc	1.21	#endif
37
38	stephd	1.1	C First call requires that we initialize everything to zero for safety
39	jmc	1.22	DO bj = myByLo(myThid), myByHi(myThid)
40			DO bi = myBxLo(myThid), myBxHi(myThid)
41			DIC_ldRec(bi,bj) = 0
42			ENDDO
43			ENDDO
44	gforget	1.18	CALL LEF_ZERO( dicwind0,myThid )
45			CALL LEF_ZERO( dicwind1,myThid )
46	stephd	1.1	CALL LEF_ZERO( atmosp0,myThid )
47			CALL LEF_ZERO( atmosp1,myThid )
48			CALL LEF_ZERO( silica0,myThid )
49			CALL LEF_ZERO( silica1,myThid )
50			CALL LEF_ZERO( ice0,myThid )
51			CALL LEF_ZERO( ice1,myThid )
52	stephd	1.4	#ifdef READ_PAR
53			CALL LEF_ZERO( par0,myThid )
54			CALL LEF_ZERO( par1,myThid )
55			#endif
56	jmc	1.24	#ifdef ALLOW_FE
57	stephd	1.1	CALL LEF_ZERO( feinput0,myThid )
58			CALL LEF_ZERO( feinput1,myThid )
59	jmc	1.24	#endif
60			#ifdef LIGHT_CHL
61			CALL LEF_ZERO( chlinput,myThid )
62			#endif
63	stephd	1.1
64	stephd	1.6	#ifdef READ_PAR
65	stephd	1.7	#ifdef USE_QSW
66	stephd	1.6	WRITE(msgBuf,'(2A)')
67			& ' DIC_INI_FORCING: You can not use READ_PAR ',
68	stephd	1.7	& ' and USE_QSW'
69	stephd	1.6	CALL PRINT_ERROR( msgBuf, myThid )
70			STOP 'ABNORMAL END: S/R DIC_INI_FORCING'
71			#endif
72			#endif
73
74	stephd	1.5	c set reasonable values to those that need at least something
75			DO bj = myByLo(myThid), myByHi(myThid)
76			DO bi = myBxLo(myThid), myBxHi(myThid)
77	jmc	1.23	DO j=1-OLy,sNy+OLy
78			DO i=1-OLx,sNx+OLx
79	dfer	1.9	WIND(i,j,bi,bj) = 5. _d 0*maskC(i,j,1,bi,bj)
80			AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
81			SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
82			FIce(i,j,bi,bj) = 0. _d 0
83	heimbach	1.20	FluxCO2(i,j,bi,bj)= 0. _d 0
84	stephd	1.5	#ifdef READ_PAR
85	dfer	1.9	PAR(i,j,bi,bj) = 100. _d 0*maskC(i,j,1,bi,bj)
86	stephd	1.5	#endif
87	jmc	1.24	#ifdef LIGHT_CHL
88			C If the chlorophyll climatology is not provided, use this default value.
89			CHL(i,j,bi,bj) = 1. _d -2*maskC(i,j,1,bi,bj)
90			#endif
91	stephd	1.5	#ifdef ALLOW_FE
92	dfer	1.9	InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
93	stephd	1.5	#endif
94			ENDDO
95			ENDDO
96			ENDDO
97			ENDDO
98
99	jmc	1.22	#endif /* ALLOW_DIC */
100	stephd	1.1	RETURN
101			END