- Timestamp:
- 2021-11-28T18:59:49+01:00 (3 years ago)
- Location:
- NEMO/branches/2021/ticket2632_r14588_theta_sbcblk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/ticket2632_r14588_theta_sbcblk
- Property svn:externals
-
old new 9 9 10 10 # SETTE 11 ^/utils/CI/sette@14244 sette 11 ^/utils/CI/sette@HEAD sette 12
-
- Property svn:externals
-
NEMO/branches/2021/ticket2632_r14588_theta_sbcblk/src/TOP/PISCES/P4Z/p4zche.F90
r14086 r15548 179 179 ! 0.04°C relative to an exact computation 180 180 ! --------------------------------------------------------------------- 181 DO_3D( 1, 1, 1, 1, 1, jpk )181 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk ) 182 182 zpres = gdept(ji,jj,jk,Kmm) / 1000. 183 183 za1 = 0.04 * ( 1.0 + 0.185 * ts(ji,jj,jk,jp_tem,Kmm) + 0.035 * (salinprac(ji,jj,jk) - 35.0) ) … … 188 188 ! CHEMICAL CONSTANTS - SURFACE LAYER 189 189 ! ---------------------------------- 190 DO_2D( 1, 1, 1, 1)190 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 191 191 ! ! SET ABSOLUTE TEMPERATURE 192 192 ztkel = tempis(ji,jj,1) + 273.15 … … 195 195 ! ! LN(K0) OF SOLUBILITY OF CO2 (EQ. 12, WEISS, 1980) 196 196 ! ! AND FOR THE ATMOSPHERE FOR NON IDEAL GAS 197 zcek1 = 9 345.17/ztkel - 60.2409 + 23.3585 * LOG(zt) + zsal*(0.023517 - 0.00023656*ztkel &198 & + 0.00 47036e-4*ztkel**2)199 chemc(ji,jj,1) = EXP( zcek1 ) * 1E-6 * rhop(ji,jj,1) / 1000.! mol/(L atm)197 zcek1 = 9050.69/ztkel - 58.0931 + 22.2940 * LOG(zt) + zsal*(0.027766 - 0.00025888*ztkel & 198 & + 0.0050578e-4*ztkel**2) 199 chemc(ji,jj,1) = EXP( zcek1 ) * 1E-6 ! mol/(L atm) 200 200 chemc(ji,jj,2) = -1636.75 + 12.0408*ztkel - 0.0327957*ztkel**2 + 0.0000316528*ztkel**3 201 201 chemc(ji,jj,3) = 57.7 - 0.118*ztkel … … 204 204 ! OXYGEN SOLUBILITY - DEEP OCEAN 205 205 ! ------------------------------- 206 DO_3D( 1, 1, 1, 1, 1, jpk )206 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk ) 207 207 ztkel = tempis(ji,jj,jk) + 273.15 208 208 zsal = salinprac(ji,jj,jk) + ( 1.- tmask(ji,jj,jk) ) * 35. … … 223 223 ! CHEMICAL CONSTANTS - DEEP OCEAN 224 224 ! ------------------------------- 225 DO_3D( 1, 1, 1, 1, 1, jpk )225 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk ) 226 226 ! SET PRESSION ACCORDING TO SAUNDER (1980) 227 227 zplat = SIN ( ABS(gphit(ji,jj)*3.141592654/180.) ) … … 451 451 IF( ln_timing ) CALL timing_start('ahini_for_at') 452 452 ! 453 DO_3D( 1, 1, 1, 1, 1, jpk )453 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk ) 454 454 p_alkcb = tr(ji,jj,jk,jptal,Kbb) * 1000. / (rhop(ji,jj,jk) + rtrn) 455 455 p_dictot = tr(ji,jj,jk,jpdic,Kbb) * 1000. / (rhop(ji,jj,jk) + rtrn) … … 549 549 550 550 ! TOTAL H+ scale: conversion factor for Htot = aphscale * Hfree 551 DO_3D( 1, 1, 1, 1, 1, jpk )551 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk ) 552 552 IF (rmask(ji,jj,jk) == 1.) THEN 553 553 p_alktot = tr(ji,jj,jk,jptal,Kbb) * 1000. / (rhop(ji,jj,jk) + rtrn) … … 578 578 579 579 DO jn = 1, jp_maxniter_atgen 580 DO_3D( 1, 1, 1, 1, 1, jpk )580 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk ) 581 581 IF (rmask(ji,jj,jk) == 1.) THEN 582 582 zfact = rhop(ji,jj,jk) / 1000. + rtrn
Note: See TracChangeset
for help on using the changeset viewer.