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 3451 – NEMO

Changeset 3451


Ignore:
Timestamp:
2012-08-11T20:08:15+02:00 (12 years ago)
Author:
cetlod
Message:

banch:2012/dev_r3438_LOCEAN15_PISLOB : minor bug corrections, see ticket #972

Location:
branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    r3446 r3451  
    224224      sn_ndepo    = FLD_N( 'ndeposition',   -12     ,  'ndep'     ,  .false.   , .true.  ,   'yearly'  , ''       , ''         ) 
    225225      sn_ironsed  = FLD_N( 'ironsed'    ,   -12     ,  'bathy'    ,  .false.   , .true.  ,   'yearly'  , ''       , ''         ) 
    226       sn_hydrofe  = FLD_N( 'hydrofe'    ,   -12     ,  'epsdb'    ,  .false.   , .true.  ,   'yearly'  , ''       , ''         ) 
     226      sn_hydrofe  = FLD_N( 'hydrofe'    ,   -12     ,  'hydro'    ,  .false.   , .true.  ,   'yearly'  , ''       , ''         ) 
    227227 
    228228      REWIND( numnatp )                     ! read numnatp 
     
    483483         CALL iom_close( numhydro ) 
    484484         ! 
    485          hydrofe(:,:,:) = hydrofe(:,:,:) * hratio / cvol(:,:,:) * ryyss / fse3t(:,:,:) 
     485         hydrofe(:,:,:) = ( hydrofe(:,:,:) * hratio ) / ( cvol(:,:,:) * ryyss + rtrn ) / fse3t(:,:,:) 
     486         ! 
    486487      ENDIF 
    487488      !  
  • branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90

    r3443 r3451  
    4141   INTEGER ::  numalk  !: logical unit for talk budget 
    4242   INTEGER ::  numsil  !: logical unit for Si budget 
     43   INTEGER ::  numfer  !: logical unit for Fer budget 
    4344 
    4445   !!---------------------------------------------------------------------- 
     
    346347      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    347348      !! 
    348       REAL(wp) :: zalkbudget, zno3budget, zsilbudget 
     349      REAL(wp) :: zalkbudget, zno3budget, zsilbudget, zferbudget 
    349350      ! 
    350351      NAMELIST/nampismass/ ln_check_mass 
     
    361362         ENDIF 
    362363 
    363          IF( ln_check_mass .AND. lwp) THEN      !   Open budget file of NO3, ALK, Si 
     364         IF( ln_check_mass .AND. lwp) THEN      !   Open budget file of NO3, ALK, Si, Fer 
    364365            CALL ctl_opn( numno3, 'no3.budget' , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
    365366            CALL ctl_opn( numsil, 'sil.budget' , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
    366367            CALL ctl_opn( numalk, 'talk.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
     368            CALL ctl_opn( numfer, 'iron.budget', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea ) 
    367369         ENDIF 
    368370      ENDIF 
    369371 
    370       IF( ln_check_mass ) THEN      !   Compute the budget of NO3, ALK, Si 
     372      IF( ln_check_mass ) THEN      !   Compute the budget of NO3, ALK, Si, Fer 
    371373         zno3budget = glob_sum( (   trn(:,:,:,jpno3) + trn(:,:,:,jpnh4)  & 
    372374            &                     + trn(:,:,:,jpphy) + trn(:,:,:,jpdia)  & 
     
    381383            &                     + trn(:,:,:,jptal)                     & 
    382384            &                     + trn(:,:,:,jpcal) * 2.                ) * cvol(:,:,:)  ) 
     385         !  
     386         zferbudget = glob_sum( (   trn(:,:,:,jpfer) + trn(:,:,:,jpnfe)  & 
     387            &                     + trn(:,:,:,jpdfe) + trn(:,:,:,jpbfe)  & 
     388            &                     + trn(:,:,:,jpsfe) +                   & 
     389            &                     + trn(:,:,:,jpzoo)                     & 
     390            &                     + trn(:,:,:,jpmes) * ferat3            ) * cvol(:,:,:)  ) 
    383391 
    384392         IF( lwp ) THEN 
     
    386394            WRITE(numsil,9500) kt,  zsilbudget / areatot 
    387395            WRITE(numalk,9500) kt,  zalkbudget / areatot 
     396            WRITE(numfer,9500) kt,  zferbudget / areatot 
    388397         ENDIF 
    389398       ENDIF 
Note: See TracChangeset for help on using the changeset viewer.