- Timestamp:
- 2019-11-22T15:29:17+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p5zmeso.F90
r11536 r11949 59 59 CONTAINS 60 60 61 SUBROUTINE p5z_meso( kt, knt )61 SUBROUTINE p5z_meso( kt, knt, Kbb, Krhs ) 62 62 !!--------------------------------------------------------------------- 63 63 !! *** ROUTINE p5z_meso *** … … 67 67 !! ** Method : - ??? 68 68 !!--------------------------------------------------------------------- 69 INTEGER, INTENT(in) :: kt, knt ! ocean time step 69 INTEGER, INTENT(in) :: kt, knt ! ocean time step 70 INTEGER, INTENT(in) :: Kbb, Krhs ! time level indices 70 71 INTEGER :: ji, jj, jk 71 72 REAL(wp) :: zcompadi, zcompaph, zcompapoc, zcompaz, zcompam, zcompames … … 108 109 DO jj = 1, jpj 109 110 DO ji = 1, jpi 110 zcompam = MAX( ( tr b(ji,jj,jk,jpmes) - 1.e-9 ), 0.e0 )111 zcompam = MAX( ( tr(ji,jj,jk,jpmes,Kbb) - 1.e-9 ), 0.e0 ) 111 112 zfact = xstep * tgfunc2(ji,jj,jk) * zcompam 112 113 113 114 ! Michaelis-Menten mortality rates of mesozooplankton 114 115 ! --------------------------------------------------- 115 zrespz = resrat2 * zfact * ( tr b(ji,jj,jk,jpmes) / ( xkmort + trb(ji,jj,jk,jpmes) ) &116 zrespz = resrat2 * zfact * ( tr(ji,jj,jk,jpmes,Kbb) / ( xkmort + tr(ji,jj,jk,jpmes,Kbb) ) & 116 117 & + 3. * nitrfac(ji,jj,jk) ) 117 118 … … 119 120 ! no real reason except that it seems to be more stable and may mimic predation 120 121 ! --------------------------------------------------------------- 121 ztortz = mzrat2 * 1.e6 * zfact * tr b(ji,jj,jk,jpmes) * (1. - nitrfac(ji,jj,jk))122 ztortz = mzrat2 * 1.e6 * zfact * tr(ji,jj,jk,jpmes,Kbb) * (1. - nitrfac(ji,jj,jk)) 122 123 123 124 ! Computation of the abundance of the preys 124 125 ! A threshold can be specified in the namelist 125 126 ! -------------------------------------------- 126 zcompadi = MAX( ( tr b(ji,jj,jk,jpdia) - xthresh2dia ), 0.e0 )127 zcompaz = MAX( ( tr b(ji,jj,jk,jpzoo) - xthresh2zoo ), 0.e0 )128 zcompaph = MAX( ( tr b(ji,jj,jk,jpphy) - xthresh2phy ), 0.e0 )129 zcompapoc = MAX( ( tr b(ji,jj,jk,jppoc) - xthresh2poc ), 0.e0 )130 zcompames = MAX( ( tr b(ji,jj,jk,jpmes) - xthresh2mes ), 0.e0 )127 zcompadi = MAX( ( tr(ji,jj,jk,jpdia,Kbb) - xthresh2dia ), 0.e0 ) 128 zcompaz = MAX( ( tr(ji,jj,jk,jpzoo,Kbb) - xthresh2zoo ), 0.e0 ) 129 zcompaph = MAX( ( tr(ji,jj,jk,jpphy,Kbb) - xthresh2phy ), 0.e0 ) 130 zcompapoc = MAX( ( tr(ji,jj,jk,jppoc,Kbb) - xthresh2poc ), 0.e0 ) 131 zcompames = MAX( ( tr(ji,jj,jk,jpmes,Kbb) - xthresh2mes ), 0.e0 ) 131 132 132 133 ! Mesozooplankton grazing … … 136 137 zfoodlim = MAX( 0., zfood - MIN( 0.5 * zfood, xthresh2 ) ) 137 138 zdenom = zfoodlim / ( xkgraz2 + zfoodlim ) 138 zgraze2 = grazrat2 * xstep * tgfunc2(ji,jj,jk) * tr b(ji,jj,jk,jpmes) * (1. - nitrfac(ji,jj,jk))139 zgraze2 = grazrat2 * xstep * tgfunc2(ji,jj,jk) * tr(ji,jj,jk,jpmes,Kbb) * (1. - nitrfac(ji,jj,jk)) 139 140 140 141 ! An active switching parameterization is used here. … … 161 162 ! ------------------------------------------------------ 162 163 zgrazdc = zgraze2 * ztmp4 * zdenom 163 zgrazdn = zgrazdc * tr b(ji,jj,jk,jpndi) / ( trb(ji,jj,jk,jpdia) + rtrn)164 zgrazdp = zgrazdc * tr b(ji,jj,jk,jppdi) / ( trb(ji,jj,jk,jpdia) + rtrn)165 zgrazdf = zgrazdc * tr b(ji,jj,jk,jpdfe) / ( trb(ji,jj,jk,jpdia) + rtrn)164 zgrazdn = zgrazdc * tr(ji,jj,jk,jpndi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn) 165 zgrazdp = zgrazdc * tr(ji,jj,jk,jppdi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn) 166 zgrazdf = zgrazdc * tr(ji,jj,jk,jpdfe,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn) 166 167 zgrazz = zgraze2 * ztmp5 * zdenom 167 168 zgrazm = zgraze2 * ztmp2 * zdenom 168 169 zgraznc = zgraze2 * ztmp1 * zdenom 169 zgraznn = zgraznc * tr b(ji,jj,jk,jpnph) / ( trb(ji,jj,jk,jpphy) + rtrn)170 zgraznp = zgraznc * tr b(ji,jj,jk,jppph) / ( trb(ji,jj,jk,jpphy) + rtrn)171 zgraznf = zgraznc * tr b(ji,jj,jk,jpnfe) / ( trb(ji,jj,jk,jpphy) + rtrn)170 zgraznn = zgraznc * tr(ji,jj,jk,jpnph,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn) 171 zgraznp = zgraznc * tr(ji,jj,jk,jppph,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn) 172 zgraznf = zgraznc * tr(ji,jj,jk,jpnfe,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn) 172 173 zgrazpoc = zgraze2 * ztmp3 * zdenom 173 zgrazpon = zgrazpoc * tr b(ji,jj,jk,jppon) / ( trb(ji,jj,jk,jppoc) + rtrn)174 zgrazpop = zgrazpoc * tr b(ji,jj,jk,jppop) / ( trb(ji,jj,jk,jppoc) + rtrn)175 zgrazpof = zgrazpoc * tr b(ji,jj,jk,jpsfe) / ( trb(ji,jj,jk,jppoc) + rtrn)174 zgrazpon = zgrazpoc * tr(ji,jj,jk,jppon,Kbb) / ( tr(ji,jj,jk,jppoc,Kbb) + rtrn) 175 zgrazpop = zgrazpoc * tr(ji,jj,jk,jppop,Kbb) / ( tr(ji,jj,jk,jppoc,Kbb) + rtrn) 176 zgrazpof = zgrazpoc * tr(ji,jj,jk,jpsfe,Kbb) / ( tr(ji,jj,jk,jppoc,Kbb) + rtrn) 176 177 177 178 ! Mesozooplankton flux feeding on GOC 178 179 ! ---------------------------------- 179 180 zgrazffeg = grazflux * xstep * wsbio4(ji,jj,jk) & 180 & * tgfunc2(ji,jj,jk) * tr b(ji,jj,jk,jpgoc) * trb(ji,jj,jk,jpmes) &181 & * tgfunc2(ji,jj,jk) * tr(ji,jj,jk,jpgoc,Kbb) * tr(ji,jj,jk,jpmes,Kbb) & 181 182 & * (1. - nitrfac(ji,jj,jk)) 182 zgrazfffg = zgrazffeg * tr b(ji,jj,jk,jpbfe) / (trb(ji,jj,jk,jpgoc) + rtrn)183 zgrazffng = zgrazffeg * tr b(ji,jj,jk,jpgon) / (trb(ji,jj,jk,jpgoc) + rtrn)184 zgrazffpg = zgrazffeg * tr b(ji,jj,jk,jpgop) / (trb(ji,jj,jk,jpgoc) + rtrn)183 zgrazfffg = zgrazffeg * tr(ji,jj,jk,jpbfe,Kbb) / (tr(ji,jj,jk,jpgoc,Kbb) + rtrn) 184 zgrazffng = zgrazffeg * tr(ji,jj,jk,jpgon,Kbb) / (tr(ji,jj,jk,jpgoc,Kbb) + rtrn) 185 zgrazffpg = zgrazffeg * tr(ji,jj,jk,jpgop,Kbb) / (tr(ji,jj,jk,jpgoc,Kbb) + rtrn) 185 186 zgrazffep = grazflux * xstep * wsbio3(ji,jj,jk) & 186 & * tgfunc2(ji,jj,jk) * tr b(ji,jj,jk,jppoc) * trb(ji,jj,jk,jpmes) &187 & * tgfunc2(ji,jj,jk) * tr(ji,jj,jk,jppoc,Kbb) * tr(ji,jj,jk,jpmes,Kbb) & 187 188 & * (1. - nitrfac(ji,jj,jk)) 188 zgrazfffp = zgrazffep * tr b(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)189 zgrazffnp = zgrazffep * tr b(ji,jj,jk,jppon) / (trb(ji,jj,jk,jppoc) + rtrn)190 zgrazffpp = zgrazffep * tr b(ji,jj,jk,jppop) / (trb(ji,jj,jk,jppoc) + rtrn)189 zgrazfffp = zgrazffep * tr(ji,jj,jk,jpsfe,Kbb) / (tr(ji,jj,jk,jppoc,Kbb) + rtrn) 190 zgrazffnp = zgrazffep * tr(ji,jj,jk,jppon,Kbb) / (tr(ji,jj,jk,jppoc,Kbb) + rtrn) 191 zgrazffpp = zgrazffep * tr(ji,jj,jk,jppop,Kbb) / (tr(ji,jj,jk,jppoc,Kbb) + rtrn) 191 192 ! 192 193 zgraztotc = zgrazdc + zgrazz + zgraznc + zgrazm + zgrazpoc + zgrazffep + zgrazffeg … … 200 201 ! since they are more porous (marine snow instead of fecal pellets) 201 202 ! ---------------------------------------------------------------- 202 zratio = tr b(ji,jj,jk,jpgsi) / ( trb(ji,jj,jk,jpgoc) + rtrn )203 zratio = tr(ji,jj,jk,jpgsi,Kbb) / ( tr(ji,jj,jk,jpgoc,Kbb) + rtrn ) 203 204 zratio2 = zratio * zratio 204 205 zfracc = zproport * grazflux * xstep * wsbio4(ji,jj,jk) & 205 & * tr b(ji,jj,jk,jpgoc) * trb(ji,jj,jk,jpmes) &206 & * tr(ji,jj,jk,jpgoc,Kbb) * tr(ji,jj,jk,jpmes,Kbb) & 206 207 & * ( 0.2 + 3.8 * zratio2 / ( 1.**2 + zratio2 ) ) 207 zfracfe = zfracc * tr b(ji,jj,jk,jpbfe) / (trb(ji,jj,jk,jpgoc) + rtrn)208 zfracn = zfracc * tr b(ji,jj,jk,jpgon) / (trb(ji,jj,jk,jpgoc) + rtrn)209 zfracp = zfracc * tr b(ji,jj,jk,jpgop) / (trb(ji,jj,jk,jpgoc) + rtrn)208 zfracfe = zfracc * tr(ji,jj,jk,jpbfe,Kbb) / (tr(ji,jj,jk,jpgoc,Kbb) + rtrn) 209 zfracn = zfracc * tr(ji,jj,jk,jpgon,Kbb) / (tr(ji,jj,jk,jpgoc,Kbb) + rtrn) 210 zfracp = zfracc * tr(ji,jj,jk,jpgop,Kbb) / (tr(ji,jj,jk,jpgoc,Kbb) + rtrn) 210 211 211 212 zgrazffep = zproport * zgrazffep ; zgrazffeg = zproport * zgrazffeg … … 296 297 ! sinks 297 298 ! -------------------------------------------------------------- 298 tr a(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) + zgrarep299 tr a(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + zgraren300 tr a(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zgradoc299 tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) + zgrarep 300 tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zgraren 301 tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zgradoc 301 302 ! 302 303 IF( ln_ligand ) THEN 303 tr a(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zgradoc * ldocz304 tr(ji,jj,jk,jplgw,Krhs) = tr(ji,jj,jk,jplgw,Krhs) + zgradoc * ldocz 304 305 zz2ligprod(ji,jj,jk) = zgradoc * ldocz 305 306 ENDIF 306 307 ! 307 tr a(ji,jj,jk,jpdon) = tra(ji,jj,jk,jpdon) + zgradon308 tr a(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zgradop309 tr a(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * zgrarem310 tr a(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zgraref308 tr(ji,jj,jk,jpdon,Krhs) = tr(ji,jj,jk,jpdon,Krhs) + zgradon 309 tr(ji,jj,jk,jpdop,Krhs) = tr(ji,jj,jk,jpdop,Krhs) + zgradop 310 tr(ji,jj,jk,jpoxy,Krhs) = tr(ji,jj,jk,jpoxy,Krhs) - o2ut * zgrarem 311 tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) + zgraref 311 312 zfezoo2(ji,jj,jk) = zgraref 312 tr a(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zgrarem313 tr a(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * zgraren314 tr a(ji,jj,jk,jpmes) = tra(ji,jj,jk,jpmes) + zepsherv * zgraztotc - zrespirc &313 tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) + zgrarem 314 tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zgraren 315 tr(ji,jj,jk,jpmes,Krhs) = tr(ji,jj,jk,jpmes,Krhs) + zepsherv * zgraztotc - zrespirc & 315 316 & - ztortz - zgrazm 316 tr a(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) - zgrazdc317 tr a(ji,jj,jk,jpndi) = tra(ji,jj,jk,jpndi) - zgrazdn318 tr a(ji,jj,jk,jppdi) = tra(ji,jj,jk,jppdi) - zgrazdp319 tr a(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) - zgrazdf320 tr a(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) - zgrazz321 tr a(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) - zgraznc322 tr a(ji,jj,jk,jpnph) = tra(ji,jj,jk,jpnph) - zgraznn323 tr a(ji,jj,jk,jppph) = tra(ji,jj,jk,jppph) - zgraznp324 tr a(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zgraznf325 tr a(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) - zgraznc * trb(ji,jj,jk,jpnch) / ( trb(ji,jj,jk,jpphy) + rtrn )326 tr a(ji,jj,jk,jpdch) = tra(ji,jj,jk,jpdch) - zgrazdc * trb(ji,jj,jk,jpdch) / ( trb(ji,jj,jk,jpdia) + rtrn )327 tr a(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zgrazdc * trb(ji,jj,jk,jpdsi) / ( trb(ji,jj,jk,jpdia) + rtrn )328 tr a(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) + zgrazdc * trb(ji,jj,jk,jpdsi) / ( trb(ji,jj,jk,jpdia) + rtrn )329 330 tr a(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zgrazpoc - zgrazffep + zfracc317 tr(ji,jj,jk,jpdia,Krhs) = tr(ji,jj,jk,jpdia,Krhs) - zgrazdc 318 tr(ji,jj,jk,jpndi,Krhs) = tr(ji,jj,jk,jpndi,Krhs) - zgrazdn 319 tr(ji,jj,jk,jppdi,Krhs) = tr(ji,jj,jk,jppdi,Krhs) - zgrazdp 320 tr(ji,jj,jk,jpdfe,Krhs) = tr(ji,jj,jk,jpdfe,Krhs) - zgrazdf 321 tr(ji,jj,jk,jpzoo,Krhs) = tr(ji,jj,jk,jpzoo,Krhs) - zgrazz 322 tr(ji,jj,jk,jpphy,Krhs) = tr(ji,jj,jk,jpphy,Krhs) - zgraznc 323 tr(ji,jj,jk,jpnph,Krhs) = tr(ji,jj,jk,jpnph,Krhs) - zgraznn 324 tr(ji,jj,jk,jppph,Krhs) = tr(ji,jj,jk,jppph,Krhs) - zgraznp 325 tr(ji,jj,jk,jpnfe,Krhs) = tr(ji,jj,jk,jpnfe,Krhs) - zgraznf 326 tr(ji,jj,jk,jpnch,Krhs) = tr(ji,jj,jk,jpnch,Krhs) - zgraznc * tr(ji,jj,jk,jpnch,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn ) 327 tr(ji,jj,jk,jpdch,Krhs) = tr(ji,jj,jk,jpdch,Krhs) - zgrazdc * tr(ji,jj,jk,jpdch,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 328 tr(ji,jj,jk,jpdsi,Krhs) = tr(ji,jj,jk,jpdsi,Krhs) - zgrazdc * tr(ji,jj,jk,jpdsi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 329 tr(ji,jj,jk,jpgsi,Krhs) = tr(ji,jj,jk,jpgsi,Krhs) + zgrazdc * tr(ji,jj,jk,jpdsi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 330 331 tr(ji,jj,jk,jppoc,Krhs) = tr(ji,jj,jk,jppoc,Krhs) - zgrazpoc - zgrazffep + zfracc 331 332 prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + zfracc 332 333 conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zgrazpoc - zgrazffep 333 tr a(ji,jj,jk,jppon) = tra(ji,jj,jk,jppon) - zgrazpon - zgrazffnp + zfracn334 tr a(ji,jj,jk,jppop) = tra(ji,jj,jk,jppop) - zgrazpop - zgrazffpp + zfracp335 tr a(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) - zgrazffeg + zgrapoc - zfracc334 tr(ji,jj,jk,jppon,Krhs) = tr(ji,jj,jk,jppon,Krhs) - zgrazpon - zgrazffnp + zfracn 335 tr(ji,jj,jk,jppop,Krhs) = tr(ji,jj,jk,jppop,Krhs) - zgrazpop - zgrazffpp + zfracp 336 tr(ji,jj,jk,jpgoc,Krhs) = tr(ji,jj,jk,jpgoc,Krhs) - zgrazffeg + zgrapoc - zfracc 336 337 prodgoc(ji,jj,jk) = prodgoc(ji,jj,jk) + zgrapoc 337 338 consgoc(ji,jj,jk) = consgoc(ji,jj,jk) - zgrazffeg - zfracc 338 tr a(ji,jj,jk,jpgon) = tra(ji,jj,jk,jpgon) - zgrazffng + zgrapon - zfracn339 tr a(ji,jj,jk,jpgop) = tra(ji,jj,jk,jpgop) - zgrazffpg + zgrapop - zfracp340 tr a(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zgrazpof - zgrazfffp + zfracfe341 tr a(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) - zgrazfffg + zgrapof - zfracfe342 zfracal = tr b(ji,jj,jk,jpcal) / ( trb(ji,jj,jk,jpgoc) + rtrn )339 tr(ji,jj,jk,jpgon,Krhs) = tr(ji,jj,jk,jpgon,Krhs) - zgrazffng + zgrapon - zfracn 340 tr(ji,jj,jk,jpgop,Krhs) = tr(ji,jj,jk,jpgop,Krhs) - zgrazffpg + zgrapop - zfracp 341 tr(ji,jj,jk,jpsfe,Krhs) = tr(ji,jj,jk,jpsfe,Krhs) - zgrazpof - zgrazfffp + zfracfe 342 tr(ji,jj,jk,jpbfe,Krhs) = tr(ji,jj,jk,jpbfe,Krhs) - zgrazfffg + zgrapof - zfracfe 343 zfracal = tr(ji,jj,jk,jpcal,Kbb) / ( tr(ji,jj,jk,jpgoc,Kbb) + rtrn ) 343 344 zgrazcal = zgrazffeg * (1. - part2) * zfracal 344 345 … … 348 349 prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo) 349 350 zprcaca = part2 * zprcaca 350 tr a(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zgrazcal - zprcaca351 tr a(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + 2. * ( zgrazcal - zprcaca )352 tr a(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) - zgrazcal + zprcaca351 tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) + zgrazcal - zprcaca 352 tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + 2. * ( zgrazcal - zprcaca ) 353 tr(ji,jj,jk,jpcal,Krhs) = tr(ji,jj,jk,jpcal,Krhs) - zgrazcal + zprcaca 353 354 END DO 354 355 END DO … … 379 380 WRITE(charout, FMT="('meso')") 380 381 CALL prt_ctl_trc_info(charout) 381 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)382 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 382 383 ENDIF 383 384 !
Note: See TracChangeset
for help on using the changeset viewer.