Changeset 3294 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zrem.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zrem.F90 (modified) (20 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zrem.F90
r2773 r3294 6 6 !! History : 1.0 ! 2004 (O. Aumont) Original code 7 7 !! 2.0 ! 2007-12 (C. Ethe, G. Madec) F90 8 !! 3.4 ! 2011-06 (O. Aumont, C. Ethe) Quota model for iron 8 9 !!---------------------------------------------------------------------- 9 10 #if defined key_pisces … … 12 13 !! 'key_pisces' PISCES bio-model 13 14 !!---------------------------------------------------------------------- 14 !! p4z_rem : Compute remineralization/scavenging of organic compounds 15 !!---------------------------------------------------------------------- 16 USE trc 17 USE oce_trc ! 18 USE sms_pisces ! 19 USE prtctl_trc 20 USE p4zint 21 USE p4zopt 22 USE p4zmeso 23 USE p4zprod 24 USE p4zche 15 !! p4z_rem : Compute remineralization/scavenging of organic compounds 16 !! p4z_rem_init : Initialisation of parameters for remineralisation 17 !! p4z_rem_alloc : Allocate remineralisation variables 18 !!---------------------------------------------------------------------- 19 USE oce_trc ! shared variables between ocean and passive tracers 20 USE trc ! passive tracers common variables 21 USE sms_pisces ! PISCES Source Minus Sink variables 22 USE p4zopt ! optical model 23 USE p4zche ! chemical model 24 USE p4zprod ! Growth rate of the 2 phyto groups 25 USE p4zmeso ! Sources and sinks of mesozooplankton 26 USE p4zint ! interpolation and computation of various fields 27 USE prtctl_trc ! print control for debugging 25 28 26 29 IMPLICIT NONE … … 31 34 PUBLIC p4z_rem_alloc 32 35 33 REAL(wp), PUBLIC :: & 34 xremik = 0.3_wp , & !: 35 xremip = 0.025_wp , & !: 36 nitrif = 0.05_wp , & !: 37 xsirem = 0.015_wp , & !: 38 xlam1 = 0.005_wp , & !: 39 oxymin = 1.e-6_wp !: 40 41 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: denitr !: denitrification array 36 !! * Shared module variables 37 REAL(wp), PUBLIC :: xremik = 0.3_wp !: remineralisation rate of POC 38 REAL(wp), PUBLIC :: xremip = 0.025_wp !: remineralisation rate of DOC 39 REAL(wp), PUBLIC :: nitrif = 0.05_wp !: NH4 nitrification rate 40 REAL(wp), PUBLIC :: xsirem = 0.003_wp !: remineralisation rate of POC 41 REAL(wp), PUBLIC :: xsiremlab = 0.025_wp !: fast remineralisation rate of POC 42 REAL(wp), PUBLIC :: xsilab = 0.31_wp !: fraction of labile biogenic silica 43 REAL(wp), PUBLIC :: xlam1 = 0.005_wp !: scavenging rate of Iron 44 REAL(wp), PUBLIC :: oxymin = 1.e-6_wp !: halk saturation constant for anoxia 45 REAL(wp), PUBLIC :: ligand = 0.6E-9_wp !: ligand concentration in the ocean 46 47 48 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: denitr !: denitrification array 49 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: denitnh4 !: - - - - - 42 50 43 51 … … 59 67 !! ** Method : - ??? 60 68 !!--------------------------------------------------------------------- 61 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released62 USE wrk_nemo, ONLY: ztempbac => wrk_2d_163 USE wrk_nemo, ONLY: zdepbac => wrk_3d_2 , zolimi => wrk_3d_364 69 ! 65 70 INTEGER, INTENT(in) :: kt ! ocean time step 66 71 ! 67 72 INTEGER :: ji, jj, jk 68 REAL(wp) :: zremip, zremik , zlam1b 73 REAL(wp) :: zremip, zremik , zlam1b, zdepbac2 69 74 REAL(wp) :: zkeq , zfeequi, zsiremin, zfesatur 70 REAL(wp) :: zsatur, zsatur2, znusil 75 REAL(wp) :: zsatur, zsatur2, znusil, zdep, zfactdep 71 76 REAL(wp) :: zbactfer, zorem, zorem2, zofer 72 REAL(wp) :: zosil, zdenom1, zscave, zaggdfe 77 REAL(wp) :: zosil, zdenom1, zscave, zaggdfe, zcoag 73 78 #if ! defined key_kriest 74 79 REAL(wp) :: zofer2, zdenom, zdenom2 … … 76 81 REAL(wp) :: zlamfac, zonitr, zstep 77 82 CHARACTER (len=25) :: charout 83 REAL(wp), POINTER, DIMENSION(:,: ) :: ztempbac 84 REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zolimi, zolimi2 78 85 !!--------------------------------------------------------------------- 79 80 IF( wrk_in_use(2, 1) .OR. wrk_in_use(3, 2,3) ) THEN 81 CALL ctl_stop('p4z_rem: requested workspace arrays unavailable') ; RETURN 82 ENDIF 86 ! 87 IF( nn_timing == 1 ) CALL timing_start('p4z_rem') 88 ! 89 ! Allocate temporary workspace 90 CALL wrk_alloc( jpi, jpj, ztempbac ) 91 CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zolimi, zolimi2 ) 83 92 84 93 ! Initialisation of temprary arrys 85 94 zdepbac (:,:,:) = 0._wp 86 95 zolimi (:,:,:) = 0._wp 96 zolimi2 (:,:,:) = 0._wp 87 97 ztempbac(:,:) = 0._wp 88 98 … … 93 103 DO jj = 1, jpj 94 104 DO ji = 1, jpi 95 IF( fsdept(ji,jj,jk) < 120. ) THEN 105 zdep = MAX( hmld(ji,jj), heup(ji,jj) ) 106 IF( fsdept(ji,jj,jk) < zdep ) THEN 96 107 zdepbac(ji,jj,jk) = MIN( 0.7 * ( trn(ji,jj,jk,jpzoo) + 2.* trn(ji,jj,jk,jpmes) ), 4.e-6 ) 97 108 ztempbac(ji,jj) = zdepbac(ji,jj,jk) 98 109 ELSE 99 zdepbac(ji,jj,jk) = MIN( 1., 120./ fsdept(ji,jj,jk) ) * ztempbac(ji,jj)110 zdepbac(ji,jj,jk) = MIN( 1., zdep / fsdept(ji,jj,jk) ) * ztempbac(ji,jj) 100 111 ENDIF 101 112 END DO … … 117 128 DO jj = 1, jpj 118 129 DO ji = 1, jpi 130 zstep = xstep 119 131 # if defined key_degrad 120 zstep = xstep * facvol(ji,jj,jk) 121 # else 122 zstep = xstep 132 zstep = zstep * facvol(ji,jj,jk) 123 133 # endif 124 134 ! DOC ammonification. Depends on depth, phytoplankton biomass … … 126 136 ! of the bacterial activity. 127 137 zremik = xremik * zstep / 1.e-6 * xlimbac(ji,jj,jk) * zdepbac(ji,jj,jk) 128 zremik = MAX( zremik, 5.5e-4 * xstep ) 129 138 zremik = MAX( zremik, 2.e-4 * xstep ) 130 139 ! Ammonification in oxic waters with oxygen consumption 131 140 ! ----------------------------------------------------- 132 zolimi(ji,jj,jk) = MIN( ( trn(ji,jj,jk,jpoxy) - rtrn ) / o2ut, & 133 & zremik * ( 1.- nitrfac(ji,jj,jk) ) * trn(ji,jj,jk,jpdoc) ) 134 141 zolimi (ji,jj,jk) = zremik * ( 1.- nitrfac(ji,jj,jk) ) * trn(ji,jj,jk,jpdoc) 142 zolimi2(ji,jj,jk) = MIN( ( trn(ji,jj,jk,jpoxy) - rtrn ) / o2ut, zolimi(ji,jj,jk) ) 135 143 ! Ammonification in suboxic waters with denitrification 136 144 ! ------------------------------------------------------- 137 denitr(ji,jj,jk) = MIN( ( trn(ji,jj,jk,jpno3) - rtrn ) / rdenit, &145 denitr(ji,jj,jk) = MIN( ( trn(ji,jj,jk,jpno3) - rtrn ) / rdenit, & 138 146 & zremik * nitrfac(ji,jj,jk) * trn(ji,jj,jk,jpdoc) ) 139 END DO 140 END DO 141 END DO 142 143 DO jk = 1, jpkm1 144 DO jj = 1, jpj 145 DO ji = 1, jpi 147 ! 146 148 zolimi (ji,jj,jk) = MAX( 0.e0, zolimi (ji,jj,jk) ) 149 zolimi2(ji,jj,jk) = MAX( 0.e0, zolimi2(ji,jj,jk) ) 147 150 denitr (ji,jj,jk) = MAX( 0.e0, denitr (ji,jj,jk) ) 148 END DO 149 END DO 150 END DO 151 152 DO jk = 1, jpkm1 153 DO jj = 1, jpj 154 DO ji = 1, jpi 151 ! 152 END DO 153 END DO 154 END DO 155 156 157 DO jk = 1, jpkm1 158 DO jj = 1, jpj 159 DO ji = 1, jpi 160 zstep = xstep 155 161 # if defined key_degrad 156 zstep = xstep * facvol(ji,jj,jk) 157 # else 158 zstep = xstep 162 zstep = zstep * facvol(ji,jj,jk) 159 163 # endif 160 164 ! NH4 nitrification to NO3. Ceased for oxygen concentrations 161 165 ! below 2 umol/L. Inhibited at strong light 162 166 ! ---------------------------------------------------------- 163 zonitr = 164 167 zonitr =nitrif * zstep * trn(ji,jj,jk,jpnh4) / ( 1.+ emoy(ji,jj,jk) ) * ( 1.- nitrfac(ji,jj,jk) ) 168 denitnh4(ji,jj,jk) = nitrif * zstep * trn(ji,jj,jk,jpnh4) * nitrfac(ji,jj,jk) 165 169 ! Update of the tracers trends 166 170 ! ---------------------------- 167 168 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) - zonitr 169 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zonitr 171 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) - zonitr - denitnh4(ji,jj,jk) 172 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zonitr - rdenita * denitnh4(ji,jj,jk) 170 173 tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2nit * zonitr 171 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - rno3 * zonitr 172 174 tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2 * rno3 * zonitr + rno3 * ( rdenita - 1. ) * denitnh4(ji,jj,jk) 173 175 END DO 174 176 END DO … … 189 191 ! studies (especially at Papa) have shown this uptake to be significant 190 192 ! ---------------------------------------------------------- 191 z bactfer = 15.e-6 * rfact2 * 4.* 0.4 * prmax(ji,jj,jk) &192 & * ( xlimphy(ji,jj,jk) * zdepbac(ji,jj,jk))&193 & * ( xlimphy(ji,jj,jk) * zdepbac(ji,jj,jk))&194 & / ( xkgraz2 + zdepbac(ji,jj,jk) )&195 & 193 zdepbac2 = zdepbac(ji,jj,jk) * zdepbac(ji,jj,jk) 194 zbactfer = 20.e-6 * rfact2 * prmax(ji,jj,jk) & 195 & * trn(ji,jj,jk,jpfer) / ( 5E-10 + trn(ji,jj,jk,jpfer) ) & 196 & * zdepbac2 / ( xkgraz2 + zdepbac(ji,jj,jk) ) & 197 & * ( 0.5 + SIGN( 0.5, trn(ji,jj,jk,jpfer) -2.e-11 ) ) 196 198 197 199 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer … … 214 216 DO jj = 1, jpj 215 217 DO ji = 1, jpi 218 zstep = xstep 216 219 # if defined key_degrad 217 zstep = xstep * facvol(ji,jj,jk) 218 # else 219 zstep = xstep 220 zstep = zstep * facvol(ji,jj,jk) 220 221 # endif 221 222 ! POC disaggregation by turbulence and bacterial activity. 222 223 ! ------------------------------------------------------------- 223 zremip = xremip * zstep * tgfunc(ji,jj,jk) * ( 1.- 0. 5* nitrfac(ji,jj,jk) )224 zremip = xremip * zstep * tgfunc(ji,jj,jk) * ( 1.- 0.7 * nitrfac(ji,jj,jk) ) 224 225 225 226 ! POC disaggregation rate is reduced in anoxic zone as shown by … … 266 267 DO jj = 1, jpj 267 268 DO ji = 1, jpi 269 zstep = xstep 268 270 # if defined key_degrad 269 zstep = xstep * facvol(ji,jj,jk) 270 # else 271 zstep = xstep 271 zstep = zstep * facvol(ji,jj,jk) 272 272 # endif 273 273 ! Remineralization rate of BSi depedant on T and saturation 274 274 ! --------------------------------------------------------- 275 zsatur = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 276 zsatur = MAX( rtrn, zsatur ) 277 zsatur2 = zsatur * ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**4 278 znusil = 0.225 * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2**9 279 zsiremin = xsirem * zstep * znusil 280 zosil = zsiremin * trn(ji,jj,jk,jpdsi) 281 275 zsatur = ( sio3eq(ji,jj,jk) - trn(ji,jj,jk,jpsil) ) / ( sio3eq(ji,jj,jk) + rtrn ) 276 zsatur = MAX( rtrn, zsatur ) 277 zsatur2 = zsatur * ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**4 278 znusil = 0.225 * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2**9.25 279 zdep = MAX( hmld(ji,jj), heup(ji,jj) ) 280 zdep = MAX( 0., fsdept(ji,jj,jk) - zdep ) 281 zfactdep = xsilab * EXP(-( xsiremlab - xsirem ) * zdep / wsbio2 ) 282 zsiremin = ( xsiremlab * zfactdep + xsirem * ( 1. - zfactdep ) ) * zstep * znusil 283 zosil = zsiremin * trn(ji,jj,jk,jpdsi) 284 ! 282 285 tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zosil 283 286 tra(ji,jj,jk,jpsil) = tra(ji,jj,jk,jpsil) + zosil … … 293 296 ENDIF 294 297 295 zfesatur = 0.6e-9298 zfesatur = ligand 296 299 !CDIR NOVERRCHK 297 300 DO jk = 1, jpkm1 … … 300 303 !CDIR NOVERRCHK 301 304 DO ji = 1, jpi 305 zstep = xstep 302 306 # if defined key_degrad 303 zstep = xstep * facvol(ji,jj,jk) 304 # else 305 zstep = xstep 307 zstep = zstep * facvol(ji,jj,jk) 306 308 # endif 307 309 ! Compute de different ratios for scavenging of iron … … 312 314 & ( trn(ji,jj,jk,jppoc) + trn(ji,jj,jk,jpdsi) + trn(ji,jj,jk,jpcal) + rtrn ) 313 315 #else 314 zdenom = 1. / ( trn(ji,jj,jk,jppoc) + trn(ji,jj,jk,jpgoc) & 315 & + trn(ji,jj,jk,jpdsi) + trn(ji,jj,jk,jpcal) + rtrn ) 316 316 zdenom = 1. / ( trn(ji,jj,jk,jppoc) + trn(ji,jj,jk,jpgoc) + trn(ji,jj,jk,jpdsi) + trn(ji,jj,jk,jpcal) + rtrn ) 317 317 zdenom1 = trn(ji,jj,jk,jppoc) * zdenom 318 318 zdenom2 = trn(ji,jj,jk,jpgoc) * zdenom … … 337 337 ! Increased scavenging for very high iron concentrations 338 338 ! found near the coasts due to increased lithogenic particles 339 ! and let s say itunknown processes (precipitation, ...)339 ! and let say it is unknown processes (precipitation, ...) 340 340 ! ----------------------------------------------------------- 341 zlam1b = xlam1 * MAX( 0.e0, ( trn(ji,jj,jk,jpfer) * 1.e9 - 1. ) ) 342 zcoag = zfeequi * zlam1b * zstep 341 343 zlamfac = MAX( 0.e0, ( gphit(ji,jj) + 55.) / 30. ) 342 344 zlamfac = MIN( 1. , zlamfac ) 345 zdep = MIN(1., 1000. / fsdept(ji,jj,jk) ) 343 346 #if ! defined key_kriest 344 347 zlam1b = ( 80.* ( trn(ji,jj,jk,jpdoc) + 35.e-6 ) & 345 & + 698.* trn(ji,jj,jk,jppoc) + 1.05e4 * trn(ji,jj,jk,jpgoc) ) & 346 & * xdiss(ji,jj,jk) + 1E-4 * (1.-zlamfac) & 347 & + xlam1 * MAX( 0.e0, ( trn(ji,jj,jk,jpfer) * 1.e9 - 1.) ) 348 #else 349 zlam1b = ( 80.* (trn(ji,jj,jk,jpdoc) + 35E-6) & 348 & + 698.* trn(ji,jj,jk,jppoc) + 1.05e4 * trn(ji,jj,jk,jpgoc) ) & 349 & * xdiss(ji,jj,jk) + 1E-4 * ( 1. - zlamfac ) * zdep 350 #else 351 zlam1b = ( 80.* (trn(ji,jj,jk,jpdoc) + 35E-6) & 350 352 & + 698.* trn(ji,jj,jk,jppoc) ) & 351 & * xdiss(ji,jj,jk) + 1E-4 * (1.-zlamfac) & 352 & + xlam1 * MAX( 0.e0, ( trn(ji,jj,jk,jpfer) * 1.e9 - 1.) ) 353 #endif 354 353 & * xdiss(ji,jj,jk) + 1E-4 * ( 1. - zlamfac ) * zdep 354 #endif 355 355 zaggdfe = zlam1b * zstep * 0.5 * ( trn(ji,jj,jk,jpfer) - zfeequi ) 356 357 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zscave - zaggdfe 358 356 tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zscave - zaggdfe - zcoag 359 357 #if defined key_kriest 360 358 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zscave * zdenom1 … … 378 376 379 377 DO jk = 1, jpkm1 380 tra(:,:,jk,jppo4) = tra(:,:,jk,jppo4) + zolimi (:,:,jk) + denitr(:,:,jk)381 tra(:,:,jk,jpnh4) = tra(:,:,jk,jpnh4) + zolimi (:,:,jk) + denitr(:,:,jk)382 tra(:,:,jk,jpno3) = tra(:,:,jk,jpno3) - denitr (:,:,jk) * rdenit383 tra(:,:,jk,jpdoc) = tra(:,:,jk,jpdoc) - zolimi (:,:,jk) - denitr(:,:,jk)384 tra(:,:,jk,jpoxy) = tra(:,:,jk,jpoxy) - zolimi (:,:,jk) * o2ut385 tra(:,:,jk,jpdic) = tra(:,:,jk,jpdic) + zolimi (:,:,jk) + denitr(:,:,jk)386 tra(:,:,jk,jptal) = tra(:,:,jk,jptal) + denitr(:,:,jk) * rno3 * rdenit378 tra(:,:,jk,jppo4) = tra(:,:,jk,jppo4) + zolimi (:,:,jk) + denitr(:,:,jk) 379 tra(:,:,jk,jpnh4) = tra(:,:,jk,jpnh4) + zolimi (:,:,jk) + denitr(:,:,jk) 380 tra(:,:,jk,jpno3) = tra(:,:,jk,jpno3) - denitr (:,:,jk) * rdenit 381 tra(:,:,jk,jpdoc) = tra(:,:,jk,jpdoc) - zolimi (:,:,jk) - denitr(:,:,jk) 382 tra(:,:,jk,jpoxy) = tra(:,:,jk,jpoxy) - zolimi2(:,:,jk) * o2ut 383 tra(:,:,jk,jpdic) = tra(:,:,jk,jpdic) + zolimi (:,:,jk) + denitr(:,:,jk) 384 tra(:,:,jk,jptal) = tra(:,:,jk,jptal) + rno3 * ( zolimi(:,:,jk) + ( rdenit + 1.) * denitr(:,:,jk) ) 387 385 END DO 388 386 … … 393 391 ENDIF 394 392 ! 395 IF( wrk_not_released(2, 1) .OR. & 396 wrk_not_released(3, 2,3) ) CALL ctl_stop('p4z_rem: failed to release workspace arrays') 393 CALL wrk_dealloc( jpi, jpj, ztempbac ) 394 CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zolimi, zolimi2 ) 395 ! 396 IF( nn_timing == 1 ) CALL timing_stop('p4z_rem') 397 397 ! 398 398 END SUBROUTINE p4z_rem … … 411 411 !! 412 412 !!---------------------------------------------------------------------- 413 NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, x lam1, oxymin414 !!----------------------------------------------------------------------415 416 REWIND( numnat ) ! read numnat417 READ ( numnat , nampisrem )413 NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xsiremlab, xsilab, & 414 & xlam1, oxymin, ligand 415 416 REWIND( numnatp ) ! read numnatp 417 READ ( numnatp, nampisrem ) 418 418 419 419 IF(lwp) THEN ! control print … … 424 424 WRITE(numout,*) ' remineralization rate of DOC xremik =', xremik 425 425 WRITE(numout,*) ' remineralization rate of Si xsirem =', xsirem 426 WRITE(numout,*) ' fast remineralization rate of Si xsiremlab =', xsiremlab 427 WRITE(numout,*) ' fraction of labile biogenic silica xsilab =', xsilab 426 428 WRITE(numout,*) ' scavenging rate of Iron xlam1 =', xlam1 427 429 WRITE(numout,*) ' NH4 nitrification rate nitrif =', nitrif 428 430 WRITE(numout,*) ' halk saturation constant for anoxia oxymin =', oxymin 431 WRITE(numout,*) ' ligand concentration in the ocean ligand =', ligand 429 432 ENDIF 430 433 ! 431 nitrfac(:,:,:) = 0._wp 432 denitr (:,:,:) = 0._wp 434 nitrfac (:,:,:) = 0._wp 435 denitr (:,:,:) = 0._wp 436 denitnh4(:,:,:) = 0._wp 433 437 ! 434 438 END SUBROUTINE p4z_rem_init … … 439 443 !! *** ROUTINE p4z_rem_alloc *** 440 444 !!---------------------------------------------------------------------- 441 ALLOCATE( denitr(jpi,jpj,jpk), STAT=p4z_rem_alloc )445 ALLOCATE( denitr(jpi,jpj,jpk), denitnh4(jpi,jpj,jpk), STAT=p4z_rem_alloc ) 442 446 ! 443 447 IF( p4z_rem_alloc /= 0 ) CALL ctl_warn('p4z_rem_alloc: failed to allocate arrays')
Note: See TracChangeset
for help on using the changeset viewer.
