Changeset 5611 for trunk/NEMOGCM
- Timestamp:
- 2015-07-17T19:42:59+02:00 (9 years ago)
- Location:
- trunk/NEMOGCM/NEMO
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r5109 r5611 329 329 ! 330 330 ! 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) 332 333 en(:,:,2) = MAX(en(:,:,2), rn_emin ) 333 334 z_elem_a(:,:,2) = 0._wp … … 350 351 z_elem_a(:,:,2) = 0._wp 351 352 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) 353 355 354 356 en(:,:,2) = en(:,:,2) + zflxs(:,:)/fse3w(:,:,2) -
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90
r5385 r5611 58 58 59 59 !--- 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 62 64 REAL(wp), DIMENSION(2) :: zrs ! ice-ocean salinity ratio, 1 - global, 2- Baltic 63 65 REAL(wp) :: zsice_bal ! prescribed ice salinity in the Baltic … … 80 82 ! fluxes 81 83 82 !--- Global case83 IF ( cn_trc_o(jpdic) == 'GL ' ) trc_o(:,:,jpdic) = 1.99e-3_wp84 IF ( cn_trc_o(jpdoc) == 'GL ' ) trc_o(:,:,jpdoc) = 2.04e-5_wp85 IF ( cn_trc_o(jptal) == 'GL ' ) trc_o(:,:,jptal) = 2.31e-3_wp86 IF ( cn_trc_o(jpoxy) == 'GL ' ) trc_o(:,:,jpoxy) = 2.47e-4_wp87 IF ( cn_trc_o(jpcal) == 'GL ' ) trc_o(:,:,jpcal) = 1.04e-8_wp88 IF ( cn_trc_o(jppo4) == 'GL ' ) trc_o(:,:,jppo4) = 5.77e-7_wp / po4r89 IF ( cn_trc_o(jppoc) == 'GL ' ) trc_o(:,:,jppoc) = 1.27e-6_wp90 # if ! defined key_kriest 91 IF ( cn_trc_o(jpgoc) == 'GL ' ) trc_o(:,:,jpgoc) = 5.23e-8_wp92 IF ( cn_trc_o(jpbfe) == 'GL ' ) trc_o(:,:,jpbfe) = 9.84e-13_wp93 # else 94 IF ( cn_trc_o(jpnum) == 'GL ' ) trc_o(:,:,jpnum) = 0. ! could not get this value since did not use it95 # endif 96 IF ( cn_trc_o(jpsil) == 'GL ' ) trc_o(:,:,jpsil) = 7.36e-6_wp97 IF ( cn_trc_o(jpdsi) == 'GL ' ) trc_o(:,:,jpdsi) = 1.07e-7_wp98 IF ( cn_trc_o(jpgsi) == 'GL ' ) trc_o(:,:,jpgsi) = 1.53e-8_wp99 IF ( cn_trc_o(jpphy) == 'GL ' ) trc_o(:,:,jpphy) = 9.57e-8_wp100 IF ( cn_trc_o(jpdia) == 'GL ' ) trc_o(:,:,jpdia) = 4.24e-7_wp101 IF ( cn_trc_o(jpzoo) == 'GL ' ) trc_o(:,:,jpzoo) = 6.07e-7_wp102 IF ( cn_trc_o(jpmes) == 'GL ' ) trc_o(:,:,jpmes) = 3.44e-7_wp103 IF ( cn_trc_o(jpfer) == 'GL ' ) trc_o(:,:,jpfer) = 4.06e-10_wp104 IF ( cn_trc_o(jpsfe) == 'GL ' ) trc_o(:,:,jpsfe) = 2.51e-11_wp105 IF ( cn_trc_o(jpdfe) == 'GL ' ) trc_o(:,:,jpdfe) = 6.57e-12_wp106 IF ( cn_trc_o(jpnfe) == 'GL ' ) trc_o(:,:,jpnfe) = 1.76e-11_wp107 IF ( cn_trc_o(jpnch) == 'GL ' ) trc_o(:,:,jpnch) = 1.67e-7_wp108 IF ( cn_trc_o(jpdch) == 'GL ' ) trc_o(:,:,jpdch) = 1.02e-7_wp109 IF ( cn_trc_o(jpno3) == 'GL ' ) trc_o(:,:,jpno3) = 5.79e-6_wp / rno3110 IF ( cn_trc_o(jpnh4) == 'GL ' ) trc_o(:,:,jpnh4) = 3.22e-7_wp / rno384 !--- 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 111 113 112 114 !--- 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 ; ENDIF114 IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdoc) = 6.00e-6_wp ; END WHERE ; ENDIF115 IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jptal) = 2.13e-3_wp ; END WHERE ; ENDIF116 IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpoxy) = 3.65e-4_wp ; END WHERE ; ENDIF117 IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpcal) = 1.50e-9_wp ; END WHERE ; ENDIF118 IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppo4) = 4.09e-7_wp / po4r ; END WHERE ; ENDIF119 IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppoc) = 4.05e-7_wp ; END WHERE ; ENDIF120 # 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 ; ENDIF122 IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpbfe) = 7.03e-13_wp ; END WHERE ; ENDIF123 # else 124 IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnum) = 0.00e-00_wp ; END WHERE ; ENDIF125 # endif 126 IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsil) = 6.87e-6_wp ; END WHERE ; ENDIF127 IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdsi) = 1.73e-7_wp ; END WHERE ; ENDIF128 IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgsi) = 7.93e-9_wp ; END WHERE ; ENDIF129 IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpphy) = 5.25e-7_wp ; END WHERE ; ENDIF130 IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdia) = 7.75e-7_wp ; END WHERE ; ENDIF131 IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpzoo) = 3.34e-7_wp ; END WHERE ; ENDIF132 IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpmes) = 2.49e-7_wp ; END WHERE ; ENDIF133 IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpfer) = 1.43e-9_wp ; END WHERE ; ENDIF134 IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsfe) = 2.21e-11_wp ; END WHERE ; ENDIF135 IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdfe) = 2.04e-11_wp ; END WHERE ; ENDIF136 IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnfe) = 1.75e-11_wp ; END WHERE ; ENDIF137 IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnch) = 1.46e-07_wp ; END WHERE ; ENDIF138 IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdch) = 2.36e-07_wp ; END WHERE ; ENDIF139 IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpno3) = 3.51e-06_wp / rno3 ; END WHERE ; ENDIF140 IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnh4) = 6.15e-08_wp / rno3 ; END WHERE ; ENDIF115 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 141 143 142 144 !--- 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 ; ENDIF144 IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpdoc) = 7.02e-6_wp ; END WHERE ; ENDIF145 IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jptal) = 2.37e-3_wp ; END WHERE ; ENDIF146 IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpoxy) = 3.42e-4_wp ; END WHERE ; ENDIF147 IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpcal) = 3.17e-9_wp ; END WHERE ; ENDIF148 IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jppo4) = 1.88e-6_wp / po4r ; END WHERE ; ENDIF149 IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jppoc) = 1.13e-6_wp ; END WHERE ; ENDIF150 # 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 ; ENDIF152 IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpbfe) = 5.63e-13_wp ; END WHERE ; ENDIF153 # else 154 IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpnum) = 0.00e-00_wp ; END WHERE ; ENDIF155 # endif 156 IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpsil) = 4.96e-5_wp ; END WHERE ; ENDIF157 IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpdsi) = 5.63e-7_wp ; END WHERE ; ENDIF158 IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpgsi) = 5.35e-8_wp ; END WHERE ; ENDIF159 IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpphy) = 8.10e-7_wp ; END WHERE ; ENDIF160 IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpdia) = 5.77e-7_wp ; END WHERE ; ENDIF161 IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpzoo) = 6.68e-7_wp ; END WHERE ; ENDIF162 IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpmes) = 3.55e-7_wp ; END WHERE ; ENDIF163 IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpfer) = 1.62e-10_wp ; END WHERE ; ENDIF164 IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpsfe) = 2.29e-11_wp ; END WHERE ; ENDIF165 IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpdfe) = 8.75e-12_wp ; END WHERE ; ENDIF166 IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpnfe) = 1.48e-11_wp ; END WHERE ; ENDIF167 IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpnch) = 2.02e-7_wp ; END WHERE ; ENDIF168 IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpdch) = 1.60e-7_wp ; END WHERE ; ENDIF169 IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpno3) = 2.64e-5_wp / rno3 ; END WHERE ; ENDIF170 IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) < 00._wp ) ; trc_o(:,:,jpnh4) = 3.39e-7_wp / rno3 ; END WHERE ; ENDIF145 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 171 173 172 174 !--- 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 206 219 207 220 !----------------------------- … … 217 230 218 231 DO jn = jp_pcs0, jp_pcs1 219 IF 220 IF 221 IF 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 222 235 END DO 223 236 … … 227 240 DO jn = jp_pcs0, jp_pcs1 228 241 !-- 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) 231 243 trc_i(:,:,jn) = zratio(jn,1) * trc_o(:,:,jn) 232 ELSE !! prescribed concentration244 ELSE ! prescribed concentration 233 245 trc_i(:,:,jn) = trc_ice_prescr(jn) 234 246 ENDIF 235 247 236 248 !-- 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) 240 251 WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND. & 241 252 54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 242 253 trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn) 243 254 END WHERE 244 ELSE ! !prescribed tracer concentration in ice255 ELSE ! prescribed tracer concentration in ice 245 256 WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND. & 246 257 54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp )
Note: See TracChangeset
for help on using the changeset viewer.