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 7646 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90 – NEMO

Ignore:
Timestamp:
2017-02-06T10:25:03+01:00 (7 years ago)
Author:
timgraham
Message:

Merge of dev_merge_2016 into trunk. UPDATE TO ARCHFILES NEEDED for XIOS2.
LIM_SRC_s/limrhg.F90 to follow in next commit due to change of kind (I'm unable to do it in this commit).
Merged using the following steps:

1) svn merge --reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk .
2) Resolve minor conflicts in sette.sh and namelist_cfg for ORCA2LIM3 (due to a change in trunk after branch was created)
3) svn commit
4) svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
5) svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2016/dev_merge_2016 .
6) At this stage I checked out a clean copy of the branch to compare against what is about to be committed to the trunk.
6) svn commit #Commit code to the trunk

In this commit I have also reverted a change to Fcheck_archfile.sh which was causing problems on the Paris machine.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90

    r5725 r7646  
    55   !!====================================================================== 
    66   !! History :  3.5  ! 2013    (M. Vancoppenolle, O. Aumont, G. Madec), original code 
    7    !! Comment ! probably not properly done when the second particle export 
    8    !! scheme (kriest) is used 
    9    !!---------------------------------------------------------------------- 
    10 #if defined key_pisces || defined key_pisces_reduced 
    11    !!---------------------------------------------------------------------- 
    12    !!   'key_pisces'                                       PISCES bio-model 
    137   !!---------------------------------------------------------------------- 
    148   !! trc_ice_pisces   : PISCES fake sea ice model setting 
     
    1812   USE oce_trc         ! Shared variables between ocean and passive tracers 
    1913   USE trc             ! Passive tracers common variables  
    20    USE phycst          ! Ocean physics parameters 
    2114   USE sms_pisces      ! PISCES Source Minus Sink variables 
    2215   USE in_out_manager 
     
    3730      !!---------------------------------------------------------------------- 
    3831 
    39       IF( lk_p4z ) THEN  ;   CALL p4z_ice_ini   !  PISCES 
    40       ELSE               ;   CALL p2z_ice_ini   !  LOBSTER 
     32      IF( ln_p4z .OR. ln_p5z ) THEN  ;   CALL p4z_ice_ini   !  PISCES 
     33      ELSE                           ;   CALL p2z_ice_ini   !  LOBSTER 
    4134      ENDIF 
    4235 
     
    4538 
    4639   SUBROUTINE p4z_ice_ini 
    47  
    48 #if defined key_pisces  
    4940      !!---------------------------------------------------------------------- 
    5041      !!                   ***  ROUTINE p4z_ice_ini *** 
     
    7566 
    7667                                        !--- Dummy variables 
    77       REAL(wp), DIMENSION(jp_pisces,2) :: zratio  ! effective ice-ocean tracer cc ratio 
    78       REAL(wp), DIMENSION(jp_pisces,4) :: zpisc   ! prescribes concentration  
     68      REAL(wp), DIMENSION(jpmaxtrc,2) :: zratio  ! effective ice-ocean tracer cc ratio 
     69      REAL(wp), DIMENSION(jpmaxtrc,4) :: zpisc   ! prescribes concentration  
    7970      !                                            !  1:global, 2:Arctic, 3:Antarctic, 4:Baltic 
    8071 
     
    10798      zpisc(jppo4,1) =  5.77e-7_wp / po4r  
    10899      zpisc(jppoc,1) =  1.27e-6_wp   
    109 #  if ! defined key_kriest 
    110100      zpisc(jpgoc,1) =  5.23e-8_wp   
    111101      zpisc(jpbfe,1) =  9.84e-13_wp  
    112 #  else 
    113       zpisc(jpnum,1) = 0. ! could not get this value since did not use it 
    114 #  endif 
    115102      zpisc(jpsil,1) =  7.36e-6_wp   
    116103      zpisc(jpdsi,1) =  1.07e-7_wp  
     
    129116      zpisc(jpnh4,1) =  3.22e-7_wp / rno3 
    130117 
     118      ! ln_p5z 
     119      zpisc(jppic,1) =  9.57e-8_wp 
     120      zpisc(jpnpi,1) =  9.57e-8_wp 
     121      zpisc(jpppi,1) =  9.57e-8_wp 
     122      zpisc(jppfe,1) =  1.76e-11_wp 
     123      zpisc(jppch,1) =  1.67e-7_wp 
     124      zpisc(jpnph,1) =  9.57e-8_wp 
     125      zpisc(jppph,1) =  9.57e-8_wp 
     126      zpisc(jpndi,1) =  4.24e-7_wp 
     127      zpisc(jppdi,1) =  4.24e-7_wp 
     128      zpisc(jppon,1) =  9.57e-8_wp 
     129      zpisc(jppop,1) =  9.57e-8_wp 
     130      zpisc(jpdon,1) =  2.04e-5_wp 
     131      zpisc(jpdop,1) =  2.04e-5_wp 
     132      zpisc(jpgon,1) =  5.23e-8_wp 
     133      zpisc(jpgop,1) =  5.23e-8_wp 
     134 
    131135      !--- Arctic specificities (dissolved inorganic & DOM) 
    132136      zpisc(jpdic,2) =  1.98e-3_wp  
     
    137141      zpisc(jppo4,2) =  4.09e-7_wp / po4r  
    138142      zpisc(jppoc,2) =  4.05e-7_wp   
    139 #  if ! defined key_kriest 
    140143      zpisc(jpgoc,2) =  2.84e-8_wp   
    141144      zpisc(jpbfe,2) =  7.03e-13_wp  
    142 #  else 
    143       zpisc(jpnum,2) =  0.00e-00_wp  
    144 #  endif 
    145145      zpisc(jpsil,2) =  6.87e-6_wp   
    146146      zpisc(jpdsi,2) =  1.73e-7_wp  
     
    159159      zpisc(jpnh4,2) =  6.15e-08_wp / rno3  
    160160 
     161      ! ln_p5z 
     162      zpisc(jppic,2) =  5.25e-7_wp 
     163      zpisc(jpnpi,2) =  5.25e-7_wp 
     164      zpisc(jpppi,2) =  5.25e-7_wp 
     165      zpisc(jppfe,2) =  1.75e-11_wp 
     166      zpisc(jppch,2) =  1.46e-07_wp 
     167      zpisc(jpnph,2) =  5.25e-7_wp 
     168      zpisc(jppph,2) =  5.25e-7_wp 
     169      zpisc(jpndi,2) =  7.75e-7_wp 
     170      zpisc(jppdi,2) =  7.75e-7_wp 
     171      zpisc(jppon,2) =  4.05e-7_wp 
     172      zpisc(jppop,2) =  4.05e-7_wp 
     173      zpisc(jpdon,2) =  6.00e-6_wp 
     174      zpisc(jpdop,2) =  6.00e-6_wp 
     175      zpisc(jpgon,2) =  2.84e-8_wp 
     176      zpisc(jpgop,2) =  2.84e-8_wp 
     177 
    161178      !--- Antarctic specificities (dissolved inorganic & DOM) 
    162179      zpisc(jpdic,3) =  2.20e-3_wp   
     
    167184      zpisc(jppo4,3) =  1.88e-6_wp / po4r   
    168185      zpisc(jppoc,3) =  1.13e-6_wp   
    169 #  if ! defined key_kriest 
    170186      zpisc(jpgoc,3) =  2.89e-8_wp   
    171187      zpisc(jpbfe,3) =  5.63e-13_wp  
    172 #  else 
    173       zpisc(jpnum,3) =  0.00e-00_wp  
    174 #  endif 
    175188      zpisc(jpsil,3) =  4.96e-5_wp   
    176189      zpisc(jpdsi,3) =  5.63e-7_wp  
     
    189202      zpisc(jpnh4,3) =  3.39e-7_wp / rno3   
    190203 
     204      ! ln_p5z 
     205      zpisc(jppic,3) =  8.10e-7_wp 
     206      zpisc(jpnpi,3) =  8.10e-7_wp 
     207      zpisc(jpppi,3) =  8.10e-7_wp  
     208      zpisc(jppfe,3) =  1.48e-11_wp 
     209      zpisc(jppch,3) =  2.02e-7_wp 
     210      zpisc(jpnph,3) =  9.57e-8_wp 
     211      zpisc(jppph,3) =  9.57e-8_wp 
     212      zpisc(jpndi,3) =  5.77e-7_wp 
     213      zpisc(jppdi,3) =  5.77e-7_wp 
     214      zpisc(jppon,3) =  1.13e-6_wp 
     215      zpisc(jppop,3) =  1.13e-6_wp 
     216      zpisc(jpdon,3) =  7.02e-6_wp 
     217      zpisc(jpdop,3) =  7.02e-6_wp 
     218      zpisc(jpgon,3) =  2.89e-8_wp 
     219      zpisc(jpgop,3) =  2.89e-8_wp 
     220 
     221 
    191222      !--- Baltic Sea particular case for ORCA configurations 
    192223      zpisc(jpdic,4) = 1.14e-3_wp 
     
    197228      zpisc(jppo4,4) = 2.85e-9_wp / po4r 
    198229      zpisc(jppoc,4) = 4.84e-7_wp 
    199 #  if ! defined key_kriest 
    200230      zpisc(jpgoc,4) = 1.05e-8_wp 
    201231      zpisc(jpbfe,4) = 4.97e-13_wp 
    202 #  else 
    203       zpisc(jpnum,4) = 0. ! could not get this value 
    204 #  endif 
    205232      zpisc(jpsil,4) = 4.91e-5_wp 
    206233      zpisc(jpdsi,4) = 3.25e-7_wp 
     
    218245      zpisc(jpno3,4) = 5.36e-5_wp / rno3 
    219246      zpisc(jpnh4,4) = 7.18e-7_wp / rno3 
     247 
     248      ! ln_p5z 
     249      zpisc(jppic,4) =  6.64e-7_wp 
     250      zpisc(jpnpi,4) =  6.64e-7_wp 
     251      zpisc(jpppi,4) =  6.64e-7_wp 
     252      zpisc(jppfe,4) =  3.89e-11_wp 
     253      zpisc(jppch,4) =  1.17e-7_wp 
     254      zpisc(jpnph,4) =  6.64e-7_wp 
     255      zpisc(jppph,4) =  6.64e-7_wp 
     256      zpisc(jpndi,4) =  3.41e-7_wp 
     257      zpisc(jppdi,4) =  3.41e-7_wp 
     258      zpisc(jppon,4) =  4.84e-7_wp 
     259      zpisc(jppop,4) =  4.84e-7_wp 
     260      zpisc(jpdon,4) =  1.06e-5_wp 
     261      zpisc(jpdop,4) =  1.06e-5_wp 
     262      zpisc(jpgon,4) =  1.05e-8_wp 
     263      zpisc(jpgop,4) =  1.05e-8_wp 
    220264  
    221265      DO jn = jp_pcs0, jp_pcs1 
     
    225269            WHERE( gphit(:,:) <  0._wp ) ; trc_o(:,:,jn) = zpisc(jn,3) ; END WHERE ! Antarctic  
    226270         ENDIF 
    227          IF( cp_cfg == "orca" ) THEN     !  Baltic Sea particular case for ORCA configurations 
     271         IF( cn_cfg == "orca" ) THEN     !  Baltic Sea particular case for ORCA configurations 
    228272             WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    229273                    54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
     
    264308        
    265309         !-- Baltic 
    266          IF( cp_cfg == "orca" ) THEN  ! Baltic treated seperately for ORCA configs 
     310         IF( cn_cfg == "orca" ) THEN  ! Baltic treated seperately for ORCA configs 
    267311            IF ( trc_ice_ratio(jn) >= - 1._wp ) THEN ! no prescribed conc. ; typically everything but iron)  
    268312               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
     
    279323      ! 
    280324      END DO ! jn 
    281 #endif 
    282  
     325      ! 
    283326   END SUBROUTINE p4z_ice_ini 
    284327 
    285328   SUBROUTINE p2z_ice_ini 
    286 #if defined key_pisces_reduced  
    287329      !!---------------------------------------------------------------------- 
    288330      !!                   ***  ROUTINE p2z_ice_ini *** 
     
    290332      !! ** Purpose :   Initialisation of the LOBSTER biochemical model 
    291333      !!---------------------------------------------------------------------- 
    292 #endif 
    293334   END SUBROUTINE p2z_ice_ini 
    294335 
    295  
    296 #else 
    297    !!---------------------------------------------------------------------- 
    298    !!   Dummy module                            No PISCES biochemical model 
    299    !!---------------------------------------------------------------------- 
    300 CONTAINS 
    301    SUBROUTINE trc_ice_ini_pisces         ! Empty routine 
    302    END SUBROUTINE trc_ice_ini_pisces 
    303 #endif 
    304336 
    305337   !!====================================================================== 
Note: See TracChangeset for help on using the changeset viewer.