/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_init_vari.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_init_vari.F

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

revision 1.7 by stephd, Tue Oct 23 16:39:32 2012 UTC revision 1.12 by jahn, Fri Jun 21 12:49:47 2013 UTC
# Line 106  c initialize total number of functional Line 106  c initialize total number of functional
106    
107  c reduce amount of diaz  c reduce amount of diaz
108  #ifdef ALLOW_DIAZ  #ifdef ALLOW_DIAZ
109           IF (nIter0.EQ.PTRACERS_Iter0) THEN
110          do np = 1, npmax          do np = 1, npmax
111           if (diazotroph(np) .eq. 1. _d 0) then           if (diazotroph(np) .eq. 1. _d 0) then
112           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
# Line 122  c reduce amount of diaz Line 123  c reduce amount of diaz
123           ENDDO           ENDDO
124           endif           endif
125          enddo          enddo
126           ENDIF
127  #endif  #endif
128    
129          _BEGIN_MASTER( myThid )          _BEGIN_MASTER( myThid )
# Line 189  c  no temp/nutrient limitation Line 191  c  no temp/nutrient limitation
191  #endif  #endif
192          _END_MASTER( myThid )          _END_MASTER( myThid )
193    
194  #ifdef DYNAMIC_CHL  #ifdef GEIDER
195  C this initializes fields...  C this initializes fields...
196          call MONOD_CHECK_CHL(myThid)          call MONOD_CHECK_CHL(myThid)
197  #endif  #endif
# Line 323  c       myProcId and myThid Line 325  c       myProcId and myThid
325          ENDDO          ENDDO
326         ENDDO         ENDDO
327        ENDDO        ENDDO
328    #endif
329    
330        IF ( .NOT. ( startTime .EQ. baseTime .AND.  nIter0 .EQ. 0        IF ( .NOT. ( startTime .EQ. baseTime .AND.  nIter0 .EQ. 0
331       &     .AND. pickupSuff .EQ. ' ') ) THEN       &     .AND. pickupSuff .EQ. ' ') ) THEN
332  COJ should probably initialize from a file when nIter0 .EQ. 0  COJ should probably initialize from a file when nIter0 .EQ. 0
333          CALL DARWIN_READ_PICKUP( nIter0, myThid )          CALL DARWIN_READ_PICKUP( nIter0, myThid )
334        ENDIF        ENDIF
 #endif  
335  c  c
336  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
337  c set arrays to zero if first timestep  c set arrays to zero if first timestep
# Line 462  c         CALL TIMEAVE_RESET(SURave,   1 Line 465  c         CALL TIMEAVE_RESET(SURave,   1
465  #ifdef CHECK_CONS  #ifdef CHECK_CONS
466         IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN         IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN
467  coj find unused units for darwin_cons output  coj find unused units for darwin_cons output
468          CALL MDSFINDUNIT( DAR_cons_unit1, mythid )          CALL MDSFINDUNIT( DAR_cons_unitP, mythid )
469          open(DAR_cons_unit1,file='darwin_cons_P.txt',status='unknown')          open(DAR_cons_unitP,file='darwin_cons_P.txt',status='unknown')
470          CALL MDSFINDUNIT( DAR_cons_unit2, mythid )          CALL MDSFINDUNIT( DAR_cons_unitN, mythid )
471          open(DAR_cons_unit2,file='darwin_cons_N.txt',status='unknown')          open(DAR_cons_unitN,file='darwin_cons_N.txt',status='unknown')
472          CALL MDSFINDUNIT( DAR_cons_unit3, mythid )          CALL MDSFINDUNIT( DAR_cons_unitF, mythid )
473          open(DAR_cons_unit3,file='darwin_cons_Fe.txt',status='unknown')          open(DAR_cons_unitF,file='darwin_cons_Fe.txt',status='unknown')
474          CALL MDSFINDUNIT( DAR_cons_unit4, mythid )          CALL MDSFINDUNIT( DAR_cons_unitS, mythid )
475          open(DAR_cons_unit4,file='darwin_cons_Si.txt',status='unknown')          open(DAR_cons_unitS,file='darwin_cons_Si.txt',status='unknown')
476  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
477          CALL MDSFINDUNIT( DAR_cons_unit5, mythid )          CALL MDSFINDUNIT( DAR_cons_unitC, mythid )
478          open(DAR_cons_unit5,file='darwin_cons_C.txt',status='unknown')          open(DAR_cons_unitC,file='darwin_cons_C.txt',status='unknown')
479          CALL MDSFINDUNIT( DAR_cons_unit6, mythid )          CALL MDSFINDUNIT( DAR_cons_unitA, mythid )
480          open(DAR_cons_unit6,file='darwin_cons_A.txt',status='unknown')          open(DAR_cons_unitA,file='darwin_cons_A.txt',status='unknown')
481          CALL MDSFINDUNIT( DAR_cons_unit7, mythid )          CALL MDSFINDUNIT( DAR_cons_unitO, mythid )
482          open(DAR_cons_unit7,file='darwin_cons_O.txt',status='unknown')          open(DAR_cons_unitO,file='darwin_cons_O.txt',status='unknown')
483  #endif  #endif
484  c      myProcId and myThid  c      myProcId and myThid
485         ENDIF         ENDIF
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22