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

Changeset 15284


Ignore:
Timestamp:
2021-09-23T21:00:08+02:00 (3 years ago)
Author:
aumont
Message:

bug fixes in p4zlys and p4zsed

Location:
NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_top_cfg

    r14963 r15284  
    128128   sn_trccbc(14) = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    129129   sn_trccbc(23) = 'river.orca'      ,    120            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    130    rn_trcfac(1)  = 8.333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
    131    rn_trcfac(2)  = 8.333e+01      !  ( 1e3 /12 ) 
    132    rn_trcfac(5)  = 3.774193e+04   !  ( 1e3 / 31. * 117 ) 
     130   rn_trcfac(1)  = 8.333333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
     131   rn_trcfac(2)  = 8.333333e+01      !  ( 1e3 /12 ) 
     132   rn_trcfac(5)  = 3.774193e+03   !  ( 1e3 / 31. * 117 ) 
    133133   rn_trcfac(7)  = 3.558719e+01   !  ( 1e3 / 28.1 ) 
    134134   rn_trcfac(10) = 8.333333e+01   !  ( 1e3 / 12 

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/cfgs/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r14239 r15284  
    1 bld::tool::fppkeys   key_top key_xios 
     1bld::tool::fppkeys   key_top key_xios key_mpi2 

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/cfgs/SHARED/namelist_pisces_ref

    r14972 r15284  
    349349&nampiscal     !   parameters for Calcite chemistry 
    350350!----------------------------------------------------------------------- 
    351    kdca       =  3.       ! calcite dissolution rate constant (1/time) 
    352    nca        =  2.       ! order of dissolution reaction (dimensionless) 
     351   kdca       =  100.     ! calcite dissolution rate constant (1/time) 
     352   nca        =  4.7      ! order of dissolution reaction (dimensionless) 
    353353/ 
    354354!----------------------------------------------------------------------- 

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zlys.F90

    r15127 r15284  
    107107         ELSE 
    108108 
    109             zkd = kdca * 0.1**(nca - 0.4) 
    110             zexcess  = zexcess0**0.4 
     109            zkd = kdca * 0.2**(nca - 0.11) 
     110            zexcess  = zexcess0**0.11 
    111111            zdispot = zkd * zexcess * tr(ji,jj,jk,jpcal,Kbb) 
    112112        ENDIF 

  • NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zsed.F90

    r15127 r15284  
    295295            tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0 
    296296            tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zfact / 3.0 
     297            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zfact * 2.0 / 3.0 
    297298            tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) - zfact * 2.0 / 3.0 
    298299            tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zfact * 1.0 / 3.0 
    … …  
    321322            tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0 
    322323            tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zfact / 3.0 
     324            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zfact * 2.0 / 3.0 
    323325            ! N/P ratio of diazotrophs is supposed to be 46     
    324326            ztrpo4 = tr(ji,jj,jk,jppo4,Kbb) / ( 1E-6 + tr(ji,jj,jk,jppo4,Kbb) ) 
Note: See TracChangeset for help on using the changeset viewer.