/[MITgcm]/MITgcm_contrib/bling/pkg/bling_ini_forcing.F
ViewVC logotype

Contents of /MITgcm_contrib/bling/pkg/bling_ini_forcing.F

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


Revision 1.2 - (show annotations) (download)
Sun Feb 28 21:49:24 2016 UTC (9 years, 4 months ago) by mmazloff
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +0 -0 lines
Update to BLING version 2

1 C $Header: $
2 C $Name: $
3
4 #include "BLING_OPTIONS.h"
5
6 CBOP
7 subroutine BLING_INI_FORCING( myThid )
8
9 C ==========================================================
10 C | subroutine bling_ini_forcing
11 C | o Initialize forcing fields
12 C ==========================================================
13
14 implicit none
15
16 C === Global variables ===
17 #include "SIZE.h"
18 #include "EEPARAMS.h"
19 #include "PARAMS.h"
20 #include "GRID.h"
21 #include "BLING_VARS.h"
22 #include "BLING_LOAD.h"
23
24 C !INPUT PARAMETERS: ===================================================
25 C myThid :: thread number
26 INTEGER myThid
27 CEOP
28
29 #ifdef ALLOW_BLING
30
31 c !LOCAL VARIABLES: ===================================================
32 INTEGER bi,bj,i,j
33 #if (defined (READ_PAR) && defined (USE_QSW))
34 CHARACTER*(MAX_LEN_MBUF) msgBuf
35 #endif
36
37 C First call requires that we initialize everything to zero for safety
38 DO bj = myByLo(myThid), myByHi(myThid)
39 DO bi = myBxLo(myThid), myBxHi(myThid)
40 BLING_ldRec(bi,bj) = 0
41 ENDDO
42 ENDDO
43 CALL LEF_ZERO( dicwind0,myThid )
44 CALL LEF_ZERO( dicwind1,myThid )
45 CALL LEF_ZERO( atmosp0,myThid )
46 CALL LEF_ZERO( atmosp1,myThid )
47 CALL LEF_ZERO( silica0,myThid )
48 CALL LEF_ZERO( silica1,myThid )
49 CALL LEF_ZERO( ice0,myThid )
50 CALL LEF_ZERO( ice1,myThid )
51 #ifdef READ_PAR
52 CALL LEF_ZERO( par0,myThid )
53 CALL LEF_ZERO( par1,myThid )
54 #endif
55 CALL LEF_ZERO( feinput0,myThid )
56 CALL LEF_ZERO( feinput1,myThid )
57
58 #ifdef READ_PAR
59 #ifdef USE_QSW
60 WRITE(msgBuf,'(2A)')
61 & ' BLING_INI_FORCING: You can not use READ_PAR ',
62 & ' and USE_QSW'
63 CALL PRINT_ERROR( msgBuf, myThid )
64 STOP 'ABNORMAL END: S/R BLING_INI_FORCING'
65 #endif
66 #endif
67
68 c set reasonable values to those that need at least something
69 DO bj = myByLo(myThid), myByHi(myThid)
70 DO bi = myBxLo(myThid), myBxHi(myThid)
71 DO j=1-Oly,sNy+Oly
72 DO i=1-Olx,sNx+Olx
73 WIND(i,j,bi,bj) = 5. _d 0*maskC(i,j,1,bi,bj)
74 AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
75 SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
76 FIce(i,j,bi,bj) = 0. _d 0
77 FluxCO2(i,j,bi,bj)= 0. _d 0
78 InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
79 ENDDO
80 ENDDO
81 ENDDO
82 ENDDO
83
84 c total_atmos_carbon = 0. _d 0
85
86 #endif /* ALLOW_BLING */
87 RETURN
88 END

  ViewVC Help
Powered by ViewVC 1.1.22