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 5610 for branches/2015 – NEMO

Changeset 5610 for branches/2015


Ignore:
Timestamp:
2015-07-17T19:42:27+02:00 (9 years ago)
Author:
cetlod
Message:

bugfix to reduce lines over 132, see ticket #1564

Location:
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r5109 r5610  
    329329      !  
    330330      ! One level below 
    331       en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1 * ((zhsro(:,:)+fsdepw(:,:,2))/zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 
     331      en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1 * ((zhsro(:,:)+fsdepw(:,:,2)) & 
     332          &            / zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 
    332333      en(:,:,2) = MAX(en(:,:,2), rn_emin ) 
    333334      z_elem_a(:,:,2) = 0._wp  
     
    350351      z_elem_a(:,:,2) = 0._wp 
    351352      zkar(:,:)       = (rl_sf + (vkarmn-rl_sf)*(1.-exp(-rtrans*fsdept(:,:,1)/zhsro(:,:)) )) 
    352       zflxs(:,:)      = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 
     353      zflxs(:,:)      = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) & 
     354           &                      * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 
    353355 
    354356      en(:,:,2) = en(:,:,2) + zflxs(:,:)/fse3w(:,:,2) 
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90

    r5385 r5610  
    5858 
    5959                                        !--- Dummy variables 
    60       REAL(wp), DIMENSION(jptra,2) & 
    61                ::  zratio            ! effective ice-ocean tracer cc ratio 
     60      REAL(wp), DIMENSION(jp_pisces,2)  :: zratio  ! effective ice-ocean tracer cc ratio 
     61      REAL(wp), DIMENSION(jp_pisces,4)  :: zpisc   ! prescribes concentration  
     62      !                                            !  1:global, 2:Arctic, 3:Antarctic, 4:Baltic 
     63 
    6264      REAL(wp), DIMENSION(2) :: zrs  ! ice-ocean salinity ratio, 1 - global, 2- Baltic 
    6365      REAL(wp) :: zsice_bal          ! prescribed ice salinity in the Baltic 
     
    8082      ! fluxes 
    8183 
    82       !--- Global case  
    83       IF ( cn_trc_o(jpdic) == 'GL ' ) trc_o(:,:,jpdic) =  1.99e-3_wp  
    84       IF ( cn_trc_o(jpdoc) == 'GL ' ) trc_o(:,:,jpdoc) =  2.04e-5_wp  
    85       IF ( cn_trc_o(jptal) == 'GL ' ) trc_o(:,:,jptal) =  2.31e-3_wp  
    86       IF ( cn_trc_o(jpoxy) == 'GL ' ) trc_o(:,:,jpoxy) =  2.47e-4_wp 
    87       IF ( cn_trc_o(jpcal) == 'GL ' ) trc_o(:,:,jpcal) =  1.04e-8_wp 
    88       IF ( cn_trc_o(jppo4) == 'GL ' ) trc_o(:,:,jppo4) =  5.77e-7_wp / po4r  
    89       IF ( cn_trc_o(jppoc) == 'GL ' ) trc_o(:,:,jppoc) =  1.27e-6_wp   
    90 #  if ! defined key_kriest 
    91       IF ( cn_trc_o(jpgoc) == 'GL ' ) trc_o(:,:,jpgoc) =  5.23e-8_wp   
    92       IF ( cn_trc_o(jpbfe) == 'GL ' ) trc_o(:,:,jpbfe) =  9.84e-13_wp  
    93 #  else 
    94       IF ( cn_trc_o(jpnum) == 'GL ' ) trc_o(:,:,jpnum) = 0. ! could not get this value since did not use it 
    95 #  endif 
    96       IF ( cn_trc_o(jpsil) == 'GL ' ) trc_o(:,:,jpsil) =  7.36e-6_wp   
    97       IF ( cn_trc_o(jpdsi) == 'GL ' ) trc_o(:,:,jpdsi) =  1.07e-7_wp  
    98       IF ( cn_trc_o(jpgsi) == 'GL ' ) trc_o(:,:,jpgsi) =  1.53e-8_wp 
    99       IF ( cn_trc_o(jpphy) == 'GL ' ) trc_o(:,:,jpphy) =  9.57e-8_wp 
    100       IF ( cn_trc_o(jpdia) == 'GL ' ) trc_o(:,:,jpdia) =  4.24e-7_wp 
    101       IF ( cn_trc_o(jpzoo) == 'GL ' ) trc_o(:,:,jpzoo) =  6.07e-7_wp 
    102       IF ( cn_trc_o(jpmes) == 'GL ' ) trc_o(:,:,jpmes) =  3.44e-7_wp 
    103       IF ( cn_trc_o(jpfer) == 'GL ' ) trc_o(:,:,jpfer) =  4.06e-10_wp 
    104       IF ( cn_trc_o(jpsfe) == 'GL ' ) trc_o(:,:,jpsfe) =  2.51e-11_wp 
    105       IF ( cn_trc_o(jpdfe) == 'GL ' ) trc_o(:,:,jpdfe) =  6.57e-12_wp 
    106       IF ( cn_trc_o(jpnfe) == 'GL ' ) trc_o(:,:,jpnfe) =  1.76e-11_wp 
    107       IF ( cn_trc_o(jpnch) == 'GL ' ) trc_o(:,:,jpnch) =  1.67e-7_wp 
    108       IF ( cn_trc_o(jpdch) == 'GL ' ) trc_o(:,:,jpdch) =  1.02e-7_wp 
    109       IF ( cn_trc_o(jpno3) == 'GL ' ) trc_o(:,:,jpno3) =  5.79e-6_wp / rno3  
    110       IF ( cn_trc_o(jpnh4) == 'GL ' ) trc_o(:,:,jpnh4) =  3.22e-7_wp / rno3 
     84      !--- Global values 
     85      zpisc(jpdic,1) =  1.99e-3_wp  
     86      zpisc(jpdoc,1) =  2.04e-5_wp  
     87      zpisc(jptal,1) =  2.31e-3_wp  
     88      zpisc(jpoxy,1) =  2.47e-4_wp 
     89      zpisc(jpcal,1) =  1.04e-8_wp 
     90      zpisc(jppo4,1) =  5.77e-7_wp / po4r  
     91      zpisc(jppoc,1) =  1.27e-6_wp   
     92#  if ! defined key_kriest 
     93      zpisc(jpgoc,1) =  5.23e-8_wp   
     94      zpisc(jpbfe,1) =  9.84e-13_wp  
     95#  else 
     96      zpisc(jpnum,1) = 0. ! could not get this value since did not use it 
     97#  endif 
     98      zpisc(jpsil,1) =  7.36e-6_wp   
     99      zpisc(jpdsi,1) =  1.07e-7_wp  
     100      zpisc(jpgsi,1) =  1.53e-8_wp 
     101      zpisc(jpphy,1) =  9.57e-8_wp 
     102      zpisc(jpdia,1) =  4.24e-7_wp 
     103      zpisc(jpzoo,1) =  6.07e-7_wp 
     104      zpisc(jpmes,1) =  3.44e-7_wp 
     105      zpisc(jpfer,1) =  4.06e-10_wp 
     106      zpisc(jpsfe,1) =  2.51e-11_wp 
     107      zpisc(jpdfe,1) =  6.57e-12_wp 
     108      zpisc(jpnfe,1) =  1.76e-11_wp 
     109      zpisc(jpnch,1) =  1.67e-7_wp 
     110      zpisc(jpdch,1) =  1.02e-7_wp 
     111      zpisc(jpno3,1) =  5.79e-6_wp / rno3  
     112      zpisc(jpnh4,1) =  3.22e-7_wp / rno3 
    111113 
    112114      !--- Arctic specificities (dissolved inorganic & DOM) 
    113       IF ( cn_trc_o(jpdic) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdic) =  1.98e-3_wp  ; END WHERE ; ENDIF 
    114       IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdoc) =  6.00e-6_wp  ; END WHERE ; ENDIF 
    115       IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jptal) =  2.13e-3_wp  ; END WHERE ; ENDIF 
    116       IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpoxy) =  3.65e-4_wp  ; END WHERE ; ENDIF 
    117       IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpcal) =  1.50e-9_wp  ; END WHERE ; ENDIF 
    118       IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppo4) =  4.09e-7_wp / po4r ; END WHERE ; ENDIF 
    119       IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppoc) =  4.05e-7_wp  ; END WHERE ; ENDIF 
    120 #  if ! defined key_kriest 
    121       IF ( cn_trc_o(jpgoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgoc) =  2.84e-8_wp  ; END WHERE ; ENDIF 
    122       IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpbfe) =  7.03e-13_wp ; END WHERE ; ENDIF 
    123 #  else 
    124       IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnum) =  0.00e-00_wp ; END WHERE ; ENDIF 
    125 #  endif 
    126       IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsil) =  6.87e-6_wp  ; END WHERE ; ENDIF 
    127       IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdsi) =  1.73e-7_wp  ; END WHERE ; ENDIF 
    128       IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgsi) =  7.93e-9_wp  ; END WHERE ; ENDIF 
    129       IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpphy) =  5.25e-7_wp  ; END WHERE ; ENDIF 
    130       IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdia) =  7.75e-7_wp  ; END WHERE ; ENDIF 
    131       IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpzoo) =  3.34e-7_wp  ; END WHERE ; ENDIF 
    132       IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpmes) =  2.49e-7_wp  ; END WHERE ; ENDIF 
    133       IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpfer) =  1.43e-9_wp  ; END WHERE ; ENDIF 
    134       IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsfe) =  2.21e-11_wp ; END WHERE ; ENDIF 
    135       IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdfe) =  2.04e-11_wp ; END WHERE ; ENDIF 
    136       IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnfe) =  1.75e-11_wp ; END WHERE ; ENDIF 
    137       IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnch) =  1.46e-07_wp ; END WHERE ; ENDIF 
    138       IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdch) =  2.36e-07_wp ; END WHERE ; ENDIF 
    139       IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpno3) =  3.51e-06_wp / rno3 ; END WHERE ; ENDIF 
    140       IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnh4) =  6.15e-08_wp / rno3 ; END WHERE ; ENDIF 
     115      zpisc(jpdic,2) =  1.98e-3_wp  
     116      zpisc(jpdoc,2) =  6.00e-6_wp  
     117      zpisc(jptal,2) =  2.13e-3_wp  
     118      zpisc(jpoxy,2) =  3.65e-4_wp   
     119      zpisc(jpcal,2) =  1.50e-9_wp   
     120      zpisc(jppo4,2) =  4.09e-7_wp / po4r  
     121      zpisc(jppoc,2) =  4.05e-7_wp   
     122#  if ! defined key_kriest 
     123      zpisc(jpgoc,2) =  2.84e-8_wp   
     124      zpisc(jpbfe,2) =  7.03e-13_wp  
     125#  else 
     126      zpisc(jpnum,2) =  0.00e-00_wp  
     127#  endif 
     128      zpisc(jpsil,2) =  6.87e-6_wp   
     129      zpisc(jpdsi,2) =  1.73e-7_wp  
     130      zpisc(jpgsi,2) =  7.93e-9_wp 
     131      zpisc(jpphy,2) =  5.25e-7_wp   
     132      zpisc(jpdia,2) =  7.75e-7_wp  
     133      zpisc(jpzoo,2) =  3.34e-7_wp 
     134      zpisc(jpmes,2) =  2.49e-7_wp   
     135      zpisc(jpfer,2) =  1.43e-9_wp  
     136      zpisc(jpsfe,2) =  2.21e-11_wp  
     137      zpisc(jpdfe,2) =  2.04e-11_wp  
     138      zpisc(jpnfe,2) =  1.75e-11_wp  
     139      zpisc(jpnch,2) =  1.46e-07_wp  
     140      zpisc(jpdch,2) =  2.36e-07_wp  
     141      zpisc(jpno3,2) =  3.51e-06_wp / rno3  
     142      zpisc(jpnh4,2) =  6.15e-08_wp / rno3  
    141143 
    142144      !--- Antarctic specificities (dissolved inorganic & DOM) 
    143       IF ( cn_trc_o(jpdic) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdic) =  2.20e-3_wp  ; END WHERE ; ENDIF 
    144       IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdoc) =  7.02e-6_wp  ; END WHERE ; ENDIF 
    145       IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jptal) =  2.37e-3_wp  ; END WHERE ; ENDIF 
    146       IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpoxy) =  3.42e-4_wp  ; END WHERE ; ENDIF 
    147       IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpcal) =  3.17e-9_wp  ; END WHERE ; ENDIF 
    148       IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jppo4) =  1.88e-6_wp / po4r  ; END WHERE ; ENDIF 
    149       IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jppoc) =  1.13e-6_wp  ; END WHERE ; ENDIF 
    150 #  if ! defined key_kriest 
    151       IF ( cn_trc_o(jpgoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpgoc) =  2.89e-8_wp  ; END WHERE ; ENDIF 
    152       IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpbfe) =  5.63e-13_wp ; END WHERE ; ENDIF 
    153 #  else 
    154       IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnum) =  0.00e-00_wp ; END WHERE ; ENDIF 
    155 #  endif 
    156       IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpsil) =  4.96e-5_wp  ; END WHERE ; ENDIF 
    157       IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdsi) =  5.63e-7_wp  ; END WHERE ; ENDIF 
    158       IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpgsi) =  5.35e-8_wp  ; END WHERE ; ENDIF 
    159       IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpphy) =  8.10e-7_wp  ; END WHERE ; ENDIF 
    160       IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdia) =  5.77e-7_wp  ; END WHERE ; ENDIF 
    161       IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpzoo) =  6.68e-7_wp  ; END WHERE ; ENDIF 
    162       IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpmes) =  3.55e-7_wp  ; END WHERE ; ENDIF 
    163       IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpfer) =  1.62e-10_wp ; END WHERE ; ENDIF 
    164       IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpsfe) =  2.29e-11_wp ; END WHERE ; ENDIF 
    165       IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdfe) =  8.75e-12_wp ; END WHERE ; ENDIF 
    166       IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnfe) =  1.48e-11_wp ; END WHERE ; ENDIF 
    167       IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnch) =  2.02e-7_wp  ; END WHERE ; ENDIF 
    168       IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdch) =  1.60e-7_wp  ; END WHERE ; ENDIF 
    169       IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpno3) =  2.64e-5_wp / rno3  ; END WHERE ; ENDIF 
    170       IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnh4) =  3.39e-7_wp / rno3  ; END WHERE ; ENDIF 
     145      zpisc(jpdic,3) =  2.20e-3_wp   
     146      zpisc(jpdoc,3) =  7.02e-6_wp   
     147      zpisc(jptal,3) =  2.37e-3_wp   
     148      zpisc(jpoxy,3) =  3.42e-4_wp   
     149      zpisc(jpcal,3) =  3.17e-9_wp   
     150      zpisc(jppo4,3) =  1.88e-6_wp / po4r   
     151      zpisc(jppoc,3) =  1.13e-6_wp   
     152#  if ! defined key_kriest 
     153      zpisc(jpgoc,3) =  2.89e-8_wp   
     154      zpisc(jpbfe,3) =  5.63e-13_wp  
     155#  else 
     156      zpisc(jpnum,3) =  0.00e-00_wp  
     157#  endif 
     158      zpisc(jpsil,3) =  4.96e-5_wp   
     159      zpisc(jpdsi,3) =  5.63e-7_wp  
     160      zpisc(jpgsi,3) =  5.35e-8_wp 
     161      zpisc(jpphy,3) =  8.10e-7_wp   
     162      zpisc(jpdia,3) =  5.77e-7_wp  
     163      zpisc(jpzoo,3) =  6.68e-7_wp 
     164      zpisc(jpmes,3) =  3.55e-7_wp   
     165      zpisc(jpfer,3) =  1.62e-10_wp 
     166      zpisc(jpsfe,3) =  2.29e-11_wp  
     167      zpisc(jpdfe,3) =  8.75e-12_wp 
     168      zpisc(jpnfe,3) =  1.48e-11_wp  
     169      zpisc(jpnch,3) =  2.02e-7_wp   
     170      zpisc(jpdch,3) =  1.60e-7_wp   
     171      zpisc(jpno3,3) =  2.64e-5_wp / rno3   
     172      zpisc(jpnh4,3) =  3.39e-7_wp / rno3   
    171173 
    172174      !--- Baltic Sea particular case for ORCA configurations 
    173       IF( cp_cfg == "orca" ) THEN            ! Baltic mask 
    174          WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    175                 54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
    176          trc_o(:,:,jpdic) = 1.14e-3_wp 
    177          trc_o(:,:,jpdoc) = 1.06e-5_wp 
    178          trc_o(:,:,jptal) = 1.16e-3_wp 
    179          trc_o(:,:,jpoxy) = 3.71e-4_wp 
    180          trc_o(:,:,jpcal) = 1.51e-9_wp 
    181          trc_o(:,:,jppo4) = 2.85e-9_wp / po4r 
    182          trc_o(:,:,jppoc) = 4.84e-7_wp 
    183 #  if ! defined key_kriest 
    184          trc_o(:,:,jpgoc) = 1.05e-8_wp 
    185          trc_o(:,:,jpbfe) = 4.97e-13_wp 
    186 #  else 
    187          trc_o(:,:,jpnum) = 0. ! could not get this value 
    188 #  endif 
    189          trc_o(:,:,jpsil) = 4.91e-5_wp 
    190          trc_o(:,:,jpdsi) = 3.25e-7_wp 
    191          trc_o(:,:,jpgsi) = 1.93e-8_wp 
    192          trc_o(:,:,jpphy) = 6.64e-7_wp 
    193          trc_o(:,:,jpdia) = 3.41e-7_wp 
    194          trc_o(:,:,jpzoo) = 3.83e-7_wp 
    195          trc_o(:,:,jpmes) = 0.225e-6_wp 
    196          trc_o(:,:,jpfer) = 2.45e-9_wp 
    197          trc_o(:,:,jpsfe) = 3.89e-11_wp 
    198          trc_o(:,:,jpdfe) = 1.33e-11_wp 
    199          trc_o(:,:,jpnfe) = 2.62e-11_wp 
    200          trc_o(:,:,jpnch) = 1.17e-7_wp 
    201          trc_o(:,:,jpdch) = 9.69e-8_wp 
    202          trc_o(:,:,jpno3) = 5.36e-5_wp / rno3 
    203          trc_o(:,:,jpnh4) = 7.18e-7_wp / rno3 
    204          END WHERE 
    205       ENDIF ! cfg 
     175      zpisc(jpdic,4) = 1.14e-3_wp 
     176      zpisc(jpdoc,4) = 1.06e-5_wp 
     177      zpisc(jptal,4) = 1.16e-3_wp 
     178      zpisc(jpoxy,4) = 3.71e-4_wp 
     179      zpisc(jpcal,4) = 1.51e-9_wp 
     180      zpisc(jppo4,4) = 2.85e-9_wp / po4r 
     181      zpisc(jppoc,4) = 4.84e-7_wp 
     182#  if ! defined key_kriest 
     183      zpisc(jpgoc,4) = 1.05e-8_wp 
     184      zpisc(jpbfe,4) = 4.97e-13_wp 
     185#  else 
     186      zpisc(jpnum,4) = 0. ! could not get this value 
     187#  endif 
     188      zpisc(jpsil,4) = 4.91e-5_wp 
     189      zpisc(jpdsi,4) = 3.25e-7_wp 
     190      zpisc(jpgsi,4) = 1.93e-8_wp 
     191      zpisc(jpphy,4) = 6.64e-7_wp 
     192      zpisc(jpdia,4) = 3.41e-7_wp 
     193      zpisc(jpzoo,4) = 3.83e-7_wp 
     194      zpisc(jpmes,4) = 0.225e-6_wp 
     195      zpisc(jpfer,4) = 2.45e-9_wp 
     196      zpisc(jpsfe,4) = 3.89e-11_wp 
     197      zpisc(jpdfe,4) = 1.33e-11_wp 
     198      zpisc(jpnfe,4) = 2.62e-11_wp 
     199      zpisc(jpnch,4) = 1.17e-7_wp 
     200      zpisc(jpdch,4) = 9.69e-8_wp 
     201      zpisc(jpno3,4) = 5.36e-5_wp / rno3 
     202      zpisc(jpnh4,4) = 7.18e-7_wp / rno3 
     203  
     204      DO jn = jp_pcs0, jp_pcs1 
     205         IF( cn_trc_o(jn) == 'GL ' ) trc_o(:,:,jn) = zpisc(jn,1)  ! Global case 
     206         IF( cn_trc_o(jn) == 'AA ' ) THEN  
     207            WHERE( gphit(:,:) >= 0._wp ) ; trc_o(:,:,jn) = zpisc(jn,2) ; END WHERE ! Arctic  
     208            WHERE( gphit(:,:) <  0._wp ) ; trc_o(:,:,jn) = zpisc(jn,3) ; END WHERE ! Antarctic  
     209         ENDIF 
     210         IF( cp_cfg == "orca" ) THEN     !  Baltic Sea particular case for ORCA configurations 
     211             WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
     212                    54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
     213                    trc_o(:,:,jn) = zpisc(jn,4) 
     214            END WHERE 
     215         ENDIF  
     216      ENDDO 
     217 
     218 
    206219 
    207220      !----------------------------- 
     
    217230 
    218231      DO jn = jp_pcs0, jp_pcs1 
    219          IF ( trc_ice_ratio(jn) >= 0._wp )  zratio(jn,:) = trc_ice_ratio(jn) 
    220          IF ( trc_ice_ratio(jn) == -1._wp ) zratio(jn,:) = zrs(:) 
    221          IF ( trc_ice_ratio(jn) == -2._wp ) zratio(jn,:) = -9999.99_wp 
     232         IF( trc_ice_ratio(jn) >= 0._wp )  zratio(jn,:) = trc_ice_ratio(jn) 
     233         IF( trc_ice_ratio(jn) == -1._wp ) zratio(jn,:) = zrs(:) 
     234         IF( trc_ice_ratio(jn) == -2._wp ) zratio(jn,:) = -9999.99_wp 
    222235      END DO 
    223236 
     
    227240      DO jn = jp_pcs0, jp_pcs1 
    228241         !-- Everywhere but in the Baltic 
    229          IF ( trc_ice_ratio(jn) >= -1._wp ) THEN !! no prescribed concentration 
    230                                               !! (typically everything but iron)  
     242         IF ( trc_ice_ratio(jn) >= -1._wp ) THEN ! no prescribed conc. ; typically everything but iron)  
    231243            trc_i(:,:,jn) = zratio(jn,1) * trc_o(:,:,jn)  
    232          ELSE                                 !! prescribed concentration 
     244         ELSE                                    ! prescribed concentration 
    233245            trc_i(:,:,jn) = trc_ice_prescr(jn) 
    234246         ENDIF 
    235247        
    236248         !-- Baltic 
    237          IF( cp_cfg == "orca" ) THEN !! Baltic treated seperately for ORCA configs 
    238             IF ( trc_ice_ratio(jn) >= - 1._wp ) THEN !! no prescribed concentration 
    239                                                  !! (typically everything but iron)  
     249         IF( cp_cfg == "orca" ) THEN  ! Baltic treated seperately for ORCA configs 
     250            IF ( trc_ice_ratio(jn) >= - 1._wp ) THEN ! no prescribed conc. ; typically everything but iron)  
    240251               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    241252                      54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
    242253                     trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn)  
    243254               END WHERE 
    244             ELSE                                 !! prescribed tracer concentration in ice 
     255            ELSE                                 ! prescribed tracer concentration in ice 
    245256               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    246257                   54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
Note: See TracChangeset for help on using the changeset viewer.