Changeset 773
- Timestamp:
- 2007-12-17T17:55:32+01:00 (16 years ago)
- Location:
- branches/dev_001_GM/NEMO/TOP_SRC
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_GM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
r772 r773 19 19 USE trccfc ! CFC sms trends 20 20 USE sms ! Source Minus Sink variables 21 USE oce_trc ! ocean variables 22 USE trp_trc ! 23 21 24 USE iom 22 25 … … 179 182 DO jj = 2 , jpjm1 180 183 DO ji = 2, jpim1 181 zcoef = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,1) )* tmask(ji,jj,1) * tmask_i(ji,jj) * raass184 zcoef = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,1) * tmask(ji,jj,1) * tmask_i(ji,jj) * raass 182 185 rivpo4input = rivpo4input + rivinp(ji,jj) * zcoef 183 186 rivalkinput = rivalkinput + cotdep(ji,jj) * zcoef -
branches/dev_001_GM/NEMO/TOP_SRC/PISCES/trclsm_pisces.F90
r772 r773 1 !!---------------------------------------------------------------------- 2 !! *** trclsm.pisces.h90 *** 3 !! TOP : Definition some run parameter for PISCES biological model 4 !!---------------------------------------------------------------------- 1 MODULE trclsm_pisces 2 !!====================================================================== 3 !! *** MODULE trclsm_lobster *** 4 !! TOP : initialisation of some run parameters for PISCES bio-model 5 !!====================================================================== 5 6 !! History : - ! 1999-10 (M.A. Foujols, M. Levy) original code 6 7 !! - ! 2000-01 (L. Bopp) hamocc3, p3zd -
branches/dev_001_GM/NEMO/TOP_SRC/TRP/trcadv_tvd.F90
r772 r773 72 72 !! * Local declarations 73 73 INTEGER :: ji, jj, jk,jn ! dummy loop indices 74 REAL(wp) :: ztra ! temporary scalar75 74 76 75 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & -
branches/dev_001_GM/NEMO/TOP_SRC/TRP/trcnxt.F90
r772 r773 70 70 INTEGER, INTENT( in ) :: kt ! ocean time-step index 71 71 !! * Local declarations 72 INTEGER :: ji, jj, jk,jn ! dummy loop indices73 REAL(wp) :: zfact , ztra! temporary scalar72 INTEGER :: ji, jj, jk, jn ! dummy loop indices 73 REAL(wp) :: zfact ! temporary scalar 74 74 CHARACTER (len=22) :: charout 75 75 !!---------------------------------------------------------------------- -
branches/dev_001_GM/NEMO/TOP_SRC/TRP/trctrp.F90
r772 r773 15 15 USE trc ! ocean passive tracers variables 16 16 USE trctrp_lec ! passive tracers transport parameters 17 USE prtctl_trc ! Print control for debbuging 17 18 18 19 USE trcbbl ! bottom boundary layer (trc_bbl routine) … … 67 68 !!---------------------------------------------------------------------- 68 69 INTEGER, INTENT( in ) :: kt ! ocean time-step index 70 !! 71 CHARACTER (len=25) :: charout 69 72 !! --------------------------------------------------------------------- 70 73 -
branches/dev_001_GM/NEMO/TOP_SRC/initrc.F90
r772 r773 51 51 !! or read data or analytical formulation 52 52 !!--------------------------------------------------------------------- 53 INTEGER :: j i, jj, jk, jn ! dummy loop indices53 INTEGER :: jk, jn ! dummy loop indices 54 54 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zbt ! workspace: masked grid volume 55 55 !!--------------------------------------------------------------------- … … 64 64 END DO 65 65 #if defined key_off_degrad 66 zbt(:,:,:) = zbt(:,:,:) * facvol( ji,jj,jk) ! degrad option: reduction by facvol66 zbt(:,:,:) = zbt(:,:,:) * facvol(:,:,jk) ! degrad option: reduction by facvol 67 67 #endif 68 68 -
branches/dev_001_GM/NEMO/TOP_SRC/trcctl.F90
r772 r773 83 83 # if defined key_lobster 84 84 # include "trcctl.lobster1.h90" 85 85 86 # elif defined key_pisces 86 87 # include "trcctl.pisces.h90" 88 87 89 # elif defined key_cfc 88 90 # include "trcctl.cfc.h90" 91 89 92 # else 90 93 IF(lwp) WRITE (numout,*) ' No Source/Sink ' -
branches/dev_001_GM/NEMO/TOP_SRC/trcdtr.F90
r772 r773 41 41 !! ** Method : 42 42 !!--------------------------------------------------------------------- 43 INTEGER :: ji,jj,jk,jn 44 #if defined key_pisces 43 INTEGER :: jn 44 # if defined key_pisces 45 INTEGER :: ji, jj, jk 45 46 REAL(wp) :: alka0, oxyg0, calc0, bioma0 46 47 REAL(wp) :: silic1, po4, no3, caralk, bicarb, co3 47 # endif48 # endif 48 49 !!--------------------------------------------------------------------- 49 50 … … 52 53 IF(lwp) WRITE(numout,*) '~~~~~~~' 53 54 54 #if defined key_cfc 55 ! CFC initialisation 56 trn(:,:,:,:) = 0.e0 55 IF( lk_cfc ) THEN 56 DO jn = jp_cfc0,jp_cfc1 57 trn(:,:,:,jn) = 0.e0 ! CFC initialisation to zera 58 END DO 59 ENDIF 57 60 58 # elif defined key_pisces61 # if defined key_pisces 59 62 ! PISCES initialisation 60 63 ! --------------------- … … 75 78 trn(:,:,:,jppo4) = po4 76 79 trn(:,:,:,jppoc) = bioma0 77 # if ! defined key_kriest80 # if ! defined key_kriest 78 81 trn(:,:,:,jpgoc) = bioma0 79 82 trn(:,:,:,jpbfe) = bioma0 * 5.e-6 80 # else83 # else 81 84 trn(:,:,:,jpnum) = bioma0 / ( 6. *xkr_massp ) 82 # endif85 # endif 83 86 trn(:,:,:,jpsil) = silic1 84 87 trn(:,:,:,jpbsi) = bioma0 * 0.15 … … 119 122 IF(lwp) WRITE(numout,*) ' ' 120 123 121 # elif defined key_lobster && ( defined key_eel_r6 || defined key_eel_r2 )124 # elif defined key_lobster && ( defined key_eel_r6 || defined key_eel_r2 ) 122 125 ! LOBSTER initialisation for EEL 123 126 ! ---------------------- … … 195 198 trn(:,:,30,jpno3) = 20.01 * tmask(:,:,30) 196 199 197 # elif defined key_lobster && defined key_gyre200 # elif defined key_lobster && defined key_gyre 198 201 ! LOBSTER initialisation for GYRE 199 202 ! ---------------------- … … 216 219 END DO 217 220 218 # else221 # else 219 222 ! Default case 220 223 ! ------------ … … 223 226 END DO 224 227 225 # endif226 227 # if defined key_dtatrc228 # endif 229 230 # if defined key_dtatrc 228 231 ! Initialization of tracer from a file that may also be used for damping 229 232 CALL dta_trc( nittrc000 ) … … 231 234 IF( lutini(jn) ) trn(:,:,:,jk) = trdta(:,:,:,jn) * tmask(:,:,:) ! initialisation from file if required 232 235 END DO 233 # endif236 # endif 234 237 235 238 ! before field : … … 237 240 trb(:,:,:,:) = trn(:,:,:,:) 238 241 239 # if defined key_lobster242 # if defined key_lobster 240 243 ! initialize the POC in sediments 241 244 sedpocb(:,:) = 0.e0 242 245 sedpocn(:,:) = 0.e0 243 246 sedpoca(:,:) = 0.e0 244 # endif247 # endif 245 248 ! 246 249 END SUBROUTINE trc_dtr
Note: See TracChangeset
for help on using the changeset viewer.