New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2997 – NEMO

Changeset 2997


Ignore:
Timestamp:
2011-10-25T20:14:21+02:00 (12 years ago)
Author:
cetlod
Message:

branch dev_LOCEAN_2011 : minor corrections to ensure PISCES restartability

Location:
branches/2011/dev_LOCEAN_2011/NEMOGCM
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces

    r2987 r2997  
    4646   conc0      =  2.e-6    ! Phosphate half saturation 
    4747   conc1      =  8E-6     ! Phosphate half saturation for diatoms 
    48    conc2      =  1E-9     ! Iron half saturation for phyto 
    49    conc2m     =  3E-9     ! Max iron half saturation for phyto 
    50    conc3      =  2E-9     ! Iron half saturation for diatoms 
    51    conc3m     =  8E-9     ! Maxi iron half saturation for diatoms 
     48   conc2      =  2E-9     ! Iron half saturation for phyto 
     49   conc2m     =  4E-9     ! Max iron half saturation for phyto 
     50   conc3      =  3E-9     ! Iron half saturation for diatoms 
     51   conc3m     =  9E-9     ! Maxi iron half saturation for diatoms 
    5252   xsizedia   =  5.E-7    ! Minimum size criteria for diatoms 
    5353   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
     
    6565&nampisprod     !   parameters for phytoplankton growth 
    6666!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    67    pislope    =  2.       ! P-I slope 
     67   pislope    =  3.       ! P-I slope 
    6868   pislope2   =  2.       ! P-I slope  for diatoms 
    6969   excret     =  0.05     ! excretion ratio of phytoplankton 
     
    9191!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    9292   part2      =  0.75    ! part of calcite not dissolved in mesozoo guts 
    93    grazrat2   =  0.9      ! maximal mesozoo grazing rate 
     93   grazrat2   =  0.7      ! maximal mesozoo grazing rate 
    9494   resrat2    =  0.005    ! exsudation rate of mesozooplankton 
    95    mzrat2     =  0.04     ! mesozooplankton mortality rate 
     95   mzrat2     =  0.03     ! mesozooplankton mortality rate 
    9696   xprefc     =  1.       ! zoo preference for phyto 
    9797   xprefp     =  0.3      ! zoo preference for POC 
     
    104104   xthresh2   =  0.       ! Food threshold for grazing 
    105105   xkgraz2    =  20.E-6   ! half sturation constant for meso grazing 
    106    epsher2    =  0.3      ! Efficicency of Mesozoo growth 
     106   epsher2    =  0.33     ! Efficicency of Mesozoo growth 
    107107   sigma2     =  0.6      ! Fraction of mesozoo excretion as DOM 
    108108   unass2     =  0.3      ! non assimilated fraction of P by mesozoo 
     
    116116   resrat     =  0.03     ! exsudation rate of zooplankton 
    117117   mzrat      =  0.0      ! zooplankton mortality rate 
    118    xpref2c    =  0.2      ! Microzoo preference for POM 
     118   xpref2c    =  0.1      ! Microzoo preference for POM 
    119119   xpref2p    =  1.       ! Microzoo preference for Nanophyto 
    120120   xpref2d    =  0.6      ! Microzoo preference for Diatoms 
     
    124124   xthresh    =  0.       ! Food threshold for feeding 
    125125   xkgraz     =  20.E-6   ! half sturation constant for grazing 
    126    epsher     =  0.3      ! Efficiency of microzoo growth 
     126   epsher     =  0.33     ! Efficiency of microzoo growth 
    127127   sigma1     =  0.6      ! Fraction of microzoo excretion as DOM 
    128128   unass      =  0.3      ! non assimilated fraction of phyto by zoo 
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_pisces

    r2987 r2997  
    4646   conc0      =  2.e-6    ! Phosphate half saturation 
    4747   conc1      =  8E-6     ! Phosphate half saturation for diatoms 
    48    conc2      =  1E-9     ! Iron half saturation for phyto 
    49    conc2m     =  3E-9     ! Max iron half saturation for phyto 
    50    conc3      =  2E-9     ! Iron half saturation for diatoms 
    51    conc3m     =  8E-9     ! Maxi iron half saturation for diatoms 
     48   conc2      =  2E-9     ! Iron half saturation for phyto 
     49   conc2m     =  4E-9     ! Max iron half saturation for phyto 
     50   conc3      =  3E-9     ! Iron half saturation for diatoms 
     51   conc3m     =  9E-9     ! Maxi iron half saturation for diatoms 
    5252   xsizedia   =  5.E-7    ! Minimum size criteria for diatoms 
    5353   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
     
    6565&nampisprod     !   parameters for phytoplankton growth 
    6666!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    67    pislope    =  2.       ! P-I slope 
     67   pislope    =  3.       ! P-I slope 
    6868   pislope2   =  2.       ! P-I slope  for diatoms 
    6969   excret     =  0.05     ! excretion ratio of phytoplankton 
     
    9191!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    9292   part2      =  0.75    ! part of calcite not dissolved in mesozoo guts 
    93    grazrat2   =  0.9      ! maximal mesozoo grazing rate 
     93   grazrat2   =  0.7      ! maximal mesozoo grazing rate 
    9494   resrat2    =  0.005    ! exsudation rate of mesozooplankton 
    95    mzrat2     =  0.04     ! mesozooplankton mortality rate 
     95   mzrat2     =  0.03     ! mesozooplankton mortality rate 
    9696   xprefc     =  1.       ! zoo preference for phyto 
    9797   xprefp     =  0.3      ! zoo preference for POC 
     
    104104   xthresh2   =  0.       ! Food threshold for grazing 
    105105   xkgraz2    =  20.E-6   ! half sturation constant for meso grazing 
    106    epsher2    =  0.3      ! Efficicency of Mesozoo growth 
     106   epsher2    =  0.33     ! Efficicency of Mesozoo growth 
    107107   sigma2     =  0.6      ! Fraction of mesozoo excretion as DOM 
    108108   unass2     =  0.3      ! non assimilated fraction of P by mesozoo 
     
    116116   resrat     =  0.03     ! exsudation rate of zooplankton 
    117117   mzrat      =  0.0      ! zooplankton mortality rate 
    118    xpref2c    =  0.2      ! Microzoo preference for POM 
     118   xpref2c    =  0.1      ! Microzoo preference for POM 
    119119   xpref2p    =  1.       ! Microzoo preference for Nanophyto 
    120120   xpref2d    =  0.6      ! Microzoo preference for Diatoms 
     
    124124   xthresh    =  0.       ! Food threshold for feeding 
    125125   xkgraz     =  20.E-6   ! half sturation constant for grazing 
    126    epsher     =  0.3      ! Efficiency of microzoo growth 
     126   epsher     =  0.33     ! Efficiency of microzoo growth 
    127127   sigma1     =  0.6      ! Fraction of microzoo excretion as DOM 
    128128   unass      =  0.3      ! non assimilated fraction of phyto by zoo 
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zsink.F90

    r2977 r2997  
    457457         DO jj = 1, jpj 
    458458            DO ji = 1,jpi 
    459                zmax  = MAX( heup(ji,jj), hmld(ji,jj) ) 
    460                zfact = MAX( 0., fsdepw(ji,jj,jk+1) - zmax ) / 5000._wp 
     459      !         zmax  = MAX( heup(ji,jj), hmld(ji,jj) ) 
     460      !         zfact = MAX( 0., fsdepw(ji,jj,jk+1) - zmax ) / 5000._wp 
     461               zmax = hmld(ji,jj) 
     462               zfact = MAX( 0., fsdepw(ji,jj,jk+1) - zmax ) / 4000._wp 
    461463               wsbio4(ji,jj,jk) = wsbio2 + ( 200.- wsbio2 ) * zfact 
    462464            END DO 
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/trc.F90

    r2977 r2997  
    2323   !! parameters for the control of passive tracers 
    2424   !! -------------------------------------------------- 
    25    INTEGER, PUBLIC                                                 ::   numnat        !: the number of the passive tracer NAMELIST 
     25   INTEGER, PUBLIC                                                 ::   numnat        !: logicla unit for the passive tracer NAMELIST 
     26   INTEGER, PUBLIC                                                 ::   numstr        !: logical unit for tracer statistics 
    2627 
    2728   !! passive tracers fields (before,now,after) 
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r2977 r2997  
    9494      IF( lk_my_trc  ) THEN   ;   CALL trc_ini_my_trc       ! MY_TRC  tracers 
    9595      ELSE                    ;   IF(lwp) WRITE(numout,*) '          MY_TRC not used' 
     96      ENDIF 
     97 
     98      IF( lwp ) THEN 
     99         ! 
     100         CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
     101         ! 
    96102      ENDIF 
    97103 
     
    153159         WRITE(numout,*) '          *** Total volume of ocean                = ', areatot 
    154160         WRITE(numout,*) '          *** Total inital content of all tracers ' 
     161         WRITE(numout,*) 
    155162         DO jn = 1, jptra 
    156             WRITE(numout,*) ' tracer nb : ', jn, '  name : ', ctrcnm(jn), ' initial content :', trai(jn) 
     163            WRITE(numout,9000) jn, TRIM( ctrcnm(jn) ), trai(jn) 
    157164         ENDDO 
    158165         WRITE(numout,*) 
    159166      ENDIF 
    160  
     167      WRITE(numout,*) 
    161168      IF(ln_ctl) THEN            ! print mean trends (used for debugging) 
    162169         CALL prt_ctl_trc_init 
     
    165172         CALL prt_ctl_trc( tab4d=trn, mask=tmask, clinfo=ctrcnm ) 
    166173      ENDIF 
     1749000  FORMAT(' tracer nb : ',i2,'      name :',a10,'      initial content :',e18.10) 
    167175      ! 
    168176   END SUBROUTINE trc_init 
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r2979 r2997  
    303303         zmean  = ztraf / areatot 
    304304         zdrift = ( ( ztraf - trai(jn) ) / ( trai(jn) + 1.e-12 )  ) * 100._wp 
    305          IF(lwp) WRITE(numout,*) ' tracer nb : ', jn,'   ', TRIM( ctrcnm(jn) ) , & 
    306             &    ' mean = ', zmean, ' min = ', zmin, ' max = ', zmax, ' drift = ', zdrift, ' %' 
     305         IF(lwp) WRITE(numout,9000) jn, TRIM( ctrcnm(jn) ), zmean, zmin, zmax, zdrift 
    307306      END DO 
    308307      WRITE(numout,*)  
    309        
     3089000  FORMAT(' tracer nb :',i2,'    name :',a10,'    mean :',e18.10,'    min :',e18.10, & 
     309      &      '    max :',e18.10,'    drift :',e18.10, ' %') 
     310      ! 
    310311   END SUBROUTINE trc_rst_stat 
    311312 
  • branches/2011/dev_LOCEAN_2011/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

    r2977 r2997  
    4747      !!              Update the passive tracers 
    4848      !!------------------------------------------------------------------- 
    49       INTEGER, INTENT( in ) ::  kt  ! ocean time-step index 
    50       INTEGER               ::  jk  ! 
     49      INTEGER, INTENT( in ) ::  kt      ! ocean time-step index 
     50      INTEGER               ::  jk, jn  ! dummy loop indices 
     51      REAL(wp)              ::  ztrai 
    5152      CHARACTER (len=25)    ::  charout 
    5253      !!------------------------------------------------------------------- 
     
    9091         ! 
    9192      ENDIF 
     93      ! 
     94      ztrai = 0._wp                                                   !  content of all tracers 
     95      DO jn = 1, jptra 
     96         ztrai = ztrai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:)   ) 
     97      END DO 
     98      IF( lwp ) WRITE(numstr,9300) kt,  ztrai / areatot 
     999300  FORMAT(' iter :',i4,'    global content :',e18.10) 
    92100 
    93101   END SUBROUTINE trc_stp 
Note: See TracChangeset for help on using the changeset viewer.