/[MITgcm]/MITgcm/pkg/bling/bling_init_varia.F
ViewVC logotype

Diff of /MITgcm/pkg/bling/bling_init_varia.F

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

revision 1.2 by mmazloff, Mon Sep 12 20:00:28 2016 UTC revision 1.3 by mmazloff, Thu Oct 27 17:55:31 2016 UTC
# Line 18  C     === Global variables === Line 18  C     === Global variables ===
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19  #include "PARAMS.h"  #include "PARAMS.h"
20  #include "GRID.h"  #include "GRID.h"
21    #include "PTRACERS_SIZE.h"
22    #include "PTRACERS_PARAMS.h"
23  #include "BLING_VARS.h"  #include "BLING_VARS.h"
24  #ifdef ALLOW_BLING_COST  #ifdef ALLOW_BLING_COST
25  # include "BLING_COST.h"  # include "BLING_COST.h"
# Line 53  C--   Initialise global variables Line 55  C--   Initialise global variables
55         ENDDO         ENDDO
56        ENDDO        ENDDO
57    
58    C     Now read phyto plankton abundance ICs and exchange
59          IF (nIter0.EQ.PTRACERS_Iter0) THEN
60            IF ( bling_psmFile .NE. ' ' ) THEN
61             CALL READ_FLD_XYZ_RL(bling_psmFile,' ',
62         &        phyto_sm(1-OLx,1-OLy,1,1,1),0,myThid)
63             _EXCH_XYZ_RL(phyto_sm(1-OLx,1-OLy,1,1,1),myThid)
64            ENDIF
65            IF (bling_plgFile .NE. ' ' ) THEN
66             CALL READ_FLD_XYZ_RL(bling_plgFile,' ',
67         &        phyto_lg(1-OLx,1-OLy,1,1,1),0,myThid)
68             _EXCH_XYZ_RL(phyto_lg(1-OLx,1-OLy,1,1,1),myThid)
69            ENDIF
70            IF ( bling_PdiazFile .NE. ' ' ) THEN
71             CALL READ_FLD_XYZ_RL(bling_pdiazFile,' ',
72         &        phyto_diaz(1-OLx,1-OLy,1,1,1),0,myThid)
73             _EXCH_XYZ_RL(phyto_diaz(1-OLx,1-OLy,1,1,1),myThid)
74            ENDIF
75          ENDIF
76    C     Apply mask
77           DO bj = myByLo(myThid), myByHi(myThid)
78            DO bi = myBxLo(myThid), myBxHi(myThid)
79             DO k=1,Nr
80              DO j=1-OLy,sNy+OLy
81               DO i=1-OLx,sNx+OLx
82                IF (maskC(i,j,k,bi,bj).EQ.0.) THEN
83                 phyto_sm(i,j,k,bi,bj)   = 0. _d 0
84                 phyto_lg(i,j,k,bi,bj)   = 0. _d 0
85                 phyto_diaz(i,j,k,bi,bj) = 0. _d 0
86                ENDIF
87               ENDDO
88              ENDDO
89             ENDDO
90            ENDDO
91           ENDDO
92    
93    C      IF ( nIter0.GT.PTRACERS_Iter0 .OR.
94    C     &    (nIter0.EQ.PTRACERS_Iter0 .AND. pickupSuff.NE.' ')
95    C IF ABOVE TRUE phyto READ FROM PICKUP AND OVERWRITES IC READ
96    C PICKUP CALLED FROM bling_carbonate_init.F
97    C
98    
99  # ifdef ALLOW_BLING_COST  # ifdef ALLOW_BLING_COST
100        totcost =  0. _d 0        totcost =  0. _d 0
101  # endif  # endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22