/[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.1 by jahn, Wed Apr 13 18:56:25 2011 UTC revision 1.6 by jahn, Thu Aug 23 21:48:25 2012 UTC
# Line 88  c test.................... Line 88  c test....................
88  c set up ecosystem coefficients  c set up ecosystem coefficients
89  c  c
90  c seed randomization  c seed randomization
91            _BEGIN_MASTER( myThid )
92    
93          CALL DARWIN_RANDOM_INIT(darwin_seed, myThid)          CALL DARWIN_RANDOM_INIT(darwin_seed, myThid)
94    
95  c initialize total number of functional groups tried  c initialize total number of functional groups tried
# Line 100  c initialize total number of functional Line 102  c initialize total number of functional
102  #endif  #endif
103           end do           end do
104    
105            _END_MASTER( myThid )
106    
107  c reduce amount of diaz  c reduce amount of diaz
108  #ifdef ALLOW_DIAZ  #ifdef ALLOW_DIAZ
109          do np = 1, npmax          do np = 1, npmax
# Line 120  c reduce amount of diaz Line 124  c reduce amount of diaz
124          enddo          enddo
125  #endif  #endif
126    
127            _BEGIN_MASTER( myThid )
128    
129  c initialize zooplankton  c initialize zooplankton
130           call MONOD_GENERATE_ZOO(MyThid)           call MONOD_GENERATE_ZOO(MyThid)
131    
# Line 163  c find mean Line 169  c find mean
169  #endif  #endif
170  #ifdef DYNAMIC_CHL  #ifdef DYNAMIC_CHL
171  c check Chl fields are reasonable  c check Chl fields are reasonable
172    #ifndef WAVEBANDS
173             do np = 1,npmax
174    c C:CHl minimum: chosen to be Chl:C at high light (2000uEin/m2/s) and
175    c  no temp/nutrient limitation
176              chl2cmin(np)=chl2cmax(np)/
177         &            (1+(chl2cmax(np)*alphachl(np)*2000. _d 0)/
178         &               (2*pcmax(np)))
179              chl2cmin(np)=0. _d 0
180             enddo
181    #else
182             do np=1,npmax             do np=1,npmax
183               chl2cmin(np)=chl2cmax(np)/               chl2cmin(np)=chl2cmax(np)/
184       &            (1+(chl2cmax(np)* alpha_mean(np) *2000. _d 0)/       &            (1+(chl2cmax(np)* alpha_mean(np) *2000. _d 0)/
185       &               (2*pcmax(np)))       &               (2*pcmax(np)))
186             enddo             enddo
            call MONOD_CHECK_CHL(myThid)  
187  #endif  #endif
188  #endif  #endif
189  c ANNA endif  #endif
190            _END_MASTER( myThid )
191    
192    #ifdef DYNAMIC_CHL
193    C this initializes fields...
194            call MONOD_CHECK_CHL(myThid)
195    #endif
196    
197           _BEGIN_MASTER(myThid)           IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN
198  c write out initial phyto characteristics  c write out initial phyto characteristics
199  #ifndef GEIDER  #ifndef GEIDER
200           CALL MDSFINDUNIT( IniUnit1, mythid )           CALL MDSFINDUNIT( IniUnit1, mythid )
# Line 268  c ANNA outputs mQyield as 10^(4) mmol C Line 289  c ANNA outputs mQyield as 10^(4) mmol C
289           close(IniUnit1)           close(IniUnit1)
290   110     format(3f4.0,f6.2,4f4.0,f5.1,4f7.3,2e11.2,2f9.4,f6.1,2i5)   110     format(3f4.0,f6.2,4f4.0,f5.1,4f7.3,2e11.2,2f9.4,f6.1,2i5)
291   111     format(3f4.0,f6.2,4f4.0,f5.1,4f7.3,2f6.1,2f9.4,f6.1,2i5)   111     format(3f4.0,f6.2,4f4.0,f5.1,4f7.3,2f6.1,2f9.4,f6.1,2i5)
292    c       myProcId and myThid
293            ENDIF
294    
295           CALL LEF_ZERO( fice,myThid )           CALL LEF_ZERO( fice,myThid )
296           CALL LEF_ZERO( inputFe,myThid )           CALL LEF_ZERO( inputFe,myThid )
# Line 327  c set arrays to zero if first timestep Line 350  c set arrays to zero if first timestep
350       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
351             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,i),             CALL TIMEAVE_RESET(Euave(1-OLx,1-OLy,1,1,1,i),
352       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
353               CALL TIMEAVE_RESET(Estave(1-OLx,1-OLy,1,1,1,i),
354         &                                    Nr,bi,bj,myThid)
355             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,i),             CALL TIMEAVE_RESET(Eutave(1-OLx,1-OLy,1,1,1,i),
356       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
357            enddo            enddo
358  #endif  #endif
359    #ifdef DAR_DIAG_IRR_AMPS
360              do i=1,tlam
361               CALL TIMEAVE_RESET(amp1ave(1-OLx,1-OLy,1,1,1,i),
362         &                                    Nr,bi,bj,myThid)
363               CALL TIMEAVE_RESET(amp2ave(1-OLx,1-OLy,1,1,1,i),
364         &                                    Nr,bi,bj,myThid)
365              enddo
366    #endif
367  #ifdef DAR_DIAG_ABSORP  #ifdef DAR_DIAG_ABSORP
368            do i=1,tlam            do i=1,tlam
369             CALL TIMEAVE_RESET(aave(1-OLx,1-OLy,1,1,1,i),             CALL TIMEAVE_RESET(aave(1-OLx,1-OLy,1,1,1,i),
# Line 355  c set arrays to zero if first timestep Line 388  c set arrays to zero if first timestep
388       &                                    Nr,bi,bj,myThid)       &                                    Nr,bi,bj,myThid)
389            enddo            enddo
390  #endif  #endif
391    #ifdef DAR_RADTRANS
392              CALL TIMEAVE_RESET(rmudave(1-OLx,1-OLy,1,1),
393         &                                   1,bi,bj,myThid)
394    #endif
395  c ANNA_TAVE  c ANNA_TAVE
396  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
397            do np=1,npmax            do np=1,npmax
# Line 399  c         CALL TIMEAVE_RESET(SURave,   1 Line 436  c         CALL TIMEAVE_RESET(SURave,   1
436  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
437    
438  #ifdef CHECK_CONS  #ifdef CHECK_CONS
439           IF ( myProcId.EQ.0 .AND. myThid.EQ.1 ) THEN
440  coj find unused units for darwin_cons output  coj find unused units for darwin_cons output
441          CALL MDSFINDUNIT( DAR_cons_unit1, mythid )          CALL MDSFINDUNIT( DAR_cons_unit1, mythid )
442          open(DAR_cons_unit1,file='darwin_cons_P.txt',status='unknown')          open(DAR_cons_unit1,file='darwin_cons_P.txt',status='unknown')
# Line 416  coj find unused units for darwin_cons ou Line 454  coj find unused units for darwin_cons ou
454          CALL MDSFINDUNIT( DAR_cons_unit7, mythid )          CALL MDSFINDUNIT( DAR_cons_unit7, mythid )
455          open(DAR_cons_unit7,file='darwin_cons_O.txt',status='unknown')          open(DAR_cons_unit7,file='darwin_cons_O.txt',status='unknown')
456  #endif  #endif
457    c      myProcId and myThid
458           ENDIF
459  #endif  #endif
          _END_MASTER(myThid)  
460    
461  c test....................  c test....................
462  c       write(6,*)'finishing darwin_init_vari '  c       write(6,*)'finishing darwin_init_vari '
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22