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

Changeset 15287


Ignore:
Timestamp:
2021-09-24T13:11:02+02:00 (15 months ago)
Author:
cetlod
Message:

trunk : bugfixes in PISCES

Location:
NEMO/trunk
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_top_cfg

    r14608 r15287  
    114114   sn_trcsbc(14) = 'dust.orca.new'   ,       -1          , 'dustfer'     ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    115115   sn_trcsbc(23) = 'ndeposition.orca',      -12          , 'ndep'        ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    116    rn_trsfac(5)  = 8.264e-02   !  (  0.021 / 31. * 122 ) 
    117    rn_trsfac(7)  = 3.313e-01     !  ( 8.8   / 28.1 ) 
    118    rn_trsfac(14) = 6.266e-04   !  (  0.035 / 55.85 ) 
    119    rn_trsfac(23) =  5.4464e-01  !  ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.625/14 ) 
    120    rn_sbc_time   =  1.          !  Time scaling factor for SBC and CBC data (seconds in a day) 
     116   rn_trsfac(5)  = 7.9258065e-02    !  (  0.021 / 31. * 117 ) 
     117   rn_trsfac(7)  = 3.1316726e-01    !  ( 8.8   / 28.1 ) 
     118   rn_trsfac(14) = 6.2667860e-04    !  (  0.035 / 55.85 ) 
     119   rn_trsfac(23) = 5.2232143e-01    !  ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.3125/14 ) 
     120   rn_sbc_time   = 1.               !  Time scaling factor for SBC and CBC data (seconds in a day) 
    121121   ! 
    122    sn_trccbc(1)  = 'river.orca'      ,    -12            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    123    sn_trccbc(2)  = 'river.orca'      ,    -12            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    124    sn_trccbc(5)  = 'river.orca'      ,    -12            , 'riverdip'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    125    sn_trccbc(7)  = 'river.orca'      ,    -12            , 'riverdsi'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    126    sn_trccbc(10) = 'river.orca'      ,    -12            , 'riverdoc'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    127    sn_trccbc(14) = 'river.orca'      ,    -12            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    128    sn_trccbc(23) = 'river.orca'      ,    -12            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    129    rn_trcfac(1)  = 8.333e+01   !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
    130    rn_trcfac(2)  = 8.333e+01   !  ( 1e3 /12 ) 
    131    rn_trcfac(5)  = 3.935e+04   !  ( 1e3 / 31. * 122 ) 
    132    rn_trcfac(7)  = 3.588e+01   !  ( 1e3 / 28.1 ) 
    133    rn_trcfac(10) = 8.333e+01   !  ( 1e3 / 12 
    134    rn_trcfac(14) = 4.166e-03   !  ( 1e3 / 12 * 5e-5 ) 
    135    rn_trcfac(23) = 5.446e+02   !  (  1e3 / 14 * 7.625 ) 
    136    rn_cbc_time   = 3.1536e+7   !  Time scaling factor for CBC data (seconds in a year) 
     122   sn_trccbc(1)  = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     123   sn_trccbc(2)  = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     124   sn_trccbc(5)  = 'river.orca'      ,    120            , 'riverdip'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     125   sn_trccbc(7)  = 'river.orca'      ,    120            , 'riverdsi'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     126   sn_trccbc(10) = 'river.orca'      ,    120            , 'riverdoc'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     127   sn_trccbc(14) = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     128   sn_trccbc(23) = 'river.orca'      ,    120            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     129   rn_trcfac(1)  = 8.333333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
     130   rn_trcfac(2)  = 8.333333e+01      !  ( 1e3 /12 ) 
     131   rn_trcfac(5)  = 3.774193e+03   !  ( 1e3 / 31. * 117 ) 
     132   rn_trcfac(7)  = 3.558719e+01   !  ( 1e3 / 28.1 ) 
     133   rn_trcfac(10) = 8.333333e+01   !  ( 1e3 / 12 
     134   rn_trcfac(14) = 4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
     135   rn_trcfac(23) = 5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
     136   rn_cbc_time   = 3.1536e+7      !  Time scaling factor for CBC data (seconds in a year) 
    137137/ 
    138138!---------------------------------------------------------------------- 
  • NEMO/trunk/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_top_cfg

    r14608 r15287  
    114114   sn_trcsbc(14) = 'dust.orca.new'   ,       -1          , 'dustfer'     ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    115115   sn_trcsbc(23) = 'ndeposition.orca',      -12          , 'ndep'        ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    116    rn_trsfac(5)  = 8.264e-02   !  (  0.021 / 31. * 122 ) 
    117    rn_trsfac(7)  = 3.313e-01     !  ( 8.8   / 28.1 ) 
    118    rn_trsfac(14) = 6.266e-04   !  (  0.035 / 55.85 ) 
    119    rn_trsfac(23) =  5.4464e-01  !  ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.625/14 ) 
    120    rn_sbc_time   =  1.          !  Time scaling factor for SBC and CBC data (seconds in a day) 
     116   rn_trsfac(5)  = 7.9258065e-02    !  (  0.021 / 31. * 117 ) 
     117   rn_trsfac(7)  = 3.1316726e-01    !  ( 8.8   / 28.1 ) 
     118   rn_trsfac(14) = 6.2667860e-04    !  (  0.035 / 55.85 ) 
     119   rn_trsfac(23) = 5.2232143e-01    !  ( From kgN m-2 s-1 to molC l-1 ====> zfact = 7.3125/14 ) 
     120   rn_sbc_time   = 1.               !  Time scaling factor for SBC and CBC data (seconds in a day) 
    121121   ! 
    122    sn_trccbc(1)  = 'river.orca'      ,    -12            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    123    sn_trccbc(2)  = 'river.orca'      ,    -12            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    124    sn_trccbc(5)  = 'river.orca'      ,    -12            , 'riverdip'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    125    sn_trccbc(7)  = 'river.orca'      ,    -12            , 'riverdsi'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    126    sn_trccbc(10) = 'river.orca'      ,    -12            , 'riverdoc'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    127    sn_trccbc(14) = 'river.orca'      ,    -12            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    128    sn_trccbc(23) = 'river.orca'      ,    -12            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    129    rn_trcfac(1)  = 8.333e+01   !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
    130    rn_trcfac(2)  = 8.333e+01   !  ( 1e3 /12 ) 
    131    rn_trcfac(5)  = 3.935e+04   !  ( 1e3 / 31. * 122 ) 
    132    rn_trcfac(7)  = 3.588e+01   !  ( 1e3 / 28.1 ) 
    133    rn_trcfac(10) = 8.333e+01   !  ( 1e3 / 12 
    134    rn_trcfac(14) = 4.166e-03   !  ( 1e3 / 12 * 5e-5 ) 
    135    rn_trcfac(23) = 5.446e+02   !  (  1e3 / 14 * 7.625 ) 
    136    rn_cbc_time   = 3.1536e+7   !  Time scaling factor for CBC data (seconds in a year) 
     122   sn_trccbc(1)  = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     123   sn_trccbc(2)  = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     124   sn_trccbc(5)  = 'river.orca'      ,    120            , 'riverdip'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     125   sn_trccbc(7)  = 'river.orca'      ,    120            , 'riverdsi'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     126   sn_trccbc(10) = 'river.orca'      ,    120            , 'riverdoc'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     127   sn_trccbc(14) = 'river.orca'      ,    120            , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     128   sn_trccbc(23) = 'river.orca'      ,    120            , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     129   rn_trcfac(1)  = 8.333333e+01      !  ( data in Mg/m2/yr : 1e3/12/ryyss) 
     130   rn_trcfac(2)  = 8.333333e+01      !  ( 1e3 /12 ) 
     131   rn_trcfac(5)  = 3.774193e+03   !  ( 1e3 / 31. * 117 ) 
     132   rn_trcfac(7)  = 3.558719e+01   !  ( 1e3 / 28.1 ) 
     133   rn_trcfac(10) = 8.333333e+01   !  ( 1e3 / 12 
     134   rn_trcfac(14) = 4.166667e-03   !  ( 1e3 / 12 * 5e-5 ) 
     135   rn_trcfac(23) = 5.223214e+02   !  (  1e3 / 14 * 7.3125 ) 
     136   rn_cbc_time   = 3.1536e+7      !  Time scaling factor for CBC data (seconds in a year) 
    137137/ 
    138138!---------------------------------------------------------------------- 
  • NEMO/trunk/cfgs/SHARED/namelist_pisces_ref

    r12845 r15287  
    339339&nampiscal     !   parameters for Calcite chemistry 
    340340!----------------------------------------------------------------------- 
    341    kdca       =  6.       ! calcite dissolution rate constant (1/time) 
    342    nca        =  1.       ! order of dissolution reaction (dimensionless) 
     341   kdca       =  100.     ! calcite dissolution rate constant (1/time) 
     342   nca        =  4.7      ! order of dissolution reaction (dimensionless) 
    343343/ 
    344344!----------------------------------------------------------------------- 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zlys.F90

    r15090 r15287  
    6060      INTEGER  ::   ji, jj, jk, jn 
    6161      REAL(wp) ::   zdispot, zfact, zcalcon 
    62       REAL(wp) ::   zomegaca, zexcess, zexcess0 
     62      REAL(wp) ::   zomegaca, zexcess, zexcess0, zkd 
    6363      CHARACTER (len=25) ::   charout 
    6464      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zco3, zcaldiss, zhinit, zhi, zco3sat 
     
    9898         excess(ji,jj,jk) = 1._wp - zomegaca 
    9999         zexcess0 = MAX( 0., excess(ji,jj,jk) ) 
    100          zexcess  = zexcess0**nca 
    101100 
    102          ! AMOUNT CACO3 (12C) THAT RE-ENTERS SOLUTION 
    103          !       (ACCORDING TO THIS FORMULATION ALSO SOME PARTICULATE 
    104          !       CACO3 GETS DISSOLVED EVEN IN THE CASE OF OVERSATURATION) 
    105          zdispot = kdca * zexcess * tr(ji,jj,jk,jpcal,Kbb) 
     101         IF( zomegaca < 0.8 ) THEN 
     102            zexcess = zexcess0**nca 
     103            ! AMOUNT CACO3 THAT RE-ENTERS SOLUTION 
     104            zdispot = kdca * zexcess * tr(ji,jj,jk,jpcal,Kbb) 
     105         ELSE 
     106            zkd = kdca * 0.2**(nca - 0.11) 
     107            zexcess = zexcess0**0.11 
     108            zdispot = zkd * zexcess * tr(ji,jj,jk,jpcal,Kbb) 
     109        ENDIF 
     110 
    106111        !  CHANGE OF [CO3--] , [ALK], PARTICULATE [CACO3], 
    107112        !       AND [SUM(CO2)] DUE TO CACO3 DISSOLUTION/PRECIPITATION 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zsed.F90

    r15090 r15287  
    265265            tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0 
    266266            tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zfact / 3.0 
     267            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zfact * 2.0 / 3.0             
    267268            tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) - zfact * 2.0 / 3.0 
    268269            tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zfact * 1.0 / 3.0 
     
    282283            tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0 
    283284            tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zfact / 3.0 
     285            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zfact * 2.0 / 3.0             
    284286            tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) - 16.0 / 46.0 * zfact * ( 1.0 - 1.0 / 3.0 ) & 
    285287            &                     * ztrpo4(ji,jj,jk) / (ztrpo4(ji,jj,jk) + ztrdop(ji,jj,jk) + rtrn) 
Note: See TracChangeset for help on using the changeset viewer.